Abstract
The boundary truncation of the kernel function affects the numerical accuracy and calculation stability of the smooth particle hydrodynamics (SPH) method and has been one of the key research fields for this method. In this paper, an efficient and stable boundary processing method for the SPH method was introduced by adopting an improved boundary interpolation method (i.e., the improved Shepard method) which needs only the sum of direct accumulation for fixed-boundary particles to improve the numerical stability and computational efficiency of the fixed ghost particle method. The improvement effect of the method was demonstrated by comparing it with different interpolation methods using the cases of still water, a wave generated by dam-breaking, and a solitary wave attacking problem with fixed walls and a moveable wall. The results showed that the new boundary processing method for SPH can help remarkably improve the efficiency of calculation and reduce the oscillations of pressure when simulating various flows.
Highlights
The smoothed particle hydrodynamics (SPH) method is a type of Lagrange meshless method for obtaining approximate numerical solutions for equations of fluid dynamics by replacing the fluid with a set of particles
Fluid is considered to be incompressible if its density variation is less than 1% (Monaghan, 2005, [4]), and the pressure varies according to the equation of state (EOS) rather than solving the Poisson equation
The solitary wave impact on a fixed wall with different inclination angles and a moveable wall with a spring-mass system was adopted with different interpolation methods
Summary
The smoothed particle hydrodynamics (SPH) method is a type of Lagrange meshless method for obtaining approximate numerical solutions for equations of fluid dynamics by replacing the fluid with a set of particles. This method was first proposed by Lucy (1977, [1]) and Monaghan and Gingold (1977, [2]) who applied it to astrophysical problems. In this method, fluid is considered to be incompressible if its density variation is less than 1% (Monaghan, 2005, [4]), and the pressure varies according to the equation of state (EOS) rather than solving the Poisson equation. The WCSPH method has become increasingly generalized to handle the pressure stability of flow field with diffusive terms (Antuono, 2012, [8]) and the turbulence model (Padova, 2017, [9]; 2018, [10])
Published Version (Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have