Abstract:As a flexible Lagrangian particle method, smoothed particle hydrodynamics (SPH) can easily capture large interface deformation, breaking, merging and splashing, which is quite suitable for the simulation of water surface waves. The classic weakly compressible smoothed particle hydrodynamics (WCSPH) algorithm for incompressible fluid flow usually needs small time steps to ensure numerical stability, which makes it computational time consuming. The major reason is that the pressure is linked with the density deviation which should be zero everywhere in incompressible fluid flows. Generally, the equation of state for a weakly compressible fluid is adopted to simulate incompressible fluid. To enforce the incompressibility, the stiffness of fluid is set to be large enough to keep the density deviation less than an allowable error. That will lead to a problem that a tiny error of computed density will cause a large pressure error and further numerical instability. In order to obtain rational pressure distribution and good computational efficiency as well, a new algorithm is proposed considering that flow field is globally incompressible from the view of spatial averaged flow field and weakly compressible locally due to the computational error from the view of particles. The basic conception of the proposed modified SPH is to split the pressure into a global part and local fluctuating part. The global part of pressure or the spatial averaged pressure on the resolution of background mesh is obtained by solving the pressure Poisson equation based on the velocity divergence free condition. The local fluctuating pressure is related to the local particle density variation by solving the local pressure Poisson equation, which is applied to maintain particles’ spacing in relative equilibrium. Based on this idea, formulas for the spatial averaged pressure and local fluctuating pressure are derived. The spatial averaged pressure is used as a primary estimation of the pressure gradient force and then the local fluctuating pressure is added as a correction of local density deviation. The proposed algorithm is verified by simulating the classic dam broken problem and then applied to simulate water wave breaking. A comparison of the particle configurations with previous studies at different times is shown. Good agreement between the present results and the previous results can be found. A discussion on computational efficiency is also carried out. The result shows that the proposed algorithm has better stability and can greatly reduce the computational time cost.

A Corrected Solid Boundary Treatment Method for Smoothed Particle Hydrodynamics

On methods of processing solid boundary condition in SPH simulation

A Method to Improve First Order Approximation of Smoothed Particle Hydrodynamics

On the Treatment of Solid Boundary in Smoothed Particle Hydrodynamics

A New Solid Boundary Treatment Algorithm for Smoothed Particle Hydrodynamics

An improved SPH method for modeling liquid sloshing dynamics

Smoothed Particle Hydrodynamics Modeling of Free Surface Flow

Simulations of Droplet Spreading and Solidification Using an Improved SPH Model

A moving least square immersed boundary method for SPH with thin-walled rigid structures

Pressure corrected SPH for fluid animation

A moving least square immersed boundary method for SPH with thin-walled structures

A Modified SPH Model for Simulating Water Surface Waves

Robust Solid Boundary Treatment for Compressible Smoothed Particle Hydrodynamics

Improved Smoothed Particle Hydrodynamics with Rans for Free-Surface Flow Problems

Semi‐analytical Solid Boundary Conditions for Free Surface Flows

Particle Distribution and Numerical Stability in Smoothed Particle Hydrodynamics Method

Meshfree Modeling of a Fluid‐particle Two‐phase Flow with an Improved SPH Method

Material Point Method and Smoothed Particle Hydrodynamics Simulations of Fluid Flow Problems: a Comparative Study

Application of Smoothed Particle Hydrodynamics Method to Free Surface and Solidification Problems

AN IMPROVED GHOST-CELL IMMERSED BOUNDARY METHOD FOR SOLVING SUPERSONIC FLOW PROBLEMS

Dynamical pressure boundary condition for weakly-compressible smoothed particle hydrodynamics