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 Coupled Phase-Field and Volume-of-fluid Method for Accurate Representation of Limiting Water Wave Deformation

Numerical Simulation of Two-Dimensional Dam Failure and Free-Side Deformation Flow Studies

A CIP-based Numerical Simulation of Wave Interaction with a Fluid-Filled Membrane Submerged Breakwater

Numerical Study on the Wave Dissipating Performance of a Submerged Heaving Plate Breakwater

High-fidelity Simulation of Regular Waves Based on Multi-Moment Finite Volume Formulation and THINC Method

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

Numerical Simulation of Wave Interaction with Coastal Structures Using a CIP-based Method

A CIP-based Numerical Simulation of Freak Wave Impact on a Floating Body

Numerical Simulation of the Interaction Between Solitary Waves and Underwater Barriers Using a VPM–THINC/QQ-Coupled Model

Coupled Computational Fluid Dynamics-Discrete Element Method for Simulating the Interactions Between Ship-Induced Waves and Riprap on Restricted Waterway Banks

Development of a Volume of Fluid Method for Computing Interfacial Incompressible Fluid Flows

Impact of Multiphase Flow Simulation of Breaking Waves on a Vertical Seawall with a Recurved Parapet

Numerical Study of Violent Impact Flow Using a CIP-Based Model

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

Applications of Improved 2-D Numerical Wave Tanks in Wave Breaking

Enhanced Ocean Wave Modeling by Including Effect of Breaking under Both Deep- and Shallow-Water Conditions

Numerical Simulation of Fluid–structure Interaction Using a Combined Volume of Fluid and Immersed Boundary Method

Numerical Study of Air Compressibility Effects in Breaking Wave Impacts Using a CIP-based Model

A Modified SPH Model for Simulating Water Surface Waves

Two-Dimensional Free-Surface Flow Modeling for Wave-Structure Interactions and Induced Motions of Floating Bodies

Experimental and Numerical Study of Regular Waves Past a Submerged Breakwater