An accurate finite-difference method for ablation-type Stefan problems

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:236

Article

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

[1] Univ Limerick, Dept Math & Stat, MACSI, Limerick, Ireland

关键词: Ablation;    Stefan problem;    Keller box scheme;    Boundary immobilization;    Starting solutions;   

DOI  :  10.1016/j.cam.2012.05.011

来源: Elsevier

【 摘 要 】

A recently derived numerical algorithm for one-dimensional time-dependent Stefan problems is extended for the purpose of solving one-phase ablation-type moving boundary problems; 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 in order to preserve second-order accuracy in both time and space. The analysis also determines that the ablation front initially moves as the time raised to the power 3/2; hence, it evolves considerably more slowly than the phase-change front in the classical Stefan problem with isothermal cooling. (C) 2012 Elsevier B.V. All rights reserved.

【 授权许可】

Free