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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:439

Article

Elarif, Ali¹  Faugeras, Blaise¹  Rapetti, Francesca¹ 

[1] Univ Cote Azur, INRIA, CNRS, LJAD, Parc Valrose, F-06108 Nice 02, France

关键词: Tokamak;    Equilibrium;    Reduced Hsieh-Clough-Tocher finite element;    Newton method;    Inverse problem;    Geometrical coefficients;   

DOI  :  10.1016/j.jcp.2021.110388

来源: Elsevier

【 摘 要 】

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.

【 授权许可】

Free