Tokamak plasma high field side response to an n = 3 magnetic perturbation: A comparison of 3D equilibrium solutions from seven different codes

A. Reiman, N. M. Ferraro, A. Turnbull, J. K. Park, Antoine Cerfon, T. E. Evans, M. J. Lanctot, E. A. Lazarus, Y. Liu, G. McFadden, D. Monticello, Y. Suzuki

Research output: Contribution to journalArticle

Abstract

In comparing equilibrium solutions for a DIII-D shot that is amenable to analysis by both stellarator and tokamak three-dimensional (3D) equilibrium codes, a significant disagreement has been seen between solutions of the VMEC stellarator equilibrium code and solutions of tokamak perturbative 3D equilibrium codes. The source of that disagreement has been investigated, and that investigation has led to new insights into the domain of validity of the different equilibrium calculations, and to a finding that the manner in which localized screening currents at low order rational surfaces are handled can affect global properties of the equilibrium solution. The perturbative treatment has been found to break down at surprisingly small perturbation amplitudes due to overlap of the calculated perturbed flux surfaces, and that treatment is not valid in the pedestal region of the DIII-D shot studied. The perturbative treatment is valid, however, further into the interior of the plasma, and flux surface overlap does not account for the disagreement investigated here. Calculated equilibrium solutions for simple model cases and comparison of the 3D equilibrium solutions with those of other codes indicate that the disagreement arises from a difference in handling of localized currents at low order rational surfaces, with such currents being absent in VMEC and present in the perturbative codes. The significant differences in the global equilibrium solutions associated with the presence or absence of very localized screening currents at rational surfaces suggests that it may be possible to extract information about localized currents from appropriate measurements of global equilibrium plasma properties. That would require improved diagnostic capability on the high field side of the tokamak plasma, a region difficult to access with diagnostics.

Original languageEnglish (US)
Article number063026
JournalNuclear Fusion
Volume55
Issue number6
DOIs
StatePublished - Jun 1 2015

Fingerprint

perturbation
stellarators
shot
screening
plasma equilibrium
breakdown

Keywords

  • 3D equilibrium
  • ELM
  • magnetic perturbation

ASJC Scopus subject areas

  • Condensed Matter Physics
  • Nuclear and High Energy Physics

Cite this

Tokamak plasma high field side response to an n = 3 magnetic perturbation : A comparison of 3D equilibrium solutions from seven different codes. / Reiman, A.; Ferraro, N. M.; Turnbull, A.; Park, J. K.; Cerfon, Antoine; Evans, T. E.; Lanctot, M. J.; Lazarus, E. A.; Liu, Y.; McFadden, G.; Monticello, D.; Suzuki, Y.

In: Nuclear Fusion, Vol. 55, No. 6, 063026, 01.06.2015.

Research output: Contribution to journalArticle

Reiman, A, Ferraro, NM, Turnbull, A, Park, JK, Cerfon, A, Evans, TE, Lanctot, MJ, Lazarus, EA, Liu, Y, McFadden, G, Monticello, D & Suzuki, Y 2015, 'Tokamak plasma high field side response to an n = 3 magnetic perturbation: A comparison of 3D equilibrium solutions from seven different codes', Nuclear Fusion, vol. 55, no. 6, 063026. https://doi.org/10.1088/0029-5515/55/6/063026
Reiman, A. ; Ferraro, N. M. ; Turnbull, A. ; Park, J. K. ; Cerfon, Antoine ; Evans, T. E. ; Lanctot, M. J. ; Lazarus, E. A. ; Liu, Y. ; McFadden, G. ; Monticello, D. ; Suzuki, Y. / Tokamak plasma high field side response to an n = 3 magnetic perturbation : A comparison of 3D equilibrium solutions from seven different codes. In: Nuclear Fusion. 2015 ; Vol. 55, No. 6.
@article{d87f300f05b741cdbe502595afba64ed,
title = "Tokamak plasma high field side response to an n = 3 magnetic perturbation: A comparison of 3D equilibrium solutions from seven different codes",
abstract = "In comparing equilibrium solutions for a DIII-D shot that is amenable to analysis by both stellarator and tokamak three-dimensional (3D) equilibrium codes, a significant disagreement has been seen between solutions of the VMEC stellarator equilibrium code and solutions of tokamak perturbative 3D equilibrium codes. The source of that disagreement has been investigated, and that investigation has led to new insights into the domain of validity of the different equilibrium calculations, and to a finding that the manner in which localized screening currents at low order rational surfaces are handled can affect global properties of the equilibrium solution. The perturbative treatment has been found to break down at surprisingly small perturbation amplitudes due to overlap of the calculated perturbed flux surfaces, and that treatment is not valid in the pedestal region of the DIII-D shot studied. The perturbative treatment is valid, however, further into the interior of the plasma, and flux surface overlap does not account for the disagreement investigated here. Calculated equilibrium solutions for simple model cases and comparison of the 3D equilibrium solutions with those of other codes indicate that the disagreement arises from a difference in handling of localized currents at low order rational surfaces, with such currents being absent in VMEC and present in the perturbative codes. The significant differences in the global equilibrium solutions associated with the presence or absence of very localized screening currents at rational surfaces suggests that it may be possible to extract information about localized currents from appropriate measurements of global equilibrium plasma properties. That would require improved diagnostic capability on the high field side of the tokamak plasma, a region difficult to access with diagnostics.",
keywords = "3D equilibrium, ELM, magnetic perturbation",
author = "A. Reiman and Ferraro, {N. M.} and A. Turnbull and Park, {J. K.} and Antoine Cerfon and Evans, {T. E.} and Lanctot, {M. J.} and Lazarus, {E. A.} and Y. Liu and G. McFadden and D. Monticello and Y. Suzuki",
year = "2015",
month = "6",
day = "1",
doi = "10.1088/0029-5515/55/6/063026",
language = "English (US)",
volume = "55",
journal = "Nuclear Fusion",
issn = "0029-5515",
publisher = "IOP Publishing Ltd.",
number = "6",

}

TY - JOUR

T1 - Tokamak plasma high field side response to an n = 3 magnetic perturbation

T2 - A comparison of 3D equilibrium solutions from seven different codes

AU - Reiman, A.

AU - Ferraro, N. M.

AU - Turnbull, A.

AU - Park, J. K.

AU - Cerfon, Antoine

AU - Evans, T. E.

AU - Lanctot, M. J.

AU - Lazarus, E. A.

AU - Liu, Y.

AU - McFadden, G.

AU - Monticello, D.

AU - Suzuki, Y.

PY - 2015/6/1

Y1 - 2015/6/1

N2 - In comparing equilibrium solutions for a DIII-D shot that is amenable to analysis by both stellarator and tokamak three-dimensional (3D) equilibrium codes, a significant disagreement has been seen between solutions of the VMEC stellarator equilibrium code and solutions of tokamak perturbative 3D equilibrium codes. The source of that disagreement has been investigated, and that investigation has led to new insights into the domain of validity of the different equilibrium calculations, and to a finding that the manner in which localized screening currents at low order rational surfaces are handled can affect global properties of the equilibrium solution. The perturbative treatment has been found to break down at surprisingly small perturbation amplitudes due to overlap of the calculated perturbed flux surfaces, and that treatment is not valid in the pedestal region of the DIII-D shot studied. The perturbative treatment is valid, however, further into the interior of the plasma, and flux surface overlap does not account for the disagreement investigated here. Calculated equilibrium solutions for simple model cases and comparison of the 3D equilibrium solutions with those of other codes indicate that the disagreement arises from a difference in handling of localized currents at low order rational surfaces, with such currents being absent in VMEC and present in the perturbative codes. The significant differences in the global equilibrium solutions associated with the presence or absence of very localized screening currents at rational surfaces suggests that it may be possible to extract information about localized currents from appropriate measurements of global equilibrium plasma properties. That would require improved diagnostic capability on the high field side of the tokamak plasma, a region difficult to access with diagnostics.

AB - In comparing equilibrium solutions for a DIII-D shot that is amenable to analysis by both stellarator and tokamak three-dimensional (3D) equilibrium codes, a significant disagreement has been seen between solutions of the VMEC stellarator equilibrium code and solutions of tokamak perturbative 3D equilibrium codes. The source of that disagreement has been investigated, and that investigation has led to new insights into the domain of validity of the different equilibrium calculations, and to a finding that the manner in which localized screening currents at low order rational surfaces are handled can affect global properties of the equilibrium solution. The perturbative treatment has been found to break down at surprisingly small perturbation amplitudes due to overlap of the calculated perturbed flux surfaces, and that treatment is not valid in the pedestal region of the DIII-D shot studied. The perturbative treatment is valid, however, further into the interior of the plasma, and flux surface overlap does not account for the disagreement investigated here. Calculated equilibrium solutions for simple model cases and comparison of the 3D equilibrium solutions with those of other codes indicate that the disagreement arises from a difference in handling of localized currents at low order rational surfaces, with such currents being absent in VMEC and present in the perturbative codes. The significant differences in the global equilibrium solutions associated with the presence or absence of very localized screening currents at rational surfaces suggests that it may be possible to extract information about localized currents from appropriate measurements of global equilibrium plasma properties. That would require improved diagnostic capability on the high field side of the tokamak plasma, a region difficult to access with diagnostics.

KW - 3D equilibrium

KW - ELM

KW - magnetic perturbation

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

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

U2 - 10.1088/0029-5515/55/6/063026

DO - 10.1088/0029-5515/55/6/063026

M3 - Article

VL - 55

JO - Nuclear Fusion

JF - Nuclear Fusion

SN - 0029-5515

IS - 6

M1 - 063026

ER -