Scaling theory of melting with natural convection in an enclosure

This study identifies the most basic scales and regimes of the phenomenon of melting with natural convection in an enclosure heated from the side. In the first part of the study the method of scale analysis is used to show that the phenomenon consists of a sequence of four regimes: (a) the pure conduction regime, (b) the mixed regime in which the upper portion of the liquid gap is ruled by convection and the lower portion by conduction, (c) the convection regime and, finally, (d) the last or ‘shrinking solid’ regime. For the first three regimes the scaling theory predicts a Nusselt number vs time curve that has features similar to a van der Waals isotherm, in particular, a clear Nu minimum of order Ra14 at a time SteFo of order Ra−12, where Ste is the liquid superheat Stefan number and Fo the Fourier number based on H. The corresponding average melting front position has an inflexion point at a time of order Ra−12. The theory shows further that during the fourth regime the solid disappears during a SteFo time interval of order Ra−14. The second part of the study consists of numerical experiments the purpose of which is to verify the correctness of the theory constructed in the first part. The numerical simulations are based on the quasi-stationary front approximation and the quasi-steady natural convection assumption. The parametric domain covered by these simulations is 0 ⩽ Ra ⩽ 108, 0 < SteFo < 0.2, Pr = 50 and H/L = 1, where L is the horizontal dimension of the enclosure and Ra the Rayleigh number based on H. Closed form correlations for both Nu and the melting front location time functions are developed by combining the theoretical and numerical conclusions of the study.

[1]  Adrian Bejan,et al.  BUCKLING FLOWS: A NEW FRONTIER IN FLUID MECHANICS , 1987 .

[2]  M. Okada,et al.  Analysis of heat transfer during melting from a vertical wall , 1984 .

[3]  G. Batchelor,et al.  Heat transfer by free convection across a closed cavity between vertical boundaries at different temperatures , 1954 .

[4]  Clifford Goodman,et al.  American Society of Mechanical Engineers , 1988 .

[5]  H. Beer,et al.  Experimental investigation of melting heat transfer with regard to different geometric arrangements , 1984 .

[6]  Adrian Bejan,et al.  Note on Gill's solution for free convection in a vertical enclosure , 1979, Journal of Fluid Mechanics.

[7]  Ashok J. Gadgil,et al.  Analysis of Two-Dimensional Melting in Rectangular Enclosures in Presence of Convection , 1984 .

[8]  A. Bejan,et al.  Natural convection in a differentially heated corner region , 1985 .

[9]  R. Viskanta,et al.  Photographic observation of the solid-liquid interface motion during melting of a solid heated from an isothermal vertical wall , 1978 .

[10]  A. Bejan Convection Heat Transfer , 1984 .

[11]  D. E. Cormack,et al.  Natural convection in a shallow cavity with differentially heated end walls. Part 1. Asymptotic theory , 1974, Journal of Fluid Mechanics.

[12]  B. W. Webb,et al.  On the characteristic length scale for correlating melting heat transfer data , 1985 .

[13]  A. Chapman,et al.  Natural-convection on a finite-size horizontal plate , 1969 .

[14]  J. Lienhard,et al.  On the commonality of equations for natural convection from immersed bodies , 1973 .

[15]  Dominique Gobin,et al.  Melting in Rectangular Enclosures: Experiments and Numerical Simulations , 1985 .

[16]  Ching-Jenq Ho,et al.  Heat Transfer During Melting From an Isothermal Vertical Wall , 1984 .