A fourth-order kernel for improving numerical accuracy and stability in Eulerian and total Lagrangian SPH

TL;DR

This paper introduces a new fourth-order truncated Laguerre-Gauss kernel for SPH that enhances numerical accuracy and stability without increasing computational cost, outperforming traditional kernels like Wendland.

Contribution

It proposes a novel fourth-order truncated Laguerre-Gauss kernel that reduces errors and avoids pair-instability, improving SPH accuracy and stability.

Findings

Significantly improved numerical accuracy in benchmark tests

Reduced relaxation residue compared to existing kernels

Maintained computational efficiency equivalent to Wendland kernel

Abstract

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.