A new implicit surface tension implementation for interfacial flows

A new implementation of surface tension effects in interfacial flow codes is proposed which is both fully implicit in space, that is the interface never has to be reconstructed, and also semi-implicit in time, with semi-implicit referring to the time integration of the surface tension forces. The main idea is to combine two previously separate techniques to yield a new expression for the capillary forces. The first is the continuum surface force (CSF) method, which is used to regularize the discontinuous surface tension force term. The regularization can be elegantly implemented with the use of distance functions, which makes the level set method a suitable choice for the interface-tracking algorithm. The second is to use a finite element discretization together with the Laplace–Beltrami operator, which enables simple reformulation of the surface tension term into its semi-implicit equivalent. The performance of the new method is benchmarked against standard explicit methods, where it is shown that the new method is significantly more robust for the chosen test problems when the time steps exceed the numerical capillary time step restriction. Some improvements are also found in the average number of nonlinear iterations and linear multigrid steps taken while solving the momentum equations. Copyright © 2005 John Wiley & Sons, Ltd.

[1]  P. Minev,et al.  A finite element technique for multifluid incompressible flow using Eulerian grids , 2003 .

[2]  S. Turek,et al.  FEATFLOW - Finite element software for the incompressible Navier-Stokes equations - User Manual Rele , 1998 .

[3]  Stefan Turek,et al.  THE EIKONAL EQUATION: NUMERICAL EFFICIENCY VS. ALGORITHMIC COMPLEXITY ON QUADRILATERAL GRIDS , 2005 .

[4]  D. Kuzmin,et al.  High-resolution FEM-TVD schemes based on a fully multidimensional flux limiter , 2004 .

[5]  G. Dziuk,et al.  An algorithm for evolutionary surfaces , 1990 .

[6]  Eberhard Bänsch,et al.  Finite element discretization of the Navier–Stokes equations with a free capillary surface , 2001, Numerische Mathematik.

[7]  E. Bänsch,et al.  Uniaxial extensional flows in liquid bridges , 2004, Journal of Fluid Mechanics.

[8]  C. Peskin Numerical analysis of blood flow in the heart , 1977 .

[9]  J. Hochstein,et al.  An implicit surface tension model , 1996 .

[10]  L. Tobiska,et al.  FINITE ELEMENT SIMULATION OF A DROPLET IMPINGING A HORIZONTAL SURFACE , 2005 .

[11]  A. Reusken,et al.  A finite element based level set method for two-phase incompressible flows , 2006 .

[12]  I. Holopainen Riemannian Geometry , 1927, Nature.

[13]  J. Z. Zhu,et al.  The superconvergent patch recovery and a posteriori error estimates. Part 1: The recovery technique , 1992 .

[14]  A. Smolianski Numerical Modeling of Two-Fluid Interfacial Flows , 2001 .

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

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

[17]  R. Clift,et al.  Bubbles, Drops, and Particles , 1978 .