Skip to Main content Skip to Navigation

A finite element method with overlapping meshes for free-boundary axisymmetric plasma equilibria in realistic geometries

Holger Heumann 1, 2 Francesca Rapetti 1
2 CASTOR - Control, Analysis and Simulations for TOkamak Research
CRISAM - Inria Sophia Antipolis - Méditerranée , JAD - Laboratoire Jean Alexandre Dieudonné
Abstract : Existing finite element implementations for the computation of free-boundary axisymmetric plasma equilibria approximate the unknown poloidal flux function by standard {lowest order} continuous finite elements with discontinuous gradients. \hh{As a consequence,} the location of critical points of the poloidal flux, that are of paramount importance in tokamak engineering, is constrained to nodes of the mesh \hh{leading} to undesired jumps in transient problems. Moreover, recent numerical results for the self-consistent coupling of equilibrium with resistive diffusion and transport suggest the necessity of higher regularity when approximating the flux map. In this work we propose a mortar element method that employs two overlapping meshes. One mesh with Cartesian quadrilaterals covers the vacuum \rfA{chamber domain accessible by the plasma} and one mesh with triangles discretizes the region outside. The two meshes overlap in a narrow region. 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 \hh{outside this region}. The continuity of the numerical solution in the region of overlap is weakly enforced by a mortar-like \hh{mapping}.
Document type :
Complete list of metadata

Cited literature [58 references]  Display  Hide  Download
Contributor : Holger Heumann Connect in order to contact the contributor
Submitted on : Tuesday, November 15, 2016 - 6:32:57 PM
Last modification on : Thursday, January 20, 2022 - 5:32:37 PM
Long-term archiving on: : Thursday, March 16, 2017 - 1:38:28 PM


Files produced by the author(s)


  • HAL Id : hal-01322816, version 3


Holger Heumann, Francesca Rapetti. A finite element method with overlapping meshes for free-boundary axisymmetric plasma equilibria in realistic geometries. [Research Report] RR-8916, Inria Sophia Antipolis. 2016. ⟨hal-01322816v3⟩



Les métriques sont temporairement indisponibles