A front-tracking immersed-boundary framework for simulating Lagrangian melting problems

TL;DR

This paper introduces a novel front-tracking immersed-boundary numerical framework for simulating Lagrangian melting problems, incorporating a dynamic remeshing technique to accurately track volume changes during phase transition.

Contribution

It extends the moving least squares immersed-boundary method to include melting by implementing the Stefan condition and introduces a scalable, volume-conserving remeshing procedure for Lagrangian solids.

Findings

Successfully simulates melting with accurate volume tracking.

Develops a dynamic remeshing method with negligible computational cost.

Ensures volume conservation to machine precision during remeshing.

Abstract

In so-called Lagrangian melting problems, a solid immersed in a fluid medium is free to rotate and translate in tandem with its phase-change from solid to liquid. Such configurations may be classified as a fluid-solid interaction (FSI) problem coupled to phase-change. Our present work proposes a numerical method capable of simulating these Lagrangian melting problems and adopts a front-tracking immersed-boundary (IB) method. We use the moving least squares IB framework, a well-established method for simulating a diverse range of FSI problems and extend this framework to accommodate melting by additionally imposing the Stefan condition at the interface. In the spirit of canonical front-tracking methods, the immersed solid is represented by a discrete triangulated mesh which is separate from the Eulerian mesh in which the governing flow equations are solved. A known requirement for these methods is the need for comparable Eulerian and Lagrangian grid spacings to stabilise interpolation and spreading operations between the two grids. For a melting object, this requirement is inevitably violated unless interventional remeshing is introduced. Our work therefore presents a novel dynamic remeshing procedure to overcome this. The remeshing is based on a gradual coarsening of the triangulated Lagrangian mesh and amounts to a negligible computational burden per timestep owing to the incremental and local nature of its operations, making it a scalable approach. Moreover, the coarsening is coupled to a volume-conserving smoothing procedure detailed by Kuprat et al. (2001), ensuring a zero net volume change in the remeshing step to machine precision. This added feature makes our present method highly specialised to the study of melting problems, where precise measurements of the melting solid's volume is often the primary predictive quantity of interest.