# Stefan problem

In mathematics and its applications, particularly to phase transitions in matter, a Stefan problem is a particular kind of boundary value problem for a system of partial differential equations (PDE), in which the boundary between the phases can move with time. The classical Stefan problem aims to describe the evolution of the boundary between two phases of a material undergoing a phase change, for example the melting of a solid, such as ice to water. This is accomplished by solving heat equations in both regions, subject to given boundary and initial conditions. At the interface between the phases (in the classical problem) the temperature is set to the phase change temperature. To close the mathematical system a further equation, the Stefan condition, is required. This is an energy balance which defines the position of the moving interface. Note that this evolving boundary is an unknown (hyper-)surface; hence, Stefan problems are examples of free boundary problems.

Analogous problems occur, for example, in the study of porous media flow, mathematical finance and crystal growth from monomer solutions.[1]

## Historical note

The problem is named after Josef Stefan (Jožef Stefan), the Slovenian physicist who introduced the general class of such problems around 1890 in a series of four papers concerning the freezing of the ground and the formation of sea ice.[2] However, some 60 years earlier, in 1831, an equivalent problem, concerning the formation of the Earth's crust, had been studied by Lamé and Clapeyron. Stefan's problem admits a similarity solution, this is often termed the Neumann solution, which was allegedly presented in a series of lectures in the early 1860s.

A comprehensive description of the history of Stefan problems may be found in Rubinstein.[3]

## Premises to the mathematical description

From a mathematical point of view, the phases are merely regions in which the solutions of the underlying PDE are continuous and differentiable up to the order of the PDE. In physical problems such solutions represent properties of the medium for each phase. The moving boundaries (or interfaces) are infinitesimally thin surfaces that separate adjacent phases; therefore, the solutions of the underlying PDE and its derivatives may suffer discontinuities across interfaces.

The underlying PDEs are not valid at the phase change interfaces; therefore, an additional condition—the Stefan condition—is needed to obtain closure. The Stefan condition expresses the local velocity of a moving boundary, as a function of quantities evaluated at either side of the phase boundary, and is usually derived from a physical constraint. In problems of heat transfer with phase change, for instance, conservation of energy dictates that the discontinuity of heat flux at the boundary must be accounted for by the rate of latent heat release (which is proportional to the local velocity of the interface).

## Mathematical formulation

### The one-dimensional one-phase Stefan problem

The one-phase Stefan problem is based on an assumption that one of the material phases may be neglected. Typically this is achieved by assuming that a phase is at the phase change temperature and hence any variation from this leads to a change of phase. This is a mathematically convenient approximation, which simplifies analysis whilst still demonstrating the essential ideas behind the process. A further standard simplification is to work in non-dimensional format, such that the temperature at the interface may be set to zero and far-field values to +1 or -1.

Consider a semi-infinite one-dimensional block of ice initially at melting temperature u ≡ 0 for x ∈ [0, +∞). The most well-known form of Stefan problem involves melting via an imposed constant temperature at the left hand boundary, leaving a region [0, s(t)] occupied by water. The melted depth, denoted by s(t), is an unknown function of time. The Stefan problem is defined by

{\displaystyle {\begin{aligned}{\frac {\partial u}{\partial t}}&={\frac {\partial ^{2}u}{\partial x^{2}}}&&{\text{in }}\{(x,t):00\},&&{\text{the heat equation}},\\[6pt]u(0,t)&=1,&&t>0,&&{\text{fixed temperature, above the melt temperature, on the left boundary}},&&\\[6pt]u{\big (}s(t),t{\big )}&=0,&&t>0,&&{\text{the interface at the melting temperature}},\\[6pt]\beta {\frac {\mathrm {d} s}{\mathrm {d} t}}&=-{\frac {\partial u}{\partial x}}{\big (}s(t),t{\big )},&&t>0,&&{\text{Stefan condition}},\\[6pt]u(x,0)&=0,&&x\geq 0,&&{\text{initial temperature distribution}},\\[6pt]s(0)&=0,&&&&{\text{initial depth of the melted ice block}}.\end{aligned}}}
where β is the Stefan number, the ratio of latent to specific sensible heat (where specific indicates it is divided by the mass). Note this definition follows naturally from the non-dimensionalisation and is used in many texts [4][5] however it may also be defined as the inverse of this (for example in the Wikipedia entry, Stefan number).
The Neumann solution, obtained by using self-similar variables, indicates that the position of the boundary is given by ${\textstyle s=2\lambda {\sqrt {t}}}$  where λ satisfies the transcendental equation
{\displaystyle {\begin{aligned}\beta \lambda ={\frac {1}{\sqrt {\pi }}}{\frac {{\text{exp}}(-\lambda ^{2})}{{\text{erf}}(\lambda )}}.\end{aligned}}}

The temperature in the liquid is then given by
${\displaystyle T=1-{\frac {{\text{erf(}}x/(2{\sqrt {t}}))}{{\text{erf}}(\lambda )}}.}$

## Applications

Apart from modelling melting of solids, Stefan problem is also used as a model for the asymptotic behaviour (in time) of more complex problems. For example, Pego[6] uses matched asymptotic expansions to prove that Cahn-Hilliard solutions for phase separation problems behave as solutions to a nonlinear Stefan problem at an intermediate time scale. Additionally, the solution of the Cahn–Hilliard equation for a binary mixture is reasonably comparable with the solution of a Stefan problem.[7] In this comparison, the Stefan problem was solved using a front-tracking, moving-mesh method with homogeneous Neumann boundary conditions at the outer boundary. Also, Stefan problems can be applied to describe phase transformations.[8]

The Stefan problem also has a rich inverse theory; in such problems, the meting depth (or curve or hypersurface) s is the known datum and the problem is to find u or f.[9]

## Advanced forms of Stefan problem

The classical Stefan problem deals with stationary materials with constant thermophysical properties (usually irrespective of phase), a constant phase change temperature and, in the example above, an instantaneous switch from the initial temperature to a distinct value at the boundary. In practice thermal properties may vary and specifically always do when the phase changes. The jump in density at phase change induces a fluid motion: the resultant kinetic energy does not figure in the standard energy balance. With an instantaneous temperature switch the initial fluid velocity is infinite, resulting in an initial infinite kinetic energy. In fact the liquid layer is often in motion, thus requiring advection or convection terms in the heat equation. The melt temperature may vary with size, curvature or speed of the interface. It is impossible to instantaneously switch temperatures and then difficult to maintain an exact fixed boundary temperature. Further, at the nanoscale the temperature may not even follow Fourier's law.

A number of these issues have been tackled in recent years for a variety of physical applications. In the solidification of supercooled melts an analysis where the phase change temperature depends on the interface velocity may be found in Font et al.[10] Nanoscale solidification, with variable phase change temperature and energy/density effects are modelled in.[11][12] Solidification with flow in a channel has been studied, in the context of lava[13] and microchannels,[14] or with a free surface in the context of water freezing over an ice layer.[15][16] A general model including different properties in each phase, variable phase change temperature and heat equations based on either Fourier's law or the Guyer-Krumhansl equation is analysed in.[17]

## Notes

1. ^ Applied partial differential equations. Ockendon, J. R. (Rev. ed.). Oxford: Oxford University Press. 2003. ISBN 0-19-852770-5. OCLC 52486357.CS1 maint: others (link)
2. ^ (Vuik 1993, p. 157).
3. ^ RUBINSTEIN, L. I. (2016). STEFAN PROBLEM. [Place of publication not identified]: American Mathematical Society. ISBN 978-1-4704-2850-1. OCLC 973324855.
4. ^ Davis, Stephen H., 1939-. Theory of solidification. Cambridge. ISBN 978-0-511-01924-1. OCLC 232161077.CS1 maint: multiple names: authors list (link)
5. ^ Fowler, A. C. (Andrew Cadle), 1953- (1997). Mathematical models in the applied sciences. Cambridge: Cambridge University Press. ISBN 0-521-46140-5. OCLC 36621805.CS1 maint: multiple names: authors list (link)
6. ^ R. L. Pego. (1989). Front Migration in the Nonlinear Cahn-Hilliard Equation. Proc. R. Soc. Lond. A.,422:261–278.
7. ^ Vermolen, F. J.; Gharasoo, M. G.; Zitha, P. L. J.; Bruining, J. (2009). "Numerical Solutions of Some Diffuse Interface Problems: The Cahn–Hilliard Equation and the Model of Thomas and Windle". International Journal for Multiscale Computational Engineering. 7 (6): 523–543. doi:10.1615/IntJMultCompEng.v7.i6.40.
8. ^ Alvarenga HD, Van de Putter T, Van Steenberge N, Sietsma J, Terryn H (Apr 2009). "Influence of Carbide Morphology and Microstructure on the Kinetics of Superficial Decarburization of C-Mn Steels". Metallurgical and Materials Transactions A. 46: 123–133. Bibcode:2015MMTA...46..123A. doi:10.1007/s11661-014-2600-y. S2CID 136871961.
9. ^ (Kirsch 1996).
10. ^ Font, F.; Mitchell, S. L.; Myers, T. G. (2013-07-01). "One-dimensional solidification of supercooled melts". International Journal of Heat and Mass Transfer. 62: 411–421. doi:10.1016/j.ijheatmasstransfer.2013.02.070. ISSN 0017-9310.
11. ^ Myers, T. G. (2016-08-01). "Mathematical modelling of phase change at the nanoscale". International Communications in Heat and Mass Transfer. 76: 59–62. doi:10.1016/j.icheatmasstransfer.2016.05.005. ISSN 0735-1933.
12. ^ Font, F.; Myers, T. G.; Mitchell, S. L. (February 2015). "A mathematical model for nanoparticle melting with density change". Microfluidics and Nanofluidics. 18 (2): 233–243. doi:10.1007/s10404-014-1423-x. ISSN 1613-4982. S2CID 54087370.
13. ^ Lister, J.R. (1994). "The solidification of buoyancy-driven flow in a flexible-walled channel. Part 1. Constant-volume release". Journal of Fluid Mechanics. 272: 21–44. Bibcode:1994JFM...272...21L. doi:10.1017/S0022112094004362.
14. ^ Myers, T. G.; Low, J. (October 2011). "An approximate mathematical model for solidification of a flowing liquid in a microchannel". Microfluidics and Nanofluidics. 11 (4): 417–428. doi:10.1007/s10404-011-0807-4. ISSN 1613-4982. S2CID 97060677.
15. ^ Myers, T. G.; Charpin, J. P. F.; Chapman, S. J. (August 2002). "The flow and solidification of a thin fluid film on an arbitrary three-dimensional surface". Physics of Fluids. 14 (8): 2788–2803. Bibcode:2002PhFl...14.2788M. doi:10.1063/1.1488599. hdl:2117/102903. ISSN 1070-6631.
16. ^ Myers, T.G.; Charpin, J.P.F. (December 2004). "A mathematical model for atmospheric ice accretion and water flow on a cold surface". International Journal of Heat and Mass Transfer. 47 (25): 5483–5500. doi:10.1016/j.ijheatmasstransfer.2004.06.037.
17. ^ Myers, T. G.; Hennessy, M. G.; Calvo-Schwarzwälder, M. (2020-03-01). "The Stefan problem with variable thermophysical properties and phase change temperature". International Journal of Heat and Mass Transfer. 149: 118975. arXiv:1904.05698. doi:10.1016/j.ijheatmasstransfer.2019.118975. ISSN 0017-9310. S2CID 115147121.