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

Carl Dunlea,Ivan Khalzov
2024-04-26
Abstract:The DELiTE (Differential Equations on Linear Triangular Elements) framework was developed for spatial discretisation of partial differential equations on an unstructured triangular grid in axisymmetric geometry. The framework is based on discrete differential operators in matrix form, which are derived using linear finite elements and mimic some of the properties of their continuous counterparts. A single-fluid two-temperature MHD code is implemented in this framework. The inherent properties of the operators are used in the code to ensure global conservation of energy, particle count, toroidal flux, and angular momentum. The code was applied to study a novel experiment in which a compact torus (CT), produced with a magnetized Marshall gun, is magnetically levitated off an insulating wall and then magnetically compressed through the action of currents in the levitation/compression coils located outside the wall. We present numerical models for CT formation, levitation, and magnetic compression, and comparisons between simulated and experimental diagnostics.
Numerical Analysis,Plasma Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a globally conservative finite - element magnetohydrodynamics (MHD) code and apply it to the study of the formation, levitation, and magnetic compression of compact tori (CT). Specifically, the authors aim to ensure the global conservation of physical quantities such as mass, energy, toroidal flux, and angular momentum through numerical simulation, thereby improving the physical fidelity of the numerical model. ### Main problems: 1. **Developing a globally conservative numerical method**: To ensure the reliability of numerical solutions in complex physical systems, the authors developed a framework of discrete differential operators (DELiTE) based on the linear finite - element method to achieve global conservation of mass, energy, toroidal flux, and angular momentum. 2. **Application to compact torus experiments**: This code was used to simulate the formation, levitation, and magnetic compression processes of compact tori (CT). A compact torus is a plasma structure generated by a magnetized Marshall gun, which can levitate on an insulating wall and be magnetically compressed by an external coil current. The authors hope to verify the experimental results through numerical simulation and further understand these physical processes. ### Specific objectives: - **Ensure global conservation**: Design discrete differential operators so that the numerical solutions maintain global conservation of mass, energy, toroidal flux, and angular momentum within the computational domain. - **Simulate compact torus experiments**: Use the developed code to perform numerical simulations of the formation, levitation, and magnetic compression processes of compact tori and compare them with experimental data to verify the validity of the model. ### Key contributions: - **Design of discrete differential operators**: Proposed discrete differential operators based on the linear finite - element method, ensuring the global conservation property of numerical solutions. - **Numerical simulation of compact torus experiments**: Verified the formation, levitation, and magnetic compression processes of compact tori through numerical simulation and demonstrated good agreement with experimental data. In summary, the main purpose of this paper is to improve the physical fidelity of magnetohydrodynamics simulations by developing a globally conservative numerical method and applying it to the study of compact torus experiments.