New Mass-Conserving Algorithm for Level Set Redistancing on Unstructured Meshes

The level set method is becoming increasingly popular for the simulation of several problems that involve interfaces. The level set function is advected by some velocity field, with the zero-level set of the function defining the position of the interface. The advection distorts the initial shape of the level set function, which needs to be re-initialized to a smooth function preserving the position of the zero-level set. Many algorithms re-initialize the level set function to (some approximation of) the signed distance from the interface. Efficient algorithms for level set redistancing on Cartesian meshes have become available over the last years, but unstructured meshes have received little attention. This presentation concerns algorithms for construction of a distance function from the zero-level set, in such a way that mass is conserved on arbitrary unstructured meshes. The algorithm is consistent with the hyperbolic character of the distance equation (∥d∥=1) and can be localized on a narrow band close to the interface, saving computing effort. The mass-correction step is weighted according to local mass differences, an improvement over usual global rebalancing techniques.

[1]  Tayfun E. Tezduyar,et al.  Moving‐interface computations with the edge‐tracked interface locator technique (ETILT) , 2005 .

[2]  Tayfun E. Tezduyar,et al.  Stabilized-finite-element/interface-capturing technique for parallel computation of unsteady flows with interfaces , 2000 .

[3]  S. Osher,et al.  A Level Set Formulation of Eulerian Interface Capturing Methods for Incompressible Fluid Flows , 1996 .

[4]  Mark Sussman,et al.  An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow , 1999, SIAM J. Sci. Comput..

[5]  Marek Behr,et al.  Enhanced-Discretization Interface-Capturing Technique (EDICT) for computation of unsteady flows with interfaces , 1998 .

[6]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[7]  David L. Chopp,et al.  Some Improvements of the Fast Marching Method , 2001, SIAM J. Sci. Comput..

[8]  R. Codina,et al.  A numerical model to track two‐fluid interfaces based on a stabilized finite element method and the level set technique , 2002 .

[9]  J. Sethian.  Advances in fast marching and level set methods for propagating interfaces , 1999 .

[10]  James A. Sethian,et al.  Fast marching methods and level set methods for propagating interfaces , 1998 .

[11]  Ian M. Mitchell,et al.  A hybrid particle level set method for improved interface capturing , 2002 .

[12]  J. Sethian,et al.  Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains , 1998 .

[13]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[14]  Djamel Lakehal,et al.  Interface tracking towards the direct simulation of heat and mass transfer in multiphase flows , 2002 .

[15]  D. Chopp.  Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[16]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .