Abstract:Force-free plasmas are a good approximation where the plasma pressure is tiny compared with the magnetic pressure, which is the case during the cold vertical displacement event (VDE) of a major disruption in a tokamak. On time scales long compared with the transit time of Alfven waves, the evolution of a force-free plasma is most efficiently described by the quasi-static magnetohydrodynamic (MHD) model, which ignores the plasma inertia. Here we consider a regularized quasi-static MHD model for force-free plasmas in tokamak disruptions and propose a mimetic finite difference (MFD) algorithm. The full geometry of an ITER-like tokamak reactor is treated, with a blanket module region, a vacuum vessel region, and the plasma region. Specifically, we develop a parallel, fully implicit, and scalable MFD solver based on PETSc and its DMStag data structure for the discretization of the five-field quasi-static perpendicular plasma dynamics model on a 3D structured mesh. The MFD spatial discretization is coupled with a fully implicit DIRK scheme. The algorithm exactly preserves the divergence-free condition of the magnetic field under the resistive Ohm's law. The preconditioner employed is a four-level fieldsplit preconditioner, which is created by combining separate preconditioners for individual fields, that calls multigrid or direct solvers for sub-blocks or exact factorization on the separate fields. The numerical results confirm the divergence-free constraint is strongly satisfied and demonstrate the performance of the fieldsplit preconditioner and overall algorithm. The simulation of ITER VDE cases over the actual plasma current diffusion time is also presented.

Flux-aligned Quad Mesh Generation in Magnetohydrodynamic Simulation.

Magnetic mesh generation and field line reconstruction for scrape-off layer and divertor modeling in stellarators

A 2-D Axisymmetric Magneto-Hydrodynamic Model of a DC Arc Plasma Torch and Its Solution Methodology

Thermal Mixing Enhancement of Liquid Metal MHD Free-Surface Flow by Optimizing Vortex Generator Arrays

Meshless Methods for Magnetohydrodynamics with Vector Potential

Unstructured mesh tools for magnetically confined fusion system simulations

Advances in Direct Numerical Simulation for MHD Modeling of Free Surface Flows

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

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

Asymptotic quasisymmetric high-beta three-dimensional magnetohydrodynamic equilibria near axisymmetry

Lattice Boltzmann Methods for Magnetohydrodynamic Flows in Fusion Applications

Two-fluid and magnetohydrodynamic modelling of magnetic reconnection in the MAST spherical tokamak and the solar corona

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

Nonlinear extended magnetohydrodynamics simulation using high-order finite elements

An Adaptive Moving Mesh Method for Two-Dimensional Ideal Magnetohydrodynamics

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

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

Mesh generation for confined fusion plasma simulation

A High-Order Gas-Kinetic Method for Multidimensional Ideal Magnetohydrodynamics

Numerical simulation of turbulent flow in FlowCube platform under a strong magnetic field

Simplicial Approximation of Deforming 3D Spaces for Visualizing Fusion Plasma Simulation Data