On the numerical solution of two-phase Stefan problems with heat-flux boundary conditions

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:264

Article

Mitchell, S. L.¹  Vynnycky, M.¹ 

[1] Univ Limerick, Dept Math & Stat, Math Applicat Consortium Sci & Ind MACSI, Limerick, Ireland

关键词: Stefan problem;    Keller box scheme;    Boundary immobilization;    Starting solutions;    Two-phase;   

DOI  :  10.1016/j.cam.2014.01.003

来源: Elsevier

【 摘 要 】

A recently derived numerical algorithm for one-dimensional one-phase Stefan problems is extended for the purpose of two-phase moving boundary problems in which the second phase first appears only after a finite delay time; this can occur if the phase change is caused by a heat-flux boundary condition. In tandem with the Keller box finite-difference scheme, the so-called boundary immobilization method is used. An important component of the work is the use of variable transformations that must be built into the numerical algorithm to resolve the boundary-condition discontinuity that is associated with the onset of phase change. This allows the delay time until solidification begins to be determined, and gives second-order accuracy in both time and space. (C) 2014 Elsevier B.V. All rights reserved.

【 授权许可】

Free