2D Fourier finite element formulation for magnetostatics in curvilinear coordinates with a symmetry direction

Albert, Christopher G.; Bíró, Oszkár; Lainer, Patrick

doi:10.1016/j.cpc.2022.108401

Cited by 1 publication

(2 citation statements)

References 14 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…On this mesh, we obtain 𝛿𝑩 p from 𝛿𝑱 via the finite element method presented in [17], thus realizing the operator M in (1). As noted in [15], we extend the computational domain with a current-free region outside the separatrix, and impose a homogenenous Dirichlet boundary condition on this extended domain 𝛺.…”

Section: Geometrical Considerationsmentioning

confidence: 99%

“…Here, the implementation of ideal MHD plasma response currents is presented that includes localized shielding currents near resonant surfaces. This permits the self-consistent solution of the linear ideal MHD plasma response to resonant 3D perturbations with iterative coupling to a 3D field solver for Ampère's law [17]. The resulting code MEPHIT may act as a drop-in replacement for kinetics in ideal regions, as an alternative geometrical correction for KiLCA, and for benchmarking in the ongoing effort on 3D kinetic equilibria.…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Iterative approach to linear ideal MHD modeling of plasma response to 3D magnetic perturbations in tokamaks

Lainer

Kasilov

Kernbichler

et al. 2022

J. Phys.: Conf. Ser.

View full text Add to dashboard Cite

The class of plasma instabilities known as edge-localized modes (ELMs) is of special concern in tokamaks operating in high-confinement mode, such as ASDEX Upgrade and ITER. One strategy for ELM mitigation is the application of resonant magnetic perturbations (RMPs) via external coils. Kinetic modeling accurately describes the plasma response to these RMPs ab initio, particularly the parallel shielding currents at resonant surfaces. Away from resonant surfaces, ideal magnetohydrodynamics (iMHD) is expected to yield sufficiently accurate results, providing a computationally less expensive option that could complement kinetic modeling. The code MEPHIT has been developed to solve the linearized iMHD equations in a way that is compatible with iterative kinetic modeling approaches. We consider an axisymmetric iMHD equilibrium in realistic tokamak geometry under the influence of a quasi-static non-axisymmetric external perturbation from ELM mitigation coils. The plasma responds to this external magnetic perturbation with a current perturbation, which in turn produces a magnetic field perturbation. The resulting fixed-point equation can be solved in a self-consistent manner by preconditioned iterations in which Ampère’s equation and the magnetic differential equations for pressure and current are solved in alternation until convergence is reached. After expansion in toroidal Fourier harmonics, these equations are solved on a triangular mesh in the poloidal plane using finite elements. These results are then benchmarked against established codes.

show abstract

Section: Geometrical Considerationsmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%