Abstract:Equilibriums in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad--Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad--Shafranov solver, since it is a fully nonlinear optimization problem that simultaneouly solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad--Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two non-trivial constraints, which correspond to the nonlinear finite element discretization of the Grad--Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. It is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

A new flux coordinates-based solver for fixed-boundary tokamak equilibrium with toroidal flow

Kinetic Equilibrium Reconstruction on EAST Tokamak

Dynamic Optimization Of Tokamak Plasmas Via Control Parameterization And The Time-Scaling Transformation

Poloidal Magnetic Flux Profile Control in Tokamaks Via Normalized Coprime Factorization Robust Control

Extremum-Seeking-Based Receding-Horizon Optimal Control of Plasma Current Profile in the DIII-D Tokamak

Closed-Loop Tracking Control Of Poloidal Magnetic Flux Profile In Tokamaks

Gyrokinetic Particle Simulation of Microturbulence for General Magnetic Geometry and Experimental Profiles

Solving the Grad–Shafranov equation using spectral elements for tokamak equilibrium with toroidal rotation

Free-boundary Plasma Equilibria with Toroidal Plasma Flows

FEQIS: A free-boundary equilibrium solver for integrated modeling of tokamak plasmas

Plasma Equilibrium Calculation in J-TEXT Tokamak

An adaptive Newton-based free-boundary Grad-Shafranov solver

Reconstruction of the equilibrium of the plasma in a Tokamak and identification of the current density profile in real time

FLIPEC, an Ideal MHD free-boundary axisymmetric equilibrium solver in the presence of macroscopic flows.

Vacuum Solution for Solov’ev’s Equilibrium Configuration in Tokamaks

An Efficient Technique for Magnetic Equilibrium Reconstruction with Q Profile Constraints and Its Application on the East Tokamak

ATEQ: Adaptive Toroidal Equilibrium code

Free-boundary, High-Beta Equilibrium in a Large Aspect Ratio Tokamak with Nearly Circular Plasma Boundary

Zonal Flows Driven by the Turbulent Energy Flux and the Turbulent Toroidal Reynolds Stress in a Magnetic Fusion Torus

Corrigendum: The residual zonal flow in tokamak plasmas rotating toroidally at sonic speed (2014 Nucl. Fusion 54 042002)

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices