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 language | English (US) |
---|---|
Pages (from-to) | 72-88 |
Number of pages | 17 |
Journal | Computer Physics Communications |
Volume | 190 |
DOIs | |
State | Published - May 1 2015 |
Keywords
- 52.30.Cv
- 52.55.Fa
- Conformal mapping
- Integral equation method
- Magnetohydrodynamic equilibria
ASJC Scopus subject areas
- Hardware and Architecture
- General Physics and Astronomy