Document Type

Article

Publication Date

11-1-2017

Abstract

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.

Publication Source (Journal or Book title)

ESAIM Mathematical Modelling and Numerical Analysis

First Page

2093

Last Page

2126

Download

Share

COinS