East Asian Journal on Applied Mathematics>2013>3>2>107-119
A Stefan problem is a free boundary problem where a phase boundary moves as a function of time. In this article, we consider one-dimensional and two-dimensional enthalpy-formulated Stefan problems. The enthalpy formulation has the advantage that the governing equations stay the same, regardless of the material state (liquid or solid). Numerical solutions are obtained by implementing the Godunov method. Our simulation of the temperature distribution and interface position for the one-dimensional Stefan problem is validated against the exact solution, and the method is then applied to the two-dimensional Stefan problem with reference to cryosurgery, where extremely cold temperatures are applied to destroy cancer cells. The temperature distribution and interface position obtained provide important information to control the cryosurgery procedure.
Financed by the National Centre for Research and Development under grant No. SP/I/1/77065/10 by the strategic scientific research and experimental development program:
SYNAT - “Interdisciplinary System for Interactive Scientific and Scientific-Technical Information”.