Tokamak free-boundary plasma equilibrium computation using finite elements of class C0 and C1 within a mortar element approach

The numerical simulation of the equilibrium of the plasma in a tokamak as well as its self-consistent coupling with resistive diffusion should benefit from higher regularity of the approximation of the magnetic flux map. In this work, we propose a finite element approach on a triangular mesh of the poloidal section, that couples piece-wise linear finite elements in a region that does not contain the plasma and reduced Hsieh-Clough-Tocher finite elements elsewhere. This approach gives the flexibility to achieve easily and at low cost higher order regularity for the approximation of the flux function in the domain covered by the plasma, while preserving accurate meshing of the geometric details in the rest of the computational domain. The continuity of the numerical solution at the coupling interface is weakly enforced by mortar projection. A new technique for the computation of the geometrical coefficients is also presented. (C) 2021 Elsevier Inc. All rights reserved.

Automated summary

This study presents a numerical simulation technique that combines finite element method and triangular mesh to improve the accuracy of the equilibrium of the plasma in a tokamak and its coupling with resistive diffusion. The method achieves higher order regularity in the plasma-covered area while maintaining accuracy in meshing the geometric details in the rest of the computational domain.