An efficient and generalized consistency correction method for weakly-compressible SPH
arXiv:2304.14865 · doi:10.1016/j.cma.2023.116460
Abstract
In this paper, a new efficient and generalized consistency correction method for weakly-compressible smoothed particle hydrodynamics is proposed and successfully implemented in the simulation of violent free-surface flow exhibiting breaking and impact events for the first time. It's well known that the original kernel gradient correction (KGC) encounters numerical instability resulting from matrix inversion. The present method remedies this issue by introducing a weighted average of the KGC matrix and the identity matrix, other than directly applying KGC matrix, to achieve numerical stability meanwhile decrease numerical dissipation. To ensure momentum conservation, the correction is implemented in a particle-average pattern by rewriting the the pressure term of the Riemann solution. Furthermore, the proposed weighted KGC scheme is incorporated into the dual-criteria time-stepping framework developed by Zhang et al. (2020) \cite{22} to achieve optimized computational efficiency. A set of numerical examples in both two- and three-dimensions are investigated to demonstrate that the present method can significantly reduce numerical dissipation meanwhile exhibit a smooth pressure field for general free-surface flows.
arXiv comment44 pages 21 figures
References in corpus (7)
- Experimental investigation of dynamic pressure loads during dam break
- Smoothed Particle Hydrodynamics in Astrophysics
- SPHinXsys: an open-source multi-physics and multi-resolution library based on smoothed particle hydrodynamics
- Review on Smoothed Particle Hydrodynamics: Methodology development and recent achievement
- Dual-criteria time stepping for weakly compressible smoothed particle hydrodynamics
- An incremental-stencil WENO reconstruction for simulation of compressible two-phase flows
- Stability and accuracy of the weakly compressible SPH with particle regularization techniques