Taylor states in stellarators: A fast high-order boundary integral solver

Dhairya Malhotra, Antoine Cerfon, Lise Marie Imbert-Gérard, Michael O'Neil

Research output: Contribution to journalArticle

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.

Original languageEnglish (US)
Article number108791
JournalJournal of Computational Physics
Volume397
DOIs
StatePublished - Nov 15 2019

Fingerprint

stellarators
Boundary Integral
integral equations
Boundary integral equations
Boundary Integral Equations
Higher Order
Geometry
toroidal shells
geometry
Iterative Solver
Partition of Unity
Singular Integrals
quadratures
Integral Operator
Quadrature
Preconditioner
Integral Representation
Integral equations
unity
partitions

Keywords

  • Generalized Debye sources
  • Laplace-Beltrami
  • Plasma equilibria
  • Stellarator
  • Taylor state

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Cite this

Taylor states in stellarators : A fast high-order boundary integral solver. / Malhotra, Dhairya; Cerfon, Antoine; Imbert-Gérard, Lise Marie; O'Neil, Michael.

In: Journal of Computational Physics, Vol. 397, 108791, 15.11.2019.

Research output: Contribution to journalArticle

@article{62ec375de39f488a82818abe1496ee09,
title = "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.",
keywords = "Generalized Debye sources, Laplace-Beltrami, Plasma equilibria, Stellarator, Taylor state",
author = "Dhairya Malhotra and Antoine Cerfon and Imbert-G{\'e}rard, {Lise Marie} and Michael O'Neil",
year = "2019",
month = "11",
day = "15",
doi = "10.1016/j.jcp.2019.06.067",
language = "English (US)",
volume = "397",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Taylor states in stellarators

T2 - A fast high-order boundary integral solver

AU - Malhotra, Dhairya

AU - Cerfon, Antoine

AU - Imbert-Gérard, Lise Marie

AU - O'Neil, Michael

PY - 2019/11/15

Y1 - 2019/11/15

N2 - 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.

AB - 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.

KW - Generalized Debye sources

KW - Laplace-Beltrami

KW - Plasma equilibria

KW - Stellarator

KW - Taylor state

UR - http://www.scopus.com/inward/record.url?scp=85070217157&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85070217157&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2019.06.067

DO - 10.1016/j.jcp.2019.06.067

M3 - Article

AN - SCOPUS:85070217157

VL - 397

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

M1 - 108791

ER -