skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Verification of the ideal magnetohydrodynamic response at rational surfaces in the VMEC code

Journal Article · · Physics of Plasmas
DOI:https://doi.org/10.1063/1.4939881· OSTI ID:1254729

The VMEC nonlinear ideal MHD equilibrium code [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)] is compared against analytic linear ideal MHD theory in a screw-pinch-like configuration. The focus of such analysis is to verify the ideal MHD response at magnetic surfaces which possess magnetic transform (ι) which is resonant with spectral values of the perturbed boundary harmonics. A large aspect ratio circular cross section zero-beta equilibrium is considered. This equilibrium possess a rational surface with safety factor q = 2 at a normalized flux value of 0.5. A small resonant boundary perturbation is introduced, exciting a response at the resonant rational surface. The code is found to capture the plasma response as predicted by a newly developed analytic theory that ensures the existence of nested flux surfaces by allowing for a jump in rotational transform (ι=1/q). The VMEC code satisfactorily reproduces these theoretical results without the necessity of an explicit transform discontinuity (Δι) at the rational surface. It is found that the response across the rational surfaces depends upon both radial grid resolution and local shear (dι/dΦ, where ι is the rotational transform and Φ the enclosed toroidal flux). Calculations of an implicit Δι suggest that it does not arise due to numerical artifacts (attributed to radial finite differences in VMEC) or existence conditions for flux surfaces as predicted by linear theory (minimum values of Δι). Scans of the rotational transform profile indicate that for experimentally relevant levels of transform shear the response becomes increasing localised. Furthermore, careful examination of a large experimental tokamak equilibrium, with applied resonant fields, indicates that this shielding response is present, suggesting the phenomena is not limited to this verification exercise.

Research Organization:
Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
Sponsoring Organization:
USDOE
Grant/Contract Number:
AC02-09CH11466
OSTI ID:
1254729
Alternate ID(s):
OSTI ID: 1234781
Report Number(s):
PPPL-5205; PHPAEN
Journal Information:
Physics of Plasmas, Vol. 23, Issue 1; Related Information: A publisher's note has been published for this content:Publisher's Note: “Verification of the ideal magnetohydrodynamic response at rational surfaces in the VMEC code” [Phys. Plasmas 23, 012507 (2016)]; ISSN 1070-664X
Publisher:
American Institute of Physics (AIP)Copyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 19 works
Citation information provided by
Web of Science

References (11)

Steepest-descent moment method for three-dimensional magnetohydrodynamic equilibria journal January 1983
Hydromagnetic stability of a diffuse linear pinch journal June 1960
Existence of three-dimensional ideal-magnetohydrodynamic equilibria with current sheets journal September 2015
Comparisons of linear and nonlinear plasma response models for non-axisymmetric perturbations journal May 2013
Experimental tests of linear and nonlinear three-dimensional equilibrium models in DIII-D journal July 2015
Three-dimensional equilibrium reconstruction on the DIII-D device journal January 2015
Measurement of a helical Pfirsch–Schlüter current with reduced magnitude in HSX journal July 2013
Non-axisymmetric equilibrium reconstruction for stellarators, reversed field pinches and tokamaks journal July 2013
V3FIT: a code for three-dimensional equilibrium reconstruction journal July 2009
Computation of multi-region relaxed magnetohydrodynamic equilibria journal November 2012
Magnetic islands and singular currents at rational surfaces in three-dimensional magnetohydrodynamic equilibria journal February 2015

Cited By (5)