Cost-efficient finite-volume high-order schemes for compressible magnetohydrodynamics
Abstract
We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented up to tenth order of accuracy. The algorithmic architecture of this method is very close to that of commonly employed second-order schemes: it requires only one reconstructed value per face for each computational cell, independently of the scheme's order. This property is highly beneficial for the numerical efficiency. It results from reusing the required values already available in neighboring grid cells, in contrast to standard algorithms that require a number of reconstructions and evaluations which increases with the scheme's order of accuracy. At a given resolution, these high-order schemes present significantly less numerical dissipation than commonly employed lower-order approaches. Thus, results of comparable accuracy are achievable at a substantially coarser resolution, yielding overall performance gains. We also present a way to include physical dissipative terms: viscosity, magnetic diffusivity and cooling functions, respecting the finite-volume and constrained-transport frameworks. Benefits of this method are shown through applications in turbulent flows.
- Publication:
-
Journal of Computational Physics
- Pub Date:
- October 2024
- DOI:
- 10.1016/j.jcp.2024.113287
- arXiv:
- arXiv:2306.09856
- Bibcode:
- 2024JCoPh.51513287T
- Keywords:
-
- Partial differential equations;
- Finite-volume schemes;
- High-order methods;
- Fluid dynamics;
- Magnetohydrodynamics;
- Turbulence;
- Mathematics - Numerical Analysis;
- Physics - Computational Physics;
- Physics - Fluid Dynamics;
- Physics - Plasma Physics;
- Physics - Space Physics
- E-Print:
- Submitted to JCP. In version 1, the use of a global smoothness indicator: mass density+average of magnetic field components, suitable in supersonic flows, lead to difficulties in the convergence test (MHD vortex) and spurious oscillations for S6 and S8 in the hydro case. Version 2 uses individual smoothness indicators