Journal
JOURNAL OF COMPUTATIONAL PHYSICS
Volume 439, Issue -, Pages -Publisher
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2021.110388
Keywords
Tokamak; Equilibrium; Reduced Hsieh-Clough-Tocher finite element; Newton method; Inverse problem; Geometrical coefficients
Funding
- French National Research Agency grant SISTEM [ANR-19-CE46-0005-03]
- Euratom research and training program 2014-2018 [633053]
- Euratom research and training program 2019-2020 [633053]
Ask authors/readers for more resources
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.
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.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available