Tokamak free-boundary plasma equilibrium computations in presence of non-linear materials

Cedric Boulbe (1, 2) , Blaise Faugeras (1, 2) , Guillaume Gros (1, 2) , Francesca Rapetti (1, 2)

1 LJAD - Laboratoire Jean Alexandre Dieudonné
2 CASTOR - Control, Analysis and Simulations for TOkamak Research

We consider the axisymmetric formulation of the equilibrium problem for a hot plasma in a tokamak. We adopt a non-overlapping mortar element approach, that couples C0 piece- wise linear Lagrange finite elements in a region that does not contain the plasma and C1 piece-wise cubic reduced Hsieh–Clough–Tocher finite elements elsewhere, to approximate the magnetic flux field on a triangular mesh of the poloidal tokamak section. The inclusion of ferromagnetic parts is simplified by assuming that they fit within the axisymmetric modeling and a new formulation of the Newton algorithm for the problem solution is stated, both in the static and quasi-static evolution cases.

Domaines: Modélisation et simulation

Cedric Boulbe, Blaise Faugeras, Guillaume Gros, Francesca Rapetti. Tokamak free-boundary plasma equilibrium computations in presence of non-linear materials. Journal of Scientific Computing, 2023, 96, ⟨10.1007/s10915-023-02265-8⟩. ⟨hal-03423469v3⟩