DESC: A stellarator equilibrium solver

D. W. Dudt, E. Kolemen

Research output: Contribution to journalArticlepeer-review

13 Scopus citations


The new code DESC is presented to solve for fixed-boundary ideal magnetohydrodynamic equilibria in stellarators. The approach directly solves the equilibrium force balance as a system of nonlinear equations in the form f (x) = 0. The independent variables x represent nested magnetic flux surfaces expressed in the inverse representation with toroidal flux coordinates, and the equations f (x) quantify equilibrium force balance errors at discrete points in real space. Discretizing with global Fourier-Zernike basis functions properly treats the magnetic axis and minimizes the number of coefficients needed to describe the flux surfaces. The pseudospectral method provides great flexibility in where the errors are evaluated, and the system of equations is efficiently solved with a Newton-Raphson iteration. Equilibria are computed and compared against VMEC for both axisymmetric and non-axisymmetric examples. The results show fast convergence rates and solutions with low errors throughout the plasma volume.

Original languageEnglish (US)
Article number102513
JournalPhysics of Plasmas
Issue number10
StatePublished - Oct 1 2020

All Science Journal Classification (ASJC) codes

  • Condensed Matter Physics


Dive into the research topics of 'DESC: A stellarator equilibrium solver'. Together they form a unique fingerprint.

Cite this