A robust numerical scheme for highly compressible magnetohydrodynamics: nonlinear stability, implementation and tests
From MaRDI portal
(Redirected from Publication:543751)
Abstract: The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate that the theoretical stability properties are valid in practice with negligible compromises to accuracy. The result is a highly robust scheme with state-of-the-art efficiency. The scheme's robustness is due to entropy stability, positivity and properly discretised Powell terms. The implementation takes the form of a modification of the MHD module in the FLASH code, an adaptive mesh refinement code. We compare the new scheme with the standard FLASH implementation for MHD. Results show comparable accuracy to standard FLASH with the Roe solver, but highly improved efficiency and stability, particularly for high Mach number flows and low plasma beta. The tests include 1D shock tubes, 2D instabilities and highly supersonic, 3D turbulence. We consider turbulent flows with RMS sonic Mach numbers up to 10, typical of gas flows in the interstellar medium. We investigate both strong initial magnetic fields and magnetic field amplification by the turbulent dynamo from extremely high plasma beta. The energy spectra show a reasonable decrease in dissipation with grid refinement, and at a resolution of 512^3 grid cells we identify a narrow inertial range with the expected power-law scaling. The turbulent dynamo exhibits exponential growth of magnetic pressure, with the growth rate twice as high from solenoidal forcing than from compressive forcing. Two versions of the new scheme are presented, using relaxation-based 3-wave and 5-wave approximate Riemann solvers, respectively. The 5-wave solver is more accurate in some cases, and its computational cost is close to the 3-wave solver.
Recommendations
- A positive MUSCL-Hancock scheme for ideal magnetohydrodynamics
- scientific article; zbMATH DE number 1284337
- Entropy stable finite volume scheme for ideal compressible MHD on 2-D Cartesian meshes
- Approximate Riemann solvers and robust high-order finite volume schemes for multi-dimensional ideal MHD equations
- High resolution entropy stable scheme for ideal MHD equations
Cites work
- scientific article; zbMATH DE number 1324011 (Why is no real title available?)
- scientific article; zbMATH DE number 2174321 (Why is no real title available?)
- scientific article; zbMATH DE number 821171 (Why is no real title available?)
- scientific article; zbMATH DE number 914383 (Why is no real title available?)
- scientific article; zbMATH DE number 3105271 (Why is no real title available?)
- A method for incorporating Gauss' law into electromagnetic PIC codes
- A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics
- A multiwave approximate Riemann solver for ideal MHD based on relaxation. I: Theoretical framework
- A multiwave approximate Riemann solver for ideal MHD based on relaxation. II: Numerical implementation with 3 and 5 waves
- A positive MUSCL-Hancock scheme for ideal magnetohydrodynamics
- An HLLC Riemann solver for magneto-hydrodynamics
- An HLLC-Type Approximate Riemann Solver for Ideal Magnetohydrodynamics
- An examination of forcing in direct numerical simulations of turbulence
- An unsplit staggered mesh scheme for multidimensional magnetohydrodynamics
- An upwind differencing scheme for the equations of ideal magnetohydrodynamics
- Approximate Riemann solvers and robust high-order finite volume schemes for multi-dimensional ideal MHD equations
- Approximate Riemann solvers, parameter vectors, and difference schemes
- Entropy satisfying flux vector splittings and kinetic BGK models
- Hydromagnetic flows from accretion discs and the production of radio jets
- Hyperbolic divergence cleaning for the MHD equations
- Numerical comparison of Riemann solvers for astrophysical hydrodynamics
- On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws
- On positivity preserving finite volume schemes for Euler equations
- On the Relation Between the Upwind-Differencing Schemes of Godunov, Engquist–Osher and Roe
- Relaxation solvers for ideal MHD equations-a review
- Small-scale dynamo action during the formation of the first stars and galaxies. I: The ideal MHD limit
- Stable upwind schemes for the magnetic induction equation
- The \(\nabla \cdot B=0\) constraint in shock-capturing magnetohydrodynamics codes
- The effect of nonzero \(\bigtriangledown\cdot B\) on the numerical solution of the magnetohydrodynamic equations
- The piecewise parabolic method (PPM) for gas-dynamical simulations
Cited in
(25)- An all Mach number scheme for visco-resistive magnetically-dominated MHD flows
- A limit problem for three-dimensional ideal compressible radiation magneto-hydrodynamics
- A high performance spectral code for nonlinear MHD stability
- An efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes
- A multi well-balanced scheme for the shallow water MHD system with topography
- A 5-wave relaxation solver for the shallow water MHD system
- Entropy stable finite volume approximations for ideal magnetohydrodynamics
- Positivity-preserving finite difference WENO scheme for ten-moment equations with source term
- Systematic construction of upwind constrained transport schemes for MHD
- Modelling and entropy satisfying relaxation scheme for the nonconservative bitemperature Euler system with transverse magnetic field
- Locally divergence-free well-balanced path-conservative central-upwind schemes for rotating shallow water MHD
- A positive MUSCL-Hancock scheme for ideal magnetohydrodynamics
- A uniquely defined entropy stable matrix dissipation operator for high Mach number ideal MHD and compressible Euler simulations
- A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure
- A New Discretely Divergence-Free Positivity-Preserving High-Order Finite Volume Method for Ideal MHD Equations
- A positivity-preserving Lagrangian discontinuous Galerkin method for ideal magnetohydrodynamics equations in one-dimension
- A provably positive discontinuous Galerkin method for multidimensional ideal magnetohydrodynamics
- A semi-Lagrangian transport method for kinetic problems with application to dense-to-dilute polydisperse reacting spray flows
- The Developments and the Applications of the Numerical Algorithms in Simulating the Incompressible Magnetohydrodynamics with Complex Boundaries and Free Surfaces
- Compact third-order accurate, positive preserving and divergence-free scheme
- Provably Positive Central Discontinuous Galerkin Schemes via Geometric Quasilinearization for Ideal MHD Equations
- Numerical analysis of an artificial compression method for magnetohydrodynamic flows at low magnetic Reynolds numbers
- Geostrophic vs magneto-geostrophic adjustment and nonlinear magneto-inertia-gravity waves in rotating shallow water magnetohydrodynamics
- A new locally divergence-free path-conservative central-upwind scheme for ideal and shallow water magnetohydrodynamics
- Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes
This page was built for publication: A robust numerical scheme for highly compressible magnetohydrodynamics: nonlinear stability, implementation and tests
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q543751)
