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

Fingerprint

Conformal mapping
Magnetohydrodynamics
Elongation
conformal mapping
Fluxes
elongation
Safety factor
poloidal flux
Wave propagation
Integral equations
safety factors
Derivatives
magnetohydrodynamics
Geometry
integral equations
wave propagation
grids
shear
geometry

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)

Cite this

ECOM : A fast and accurate solver for toroidal axisymmetric MHD equilibria. / Lee, Jungpyo; Cerfon, Antoine.

In: Computer Physics Communications, Vol. 190, 01.05.2015, p. 72-88.

Research output: Contribution to journalArticle

@article{ae9831c25b3c4214959b23e932ba32f2,
title = "ECOM: A fast and accurate solver for toroidal axisymmetric MHD equilibria",
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{\"u}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.",
keywords = "52.30.Cv, 52.55.Fa, Conformal mapping, Integral equation method, Magnetohydrodynamic equilibria",
author = "Jungpyo Lee and Antoine Cerfon",
year = "2015",
month = "5",
day = "1",
doi = "10.1016/j.cpc.2015.01.015",
language = "English (US)",
volume = "190",
pages = "72--88",
journal = "Computer Physics Communications",
issn = "0010-4655",
publisher = "Elsevier",

}

TY - JOUR

T1 - ECOM

T2 - A fast and accurate solver for toroidal axisymmetric MHD equilibria

AU - Lee, Jungpyo

AU - Cerfon, Antoine

PY - 2015/5/1

Y1 - 2015/5/1

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

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

KW - 52.30.Cv

KW - 52.55.Fa

KW - Conformal mapping

KW - Integral equation method

KW - Magnetohydrodynamic equilibria

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

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

U2 - 10.1016/j.cpc.2015.01.015

DO - 10.1016/j.cpc.2015.01.015

M3 - Article

VL - 190

SP - 72

EP - 88

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

ER -