skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

This content will become publicly available on January 13, 2021

Title: Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria. Part 2. Applications

Abstract

We report that the shape gradient is a local sensitivity function defined on the surface of an object which provides the change in a characteristic quantity, or figure of merit, associated with a perturbation to the shape of the object. The shape gradient can be used for gradient-based optimization, sensitivity analysis and tolerance calculations. However, it is generally expensive to compute from finite-difference derivatives for shapes that are described by many parameters, as is the case for typical stellarator geometry. In an accompanying work (Antonsen, Paul & Landreman  J. Plasma Phys., vol. 85 (2), 2019), generalized self-adjointness relations are obtained for magnetohydrodynamic (MHD) equilibria. These describe the relation between perturbed equilibria due to changes in the rotational transform or toroidal current profiles, displacements of the plasma boundary, modifications of currents in the vacuum region or the addition of bulk forces. These are applied to efficiently compute the shape gradient of functions of MHD equilibria with an adjoint approach. In this way, the shape derivative with respect to any perturbation applied to the plasma boundary or coil shapes can be computed with only one additional MHD equilibrium solution. We demonstrate that this approach is applicable for several figures of merit of interest for stellarator configuration optimization: the magnetic well, the magnetic ripple on axis, the departure from quasisymmetry, the effective ripple in the low-collisionality$$1/\unicode[STIX]{x1D708}$$regime$$(\unicode[STIX]{x1D716}_{\text{eff}}^{3/2})$$(Nemov et al. Phys. Plasmas, vol. 6 (12), 1999, pp. 4622–4632) and several finite-collisionality neoclassical quantities. Numerical verification of this method is demonstrated for the magnetic well figure of merit with the VMEC code (Hirshman & WhitsonPhys. Fluids, vol. 26 (12), 1983, p. 3553) and for the magnetic ripple with modification of the ANIMEC code (Cooperet al. Comput. Phys. Commun., vol. 72 (1), 1992, pp. 1–13). Finally, comparisons with the direct approach demonstrate that, in order to obtain agreement within several per cent, the adjoint approach provides a factor of$$O(10^{3})$$in computational savings.

Authors:
ORCiD logo [1];  [1]; ORCiD logo [1];  [2]
  1. Univ. of Maryland, College Park, MD (United States). Inst. for Research in Electronics and Applied Physics
  2. Swiss Alps Fusion Energy (SAFE) (Switzerland)
Publication Date:
Research Org.:
Univ. of Maryland, College Park, MD (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Fusion Energy Sciences (FES) (SC-24)
OSTI Identifier:
1597686
Grant/Contract Number:  
[FG02-93ER54197; FC02-08ER54964]
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Plasma Physics
Additional Journal Information:
[ Journal Volume: 86; Journal Issue: 1]; Journal ID: ISSN 0022-3778
Publisher:
Cambridge University Press
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Citation Formats

Paul, Elizabeth J., Antonsen, Thomas, Landreman, Matt, and Cooper, W. Anthony. Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria. Part 2. Applications. United States: N. p., 2020. Web. doi:10.1017/S0022377819000916.
Paul, Elizabeth J., Antonsen, Thomas, Landreman, Matt, & Cooper, W. Anthony. Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria. Part 2. Applications. United States. doi:10.1017/S0022377819000916.
Paul, Elizabeth J., Antonsen, Thomas, Landreman, Matt, and Cooper, W. Anthony. Mon . "Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria. Part 2. Applications". United States. doi:10.1017/S0022377819000916.
@article{osti_1597686,
title = {Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria. Part 2. Applications},
author = {Paul, Elizabeth J. and Antonsen, Thomas and Landreman, Matt and Cooper, W. Anthony},
abstractNote = {We report that the shape gradient is a local sensitivity function defined on the surface of an object which provides the change in a characteristic quantity, or figure of merit, associated with a perturbation to the shape of the object. The shape gradient can be used for gradient-based optimization, sensitivity analysis and tolerance calculations. However, it is generally expensive to compute from finite-difference derivatives for shapes that are described by many parameters, as is the case for typical stellarator geometry. In an accompanying work (Antonsen, Paul & Landreman J. Plasma Phys., vol. 85 (2), 2019), generalized self-adjointness relations are obtained for magnetohydrodynamic (MHD) equilibria. These describe the relation between perturbed equilibria due to changes in the rotational transform or toroidal current profiles, displacements of the plasma boundary, modifications of currents in the vacuum region or the addition of bulk forces. These are applied to efficiently compute the shape gradient of functions of MHD equilibria with an adjoint approach. In this way, the shape derivative with respect to any perturbation applied to the plasma boundary or coil shapes can be computed with only one additional MHD equilibrium solution. We demonstrate that this approach is applicable for several figures of merit of interest for stellarator configuration optimization: the magnetic well, the magnetic ripple on axis, the departure from quasisymmetry, the effective ripple in the low-collisionality$1/\unicode[STIX]{x1D708}$regime$(\unicode[STIX]{x1D716}_{\text{eff}}^{3/2})$(Nemovet al. Phys. Plasmas, vol. 6 (12), 1999, pp. 4622–4632) and several finite-collisionality neoclassical quantities. Numerical verification of this method is demonstrated for the magnetic well figure of merit with the VMEC code (Hirshman & WhitsonPhys. Fluids, vol. 26 (12), 1983, p. 3553) and for the magnetic ripple with modification of the ANIMEC code (Cooperet al. Comput. Phys. Commun., vol. 72 (1), 1992, pp. 1–13). Finally, comparisons with the direct approach demonstrate that, in order to obtain agreement within several per cent, the adjoint approach provides a factor of$O(10^{3})$in computational savings.},
doi = {10.1017/S0022377819000916},
journal = {Journal of Plasma Physics},
number = [1],
volume = [86],
place = {United States},
year = {2020},
month = {1}
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on January 13, 2021
Publisher's Version of Record

Save / Share:

Works referenced in this record:

Anisotropic pressure bi-Maxwellian distribution function model for three-dimensional equilibria
journal, May 2006


ESTELL: A Quasi-Toroidally Symmetric Stellarator
journal, March 2013

  • Drevlak, M.; Brochard, F.; Helander, P.
  • Contributions to Plasma Physics, Vol. 53, Issue 6
  • DOI: 10.1002/ctpp.201200055

Electrostatic modification of variational principles for anisotropic plasmas
journal, January 1982


Magnetic islands and perturbed plasma equilibria
journal, July 2003

  • Nührenberg, Carolin; Boozer, Allen H.
  • Physics of Plasmas, Vol. 10, Issue 7
  • DOI: 10.1063/1.1578489

An adjoint method for neoclassical stellarator optimization
journal, August 2019

  • Paul, Elizabeth J.; Abel, Ian G.; Landreman, Matt
  • Journal of Plasma Physics, Vol. 85, Issue 5
  • DOI: 10.1017/S0022377819000527

The ∇B drift velocity of trapped particles in stellarators
journal, November 2005

  • Nemov, V. V.; Kasilov, S. V.; Kernbichler, W.
  • Physics of Plasmas, Vol. 12, Issue 11
  • DOI: 10.1063/1.2131848

3D magnetohydrodynamic equilibria with anisotropic pressure
journal, October 1992


Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria
journal, March 2019

  • Antonsen, Thomas; Paul, Elizabeth J.; Landreman, Matt
  • Journal of Plasma Physics, Vol. 85, Issue 2
  • DOI: 10.1017/S0022377819000254

Physics of the compact advanced stellarator NCSX
journal, November 2001


Neoclassical diffusion in an l = 3 stellarator
journal, January 1974


Methods for the Efficient Calculation of the (Mhd) Magnetohydrodynamic Stability Properties of Magnetically Confined Fusion Plasmas
journal, September 1990

  • Anderson, D. V.; Cooper, W. A.; Gruber, R.
  • The International Journal of Supercomputing Applications, Vol. 4, Issue 3
  • DOI: 10.1177/109434209000400305

Quasi-helically symmetric toroidal stellarators
journal, May 1988


Plasma equilibrium with rational magnetic surfaces
journal, January 1981


Three-dimensional anisotropic pressure equilibria that model balanced tangential neutral beam injection effects
journal, February 2005

  • Cooper, W. A.; Hirshman, S. P.; Yamaguchi, T.
  • Plasma Physics and Controlled Fusion, Vol. 47, Issue 3
  • DOI: 10.1088/0741-3335/47/3/011

Perturbed plasma equilibria
journal, October 2006

  • Boozer, Allen H.; Nührenberg, Carolin
  • Physics of Plasmas, Vol. 13, Issue 10
  • DOI: 10.1063/1.2353903

Quasi-helical symmetry in stellarators
journal, November 1995


Theory of plasma confinement in non-axisymmetric magnetic fields
journal, July 2014


Neoclassical transport in stellarators
journal, January 1987

  • Ho, Darwin D. -M.; Kulsrud, Russell M.
  • Physics of Fluids, Vol. 30, Issue 2
  • DOI: 10.1063/1.866395

Physics of compact stellarators
journal, May 1999

  • Hirshman, S. P.; Spong, D. A.; Whitson, J. C.
  • Physics of Plasmas, Vol. 6, Issue 5
  • DOI: 10.1063/1.873489

Steepest-descent moment method for three-dimensional magnetohydrodynamic equilibria
journal, January 1983


Properties of a new quasi-axisymmetric configuration
journal, January 2019


Differentiating the shape of stellarator coils with respect to the plasma boundary
journal, September 2018


Comparison of particle trajectories and collision operators for collisional transport in nonaxisymmetric plasmas
journal, April 2014

  • Landreman, M.; Smith, H. M.; Mollén, A.
  • Physics of Plasmas, Vol. 21, Issue 4
  • DOI: 10.1063/1.4870077

Ideal magnetohydrodynamics: Global mode analysis of three‐dimensional plasma configurations
journal, September 1993

  • Schwab, Carolin
  • Physics of Fluids B: Plasma Physics, Vol. 5, Issue 9
  • DOI: 10.1063/1.860656

Physics Design for ARIES-CS
journal, October 2008

  • Ku, L. P.; Garabedian, P. R.; Lyon, J.
  • Fusion Science and Technology, Vol. 54, Issue 3
  • DOI: 10.13182/FST08-A1899

Computing local sensitivity and tolerances for stellarator physics properties using shape gradients
journal, June 2018


The virtual-casing principle for 3D toroidal systems
journal, November 2012


Evaluation of 1/ν neoclassical transport in stellarators
journal, December 1999

  • Nemov, V. V.; Kasilov, S. V.; Kernbichler, W.
  • Physics of Plasmas, Vol. 6, Issue 12
  • DOI: 10.1063/1.873749

Physics design of a high-bbeta quasi-axisymmetric stellarator
journal, December 1999


Fast particle confinement with optimized coil currents in the W7-X stellarator
journal, April 2014


Configuration Property of the Chinese First Quasi-Axisymmetric Stellarator
journal, January 2018

  • Shimizu, Akihiro; Liu, Haifeng; Isobe, Mitsutaka
  • Plasma and Fusion Research, Vol. 13, Issue 0
  • DOI: 10.1585/pfr.13.3403123

Collisional Diffusion in Nonaxisymmetric Toroidal Systems
journal, January 1970


Physics issues of compact drift optimized stellarators
journal, June 2001


Development of a Robust Quasi-Poloidal Compact Stellarator
journal, January 2004

  • Strickler, Dennis J.; Hirshman, Steven P.; Spong, Donald A.
  • Fusion Science and Technology, Vol. 45, Issue 1
  • DOI: 10.13182/FST04-A421

Magnetic Configuration and Modular Coil Design for the Chinese First Quasi-Axisymmetric Stellarator
journal, January 2018

  • Liu, Haifeng; Shimizu, Akihiro; Isobe, Mitsutaka
  • Plasma and Fusion Research, Vol. 13, Issue 0
  • DOI: 10.1585/pfr.13.3405067