ECOM: A fast and accurate solver for toroidal axisymmetric MHD equilibria

Jungpyo Lee, Antoine Cerfon

Research output: Contribution to journalArticle

Abstract

We present ECOM (Equilibrium solver via COnformal Mapping), a fast and accurate fixed boundary solver for toroidally axisymmetric magnetohydrodynamic equilibria with or without a toroidal flow. ECOM combines conformal mapping and Fourier and integral equation methods on the unit disk to achieve exponential convergence for the poloidal flux function as well as its first and second partial derivatives. As a consequence of its high order accuracy, for dense grids and elongations comparable to or smaller than the elongation of ITER, ECOM computes key quantities such as the safety factor and the magnetic shear with higher accuracy than the finite element based code CHEASE (Lütjens et al., 1996) at equal run time. ECOM has been developed to provide equilibrium quantities and details of the flux contour geometry as inputs to stability, wave propagation and transport codes.

Original languageEnglish (US)
Pages (from-to)72-88
Number of pages17
JournalComputer Physics Communications
Volume190
DOIs
StatePublished - May 1 2015

Keywords

  • 52.30.Cv
  • 52.55.Fa
  • Conformal mapping
  • Integral equation method
  • Magnetohydrodynamic equilibria

ASJC Scopus subject areas

  • Hardware and Architecture
  • Physics and Astronomy(all)

Fingerprint Dive into the research topics of 'ECOM: A fast and accurate solver for toroidal axisymmetric MHD equilibria'. Together they form a unique fingerprint.

  • Cite this