A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations
arXiv:2305.06497 · doi:10.1016/j.jcp.2023.112493
Abstract
We present a new divergence-free and well-balanced hybrid FV/FE scheme for the incompressible viscous and resistive MHD equations on unstructured mixed-element meshes in 2 and 3 space dimensions. The equations are split into subsystems. The pressure is defined on the vertices of the primary mesh, while the velocity field and the normal components of the magnetic field are defined on an edge-based/face-based dual mesh in two and three space dimensions, respectively. This allows to account for the divergence-free conditions of the velocity field and of the magnetic field in a rather natural manner. The non-linear convective and the viscous terms are solved at the aid of an explicit FV scheme, while the magnetic field is evolved in a divergence-free manner via an explicit FV method based on a discrete form of the Stokes law in the edges/faces of each primary element. To achieve higher order of accuracy, a pw-linear polynomial is reconstructed for the magnetic field, which is guaranteed to be divergence-free via a constrained L2 projection. The pressure subsystem is solved implicitly at the aid of a classical continuous FE method in the vertices of the primary mesh. In order to maintain non-trivial stationary equilibrium solutions of the governing PDE system exactly, which are assumed to be known a priori, each step of the new algorithm takes the known equilibrium solution explicitly into account so that the method becomes exactly well-balanced. This paper includes a very thorough study of the lid-driven MHD cavity problem in the presence of different magnetic fields. We finally present long-time simulations of Soloviev equilibrium solutions in several simplified 3D tokamak configurations even on very coarse unstructured meshes that, in general, do not need to be aligned with the magnetic field lines.
arXiv comment57 pages, 33 figures, 13 tables, reference-data (supplementary electronic material) will be available after publication on the Journal web-page
References in corpus (40)
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport
- Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
- Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
- Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for One-Dimensional Flow
- Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows
- The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas
- A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows
- Space-time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting
- A pressure-based semi-implicit space-time discontinuous Galerkin method on staggered unstructured meshes for the solution of the compressible Navier-Stokes equations at all Mach numbers
- Central WENO schemes for hyperbolic conservation laws on fixed and moving unstructured meshes
- A staggered space-time discontinuous Galerkin method for the three-dimensional incompressible Navier-Stokes equations on unstructured tetrahedral meshes
- An Unstaggered Constrained Transport Method for the 3D Ideal Magnetohydrodynamic Equations
- Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations
- High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics
- High order pressure-based semi-implicit IMEX schemes for the 3D Navier-Stokes equations at all Mach numbers
- A linearized energy preserving finite element method for the dynamical incompressible magnetohydrodynamics equations
- A divergence-free semi-implicit finite volume scheme for ideal, viscous and resistive magnetohydrodynamics
- Robust Preconditioners for Incompressible MHD Models
- Arbitrary order finite volume well-balanced schemes for the Euler equations with gravity
- Well balanced Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming meshes for the Euler equations of gasdynamics with gravity
- Rarefaction Shocks, Shock Errors and Low Order of Accuracy in ZEUS
- Direct Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming unstructured meshes
- An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity
- A projection hybrid high order finite volume/finite element method for incompressible turbulent flows
- A second order positivity preserving well-balanced finite volume scheme for Euler equations with gravity for arbitrary hydrostatic equilibria
- A semi-implicit hybrid finite volume / finite element scheme for all Mach number flows on staggered unstructured meshes
- Collocation Methods for High-Order Well-Balanced Methods for Systems of Balance Laws
- High order entropy preserving ADER scheme
- A staggered semi-implicit hybrid finite volume / finite element scheme for the shallow water equations at all Froude numbers
- Design and analysis of ADER-type schemes for model advection-diffusion-reaction equations
- An Entropy Stable Nodal Discontinuous Galerkin Method for the resistive MHD Equations. Part II: Subcell Finite Volume Shock Capturing
- High-order well-balanced methods for systems of balance laws: a control-based approach
- Helicity-conservative finite element discretization for incompressible MHD systems
- An Arbitrary-Lagrangian-Eulerian hybrid finite volume/finite element method on moving unstructured meshes for the Navier-Stokes equations
- Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers
- A new thermodynamically compatible finite volume scheme for magnetohydrodynamics
- A structure-preserving finite element method for compressible ideal and resistive MHD
- A novel flexible field-aligned coordinate system for tokamak edge plasma simulation
- Well balanced finite volume schemes for shallow water equations on manifolds
- A locally field-aligned discontinuous Galerkin method for anisotropic wave equations