Hydrodynamic simulations with the Godunov SPH

We present results based on an implementation of the Godunov Smoothed Particle Hydrodynamics (GSPH), originally developed by Inutsuka (2002), in the GADGET-3 hydrodynamic code. We first review the derivation of the GSPH discretization of the equations of moment and energy conservation, starting from the convolution of these equations with the interpolating kernel. The two most important aspects of the numerical implementation of these equations are (a) the appearance of fluid velocity and pressure obtained from the solution of the Riemann problem between each pair of particles, and (b the absence of an artificial viscosity term. We carry out three different controlled hydrodynamical three-dimensional tests, namely the Sod shock tube, the development of Kelvin-Helmholtz instabilities in a shear flow test, and the "blob" test describing the evolution of a cold cloud moving against a hot wind. The results of our tests confirm and extend in a number of aspects those recently obtained by Cha (2010): (i) GSPH provides a much improved description of contact discontinuities, with respect to SPH, thus avoiding the appearance of spurious pressure forces; (ii) GSPH is able to follow the development of gas-dynamical instabilities, such as the Kevin--Helmholtz and the Rayleigh-Taylor ones; (iii) as a result, GSPH describes the development of curl structures in the shear-flow test and the dissolution of the cold cloud in the "blob" test. We also discuss in detail the effect on the performances of GSPH of changing different aspects of its implementation. The results of our tests demonstrate that GSPH is in fact a highly promising hydrodynamic scheme, also to be coupled to an N-body solver, for astrophysical and cosmological applications. [abridged]

[1]  Joseph John Monaghan,et al.  SPH and Riemann Solvers , 1997 .

[2]  Walter Dehnen,et al.  Inviscid smoothed particle hydrodynamics , 2010 .

[3]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[4]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[5]  J. Monaghan,et al.  Fundamental differences between SPH and grid methods , 2006, astro-ph/0610051.

[6]  Daniel J. Price Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH , 2007, J. Comput. Phys..

[7]  J. Monaghan,et al.  A Switch to Reduce SPH Viscosity , 1997 .

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

[9]  Volker Springel,et al.  Particle hydrodynamics with tessellation techniques , 2009, 0912.0629.

[10]  Stephan Rosswog,et al.  Astrophysical smooth particle hydrodynamics , 2009, 0903.5075.

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

[12]  H. Dejonghe,et al.  Kelvin-Helmholtz instabilities in smoothed particle hydrodynamics , 2010, 1006.1767.

[13]  V. Springel E pur si muove: Galilean-invariant cosmological hydrodynamical simulations on a moving mesh , 2009, 0901.4107.

[14]  J. Monaghan,et al.  Shock simulation by the particle method SPH , 1983 .

[15]  O. Agertz,et al.  Resolving mixing in smoothed particle hydrodynamics , 2009, 0906.0774.

[16]  V. Springel Smoothed Particle Hydrodynamics in Astrophysics , 2010, 1109.2219.

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

[18]  Nickolay Y. Gnedin,et al.  Computational Eulerian hydrodynamics and Galilean invariance , 2009, 0909.0513.

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

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

[21]  F. Vazza,et al.  Turbulent gas motions in galaxy cluster simulations: the role of smoothed particle hydrodynamics viscosity , 2005 .

[22]  Stuart L. Shapiro,et al.  Tests of Spurious Transport in Smoothed Particle Hydrodynamics , 1999 .

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

[24]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[25]  Sergei Nayakshin,et al.  Kelvin-Helmholtz instabilities with Godunov smoothed particle hydrodynamics , 2010 .

[26]  Volker Springel,et al.  Cosmological SPH simulations: The entropy equation , 2001 .

[27]  James Wadsley,et al.  On the treatment of entropy mixing in numerical cosmology , 2008 .

[28]  D. O. Astronomy,et al.  EvoL: The new Padova T-SPH parallel code for cosmological simulations - I. Basic code: gravity and hydrodynamics , 2009, 0912.1786.

[29]  V. Springel,et al.  Cosmological smoothed particle hydrodynamics simulations: the entropy equation , 2001, astro-ph/0111016.

[30]  V. Springel The Cosmological simulation code GADGET-2 , 2005, astro-ph/0505010.

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

[32]  G. Sod A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .