Hydrodynamic capabilities of an SPH code incorporating an artificial conductivity term with a gravity-based signal velocity

Abstract

This paper investigates the hydrodynamic performances of an SPH code incorporating an artificial heat conductivity term in which the adopted signal velocity is applicable when gravity is present. In accordance with previous findings it is shown that the performances of SPH to describe the development of Kelvin-Helmholtz instabilities depend strongly on the consistency of the initial condition set-up and on the leading error in the momentum equation due to incomplete kernel sampling. An error and stability analysis shows that the quartic B-spline kernel (M_5) possesses very good stability properties and we propose its use with a large neighbor number, between ~50 (2D) to ~ 100 (3D), to improve convergence in simulation results without being affected by the so-called clumping instability. SPH simulations of the blob test show that in the regime of strong supersonic flows an appropriate limiting condition, which depends on the Prandtl number, must be imposed on the artificial conductivity SPH coefficients in order to avoid an unphysical amount of heat diffusion. Results from hydrodynamic simulations that include self-gravity show profiles of hydrodynamic variables that are in much better agreement with those produced using mesh-based codes. In particular, the final levels of core entropies in cosmological simulations of galaxy clusters are consistent with those found using AMR codes. Finally, results of the Rayleigh-Taylor instability test demonstrate that in the regime of very subsonic flows the code has still several difficulties in the treatment of hydrodynamic instabilities. These problems being intrinsically due to the way in which in standard SPH gradients are calculated and not to the implementation of the artificial conductivity term.

26 pages, 15 figures, accepted for publication in A&A

References in corpus (9)

• Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
• Inviscid SPH
• Astrophysical Smooth Particle Hydrodynamics
• A test suite for quantitative comparison of hydrodynamics codes in astrophysics
• On the Origin of Cores in Simulated Galaxy Clusters
• MAGMA: a 3D, Lagrangian magnetohydrodynamics code for merger applications
• VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods

Cited by in corpus (21)

• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Simulating cosmic structure formation with the GADGET-4 code
• SPHGal: Smoothed Particle Hydrodynamics with improved accuracy for Galaxy simulations
• The Magellanic Stream System: I. Ram-pressure tails and the relics of the collision between the Magellanic Clouds
• Boosting the accuracy of SPH techniques: Newtonian and special-relativistic tests
• nIFTy galaxy cluster simulations I: dark matter & non-radiative models
• Temperature Structure of the Intra-Cluster Medium from SPH and AMR simulations
• nIFTy galaxy cluster simulations II: radiative models
• SPHYNX: an accurate density-based SPH method for astrophysical applications
• The Lagrangian hydrodynamics code MAGMA2
• On the effect of galactic outflows in cosmological simulations of disc galaxies
• A multifiltering study of turbulence in a large sample of simulated galaxy clusters
• The Kelvin-Helmholtz instability and smoothed particle hydrodynamics
• A study of cool core resiliency and entropy mixing in simulations of galaxy cluster mergers
• The role of the artificial conductivity in SPH simulations of galaxy clusters: effects on the ICM properties
• Santa Barbara Cluster Comparison Test with DISPH
• Improved Performances in Subsonic Flows of an SPH Scheme with Gradients Estimated using an Integral Approach
• Density Independent Smoothed Particle Hydrodynamics for Non-Ideal Equation of State
• Equalizing resolution in smoothed-particle hydrodynamics calculations using self-adaptive sinc kernels
• Excess entropy and energy feedback from within cluster cores up to r${}_{200}$
• Lagrangian Methods Of Cosmic Web Classification