A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics

Abstract

The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form.

Reference solution snapshots can be found at https://doi.org/10.5281/zenodo.5803759

References in corpus (7)

• Athena: A New Code for Astrophysical MHD
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
• Dust sedimentation and self-sustained Kelvin-Helmholtz turbulence in protoplanetary disk mid-planes. I. Radially symmetric simulations
• An optimisation approach for analysing nonlinear stability with transition to turbulence in fluids as an exemplar
• Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement
• An Explicit Scheme for Incorporating Ambipolar Diffusion in a Magnetohydrodynamics Code
• Conduction in low Mach number flows: Part I Linear & weakly nonlinear regimes

Cited by in corpus (60)

• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Dedalus: A Flexible Framework for Numerical Simulations with Spectral Methods
• Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
• The Black Hole Accretion Code
• A Characteristic Scale for Cold Gas
• The Impact of Star Formation Feedback on the Circumgalactic Medium
• Multiphase Gas and the Fractal Nature of Radiative Turbulent Mixing Layers
• Fornax: a Flexible Code for Multiphysics Astrophysical Simulations
• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• Filament formation in wind-cloud interactions. II. Clouds with turbulent density, velocity, and magnetic fields
• DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks
• Magnetic Fields in the Formation of the First Stars. I. Theory vs. Simulation
• Common Envelope Evolution on a Moving Mesh
• Planetary Giant Impacts: Convergence of High-Resolution Simulations using Efficient Spherical Initial Conditions and SWIFT
• The impact of magnetic fields on cold streams feeding galaxies
• From the butterfly effect to intrinsic randomness: the spontaneous growth of singular shear flows
• The Efficiency of Magnetic Field Amplification at Shocks by Turbulence
• On the Kelvin-Helmholtz instability with smooth initial conditions -- Linear theory and simulations
• Dissipative Magnetohydrodynamics for Non-Resistive Relativistic Plasmas
• Conservative finite volume scheme for first-order viscous relativistic hydrodynamics
• Benchmarking the Multi-dimensional Stellar Implicit Code MUSIC
• Mixing matters
• Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs
• Preventing pressure oscillations does not fix local linear stability issues of entropy-based split-form high-order schemes
• Extension of the sub-grid-scale gradient model for compressible magnetohydrodynamics turbulent instabilities
• Convective dynamics and disequilibrium chemistry in the atmospheres of giant planets and brown dwarfs
• The Kelvin-Helmholtz instability and smoothed particle hydrodynamics
• Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
• The moving mesh code Shadowfax
• Braginskii viscosity on an unstructured, moving mesh accelerated with super-time-stepping
• Convective heat transport in stratified atmospheres at low and high Mach number
• A Fourth-Order Adaptive Mesh Refinement Algorithm for the Multicomponent, Reacting Compressible Navier-Stokes Equations
• A Moving Mesh Hydrodynamic Solver for ChaNGa
• On the origin of intrinsic randomness of Rayleigh-Benard turbulence
• The role of physical and numerical viscosity in hydrodynamical instabilities
• Multidimensional Simulations of Core Convection
• Coupling of Damped and Growing Modes in Unstable Shear Flow
• General Relativistic Hydrodynamics on a Moving-mesh I: Static Spacetimes
• The density distribution of accreting cosmic filaments as shaped by Kelvin-Helmholtz instability
• Near-cancellation of up- and down-gradient momentum transport in forced magnetized shear-flow turbulence
• High-order Discontinuous Galerkin hydrodynamics with sub-cell shock capturing on GPUs
• Sensitivity of non-radiative cloud-wind interactions to the hydrodynamics solver
• The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes
• Nonlinear mode coupling and energetics of driven magnetized shear-flow turbulence
• Moving and Reactive Boundary Conditions in Moving-Mesh Hydrodynamics
• Multidimensional upwind hydrodynamics on unstructured meshes using Graphics Processing Units I. Two-dimensional uniform meshes
• WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
• Numerical viscosity in simulations of the two-dimensional Kelvin-Helmholtz instability
• Resistive instabilities in sinusoidal shear flows with a streamwise magnetic field
• Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows
• First-Order Viscous Relativistic Hydrodynamics on the Two-Sphere
• Mass-Conserving Implicit-Explicit Methods for Coupled Compressible Navier-Stokes Equations
• A general relativistic magnetohydrodynamics extension to mesh-less schemes in the code GIZMO
• A scalable exponential-DG approach for nonlinear conservation laws: with application to Burger and Euler equations
• A New Residual Distribution Hydrodynamics Solver for Astrophysical Simulations
• The Kelvin-Helmholtz instability and smoothed particle hydrodynamics
• Large-Scale Dynamos Driven by Shear-Flow-Induced Jets
• Pade: A code for protoplanetary disk turbulence based on Pade differencing
• Optimized smoothing kernels for SPH
• Towards Exascale Simulations of the ICM Dynamo with WENO-Wombat