A gradient-augmented level set method with an optimally local, coherent advection scheme

The level set approach represents surfaces implicitly, and advects them by evolving a level set function, which is numerically defined on an Eulerian grid. Here we present an approach that augments the level set function values by gradient information, and evolves both quantities in a fully coupled fashion. This maintains the coherence between function values and derivatives, while exploiting the extra information carried by the derivatives. The method is of comparable quality to WENO schemes, but with optimally local stencils (performing updates in time by using information from only a single adjacent grid cell). In addition, structures smaller than the grid size can be located and tracked, and the extra derivative information can be employed to obtain simple and accurate approximations to the curvature. We analyze the accuracy and the stability of the new scheme, and perform benchmark tests.

[1]  Takashi Yabe,et al.  Cubic-Interpolated Pseudo Particle (CIP) Method : Application to Nonlinear or Multi-Dimensional Problem , 1987 .

[2]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[3]  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..

[4]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[5]  Mehdi Raessi,et al.  Advecting normal vectors: A new method for calculating interface normals and curvatures when modeling two-phase flows , 2007, J. Comput. Phys..

[6]  Martin Oberlack,et al.  On symmetries and averaging of the G-equation for premixed combustion , 2001, Proceeding of Second Symposium on Turbulence and Shear Flow Phenomena.

[7]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[8]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[9]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[10]  M M Hafez,et al.  Innovative Methods for Numerical Solution of Partial Differential Equations , 2001 .

[11]  Takashi Yabe,et al.  3. Cubic Interpolated Pseudo-Particle Method (CIP) for Solving Hyperbolic-Type Equation (II. Basic Algorithm) , 1985 .

[12]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

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

[14]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[15]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[16]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[17]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme III. Upstream-centered finite-difference schemes for ideal compressible flow , 1977 .

[18]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[19]  R. D. Richtmyer,et al.  Survey of the stability of linear finite difference equations , 1956 .

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

[21]  J. Strain.  Tree Methods for Moving Interfaces , 1999 .

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

[23]  Frédéric Gibou,et al.  A second order accurate level set method on non-graded adaptive cartesian grids , 2007, J. Comput. Phys..

[24]  R. Fedkiw,et al.  The Ghost Fluid Method for Viscous Flows , 2010 .

[25]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[26]  H. Takewaki,et al.  The cubic-interpolated Pseudo particle (cip) method: application to nonlinear and multi-dimensional hyperbolic equations , 1987 .

[27]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[28]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[29]  R. LeVeque.  High-resolution conservative algorithms for advection in incompressible flow , 1996 .

[30]  R. Courant,et al.  On the solution of nonlinear hyperbolic differential equations by finite differences , 1952 .

[31]  E LorensenWilliam,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987 .

[32]  Eric F Darve,et al.  Author ' s personal copy A hybrid method for the parallel computation of Green ’ s functions , 2009 .

[33]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme I. The quest of monotonicity , 1973 .

[34]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .