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.

A Fully Implicit Parallel Solver for MHD Instabilities in a Tokamak

A mimetic finite difference based quasi-static magnetohydrodynamic solver for force-free plasmas in tokamak disruptions

A high-order finite-difference solver for direct numerical simulations of magnetohydrodynamic turbulence

An Adaptive Moving Mesh Finite Difference Scheme for Tokamak Magneto-Hydrodynamic Simulations

High-order implicit solver in conservative formulation for tokamak plasma transport equations

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

Numerical simulation of magnetohydrodynamic instability with variational principle in the tokamak

Parallel Finite Volume Simulation of the Spherical Shell Dynamo with Pseudo-Vacuum Magnetic Boundary Conditions

Nonlinear extended magnetohydrodynamics simulation using high-order finite elements

Development of a tokamak magnetohydrodynamic code with the discontinuous Galerkin and Weighted Essentially Non-Oscillatory method

On scalable liquid-metal MHD solvers for fusion breeder blanket multiphysics applications

BOUT++: A framework for parallel plasma fluid simulations

A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations

A structure-preserving semi-implicit IMEX finite volume scheme for ideal magnetohydrodynamics at all Mach and Alfvén numbers

Fully decoupled, linearized and stabilized finite volume method for the time-dependent incompressible MHD equations

Development of a gyrokinetic-MHD energetic particle simulation code Part I: MHD version

A new MHD code with adaptive mesh refinement and parallelization for astrophysics

A new flux coordinates-based solver for fixed-boundary tokamak equilibrium with toroidal flow

ATEQ: Adaptive Toroidal Equilibrium code

Implicit Solar Coronal Magnetohydrodynamic (MHD) Modeling with a Low-dissipation Hybridized AUSM-HLL Riemann Solver

Piecewise Field-Aligned Finite Element Method for Multi-Mode Nonlinear Particle Simulations in tokamak plasmas