A simple, entropy-based dissipation trigger for SPH

Abstract

Smoothed Particle Hydrodynamics (SPH) schemes need to be enhanced by dissipation mechanisms to handle shocks. Most SPH formulations rely on artificial viscosity and while this is working well in pure shocks, attention has to be payed to avoid dissipation where it is not wanted. Commonly used approaches include limiters and time-dependent dissipation parameters. The former try to distinguish shocks from other types of flows that do not require dissipation while in the latter approach the dissipation parameters are steered by some source term ("trigger") and, if not triggered, they decay to a pre-described floor value. The commonly used source terms trigger on either compression, $-\nabla \cdot \vec{v}$, or its time derive. Here we explore a novel way to trigger SPH-dissipation: based on the entropy growth rate between two time steps we identify "troubled particles" that need to have dissipation added because they are either passing through a shock wave or become noisy. Our new scheme is implemented into the Lagrangian hydrodynamics code MAGMA2 and scrutinized in a number of shock and fluid instability tests. We find excellent results in shocks and only a moderate (and desired) switch-on in instability tests, despite our conservatively chosen trigger parameters. The new scheme is robust, trivial to implement into existing SPH codes and does not add any computational overhead.

arXiv comment12 pages, 14 figures; changes in response to referee report

References in corpus (20)

• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Cosmological SPH simulations: The entropy equation
• An energy-conserving formalism for adaptive gravitational force softening in SPH and N-body codes
• Inviscid SPH
• Gasoline2: A Modern SPH Code
• Hydromagnetic turbulence in computer simulations
• Reformulation of Smoothed Particle Hydrodynamics with Riemann Solver
• SPHGal: Smoothed Particle Hydrodynamics with improved accuracy for Galaxy simulations
• SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch
• A Well-Posed Kelvin-Helmholtz Instability Test and Comparison
• MAGMA: a 3D, Lagrangian magnetohydrodynamics code for merger applications
• Astrophysical Weighted Particle Magnetohydrodynamics
• CRKSPH - A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme
• Boosting the accuracy of SPH techniques: Newtonian and special-relativistic tests
• rpSPH: a novel Smoothed Particle Hydrodynamics Algorithm
• Improving smoothed particle hydrodynamics with an integral approach to calculating gradients
• The Kelvin-Helmholtz instability and smoothed particle hydrodynamics
• Entropy-limited hydrodynamics: a novel approach to relativistic hydrodynamics
• Testing the concept of integral approach to derivatives within the smoothed particle hydrodynamics technique in astrophysical scenarios
• GenASiS: General Astrophysical Simulation System. I. Refinable Mesh and Nonrelativistic Hydrodynamics

Cited by in corpus (9)

• Sphenix: Smoothed Particle Hydrodynamics for the next generation of galaxy formation simulations
• The Lagrangian hydrodynamics code MAGMA2
• Smoothed Particle Radiation Hydrodynamics: Two-Moment method with Local Eddington Tensor Closure
• Simulating neutron star mergers with the Lagrangian Numerical Relativity code SPHINCS_BSSN
• Conservative, density-based smoothed particle hydrodynamics with improved partition of the unity and better estimation of gradients
• Robust Solid Boundary Treatment for Compressible Smoothed Particle Hydrodynamics
• The Lagrangian Numerical Relativity code SPHINCS_BSSN_v1.0
• Adaptive Compressible Smoothed Particle Hydrodynamics
• Shock-capturing particle hydrodynamics with reproducing kernels