Variable Smoothing Lengths and Energy Conservation in Smoothed Particle Hydrodynamics

Abstract

We present a new formulation of the equations of motion used in smoothed particle hydrodynamics (SPH). The spatial resolution in SPH is determined by the smoothing length, $h$, and it has become common practice for each particle to be given its own adaptive smoothing length, $h_i$. Consequently, the dynamic range that may be spatially resolved is greatly increased, but additional ($\nabla h$) terms, which account for the variability of the smoothing lengths, should appear in the equations of motion in order to satisy conservation requirements. Previous implementations of SPH have neglected these additional terms, whereas we have included them. This is achieved by defining a functional form for the $h_i$s, that depends on inter-particle distances only, and then deriving the equations of motion using a Hamiltonian formalism. A number of test calculations are presented. We find that including the $\nabla h$ terms has no detrimental effect on the ability of SPH to model known problems with reasonable accuracy. For problems where the energy conservation was rather poor, including the $\nabla h$ terms results in a dramatic improvement. In particluar non-conservation of energy during a collision between two polytropes can occur at the level of 10% when the $\nabla h$ terms are neglected. When the $\nabla h$ terms are included, then this error reduces to 0.8%.

29 pages, LaTeX. To appear in M.N.R.A.S

Cited by in corpus (43)

• GADGET: A code for collisionless and gasdynamical cosmological simulations
• Cosmological SPH simulations: The entropy equation
• Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• Smoothed Particle Hydrodynamics in Astrophysics
• Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
• Gasoline2: A Modern SPH Code
• Reformulation of Smoothed Particle Hydrodynamics with Riemann Solver
• Smoothed Particle Hydrodynamics with particle splitting, applied to self-gravitating collapse
• Astrophysical Smooth Particle Hydrodynamics
• Multiphase Smoothed-Particle Hydrodynamics
• Smoothed Particle Magnetohydrodynamics II. Variational principles and variable smoothing length terms
• Resolving mixing in Smoothed Particle Hydrodynamics
• MAGMA: a 3D, Lagrangian magnetohydrodynamics code for merger applications
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• SEREN - A new SPH code for star and planet formation simulations
• Stellar Collisions and Ultracompact X-ray Binary Formation
• Numerical solution of the non-linear Schrodinger equation using smoothed-particle hydrodynamics
• The radial structure of galaxy groups and clusters
• Tests of Spurious Transport in Smoothed Particle Hydrodynamics
• Adaptive resolution for multiphase smoothed particle hydrodynamics
• Quasars: What turns them off?
• A one-parameter family of interpolating kernels for Smoothed Particle Hydrodynamics studies
• VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods
• Hydrodynamic Simulations of Propagating Warps and Bending Waves In Accretion Discs
• Conservative, special-relativistic smoothed particle hydrodynamics
• Smoothed Particle Hydrodynamics in cosmology: a comparative study of implementations
• Density Estimators in Particle Hydrodynamics: DTFE versus regular SPH
• EvoL: The new Padova T-SPH parallel code for cosmological simulations - I. Basic code: gravity and hydrodynamics
• Conservation Laws in Smooth Particle Hydrodynamics: the DEVA Code
• Morphological Evolution of Dwarf Galaxies in the Local Group
• SPH Simulations of Accretion Disks and Narrow Rings
• Evolution of Prolate Molecular Clouds at HII Boundaries: I. Formation of fragment-core structures
• SPH simulations of Shakura-Sunyaev instability at intermediate accretion rates
• Smoothed Particle Hydrodynamics
• Ab initio simulations of accretion disks instability
• Accelerating self-gravitating hydrodynamics simulations with adaptive force updates
• An approach for solving the boundary free edge difficulties in SPH modelling: application to a viscous accretion disc in close binaries
• Formation of accretion centers in simulations of colliding uniform density H${}_2$ cores
• Precessing warped discs in close binary systems
• Dual-support Smoothed Particle Hydrodynamics
• SPH Calculations of Collisions between Main-Sequence Stars
• A Criterion for the choice of the interpolation kernel in Smoothed Particle Hydrodynamics