6533b85ffe1ef96bd12c2676

RESEARCH PRODUCT

Numerical viscosity in simulations of the two-dimensional Kelvin-Helmholtz instability

Miguel-ángel AloyM. Obergaulinger

subject

PhysicsHistoryNumerical analysisFOS: Physical sciences010103 numerical & computational mathematicsMechanicsComputational Physics (physics.comp-ph)Dissipation01 natural sciencesInstabilityComputer Science ApplicationsEducationPiecewise linear functionViscositySimple (abstract algebra)0103 physical sciencesConvergence (routing)Magnetohydrodynamic drive0101 mathematicsAstrophysics - Instrumentation and Methods for AstrophysicsInstrumentation and Methods for Astrophysics (astro-ph.IM)Physics - Computational Physics010303 astronomy & astrophysics

description

The Kelvin-Helmholtz instability serves as a simple, well-defined setup for assessing the accuracy of different numerical methods for solving the equations of hydrodynamics. We use it to extend our previous analysis of the convergence and the numerical dissipation in models of the propagation of waves and in the tearing-mode instability in magnetohydrodynamic models. To this end, we perform two-dimensional simulations with and without explicit physical viscosity at different resolutions. A comparison of the growth of the modes excited by our initial perturbations allows us to estimate the effective numerical viscosity of two spatial reconstruction schemes (fifth-order monotonicity preserving and second-order piecewise linear schemes).

https://doi.org/10.1088/1742-6596/1623/1/012018