Density Contrast SPH Interfaces

Abstract

To simulate multiple fluids realistically many important interaction effects have to be captured accurately. Smoothed Particle Hydrodynamics (SPH) has shown to be a simple, yet flexible method to cope with many fluid simulation problems in a robust way. Unfortunately, the results obtained when using SPH to simulate miscible fluids are severely affected, especially if density ratios become large. The undesirable effects reach from unphysical density and pressure variations to spurious and unnatural interface tensions, as well as severe numerical instabilities. In this work, we present a formulation based on SPH which can handle density discontinuities at interfaces between multiple fluids correctly without increasing the computational costs compared to standard SPH. The basic idea is to replace the density computation in SPH by a measure of particle densities and consequently derive new formulations for pressure and viscous forces. The new method enables the user to select the desired amount of interface tension according to the simulation problem at hand. We succeed to stably simulate multiple fluids with high density contrasts without the above described artifacts apparent in standard SPH simulations.