A conservative phase field method for solving incompressible two-phase flows

In this paper a conservative phase-field method based on the work of Sun and Beckermann [Y. Sun, C. Beckermann, Sharp interface tracking using the phase-field equation, J. Comput. Phys. 220 (2007) 626-653] for solving the two- and three-dimensional two-phase incompressible Navier-Stokes equations is proposed. The present method can preserve the total mass as the Cahn-Hilliard equation, but the calculation and implementation are much simpler than that. The dispersion-relation-preserving schemes are utilized for the advection terms while the Helmholtz smoother is applied to compute the surface-tension force term. To verify the proposed method, several benchmarks are examined and shown to have good agreements with previous results. It also shows that the satisfactions of mass conservations are guaranteed.

[1]  J. C. Martin.  An experimental study of the collapse of liquid column on a rigid horizontal plane , 1952 .

[2]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[3]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

[4]  F. Xiao,et al.  Numerical simulations of free-interface fluids by a multi-integrated moment method , 2005 .

[5]  Samuel M. Allen,et al.  Mechanisms of phase transformations within the miscibility gap of Fe-rich Fe-Al alloys , 1976 .

[6]  Peter Smereka,et al.  Axisymmetric free boundary problems , 1997, Journal of Fluid Mechanics.

[7]  Chang Shu,et al.  Diffuse interface model for incompressible two-phase flows with large density ratios , 2007, J. Comput. Phys..

[8]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

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

[10]  R. Folch,et al.  Phase-field model for Hele-Shaw flows with arbitrary viscosity contrast. I. Theoretical approach. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[11]  Martin Rieber,et al.  A numerical study on the mechanism of splashing , 1999 .

[12]  C. K. Thornhill,et al.  Part IV. An experimental study of the collapse of liquid columns on a rigid horizontal plane , 1952, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[13]  Jean Leray,et al.  Sur le mouvement d'un liquide visqueux emplissant l'espace , 1934 .

[14]  Ying Sun,et al.  Sharp interface tracking using the phase-field equation , 2007, J. Comput. Phys..

[15]  Daisuke Furihata,et al.  A stable and conservative finite difference scheme for the Cahn-Hilliard equation , 2001, Numerische Mathematik.

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

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

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

[19]  D. B. Kothe,et al.  Accuracy and Convergence of Continuum Surface Tension Models , 1998 .

[20]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy and Free Energy of a Nonuniform System. III. Nucleation in a Two‐Component Incompressible Fluid , 2013 .

[21]  G. Kreiss,et al.  A conservative level set method for two phase flow II , 2005, Journal of Computational Physics.

[22]  Kensuke Yokoi,et al.  Efficient implementation of THINC scheme: A simple and practical smoothed VOF algorithm , 2007, J. Comput. Phys..

[23]  Junseok Kim,et al.  A continuous surface tension force formulation for diffuse-interface models , 2005 .

[24]  P. Chu,et al.  A Three-Point Combined Compact Difference Scheme , 1998 .

[25]  Ranga Narayanan,et al.  Fluid Dynamics at Interfaces , 2010 .

[26]  Héctor D. Ceniceros,et al.  Computation of multiphase systems with phase field models , 2002 .

[27]  John R. Blake,et al.  Gas bubbles bursting at a free surface , 1993, Journal of Fluid Mechanics.

[28]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[29]  Kensuke Yokoi,et al.  A Numerical Method for Free-Surface Flows and Its Application to Droplet Impact on a Thin Liquid Layer , 2008, J. Sci. Comput..

[30]  Christopher R. Anderson,et al.  A VORTEX METHOD FOR FLOWS WITH SLIGHT DENSITY VARIATIONS , 1985 .

[31]  D. M. Anderson,et al.  DIFFUSE-INTERFACE METHODS IN FLUID MECHANICS , 1997 .

[32]  W. Tao,et al.  A coupled volume-of-fluid and level set ( VOSET ) method for computing incompressible two-phase flows , 2009 .

[33]  L. Quartapelle,et al.  A projection FEM for variable density incompressible flows , 2000 .

[34]  Francis H. Harlow,et al.  The Splash of a Liquid Drop , 1967 .

[35]  A. Karma,et al.  Phase-Field Simulation of Solidification , 2002 .

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

[37]  Tony W. H. Sheu,et al.  Development of a Dispersion Relation-Preserving Upwinding Scheme for Incompressible Navier–Stokes Equations on NonStaggered Grids , 2005 .

[38]  Gene H. Golub,et al.  A Fast Poisson Solver for the Finite Difference Solution of the Incompressible Navier-Stokes Equations , 1998, SIAM J. Sci. Comput..

[39]  Feng Xiao,et al.  A simple algebraic interface capturing scheme using hyperbolic tangent function , 2005 .

[40]  R. Folch,et al.  Phase-field model for Hele-Shaw flows with arbitrary viscosity contrast. II. Numerical study. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[41]  Tony W. H. Sheu,et al.  On the development of a dispersion-relation-preserving dual-compact upwind scheme for convection-diffusion equation , 2009, J. Comput. Phys..