Papers › A fourth-order kernel for improving numerical accuracy and stability in Eulerian and...
A fourth-order kernel for improving numerical accuracy and stability in Eulerian and total Lagrangian SPH
Zhentong Wang, Bo Zhang, Oskar J. Haidn, Xiangyu Hu
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
The error of smoothed particle hydrodynamics (SPH) using kernel for particle-based approximation mainly comes from smoothing and integration errors. The choice of kernels has a significant impact on the numerical accuracy, stability and computational efficiency. At present, the most popular kernels such as B-spline, truncated Gaussian (for compact support), Wendland kernels have 2nd-order smoothing error and Wendland kernel becomes mainstream in SPH community as its stability and accuracy. Due to the fact that the particle distribution after relaxation can achieve fast convergence of integration error respected to support radius, it is logical to choose kernels with higher-order smoothing error to improve the numerical accuracy. In this paper, the error of 4th-order Laguerre-Wendland kernel proposed by Litvinov et al. \cite{litvinov2015towards} is revisited and another 4th-order truncated Laguerre-Gauss kernel is further analyzed and considered to replace the widely used Wendland kernel. The proposed kernel has following three properties: One is that it avoids the pair-instability problem during the relaxation process, unlike the original truncated Gaussian kernel, and achieves much less relaxation residue than Wendland and Laguerre-Wendland kernels; One is the truncated compact support size is the same as the non-truncated compact support of Wendland kernel, which leads to both kernels' computational efficiency at the same level; Another is that the truncation error of this kernel is much less than that of Wendland kernel. Furthermore, a comprehensive set of 2D and 3D benchmark cases on Eulerian SPH for fluid dynamics and total Lagrangian SPH for solid dynamics validate the considerably improved numerical accuracy by using truncated Laguerre-Gauss kernel without introducing extra computational effort.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections