On the numerical solution of a Stefan problem with finite extinction time
2015 (English)In: Journal of Computational and Applied Mathematics, ISSN 0377-0427, Vol. 276, 98-109 p.Article in journal (Refereed) Published
In many phase-change problems of practical interest, it is important to know when a phase is depleted, a quantity referred to as the extinction time; however, there are no numerical schemes that are able to compute this with any degree of rigour or formal accuracy. In this paper, we develop such a scheme for the one-dimensional time-dependent problem of an evaporating spherical droplet. The Keller box finite-difference scheme is used, in tandem with the so-called boundary immobilization method. An important component of the work is the careful use of variable transformations that must be built into the numerical algorithm in order to preserve second-order accuracy in both time and space, in particular as regards resolving a square-root singularity in the droplet radius as the extinction time is approached.
Place, publisher, year, edition, pages
2015. Vol. 276, 98-109 p.
Evaporation, Stefan problem, Keller box scheme, Extinction time
IdentifiersURN: urn:nbn:se:kth:diva-158370DOI: 10.1016/j.cam.2014.08.023ISI: 000345478900006ScopusID: 2-s2.0-84907494982OAI: oai:DiVA.org:kth-158370DiVA: diva2:780982
QC 201501152015-01-152015-01-072015-01-15Bibliographically approved