Finite-difference methods with increased accuracy and correct initialization for one-dimensional Stefan problems. (English) Zbl 1177.80078

Summary: Although the numerical solution of one-dimensional phase-change, or Stefan, problems is well documented, a review of the most recent literature indicates that there are still unresolved issues regarding the start-up of a computation for a region that initially has zero thickness, as well as how to determine the location of the moving boundary thereafter. This paper considers the so-called boundary immobilization method for four benchmark melting problems, in tandem with three finite-difference discretization schemes. We demonstrate a combined analytical and numerical approach that eliminates completely the ad hoc treatment of the starting solution that is often used, and is numerically second-order accurate in both time and space, a point that has been consistently overlooked for this type of moving-boundary problem.


80A22 Stefan problems, phase changes, etc.
80M20 Finite difference methods applied to problems in thermodynamics and heat transfer
Full Text: DOI


[1] Caldwell, J.; Kwan, Y. Y., Numerical methods for one-dimensional Stefan problems, Commun. Numer. Meth. Eng., 20, 535-545 (2004) · Zbl 1048.65095
[2] Esen, A.; Kutluay, S., A numerical solution of the Stefan problem with a Neumann-type boundary condition by enthalpy method, App. Math. Comput., 148, 321-329 (2004) · Zbl 1034.65070
[3] Caldwell, J.; Chan, C.-C., Spherical solidification by the enthalpy method and the heat balance integral method, Appl. Math. Modell., 24, 45-53 (2000) · Zbl 0944.80001
[4] Caldwell, J.; Savović, S., Numerical solution of Stefan problem by variable space grid and boundary immobilization method, J. Math. Sci., 13, 67-79 (2002)
[5] Caldwell, J.; Kwan, Y. Y., Starting solutions for the boundary immobilization method, Commun. Numer. Meth. Eng., 21, 289-295 (2005) · Zbl 1112.80011
[6] Savović, S.; Caldwell, J., Finite-difference solution of one-dimensional Stefan problem with periodic boundary conditions, Int. J. Heat Mass Trans., 46, 2911-2916 (2003) · Zbl 1041.80004
[7] Kutluay, S.; Bahadir, A. R.; Ozdes, A., The numerical solution of one-phase classical Stefan problem, J. Comput. Appl. Math., 81, 135-144 (1997) · Zbl 0885.65102
[8] Meek, P. C.; Norbury, J., Nonlinear moving boundary problems and a Keller box scheme, SIAM J. Numer. Anal., 21, 5, 883-893 (1984) · Zbl 0558.65087
[9] Liu, F.; McElwain, D. L.S., A computationally efficient solution technique for moving-boundary problems in finite media, IMA J. Appl. Math., 59, 71-84 (1997) · Zbl 0958.80007
[10] Rizwan-uddin, A nodal method for phase change moving boundary problems, Int. J. Comput. Fluid Dyn., 11, 211-221 (1999) · Zbl 0962.76074
[11] Caldwell, J.; Chiu, C. K., Numerical solution of one-phase Stefan problems by the heat balance integral method, Part II - special small time starting procedure, Commun. Numer. Meth. Eng., 16, 585-593 (2000) · Zbl 0964.65110
[12] Caldwell, J.; Chiu, C. K., Numerical solution of one-phase Stefan problems by the heat balance integral method, Part I - cylindrical and spherical geometries, Commun. Numer. Meth. Eng., 16, 569-583 (2000) · Zbl 0964.65109
[13] Goodman, T. R., The heat-halance integral and its application to problems involving a change of phase, Trans. ASME, 80, 335-342 (1958)
[14] Mitchell, S. L.; Myers, T. G., A heat balance integral method for one-dimensional finite ablation, AIAA J. Thermophys., 22, 3, 508-514 (2008)
[15] Mitchell, S. L.; Myers, T. G., Approximate solution methods for one-dimensional solidification from an incoming fluid, Appl. Math. Comput., 202, 1, 311-5317 (2008) · Zbl 1208.76142
[16] Myers, T. G.; Mitchell, S. L.; Muchatibaya, G.; Myers, M. Y., A cubic heat balance integral method for one-dimensional melting of a finite thickness layer, Int. J. Heat Mass Trans., 50, 5305-5317 (2007) · Zbl 1140.80389
[17] Myers, T. G., Optimizing the exponent in the heat balance and refined integral methods, Int. Commun. Heat Mass Trans., 36, 2, 143-147 (2009)
[20] Rizwan-uddin, One-dimensional phase change with periodic boundary conditions, Numer. Heat Trans. A, 35, 361-372 (1999) · Zbl 0962.76074
[21] Ascher, U. M.; McLachlan, R. I., Multisymplectic box schemes and the Korteweg-de Vries equation, Appl. Numer. Math., 48, 255-269 (2004) · Zbl 1038.65138
[22] Mitchell, S. L.; Morton, K. W.; Spence, A., Analysis of box schemes for reactive flow problems, SIAM J. Sci. Comput., 27, 4, 1202-1223 (2006) · Zbl 1136.65332
[23] Strikwerda, J. C., Finite Difference Schemes and Partial Differential Equations (2004), Society for Industrial Mathematics · Zbl 1071.65118
[24] Schwerdtfeger, K.; Sato, M.; Tacke, K.-H., Stress formation in solidifying bodies. Solidification in a round continuous casting mold, Metall. Mater. Trans. B, 29B, 1057-1068 (1998)
[25] Vynnycky, M., An asymptotic model for the formation and evolution of air gaps in vertical continuous casting, Proc. Roy. Soc. A, 465, 1617-1644 (2009) · Zbl 1186.74037
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.