Simulating surface tension with smoothed particle hydrodynamics

• Published in: International journal for numerical methods in fluids Vol. 33; no. 3; pp. 333 - 353
• Main Author: Morris, Joseph P.
• Format: Journal Article
• Language: English
• Published: Chichester, UK John Wiley & Sons, Ltd 15.06.2000
• Subjects: Boundary conditions; Computer simulation; Interfaces (materials); interfacial flow; Lagrange multipliers; meshless methods; Surface tension; Viscous flow
• ISSN: 0271-2091; 1097-0363
• DOI: 10.1002/1097-0363(20000615)33:3<333::AID-FLD11>3.0.CO;2-7

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.