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.

Modeling the Dispersion of Surface Oil Slick under Breaking Waves by Using a Multiphase SPH

Oil Spill Spreading Simulation Based on an Enhanced Multi-phase Consistent Particle Method

Behavior of Oil under Breaking Waves by a Two-phase SPH Model

Fractal Simulation of Oil Spill Trajectory in the Northern Port Area of the Yueqing Bay

Study of Wave-Induced Mass Transport and Internal Mixing Based on a Two-Liquid VOF Model

Simulation of Wave Mixing of Surface Layers in the Sea Area

Research on oil boom performance based on Smoothed Particle Hydrodynamics method

Modeling oil dispersion under breaking waves. Part II: Coupling Lagrangian particle tracking with population balance model

Numerical modeling of oil spill containment by boom using SPH

Modelling of Oil Droplet Kinetics Under Breaking Waves

Modeling Oil Droplet Formation and Evolution under Breaking Waves

Multiphase CFD simulation of the nearshore spilled oil behaviors

A multi-phase SPH simulation of oil spill diffusion in seawater currents

Numerical simulation of multiphase oil behaviors in ice-covered nearshore water

A Modified SPH Model for Simulating Water Surface Waves

Simulating multi-phase sloshing flows with the SPH method

Multiphase SPH Modeling of Water Boiling on Hydrophilic and Hydrophobic Surfaces

Numerical simulation of submarine landslide tsunamis based on the smoothed particle hydrodynamics model

SPH Modeling of Multiphase Drop Dynamics

Smoothed Particle Hydrodynamics Modeling of Free Surface Flow

Numerical investigation of solitary wave breaking over a slope based on multi-phase smoothed particle hydrodynamics