A Well-Posed Kelvin-Helmholtz Instability Test and Comparison

Abstract

Recently, there has been a significant level of discussion of the correct treatment of Kelvin-Helmholtz instability in the astrophysical community. This discussion relies largely on how the KHI test is posed and analyzed. We pose a stringent test of the initial growth of the instability. The goal is to provide a rigorous methodology for verifying a code on two dimensional Kelvin-Helmholtz instability. We ran the problem in the Pencil Code, Athena, Enzo, NDSPHMHD, and Phurbas. A strict comparison, judgment, or ranking, between codes is beyond the scope of this work, though this work provides the mathematical framework needed for such a study. Nonetheless, how the test is posed circumvents the issues raised by tests starting from a sharp contact discontinuity yet it still shows the poor performance of Smoothed Particle Hydrodynamics. We then comment on the connection between this behavior to the underlying lack of zeroth-order consistency in Smoothed Particle Hydrodynamics interpolation. We comment on the tendency of some methods, particularly those with very low numerical diffusion, to produce secondary Kelvin-Helmholtz billows on similar tests. Though the lack of a fixed, physical diffusive scale in the Euler equations lies at the root of the issue, we suggest that in some methods an extra diffusion operator should be used to damp the growth of instabilities arising from grid noise. This statement applies particularly to moving-mesh tessellation codes, but also to fixed-grid Godunov schemes.

18 pages, 14 figures, ApJS accepted, sections reorganized for clarity and discussion extended

References in corpus (12)

• Athena: A New Code for Astrophysical MHD
• SPLASH: An interactive visualisation tool for Smoothed Particle Hydrodynamics simulations
• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• Smoothed Particle Hydrodynamics in Astrophysics
• Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
• Dust sedimentation and self-sustained Kelvin-Helmholtz turbulence in protoplanetary disk mid-planes. I. Radially symmetric simulations
• SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch
• Global magnetohydrodynamical models of turbulence in protoplanetary disks I. A cylindrical potential on a Cartesian grid and transport of solids
• Introducing Enzo, an AMR Cosmology Application
• Stochastic oscillations of general relativistic disks
• Phurbas: An Adaptive, Lagrangian, Meshless, Magnetohydrodynamics Code. I. Algorithm
• Phurbas: An Adaptive, Lagrangian, Meshless, Magnetohydrodynamics Code. II. Implementation and Tests

Cited by in corpus (80)

• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Enzo: An Adaptive Mesh Refinement Code for Astrophysics
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
• Global simulations of protoplanetary disks with ohmic resistivity and ambipolar diffusion
• Gasoline2: A Modern SPH Code
• THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics
• Accurate, Meshless Methods for Magneto-Hydrodynamics
• A planetesimal orbiting within the debris disc around a white dwarf star
• Images of Embedded Jovian Planet Formation At A Wide Separation Around AB Aurigae
• Convective overstability in accretion disks: 3D linear analysis and nonlinear saturation
• On the Impact of Three Dimensions in Simulations of Neutrino-Driven Core-Collapse Supernova Explosions
• A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics
• CRKSPH - A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme
• Rossby wave instability at dead zone boundaries in 3D resistive magnetohydrodynamical global models of protoplanetary disks
• New numerical solver for flows at various Mach numbers
• Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement
• Rossby wave instability does not require sharp resistivity gradients
• Viscous Kelvin-Helmholtz instabilities in highly ionised plasmas
• The impact of magnetic fields on cold streams feeding galaxies
• On shocks driven by high-mass planets in radiatively inefficient disks. II. Three-dimensional global disk simulations
• A new hydrodynamics code for Type Ia Supernovae
• White Dwarf Mergers on Adaptive Meshes I. Methodology and Code Verification
• On the Kelvin-Helmholtz instability with smooth initial conditions -- Linear theory and simulations
• SPHYNX: an accurate density-based SPH method for astrophysical applications
• The Lagrangian hydrodynamics code MAGMA2
• The Cosmological Simulation Code OpenGadget3 -- Implementation of Meshless Finite Mass
• Benchmarking the Multi-dimensional Stellar Implicit Code MUSIC
• Super massive black holes in star forming gaseous circumnuclear discs
• Numerical code for multi-component galaxies: from N-body to chemistry and magnetic fields
• Local simulations of MRI turbulence with meshless methods
• Hydrodynamic capabilities of an SPH code incorporating an artificial conductivity term with a gravity-based signal velocity
• The 3D MHD code GOEMHD3 for large-Reynolds-number astrophysical plasmas
• On shocks driven by high-mass planets in radiatively inefficient disks. I. Two-dimensional global disk simulations
• The interplay between radiation pressure and the photoelectric instability in optically thin disks of gas and dust
• Extension of the sub-grid-scale gradient model for compressible magnetohydrodynamics turbulent instabilities
• Planet-disc interaction on a freely moving mesh
• Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
• The Kelvin-Helmholtz instability and smoothed particle hydrodynamics
• The role of the artificial conductivity in SPH simulations of galaxy clusters: effects on the ICM properties
• Braginskii viscosity on an unstructured, moving mesh accelerated with super-time-stepping
• EDGE: A new approach to suppressing numerical diffusion in adaptive mesh simulations of galaxy formation
• Quantifying energetics and dissipation in magnetohydrodynamic turbulence
• Density Independent Smoothed Particle Hydrodynamics for Non-Ideal Equation of State
• Improved Performances in Subsonic Flows of an SPH Scheme with Gradients Estimated using an Integral Approach
• Molecular cloud formation as seen in synthetic Hi and molecular gas observations
• A simple, entropy-based dissipation trigger for SPH
• Simulating neutron star mergers with the Lagrangian Numerical Relativity code SPHINCS_BSSN
• The role of physical and numerical viscosity in hydrodynamical instabilities
• Simulating Hydrodynamics in Cosmology with CRK-HACC
• Smoothed particle magnetohydrodynamics with the geometric density average force expression
• Conservative, density-based smoothed particle hydrodynamics with improved partition of the unity and better estimation of gradients
• Suppression of the superfluid Kelvin-Helmholtz instability due to massive vortex cores, friction and confinement
• Equalizing resolution in smoothed-particle hydrodynamics calculations using self-adaptive sinc kernels
• Semi-conservative reduced speed of sound technique for low Mach number flows with large density variations
• High-order Discontinuous Galerkin hydrodynamics with sub-cell shock capturing on GPUs
• Sensitivity of non-radiative cloud-wind interactions to the hydrodynamics solver
• High-Frequency Voronoi Noise Reduced by Smoothed Mesh Motion
• WVTICs -- SPH initial conditions for everyone
• REMIX SPH -- improving mixing in smoothed particle hydrodynamics simulations using a generalised, material-independent approach
• Influence of adaptive mesh refinement and the hydro solver on shear-induced mass stripping in a minor-merger scenario
• WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
• Multidimensional upwind hydrodynamics on unstructured meshes using Graphics Processing Units I. Two-dimensional uniform meshes
• Axisymmetric smoothed particle magneto-hydrodynamics
• Numerical viscosity in simulations of the two-dimensional Kelvin-Helmholtz instability
• Assessing the role of the Kelvin-Helmholtz instability at the QCD cosmological phase transition
• Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows
• An unstable mode of the stratified atmosphere under the non-traditional Coriolis acceleration
• A Formulation of Consistent Particle Hydrodynamics in Strong Form
• Projective and Telescopic Projective Integration for Non-Linear Kinetic Mixtures
• Dealing with density discontinuities in planetary SPH simulations
• A New Residual Distribution Hydrodynamics Solver for Astrophysical Simulations
• An implicit-explicit solver for a two-fluid single-temperature model
• No dilute core produced in simulations of giant impacts on to Jupiter
• The Kelvin-Helmholtz instability and smoothed particle hydrodynamics
• Low mass planet migration in Hall-affected disks
• DISPATCH methods: an approximate, entropy-based Riemann solver for ideal magnetohydrodynamics
• Sharp eccentric rings in planetless hydrodynamical models of debris disks
• Shock-capturing particle hydrodynamics with reproducing kernels
• Optimized smoothing kernels for SPH