Abstract:Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with nonlinear numerical instabilities. One approach for controlling the discrete divergence is through a so-called constrained transport method, which is based on first predicting a magnetic field through a standard finite volume solver, and then correcting this field through the appropriate use of a magnetic vector potential. In this work we develop a constrained transport method for the 3D ideal MHD equations that is based on a high-resolution wave propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme developed by Rossmanith [SIAM J. Sci. Comp. 28, 1766 (2006)], and is based on the high-resolution wave propagation method of Langseth and LeVeque [J. Comp. Phys. 165, 126 (2000)]. In particular, in our extension we take great care to maintain the three most important properties of the 2D scheme: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as cell-centered; (2) we develop a high-resolution wave propagation scheme for evolving the magnetic potential; and (3) we develop a wave limiting approach that is applied during the vector potential evolution, which controls unphysical oscillations in the magnetic field. One of the key numerical difficulties that is novel to 3D is that the transport equation that must be solved for the magnetic vector potential is only weakly hyperbolic. In presenting our numerical algorithm we describe how to numerically handle this problem of weak hyperbolicity, as well as how to choose an appropriate gauge condition. The resulting scheme is applied to several numerical test cases.

Regions without flux surfaces of given class for magnetic fields in toroidal geometry

Electromagnetic Formulation of Global Gyrokinetic Particle Simulation in Toroidal Geometry

Closed-Loop Tracking Control Of Poloidal Magnetic Flux Profile In Tokamaks

Confining and escaping magnetic field lines in Tokamaks: Analysis via symplectic map

A non-twisting flux tube for local gyrokinetic simulations

Construction of an invertible mapping to boundary conforming coordinates for arbitrarily shaped toroidal domains

An Unstaggered Constrained Transport Method for the 3D Ideal Magnetohydrodynamic Equations

Generalised action-angle coordinates defined on island chains

Three Dimensional Equilibrium Solutions for A Current-Carrying Reversed-Field Pinch Plasma with A Close-Fitting Conducting Shell

A Three-Dimensional Magnetohydrodynamic Equilibrium In An Axial Coordinate With A Constant Curvature

Flux-aligned Quad Mesh Generation in Magnetohydrodynamic Simulation.

Steady Compressible 3D Euler Flows in Toroidal Volumes without Continuous Euclidean Isometries

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

Discretization form of the continuity condition at the polar axis, with application to the gyrokinetic simulation in a magnetic fusion torus

A Reduced Ideal MHD System for Nonlinear Magnetic Field Turbulence in Plasmas with Approximate Flux Surfaces

Numerical study of $δ$-function current sheets arising from resonant magnetic perturbations

Magnetohydrodynamic "cat eyes" and stabilizing effects of plasma flow

A globally conservative finite element MHD code and its application to the study of compact torus formation, levitation and magnetic compression

Zonal Flows Driven by the Turbulent Energy Flux and the Turbulent Toroidal Reynolds Stress in a Magnetic Fusion Torus

Kinetic vs magnetic chaos in toroidal plasmas: A systematic quantitative comparison

Control of Magnetic Islands in the STOR-M Tokamak Using Resonant Helical Fields