Analysis and applications of the Voronoi Implicit Interface Method

We analyze a new mathematical and numerical framework, the "Voronoi Implicit Interface Method" ("VIIM"), first introduced in Saye and Sethian (2011) R.I. Saye, J.A. Sethian, The Voronoi Implicit Interface Method for computing multiphase physics, PNAS 108 (49) (2011) 19498-19503] for tracking multiple interacting and evolving regions ("phases") whose motion is determined by complex physics (fluids, mechanics, elasticity, etc.). From a mathematical point of view, the method provides a theoretical framework for moving interface problems that involve multiple junctions, defining the motion as the formal limit of a sequence of related problems. Discretizing this theoretical framework provides a numerical methodolology which automatically handles multiple junctions, triple points and quadruple points in two dimensions, as well as triple lines, etc. in higher dimensions. Topological changes in the system occur naturally, with no surgery required. In this paper, we present the method in detail, and demonstrate several new extensions of the method to different physical phenomena, including curvature flow with surface energy densities defined on a per-phase basis, as well as multiphase fluid flow in which density, viscosity and surface tension can be defined on a per-phase basis.We test this method in a variety of ways. We perform rigorous analysis and demonstrate convergence in both two and three dimensions for a variety of evolving interface problems, including verification of von Neumann-Mullins' law in two dimensions (and its analog in three dimensions), as well as normal driven flow and curvature flow with and without constraints, demonstrating topological change and the effects of different boundary conditions. We couple the method to a second order projection method solver for incompressible fluid flow, and study the effects of membrane permeability and impermeability, large shearing torsional forces, and the effects of varying density, viscosity and surface tension on a per-phase basis. Finally, we demonstrate convergence in both space and time of a topological change in a multiphase foam.

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

[2]  Harald Garcke,et al.  A MultiPhase Field Concept: Numerical Simulations of Moving Phase Boundaries and Multiple Junctions , 1999, SIAM J. Appl. Math..

[3]  J. Sethian,et al.  The Voronoi Implicit Interface Method and Computational Challenges in Multiphase Physics , 2012 .

[4]  Baba C. Vemuri,et al.  Shape Modeling with Front Propagation: A Level Set Approach , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

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

[6]  T. Chan,et al.  A Variational Level Set Approach to Multiphase Motion , 1996 .

[7]  J. Tsitsiklis,et al.  Efficient algorithms for globally optimal trajectories , 1994, Proceedings of 1994 33rd IEEE Conference on Decision and Control.

[8]  Alexander Vladimirsky,et al.  Ordered Upwind Methods for Static Hamilton-Jacobi Equations: Theory and Algorithms , 2003, SIAM J. Numer. Anal..

[9]  P. Colella,et al.  A Conservative Adaptive Projection Method for the Variable Density Incompressible Navier-Stokes Equations , 1998 .

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

[11]  David J. Benson,et al.  Volume of fluid interface reconstruction methods for multi - material problems , 2002 .

[12]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[13]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[14]  J. Sethian,et al.  The Voronoi Implicit Interface Method for computing multiphase physics , 2011, Proceedings of the National Academy of Sciences.

[15]  P. Gennes.  The Physics Of Foams , 1999 .

[16]  Markus Bussmann,et al.  A piecewise linear approach to volume tracking a triple point , 2007 .

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

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

[19]  D. Chopp,et al.  A projection method for motion of triple junctions by level sets , 2002 .

[20]  James A. Sethian,et al.  Numerical Methods for Propagating Fronts , 1987 .

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

[22]  Selim Esedoglu,et al.  Diffusion generated motion for grain growth in two and three dimensions , 2009, J. Comput. Phys..

[23]  J. Tsitsiklis.  Efficient algorithms for globally optimal trajectories , 1995, IEEE Trans. Autom. Control..

[24]  L. Bronsard,et al.  A numerical method for tracking curve networks moving with curvature motion , 1995 .

[25]  Kenneth A. Brakke,et al.  The Surface Evolver , 1992, Exp. Math..

[26]  H. Soner,et al.  Three-phase boundary motions under constant velocities. I: The vanishing surface tension limit , 1996, Proceedings of the Royal Society of Edinburgh: Section A Mathematics.

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

[28]  P. Woodward,et al.  SLIC (Simple Line Interface Calculation) , 1976 .

[29]  David J. Srolovitz,et al.  The von Neumann relation generalized to coarsening of three-dimensional microstructures , 2007, Nature.

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

[31]  S. Osher,et al.  Motion of multiple junctions: a level set approach , 1994 .

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

[33]  A. Chorin.  Numerical solution of the Navier-Stokes equations , 1968 .

[34]  Charles S. Peskin,et al.  Numerical simulations of two-dimensional foam by the immersed boundary method , 2010, J. Comput. Phys..

[35]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[36]  J. Sethian.  Curvature and the evolution of fronts , 1985 .

[37]  Fernando Reitich,et al.  Three-phase boundary motions under constant velocities. , 1994 .