Institute of Metals Division - Freezing of Semi-Infinite Slab with Time-Dependent Surface Temperature-An Extension of Neumann's Solution

Temperature distribution as well as position of the solidified front is solved by means of "heat balance integral", for the case of freezing a slab with time-dependent surface temperature. Numerical solutions for the cases of linear variation and exponential decay for surface temperature are given. Application also includes solidification of casting. An increased interest in the heat-transfer problem involving a change of phase has been developed in recent years. Practical problems include solidification of a casting, freezing and thawing of soil, ice formation, and others. Despite the importance of the topic, a literature search reveals relatively few solutions which may be extended to practical problems. Due to mathematical complexity, which comes from the nonlinearity of governing differential equations, exact analytical solutions which have been published are limited to a few cases with simple boundary conditions. Although numerical solutions of these problems have become available after highspeed digital computers were introduced, lack of simple and general methods of programming for computation makes this approach impractical. The case of the semi-infinite slab with constant surface temperature was solved in the 19th century. In honor of Neumann, who first found an analytical solution, this kind of problem is generally referred to as the Neumann solution. As has been shown in Fig. 1, the problem to be analyzed is the freezing of a semi-infinite region (x > 0) which is initially liquid at the melting temperature with the surface (x = 0) maintained at time-dependent temperature. This analysis can be applied to the study of many casting processes where most of the heat removal is essentially unidirectional. "Heat balance integral" which has been used by Goodman1 is employed in this analysis. The concept of "heat balance integral" is to define a quantity, d(t), called the thermal layer. Beyond this thermal layer, for all practical purposes, there is no heat transfer, and the region is at the equilibrium temperature. "Heat balance integral" is obtained by integrating the conduction equation with respect to the space variable over this thermal layer d(t). Choosing some plausible temperature distribution, usually a polynomial or a finite trigonometric series in the space variable, this integral equation becomes an ordinary differential equation, with time as the independent variable, which can be solved by classical methods. The thermal layer and "heat balance integral" are equivalent, respectively, to the boundary-layer thickness and the momentum integral in hydrodynamics. A modern account of Prandtl's concept of the boundary layer and von Karman's momentum integral may be found in Ref. 2. Although the problem considered here has been limited initially to uniform melting temperature, the present method is equally applicable to the case with arbitrary temperature distribution in the liquid. STATEMENT OF THE PROBLEM The physical situation and coordinate system used in the problem to be analyzed are shown in Fig. 1. Notations used are listed in the table of nomenclature. Let us assume a semi-infinite region extending over positive x, initially at the uniform melting temperature. At the surface x = 0, a time-dependent temperature, which has a numerical value