Volume 51, Number 6, November-December 2017
|Page(s)||2093 - 2126|
|Published online||27 November 2017|
Semi-discrete error estimates and implementation of a mixed method for the Stefan problem∗
1 Department of Mathematics, Tennessee Tech University, 1 William L Jones Dr, Cookeville, TN 38505, USA.
2 Department of Mathematics, Louisiana State University, Baton Rouge, LA 70803-4918, USA.
Received: 29 July 2016
Revised: 19 January 2017
Accepted: 20 April 2017
We analyze a dual formulation and finite element method for simulating the Stefan problem with surface tension (originally presented in [C.B. Davis and S.W. Walker, Int. Free Bound. 17 (2015) 427–464]). The method uses a mixed form of the heat equation in the solid and liquid (bulk) domains, and imposes a weak formulation of the interface motion law (on the solid-liquid interface) as a constraint. The computational method uses a conforming mesh approach to accurately capture the jump conditions across the interface. Preliminary error estimates are derived, under reduced regularity assumptions, for the difference between the time semi-discrete solution and the fully discrete solution over one time step. Moreover, details of the implementation are discussed including mesh generation issues. Several simulations of interface growth (in two dimensions) are presented to illustrate the method.
Mathematics Subject Classification: 65M60 / 35K20
Key words: Stefan problem / mixed method / energy stability / error estimate; interface motion / semi-implicit scheme / re-meshing / conforming mesh
© EDP Sciences, SMAI 2017
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.