Commun. Comput. Phys.,
Notice: Undefined index: year in /var/www/html/issue/abstract/readabs.php on line 20

Notice: Undefined index: ppage in /var/www/html/issue/abstract/readabs.php on line 21

Notice: Undefined index: issue in /var/www/html/issue/abstract/readabs.php on line 23
Volume 4.

A New Stable Version of the SPH Method in Lagrangian Coordinates

A. Ferrari 1, M. Dumbser 1, E. F. Toro 1, A. Armanini 1*

1 Laboratory of Applied Mathematics, Faculty of Engineering, University of Trento, Trento 38100, Italy.

Received 15 October 2007; Accepted (in revised version) 18 January 2008
Available online 18 March 2008


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 \cite{Vila} and another is proposed in this article which mimics the classical 1D Lax-Wendroff 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 \cite{Munz}. Some well-known numerical 1D shock tube test cases \cite{Toro} 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

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: SPH, meshfree particle methods, Riemann solvers, gas dynamics in Lagrangian coordinates, Godunov type schemes.

*Corresponding author.
Email: (A. Ferrari), (M. Dumbser), (E. F. Toro), (A. Armanini)

The Global Science Journal