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

Article · Zakariae Jorti, Qi Tang, Konstantin Lipnikov, Xian-Zhu Tang

- Archive id: arxiv:2303.08337
- Published: 2023-03-15
- Venue: SSRN Electronic Journal
- DOI: https://doi.org/10.2139/ssrn.4408634
- arXiv: https://arxiv.org/abs/2303.08337
- Citations: 0
- Source URL: https://arxiv.org/abs/2303.08337
- PDF: https://arxiv.org/pdf/2303.08337
- Provenance: arxiv

## 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.

_Abstract reproduced under CC0-1.0 (arXiv descriptive metadata)._

Topics: [Magnetic Confinement Fusion](https://space-kitty.com/topics/magnetic-fusion)

---
Space Kitty Archive · https://space-kitty.com · API: https://space-kitty.com/api/v1 · MCP: https://space-kitty.com/mcp · Premium: $0.20/request via x402 (Base Sepolia (testnet)). Metadata: arXiv (CC0), OpenAlex (CC0), Europe PMC (CC BY/CC0 works), ClinicalTrials.gov.