A New Stable Version of the SPH Method in Lagrangian Coordinates

Abstract. The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach. First, we focus on the problem of the stability for two different SPH schemes, one is based on the approach of Vila [25] and another is proposed in this article which mimics the classical 1D LaxWendroff scheme. In both approaches the classical SPH artificial viscosity term is removed preserving nevertheless the linear stability of the methods, demonstrated via the von Neumann stability analysis. Moreover, the issue of the consistency for the equations of gas dynamics is analyzed. An alternative approach is proposed that consists of using Godunov-type SPH schemes in Lagrangian coordinates. This not only provides an improvement in accuracy of the numerical solutions, but also assures that the consistency condition on the gradient of the kernel function is satisfied using an equidistant distribution of particles in Lagrangian mass coordinates. Three different Riemann solvers are implemented for the first-order Godunov type SPH schemes in Lagrangian coordinates, namely the Godunov flux based on the exact Riemann solver, the Rusanov flux and a new modified Roe flux, following the work of Munz [17]. Some well-known numerical 1D shock tube test cases [22] are solved, comparing the numerical solutions of the Godunov-type SPH schemes in Lagrangian coordinates with the first-order Godunov finite volume method in Eulerian coordinates and the standard SPH scheme with Monaghan’s viscosity term. AMS subject classifications: 52B10, 65D18, 68U05, 68U07

[1]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[2]  S. A. Medin,et al.  Improvements in SPH method by means of interparticle contact algorithm and analysis of perforation tests at moderate projectile velocities , 2000 .

[3]  Leonardo Di G. Sigalotti,et al.  A shock-capturing SPH scheme based on adaptive kernel estimation , 2006, J. Comput. Phys..

[4]  J. Monaghan,et al.  A refined particle method for astrophysical problems , 1985 .

[5]  Bram van Leer,et al.  Upwind and High-Resolution Methods for Compressible Flow: From Donor Cell to Residual-Distribution Schemes , 2003 .

[6]  J. Monaghan.  SPH without a Tensile Instability , 2000 .

[7]  G. V. Bicknell.  The Equations of Motion of Particles in Smoothed Particle Hydrodynamics , 1991, SIAM J. Sci. Comput..

[8]  Wing Kam Liu,et al.  Reproducing kernel particle methods , 1995 .

[9]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[10]  K. Y. Lam,et al.  Constructing smoothing functions in smoothed particle hydrodynamics with applications , 2003 .

[11]  B. Ben Moussa,et al.  Convergence of SPH Method for Scalar Nonlinear Conservation Laws , 2000, SIAM J. Numer. Anal..

[12]  Anthony Peter Whitworth,et al.  Implementations and tests of Godunov-type particle hydrodynamics , 2003 .

[13]  Richard Courant,et al.  Supersonic Flow And Shock Waves , 1948 .

[14]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[15]  Javier Bonet,et al.  A corrected smooth particle hydrodynamics method for the simulation of debris flows , 2004 .

[16]  David H. Wagner,et al.  Equivalence of the Euler and Lagrangian equations of gas dynamics for weak solutions , 1987 .

[17]  Jean-Paul Vila,et al.  ON PARTICLE WEIGHTED METHODS AND SMOOTH PARTICLE HYDRODYNAMICS , 1999 .

[18]  M. Lastiwka,et al.  Truncation error in mesh‐free particle methods , 2006 .

[19]  Wing Kam Liu,et al.  Meshfree and particle methods and their applications , 2002 .

[20]  D. Balsara.  von Neumann stability analysis of smoothed particle hydrodynamics—suggestions for optimal algorithms , 1995 .

[21]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

[22]  Claus-Dieter Munz,et al.  On Godunov-type schemes for Lagrangian gas dynamics , 1994 .

[23]  Shu-ichiro Inutsuka.  Reformulation of Smoothed Particle Hydrodynamics with Riemann Solver , 2002 .

[24]  J. Monaghan.  On the problem of penetration in particle methods , 1989 .