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.

Sharp estimates, uniqueness and spikes condensation for superlinear free boundary problems arising in plasma physics

Sharp Non-Uniqueness for the 3D Hyperdissipative Navier-Stokes Equations: Beyond the Lions Exponent

Sharp Non-Uniqueness of Weak Solutions to 3D Magnetohydrodynamic Equations: Beyond the Lions Exponent

Stability Analysis of Two-Dimensional Ideal Flows With Applications to Viscous Fluids and Plasmas

Sharp boundary concentration for a two-dimensional nonlinear Neumann problem

Desingularization of 2D elliptic free-boundary problem with non-autonomous nonlinearity

Sharp Nonuniqueness of Solutions to Stochastic Navier–Stokes Equations

Sharp non-uniqueness for the 2D hyper-dissipative Navier-Stokes equations

Asymptotic stability of boundary layer to the multi-dimensional isentropic Euler-Poisson equations arising in plasma physics

Sharp Global Well-posedness and Scattering of the Boltzmann Equation

Logarithmic Finite-Size Correction in Non-neutral Two-Component Plasma on Sphere

Global Existence and Scattering of the Klein–Gordon–Zakharov System in Two Space Dimensions

Numerical solution of the quantum Lenard-Balescu equation for a one-component plasma

From Two-Fluid Euler-Poisson Equations to One-Fluid Euler Equations

Sharp well-posedness for the cubic NLS and mKdV in

Sharp uniform-in-time mean-field convergence for singular periodic Riesz flows

Global Solutions of the Compressible Euler-Poisson Equations for Plasma with Doping Profile for Large Initial Data of Spherical Symmetry

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

Time-Global Regularity of the Navier-Stokes System with Hyper-Dissipation--Turbulent Scenario

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

Sharp criterion of global existence for a class of nonlinear Schrödinger equations with critical exponent.