Construction of an invertible mapping to boundary conforming coordinates for arbitrarily shaped toroidal domains

Abstract

Boundary conforming coordinates are commonly used in plasma physics to describe the geometry of toroidal domains, for example, in three-dimensional magnetohydrodynamic equilibrium solvers. The magnetohydrodynamic equilibrium configuration can be approximated with an inverse map, defining nested surfaces of constant magnetic flux. For equilibrium solvers that solve for this inverse map iteratively, the initial guess for the inverse map must be well defined and invertible. Even if magnetic islands are to be included in the representation, boundary conforming coordinates can still be useful, for example to parametrize the interface surfaces in multi-region, relaxed magnetohydrodynamics or as a general-purpose, field-agnostic coordinate system in strongly shaped domains. Given a fixed boundary shape, finding a valid boundary conforming mapping can be challenging, especially for the non-convex boundaries from recent developments in stellarator optimization. In this work, we propose a new algorithm to construct such a mapping, by solving two Dirichlet-Laplace problems via a boundary integral method. We can prove that the generated harmonic map is always smooth and has a smooth inverse. Furthermore, we can find a discrete approximation of the mapping that preserves this property.

Presented at the 2024 Joint Varenna-Lausanne International Workshop on the Theory of Fusion Plasmas

References in corpus (5)

• The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas
• Computation of multi-region relaxed magnetohydrodynamic equilibria
• Mapping the space of quasisymmetric stellarators using optimized near-axis expansion
• Approach to nonlinear magnetohydrodynamic simulations in stellarator geometry
• Optimal sampling patterns for Zernike polynomials