Advecting normal vectors: A new method for calculating interface normals and curvatures when modeling two-phase flows

In simulating two-phase flows, interface normal vectors and curvatures are needed for modeling surface tension. In the traditional approach, these quantities are calculated from the spatial derivatives of a scalar function (e.g. the volume-of-fluid or the level set function) at any instant in time. The orders of accuracy of normals and curvatures calculated from these functions are studied. A new method for calculating these quantities is then presented, where the interface unit normals are advected along with whatever function represents the interface, and curvatures are calculated directly from these advected normals. To illustrate this new approach, the volume-of-fluid method is used to represent the interface and the advected normals are used for interface reconstruction. The accuracy and performance of the new method are demonstrated via test cases with prescribed velocity fields. The results are compared with those of traditional approaches.

[1]  Feng Xiao,et al.  Description of complex and sharp interface with fixed grids in incompressible and compressible fluid , 1995 .

[2]  J. Lowengrub,et al.  Evolving interfaces via gradients of geometry-dependent interior Poisson problems: application to tumor growth , 2005 .

[3]  J. López,et al.  On the reinitialization procedure in a narrow‐band locally refined level set method for interfacial flows , 2005 .

[4]  J. Sethian.  Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

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

[6]  Peter Smereka,et al.  Semi-Implicit Level Set Methods for Curvature and Surface Diffusion Motion , 2003, J. Sci. Comput..

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

[8]  Takashi Yabe,et al.  Cubic interpolated pseudo-particle method (CIP) for solving hyperbolic-type equations , 1985 .

[9]  J. A. Sethian,et al.  Fast Marching Methods , 1999, SIAM Rev..

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

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

[12]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

[13]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[14]  E. Shirani,et al.  Interface pressure calculation based on conservation of momentum for front capturing methods , 2005 .

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

[16]  D. J. Torres,et al.  On the theory and computation of surface tension: the elimination of parasitic currents through energy conservation in the second-gradient method , 2002 .

[17]  Takashi Yabe,et al.  Constructing exactly conservative scheme in a non-conservative form , 2000 .

[18]  M. Sussman.  A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles , 2003 .

[19]  P. Colella,et al.  Non-convex profile evolution in two dimensions using volume of fluids , 1997 .

[20]  M. Renardy,et al.  PROST: a parabolic reconstruction of surface tension for the volume-of-fluid method , 2002 .

[21]  Falai Chen,et al.  Fast data extrapolating , 2007 .

[22]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .

[23]  Suthee Wiri,et al.  On improving mass conservation of level set by reducing spatial discretization errors , 2005 .

[24]  S. Osher,et al.  Spatially adaptive techniques for level set methods and incompressible flow , 2006 .

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

[26]  S. Osher,et al.  A Simple Level Set Method for Solving Stefan Problems , 1997, Journal of Computational Physics.

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

[28]  L YoungsD,et al.  Time-dependent multi-material flow with large fluid distortion. , 1982 .

[29]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[30]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .

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

[32]  S. Cummins,et al.  Estimating curvature from volume fractions , 2005 .

[33]  Said I. Abdel-Khalik,et al.  Accurate representation of surface tension using the level contour reconstruction method , 2005 .

[34]  Matthew W. Williams,et al.  A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework , 2006, J. Comput. Phys..

[35]  Nasser Ashgriz,et al.  A computational method for determining curvatures , 1989 .

[36]  Nahmkeon Hur,et al.  A COUPLED LEVEL SET AND VOLUME-OF-FLUID METHOD FOR THE BUOYANCY-DRIVEN MOTION OF FLUID PARTICLES , 2002 .