Modelling supernova-driven turbulence

DOI10.1080/03091929.2019.1634705zbMATH Open1482.76070arXiv1806.01570OpenAlexW3103569200WikidataQ114640056 ScholiaQ114640056MaRDI QIDQ5065696FDOQ5065696

Authors: Frederick A. Gent, Mordecai-Mark Mac Low, Maarit J. Käpylä, Graeme R. Sarson, James F. Hollins 

Publication date: 22 March 2022

Published in: Geophysical & Astrophysical Fluid Dynamics (Search for Journal in Brave)

Abstract: High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to handling high Mach number shocks, such as those produced by supernovae (SNe). Here, we outline methods required to enable the code to efficiently and accurately model SNe, using parameters that allow stable simulation of SN-driven turbulence, in order to construct a physically realistic galactic dynamo model. These include the resolution of shocks with artificial viscosity, thermal conductivity, and mass diffusion; the correction of the mass diffusion terms; and a novel generalization of the Courant condition to include all source terms in the momentum and energy equations. We test our methods with the numerical solution of the one-dimensional (1D) Riemann shock tube (Sod, J. Comput. Phys. 1978, 27), also extended to a 1D adiabatic shock with parameters and Mach number relevant to SN shock evolution, including shocks with radiative losses. We extend our test with the three-dimensional (3D) numerical simulation of individual SN remnant evolution for a range of ambient gas densities typical of the interstellar medium and compare these to the analytical solutions of Sedov-Taylor (adiabatic) and the snowplough and Cioffi, McKee and Bertschinger (Astrophys. J. 1988, 334) results incorporating cooling and heating processes. We show that our new timestep algorithm leads to linear rather than quadratic resolution dependence as the strength of the artificial viscosity varies, because of the corresponding change in the strength of interzone gradients.

Full work available at URL: https://arxiv.org/abs/1806.01570

Recommendations

• Piecewise parabolic method on a local stencil for magnetized supersonic turbulence simulation
• Introduction
• A numerical study of supersonic turbulence
• The superN-project: current progress in modelling core collapse supernovae
• Toward TFlop simulations of supernovae

zbMATH Keywords

numerical methodsinstabilitiesartificial diffusivityhigh Mach number shockssupernova-driven turbulence

Mathematics Subject Classification ID

Turbulence (76F99) Shock waves and blast waves in fluid mechanics (76L05) Hydrodynamic and hydromagnetic problems in astronomy and astrophysics (85A30)

Cites Work

• A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws
• Title not available (Why is that?)
• The formation of a blast wave by a very intense explosion I. Theoretical discussion
• Hydromagnetic turbulence in computer simulations
• A 3D MHD model of astrophysical flows: Algorithms, tests and parallelisation
• Large-scale dynamos in turbulent convection with shear

Cited In (7)

• Reducing Turn-around Times for Supernova Simulations
• Spectral characterisation of inertial particle clustering in turbulence
• Mean fields and fluctuations in compressible magnetohydrodynamic flows
• Introduction
• Supernova remnant evolution in uniform and non-uniform media
• Title not available (Why is that?)
• The supernova fragmentation model of solar system formation

Uses Software

• ZEUS

This page was built for publication: Modelling supernova-driven turbulence

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5065696)