LAUSR

We introduce a Smoothed Particle Hydrodynamics (SPH) concept for the stabilization of the interface between two fluids. It is demonstrated that the change in the pressure gradient across the interface leads to a force imbalance. This force imbalance is attributed to the particle approximation implicit to SPH. To stabilize the interface a pressure gradient correction is proposed. In this approach the multi-fluid pressure gradients are related to the (gravitational and fluid) accelerations. This leads to a quasi-buoyancy correction for hydrostatic (stratified) flows, which is extended to non-hydrostatic flows. The result is a simple density correction which involves no parameters or coefficients. This correction is included as an extra term in the SPH momentum equation. The new concept for the stabilization of the interface is explored in five case studies and compared with other multi-fluid models. The first case is the stagnant flow in a tank: the interface remains stable up to density ratios of 1:1000 (typical for water and air) in combination with artificial wave speed ratios up to 1:4. The second and third cases are the Rayleigh-Taylor instability and the rising bubble, where a reasonable agreement between SPH and level-set models is achieved. The fourth case is an air flow across a water surface up to density ratios of 1:100, artificial wave speeds for water higher than that of air, and high air velocities. The fifth case is about the propagation of internal gravity waves up to density ratios of 1:100 and artificial wave speed ratios of 1:2. It is demonstrated that the quasi-buoyancy model may be used to stabilize the interface between two fluids up to high density ratios, with real (low) viscosities and more realistic wave speed ratios than achieved by other WCSPH multi-fluid models. Real wave speed ratios can be achieved, as long as the fluid velocities are not very high. Although the wave speeds may be artificial in many cases, correct and realistic wave speed ratios are essential in the modelling of heat transfer between two fluids (e.g. in engineering applications such as gas turbines).