Taylor States in Stellarators: A Fast High-order Boundary Integral Solver

Abstract

We present a boundary integral equation solver for computing Taylor relaxed states in non-axisymmetric solid and shell-like toroidal geometries. The computation of Taylor states in these geometries is a key element for the calculation of stepped pressure stellarator equilibria. The integral representation of the magnetic field in this work is based on the generalized Debye source formulation, and results in a well-conditioned second-kind boundary integral equation. The integral equation solver is based on a spectral discretization of the geometry and unknowns, and the computation of the associated weakly-singular integrals is performed with high-order quadrature based on a partition of unity. The resulting scheme for applying the integral operator is then coupled with an iterative solver and suitable preconditioners. Several numerical examples are provided to demonstrate the accuracy and efficiency of our method, and a direct comparison with the leading code in the field is reported.

References in corpus (3)

• Pseudo-spectral methods for the Laplace-Beltrami equation and the Hodge decomposition on surfaces of genus one
• Exact axisymmetric Taylor states for shaped plasmas
• Second-kind integral equations for the Laplace-Beltrami problem on surfaces in three dimensions

Cited by in corpus (10)

• Combined plasma-coil optimization algorithms
• Efficient high-order singular quadrature schemes in magnetic fusion
• A robust solver for elliptic PDEs in 3D complex geometries
• Fast multipole methods for evaluation of layer potentials with locally-corrected quadratures
• On the non-existence of stepped-pressure equilibria far from symmetry
• Adjoint methods for quasisymmetry of vacuum fields on a surface
• Efficient Convergent Boundary Integral Methods for Slender Bodies
• Debye source representation for type-I superconductors, I
• A fast, accurate, and easy to implement Kapur-Rokhlin quadrature scheme for singular integrals in axisymmetric geometries
• A parametrix method for elliptic surface PDEs