A simple SPH algorithm for multi‐fluid flow with high density ratios

SUMMARY In this paper, we describe an SPH algorithm for multi-fluid flow, which is efficient, simple and robust. We derive the inviscid equations of motion from a Lagrangian together with the constraint provided by the continuity equation. The viscous flow equations then follow by adding a viscous term. Rigid boundaries are simulated using boundary force particles in a manner similar to the immersed boundary method. Each fluid is approximated as weakly compressible with a speed of sound sufficiently large to guarantee that the relative density variations are typically 1%. When the SPH force interaction is between two particles of different fluids, we increase the pressure terms. This simple procedure stabilizes the interface between the fluids. The equations of motion are integrated using a time stepping rule based on a second-order symplectic integrator. When linear and angular momentum should be conserved exactly, they are conserved to within round-off errors. We test the algorithm by simulating a variety of problems involving fluids with a density ratio in the range 1–1000. The first of these is a free surface problem with no rigid boundaries. It involves the flow of an elliptical distribution with one fluid inside the other. We show that the simulations converge as the particle spacing decreases, and the results are in good agreement with the exact inviscid, incompressible theory. The second test is similar to the first but involves the nonlinear oscillation of the fluids. As in the first test, the agreement with theory is very good, and the method converges. The third test is the simulation of waves at the interface between two fluids. The method is shown to converge, and the agreement with theory is satisfactory. The fourth test is the Rayleigh–Taylor instability for a configuration considered by other authors. Key parameters are shown to converge, and the agreement with other authors is good. The fifth and final test is how well the SPH method simulates gravity currents with density ratios in the range 2–30. The results of these simulations are in very good agreement with those of other authors and in satisfactory agreement with experimental results.Copyright © 2012 John Wiley & Sons, Ltd.

[1]  A. Colagrossi,et al.  Numerical simulation of interfacial flows by smoothed particle hydrodynamics , 2003 .

[2]  Pierre Saramito,et al.  Numerical simulations of high density ratio lock-exchange flows , 2005 .

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

[4]  D. Layzer,et al.  On the Instability of Superposed Fluids in a Gravitational Field. , 1955 .

[5]  Ryan J. Lowe,et al.  The non-Boussinesq lock-exchange problem. Part 1. Theory and experiments , 2005, Journal of Fluid Mechanics.

[6]  J. Ross Macdonald,et al.  Some Simple Isothermal Equations of State , 1966 .

[7]  Paul W. Cleary,et al.  Modelling confined multi-material heat and mass flows using SPH , 1998 .

[8]  J. Monaghan.  SPH without a Tensile Instability , 2000 .

[9]  Lawrence Sirovich,et al.  Initial and Boundary Value Problems in Dissipative Gas Dynamics , 1967 .

[10]  Nikolaus A. Adams,et al.  An incompressible multi-phase SPH method , 2007, J. Comput. Phys..

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

[12]  Carl Eckart,et al.  Variation Principles of Hydrodynamics , 1960 .

[13]  Rex Britter,et al.  The propagation of intrusion fronts of high density ratios , 1993, Journal of Fluid Mechanics.

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

[15]  Lawrence Sirovich,et al.  Steady Gasdynamic Flows , 1968 .

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

[17]  Eckart Meiburg,et al.  The non-Boussinesq lock-exchange problem. Part 2. High-resolution simulations , 2005, Journal of Fluid Mechanics.

[18]  Ding Xin,et al.  On criterions for smoothed particle hydrodynamics kernels in stable field , 2005 .

[19]  David Le Touzé,et al.  An Hamiltonian interface SPH formulation for multi-fluid and free surface flows , 2009, J. Comput. Phys..

[20]  S. Cummins,et al.  An SPH Projection Method , 1999 .

[21]  A T J Hayward,et al.  Compressibility equations for liquids: a comparative study , 1967 .

[22]  P. Roache.  QUANTIFICATION OF UNCERTAINTY IN COMPUTATIONAL FLUID DYNAMICS , 1997 .

[23]  Nikolaus A. Adams,et al.  A multi-phase SPH method for macroscopic and mesoscopic flows , 2006, J. Comput. Phys..

[24]  Joe J. Monaghan,et al.  SPH particle boundary forces for arbitrary boundaries , 2009, Comput. Phys. Commun..

[25]  Holger Wendland,et al.  Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree , 1995, Adv. Comput. Math..

[26]  P. Cleary,et al.  Conduction Modelling Using Smoothed Particle Hydrodynamics , 1999 .