Simulating surface tension with smoothed particle hydrodynamics

A method for simulating two-phase flows including surface tension is presented. The approach is based upon smoothed particle hydrodynamics (SPH). The fully Lagrangian nature of SPH maintains sharp fluid–fluid interfaces without employing high-order advection schemes or explicit interface reconstruction. Several possible implementations of surface tension force are suggested and compared. The numerical stability of the method is investigated and optimal choices for numerical parameters are identified. Comparisons with a grid-based volume of fluid method for two-dimensional flows are excellent. The methods presented here apply to problems involving interfaces of arbitrary shape undergoing fragmentation and coalescence within a two-phase system and readily extend to three-dimensional problems. Boundary conditions at a solid surface, high viscosity and density ratios, and the simulation of free-surface flows are not addressed. Copyright © 2000 John Wiley & Sons, Ltd.

[1]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[2]  I. J. Schoenberg.  Contributions to the problem of approximation of equidistant data by analytic functions. Part A. On the problem of smoothing or graduation. A first class of analytic approximation formulae , 1946 .

[3]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[4]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

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

[6]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[7]  G. J. Phillips,et al.  A numerical method for three-dimensional simulations of collapsing, isothermal, magnetic gas clouds , 1985 .

[8]  L. Brookshaw,et al.  A Method of Calculating Radiative Heat Diffusion in Particle Simulations , 1985, Publications of the Astronomical Society of Australia.

[9]  David H. Sharp,et al.  Front Tracking Applied to Rayleigh–Taylor Instability , 1986 .

[10]  Gretar Tryggvason,et al.  Computations of three‐dimensional Rayleigh–Taylor instability , 1990 .

[11]  D. B. Kothe,et al.  RIPPLE: A NEW MODEL FOR INCOMPRESSIBLE FLOWS WITH FREE SURFACES , 1991 .

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

[13]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[14]  Dimos Poulikakos,et al.  Modeling of the deformation of a liquid droplet impinging upon a flat surface , 1993 .

[15]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[16]  J. Monaghan.  Simulating Free Surface Flows with SPH , 1994 .

[17]  S. Zaleski,et al.  Modelling Merging and Fragmentation in Multiphase Flows with SURFER , 1994 .

[18]  Abraham M. Lenhoff,et al.  Dynamic breakup of liquid–liquid jets , 1994 .

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

[20]  D. C. Leslie,et al.  Free surface simulations using a conservative 3D code , 1995 .

[21]  S. Attaway,et al.  Smoothed particle hydrodynamics stability analysis , 1995 .

[22]  William J. Rider,et al.  Accurate solution algorithms for incompressible multiphase flows , 1995 .

[23]  W. Rider,et al.  Stretching and tearing interface tracking methods , 1995 .

[24]  D. Balsara.  von Neumann stability analysis of smoothed particle hydrodynamics—suggestions for optimal algorithms , 1995 .

[25]  A spine-flux method for simulating free surface flows , 1995 .

[26]  Joseph Peter Morris.  A Study of the Stability Properties of SPH , 1995 .

[27]  O. Hassager,et al.  Simulation of free surfaces in 3-D with the arbitrary Lagrange-Euler method , 1995 .

[28]  L. Libersky,et al.  Smoothed Particle Hydrodynamics: Some recent improvements and applications , 1996 .

[29]  Joseph P. Morris,et al.  A Study of the Stability Properties of Smooth Particle Hydrodynamics , 1996, Publications of the Astronomical Society of Australia.

[30]  Mark A Fleming,et al.  Meshless methods: An overview and recent developments , 1996 .

[31]  D. Yue,et al.  COMPUTATION OF NONLINEAR FREE-SURFACE FLOWS , 1996 .

[32]  J. Morris,et al.  Modeling Low Reynolds Number Incompressible Flows Using SPH , 1997 .

[33]  Wei Shyy,et al.  Multiphase Dynamics in Arbitrary Geometries on Fixed Cartesian Grids , 1997 .

[34]  M. Jaeger,et al.  NUMERICAL SIMULATION OF INCOMPRESSIBLE FLOWS WITH MOVING INTERFACES , 1997 .

[35]  M. Rudman.  INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 671–691 (1997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS , 2022 .

[36]  M. Rudman.  International Journal for Numerical Methods in Fluids a Volume-tracking Method for Incompressible Multifluid Flows with Large Density Variations , 2022 .

[37]  D. B. Kothe,et al.  Accurate and robust methods for variable density incompressible flows with discontinuities , 1998 .

[38]  P. Colella,et al.  An Adaptive Level Set Approach for Incompressible Two-Phase Flows , 1997 .

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

[40]  Patrick J. Fox,et al.  A pore‐scale numerical model for flow through porous media , 1999 .

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