A robust numerical scheme for highly compressible magnetohydrodynamics: nonlinear stability, implementation and tests
From MaRDI portal
Publication:543751
DOI10.1016/J.JCP.2011.01.026zbMATH Open1316.76082DBLPjournals/jcphy/WaaganFK11arXiv1101.3007OpenAlexW2147617663WikidataQ62623661 ScholiaQ62623661MaRDI QIDQ543751FDOQ543751
C. Klingenberg, C. Federrath, Knut Waagan
Publication date: 17 June 2011
Published in: Journal of Computational Physics (Search for Journal in Brave)
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.
Full work available at URL: https://arxiv.org/abs/1101.3007
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
General questions in astronomy and astrophysics (85A04) Finite volume methods for initial value and initial-boundary value problems involving PDEs (65M08) Magnetohydrodynamics and electrohydrodynamics (76W05) Turbulence (76F99) Finite volume methods applied to problems in fluid mechanics (76M12)
Cites Work
- On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws
- Title not available (Why is that?)
- Approximate Riemann solvers, parameter vectors, and difference schemes
- Title not available (Why is that?)
- Approximate Riemann Solvers and Robust High-Order Finite Volume Schemes for Multi-Dimensional Ideal MHD Equations
- The piecewise parabolic method (PPM) for gas-dynamical simulations
- Title not available (Why is that?)
- Title not available (Why is that?)
- On the Relation Between the Upwind-Differencing Schemes of Godunov, Engquist–Osher and Roe
- The \(\nabla \cdot B=0\) constraint in shock-capturing magnetohydrodynamics codes
- On positivity preserving finite volume schemes for Euler equations
- A multiwave approximate Riemann solver for ideal MHD based on relaxation. II: Numerical implementation with 3 and 5 waves
- A method for incorporating Gauss' law into electromagnetic PIC codes
- An examination of forcing in direct numerical simulations of turbulence
- Hyperbolic divergence cleaning for the MHD equations
- An upwind differencing scheme for the equations of ideal magnetohydrodynamics
- The effect of nonzero \(\bigtriangledown\cdot B\) on the numerical solution of the magnetohydrodynamic equations
- An unsplit staggered mesh scheme for multidimensional magnetohydrodynamics
- A multiwave approximate Riemann solver for ideal MHD based on relaxation. I: Theoretical framework
- A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics
- An HLLC-Type Approximate Riemann Solver for Ideal Magnetohydrodynamics
- An HLLC Riemann solver for magneto-hydrodynamics
- Title not available (Why is that?)
- Numerical comparison of Riemann solvers for astrophysical hydrodynamics
- Stable upwind schemes for the magnetic induction equation
- A positive MUSCL-Hancock scheme for ideal magnetohydrodynamics
- Hydromagnetic flows from accretion discs and the production of radio jets
- Entropy satisfying flux vector splittings and kinetic BGK models
- Small-scale dynamo action during the formation of the first stars and galaxies
- Relaxation solvers for ideal MHD equations-a review
Cited In (23)
- A multi well-balanced scheme for the shallow water MHD system with topography
- Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes
- A new locally divergence-free path-conservative central-upwind scheme for ideal and shallow water magnetohydrodynamics
- Systematic construction of upwind constrained transport schemes for MHD
- A limit problem for three-dimensional ideal compressible radiation magneto-hydrodynamics
- Positivity-preserving finite difference WENO scheme for ten-moment equations with source term
- A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure
- Modelling and entropy satisfying relaxation scheme for the nonconservative bitemperature Euler system with transverse magnetic field
- A New Discretely Divergence-Free Positivity-Preserving High-Order Finite Volume Method for Ideal MHD Equations
- Provably Positive Central Discontinuous Galerkin Schemes via Geometric Quasilinearization for Ideal MHD Equations
- Geostrophic vs magneto-geostrophic adjustment and nonlinear magneto-inertia-gravity waves in rotating shallow water magnetohydrodynamics
- A Provably Positive Discontinuous Galerkin Method for Multidimensional Ideal Magnetohydrodynamics
- Locally divergence-free well-balanced path-conservative central-upwind schemes for rotating shallow water MHD
- Numerical analysis of an artificial compression method for magnetohydrodynamic flows at low magnetic Reynolds numbers
- Entropy stable finite volume approximations for ideal magnetohydrodynamics
- A 5-wave relaxation solver for the shallow water MHD system
- A positivity-preserving Lagrangian discontinuous Galerkin method for ideal magnetohydrodynamics equations in one-dimension
- Compact third-order accurate, positive preserving and divergence-free scheme
- The Developments and the Applications of the Numerical Algorithms in Simulating the Incompressible Magnetohydrodynamics with Complex Boundaries and Free Surfaces
- An Efficient Implementation of the Divergence Free Constraint in a Discontinuous Galerkin Method for Magnetohydrodynamics on Unstructured Meshes
- A high performance spectral code for nonlinear MHD stability
- A uniquely defined entropy stable matrix dissipation operator for high Mach number ideal MHD and compressible Euler simulations
- A semi-Lagrangian transport method for kinetic problems with application to dense-to-dilute polydisperse reacting spray flows
Uses Software
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)
