Smoothed Particle Hydrodynamics method for simulating free surface flow

Motivation and Objectives

Smoothed particle hydrodynamics (SPH) is a purely mesh-free Lagrangian method developed for astrophysical applications. Since these pioneering works, the SPH method has been successfully applied for numerical simulations of solid mechanics, fluid dynamics and fluid-structure interaction. Concerning the computation of hydrodynamic problems, the present methods either lead to violent pressure oscillations or excessive dissipation and are not able to reproduce correct physical phenomenon reliably.

Approach to Solution

We present a low-dissipation weakly-compressible SPH method for modeling free-surface flows exhibiting violent events such as impact and breaking. The key idea is to modify a Riemann solver which determines the interaction between particles by a simple limiter to decrease the intrinsic numerical dissipation. The modified Riemann solver is also extended for imposing wall boundary conditions. Numerical tests show that the method resolves free-surface flows accurately and produces smooth, accurate pressure fields.

Three-dimensional dambreak problem simulated with dp=H/30 (the total fluid particle number N=27000): free-surface profile compared with experiment.