TY - JOUR
T1 - Analysis of multidimensional conduction phase change via the enthalpy model
AU - Shamsundar, N.
AU - Sparrow, E. M.
PY - 1975/8
Y1 - 1975/8
N2 - The basis of the enthalpy model for multidimensional phase change problems in media having a distinct phase change temperature is demonstrated, and subsequent numerical applications of the model are carried out. It is shown that the mathematical representation of the enthalpy model is equivalent to the conventional conservation equations in the solid and liquid regions and at the solid-liquid interface. The model is employed in conjunction with a fully implicit finite-difference scheme to solve for solidification in a convectively cooled square container. The implicit scheme was selected because of its ability to accommodate a wide range of the Stefan number Ste. After its accuracy had been established, the solution method was used to obtain results for the local and surface- integrated heat transfer rates, boundary temperatures, solidified fraction, and interface position, all as functions of time. The results are presented with SteFo (Fo = Fourier number) as a correlating parameter, thereby facilitating their use for all Ste values in the range investigated. At low values of the Biot number, the surface-integrated heat transfer rate was relatively constant during the entire solidification period, which is a desirable characteristic for phase change thermal energy storage.
AB - The basis of the enthalpy model for multidimensional phase change problems in media having a distinct phase change temperature is demonstrated, and subsequent numerical applications of the model are carried out. It is shown that the mathematical representation of the enthalpy model is equivalent to the conventional conservation equations in the solid and liquid regions and at the solid-liquid interface. The model is employed in conjunction with a fully implicit finite-difference scheme to solve for solidification in a convectively cooled square container. The implicit scheme was selected because of its ability to accommodate a wide range of the Stefan number Ste. After its accuracy had been established, the solution method was used to obtain results for the local and surface- integrated heat transfer rates, boundary temperatures, solidified fraction, and interface position, all as functions of time. The results are presented with SteFo (Fo = Fourier number) as a correlating parameter, thereby facilitating their use for all Ste values in the range investigated. At low values of the Biot number, the surface-integrated heat transfer rate was relatively constant during the entire solidification period, which is a desirable characteristic for phase change thermal energy storage.
UR - http://www.scopus.com/inward/record.url?scp=0016543869&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=0016543869&partnerID=8YFLogxK
U2 - 10.1115/1.3450375
DO - 10.1115/1.3450375
M3 - Article
AN - SCOPUS:0016543869
VL - 97
SP - 333
EP - 340
JO - Journal of Heat Transfer
JF - Journal of Heat Transfer
SN - 0022-1481
IS - 3
ER -