A Locally Gradient-Preserving Reinitialization for Level Set Functions

The level set method commonly requires a reinitialization of the level set function due to interface motion and deformation. We extend the traditional technique for reinitializing the level set function to a method that preserves the interface gradient. The gradient of the level set function represents the stretching of the interface, which is of critical importance in many physical applications. The proposed locally gradient-preserving reinitialization (LGPR) method involves the solution of three PDEs of Hamilton–Jacobi type in succession; first the signed distance function is found using a traditional reinitialization technique, then the interface gradient is extended into the domain by a transport equation, and finally the new level set function is found by solving a generalized reinitialization equation. We prove the well-posedness of the Hamilton–Jacobi equations, with possibly discontinuous Hamiltonians, and propose numerical schemes for their solutions. A subcell resolution technique is used in the numerical solution of the transport equation to extend data away from the interface directly with high accuracy. The reinitialization technique is computationally inexpensive if the PDEs are solved only in a small band surrounding the interface. As an important application, we show how the LGPR procedure can be used to make possible the local level set approach to the Eulerian Immersed boundary method.

[1]  H. Ishii,et al.  Existence and uniqueness of solutions of Hamilton-Jacobi equations , 1986 .

[2]  Shigeaki Koike,et al.  A Beginner's Guide to the Theory of Viscosity Solutions , 2014 .

[3]  Degenerate Eikonal equations with discontinuous refraction index , 2006 .

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

[5]  A. Harten.  ENO schemes with subcell resolution , 1989 .

[6]  H. Ishii.  A simple, direct proof of uniqueness for solutions of the hamilton-jacobi equations of eikonal type , 1987 .

[7]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[8]  P. Lions,et al.  Two approximations of solutions of Hamilton-Jacobi equations , 1984 .

[9]  Joseph Teran,et al.  Peristaltic pumping and irreversibility of a Stokesian viscoelastic fluid , 2008 .

[10]  J. Sethian,et al.  A Fast Level Set Method for Propagating Interfaces , 1995 .

[11]  Stanley Osher,et al.  Fast Sweeping Algorithms for a Class of Hamilton-Jacobi Equations , 2003, SIAM J. Numer. Anal..

[12]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid , 2012 .

[13]  Wanda Strychalski,et al.  A poroelastic immersed boundary method with applications to cell biology , 2015, J. Comput. Phys..

[14]  R Malladi,et al.  Image processing via level set curvature flow. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[15]  G. Batchelor,et al.  An Introduction to Fluid Dynamics , 1968 .

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

[17]  Chohong Min,et al.  On reinitializing level set functions , 2010, J. Comput. Phys..

[18]  G. Cottet,et al.  EULERIAN FORMULATION AND LEVEL SET MODELS FOR INCOMPRESSIBLE FLUID-STRUCTURE INTERACTION , 2008 .

[19]  Hongkai Zhao,et al.  A fast sweeping method for Eikonal equations , 2004, Math. Comput..

[20]  R. Guy,et al.  Computational Challenges for Simulating Strongly Elastic Flows in Biology , 2015 .

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

[22]  G. Cottet,et al.  A level-set formulation of immersed boundary methods for fluid–structure interaction problems , 2004 .

[23]  Lisa Fauci,et al.  An actuated elastic sheet interacting with passive and active structures in a viscoelastic fluid , 2013 .

[24]  Hwan Pyo Moon,et al.  MATHEMATICAL THEORY OF MEDIAL AXIS TRANSFORM , 1997 .

[25]  M. Lai,et al.  An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity , 2000 .

[26]  C. Peskin,et al.  Implicit second-order immersed boundary methods with boundary mass , 2008 .

[27]  P. Lions.  Generalized Solutions of Hamilton-Jacobi Equations , 1982 .

[28]  C. M. Elliott,et al.  Uniqueness and error analysis for Hamilton-Jacobi equations with discontinuities , 2004 .

[29]  Stanley Osher,et al.  Computing multi-valued physical observables for the high frequency limit of symmetric hyperbolic systems , 2005 .

[30]  Huajian Gao,et al.  A Numerical Study of Electro-migration Voiding by Evolving Level Set Functions on a Fixed Cartesian Grid , 1999 .

[31]  Ricardo Cortez,et al.  Shape oscillations of a droplet in an Oldroyd-B fluid , 2011 .

[32]  Gilles Aubert,et al.  Signed distance functions and viscosity solutions of discontinuous Hamilton-Jacobi Equations , 2002 .

[33]  Soravia,et al.  Boundary value problems for Hamilton-Jacobi equations with discontinuous lagrangian , 2002 .

[34]  N. Abbott,et al.  Straining soft colloids in aqueous nematic liquid crystals , 2016, Proceedings of the National Academy of Sciences.

[35]  D. Bottino,et al.  Modeling Viscoelastic Networks and Cell Deformation in the Context of the Immersed Boundary Method , 1998 .

[36]  G. Cottet,et al.  A LEVEL SET METHOD FOR FLUID-STRUCTURE INTERACTIONS WITH IMMERSED SURFACES , 2006 .

[37]  André Lieutier,et al.  Any open bounded subset of Rn has the same homotopy type as its medial axis , 2004, Comput. Aided Des..

[38]  Robert Michael Kirby,et al.  Unconditionally stable discretizations of the immersed boundary equations , 2007, J. Comput. Phys..

[39]  Becca Thomases,et al.  Mechanisms of elastic enhancement and hindrance for finite-length undulatory swimmers in viscoelastic fluids. , 2014, Physical review letters.

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

[41]  S. Osher,et al.  A PDE-Based Fast Local Level Set Method 1 , 1998 .

[42]  David Salac,et al.  A level set projection model of lipid vesicles in general flows , 2011, J. Comput. Phys..

[43]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[44]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[45]  H. Ishii.  Hamilton-Jacobi Equations with Discontinuous Hamiltonians on Arbitrary Open Sets , 1985 .

[46]  Maurizio Falcone,et al.  An Approximation Scheme for an Eikonal Equation with Discontinuous Coefficient , 2013, SIAM J. Numer. Anal..

[47]  R. Newcomb.  VISCOSITY SOLUTIONS OF HAMILTON-JACOBI EQUATIONS , 2010 .

[48]  Pierpaolo Soravia.  Optimal control with discontinuous running cost: Eikonal equation and shape-from-shading , 2000, Proceedings of the 39th IEEE Conference on Decision and Control (Cat. No.00CH37187).

[49]  Joseph Teran,et al.  Viscoelastic fluid response can increase the speed and efficiency of a free swimmer. , 2010, Physical review letters.

[50]  Jessika Eichel,et al.  Partial Differential Equations Second Edition , 2016 .

[51]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .

[52]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .

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

[54]  Daniel N. Ostrov.  Extending viscosity solutions to Eikonal equations with discontinuous spatial dependence , 2000 .

[55]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.