A novel structure preserving semi-implicit finite volume method for viscous and resistive magnetohydrodynamics
arXiv:2012.11218 · doi:10.1002/fld.5041
Abstract
In this work we introduce a novel semi-implicit structure-preserving finite-volume/finite-difference scheme for the viscous and resistive equations of magnetohydrodynamics (MHD) based on an appropriate 3-split of the governing PDE system, which is decomposed into a first convective subsystem, a second subsystem involving the coupling of the velocity field with the magnetic field and a third subsystem involving the pressure-velocity coupling. The nonlinear convective terms are discretized explicitly, while the remaining two subsystems accounting for the Alfven waves and the magneto-acoustic waves are treated implicitly. The final algorithm is at least formally constrained only by a mild CFL stability condition depending on the velocity field of the pure hydrodynamic convection. To preserve the divergence-free constraint of the magnetic field exactly at the discrete level, a proper set of overlapping dual meshes is employed. The resulting linear algebraic systems are shown to be symmetric and therefore can be solved by means of an efficient standard matrix-free conjugate gradient algorithm. One of the peculiarities of the presented algorithm is that the magnetic field is defined on the edges of the main grid, while the electric field is on the faces. The final scheme can be regarded as a novel shock-capturing, conservative and structure preserving semi-implicit scheme for the nonlinear viscous and resistive MHD equations. Several numerical tests are presented to show the main features of our novel solver: linear-stability in the sense of Lyapunov is verified at a prescribed constant equilibrium solution; a 2nd-order of convergence is numerically estimated; shock-capturing capabilities are proven against a standard set of stringent MHD shock-problems; accuracy and robustness are verified against a nontrivial set of 2- and 3-dimensional MHD problems.
arXiv comment44 pages, 22 figures, 2 tables
References in corpus (20)
- Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
- Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for One-Dimensional Flow
- Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows
- ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement
- A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows
- High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: viscous heat-conducting fluids and elastic solids
- An Asymptotic-Preserving all-speed scheme for the Euler and Navier-Stokes equations
- Space-time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting
- An Intercomparison Between Divergence-Cleaning and Staggered Mesh Formulations for Numerical Magnetohydrodynamics
- Multi-dimensional Numerical Scheme for Resistive Relativistic MHD
- High order ADER schemes for continuum mechanics
- High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics
- Space-time adaptive ADER-DG schemes for dissipative flows: compressible Navier-Stokes and resistive MHD equations
- Efficient high order accurate staggered semi-implicit discontinuous Galerkin methods for natural convection problems
- A divergence-free semi-implicit finite volume scheme for ideal, viscous and resistive magnetohydrodynamics
- A structure-preserving staggered semi-implicit finite volume scheme for continuum mechanics
- Rarefaction Shocks, Shock Errors and Low Order of Accuracy in ZEUS
- A projection hybrid high order finite volume/finite element method for incompressible turbulent flows
- POD-Galerkin reduced order methods for combined Navier-Stokes transport equations based on a hybrid FV-FE solver
- Design and analysis of ADER-type schemes for model advection-diffusion-reaction equations