N.M. Ferraro

Comparisons of linear and nonlinear plasma response models for non-axisymmetric perturbationsa)

With the installation of non-axisymmetric coil systems on major tokamaks for the purpose of studying the prospects of ELM-free operation, understanding the plasma response to the applied fields is a crucial issue. Application of different response models, using standard tools, to DIII-D discharges with applied non-axisymmetric fields from internal coils, is shown to yield qualitatively different results. The plasma response can be treated as an initial value problem, following the system dynamically from an initial unperturbed state, or from a nearby perturbed equilibrium approach, and using both linear and nonlinear models [A. D. Turnbull, Nucl. Fusion 52, 054016 (2012)]. Criteria are discussed under which each of the approaches can yield a valid response. In the DIII-D cases studied, these criteria show a breakdown in the linear theory despite the small 10−3 relative magnitude of the applied magnetic field perturbations in this case. For nonlinear dynamical evolution simulations to reach a saturated non...

Measurement of plasma boundary displacement by n = 2 magnetic perturbations using imaging beam emission spectroscopy

Imaging beam emission spectroscopy has been used to study the displacement of the plasma boundary in ELMing H-mode discharges with a 10 Hz rotating n = 2 external magnetic field perturbation in DIII-D. The rotating magnetic field creates a helical displacement of the beam emission profile of ~2 cm on the low-field-side (LFS) midplane which rotates with the applied resonant magnetic perturbation. This shift in the beam emission profile is due primarily to a shift in the electron density profile, which is independently measured to be 1.9 cm on the LFS midplane. These boundary displacements exceed calculations for the displacement of the stable and unstable manifolds formed by the interaction of the magnetic perturbation with the divertor separatrix by a factor of 4–5, suggesting that the vacuum field model does not correctly model the effect of the magnetic perturbations even near the separatrix. The measured displacements are suggestive of a non-resonant kink response.

Role of plasma response in displacements of the tokamak edge due to applied non-axisymmetric fields

Linear, two-fluid, resistive modelling of the plasma response to applied non-axisymmetric fields shows significant displacement of edge temperature and density profiles. The calculated displacements, often of 2 cm or more in H-mode pedestals with parameters appropriate to DIII-D, are due to the helical distortions resulting from stable edge modes being driven to finite amplitude by the applied fields. In many cases, these displacements are greater in magnitude, and different in phase, than the distortions of the separatrix manifolds predicted from vacuum modelling. Comparison of these results with experimental measurements from Thomson scattering and soft x-ray imaging finds good quantitative agreement. In these experiments, the phase of the applied non-axisymmetric magnetic field was flipped or rotated in order to probe the non-axisymmetric features of the response. The poloidal structures measured by x-ray imaging show clear indications of a helical response, as opposed to simply a change in the axisymmetric transport. Inclusion of two-fluid effects and rotation are found to be important in obtaining quantitative agreement with Thomson scattering data. Modelling shows screening of islands in the H-mode pedestal, but island penetration near the top of the pedestal where the electron rotation vanishes in plasmas with co-current rotation. Enhanced transport due to these islands may provide a mechanism for maintaining the pedestal width below the stability threshold of edge-localized modes. For typical DIII-D parameters, it is shown that the linear approximation is often near or beyond the limit of validity in the H-mode edge; however, the general agreement with experimental measurements indicates that these linear results nevertheless maintain good predictive value for profile displacements.

Integrated modeling applications for tokamak experiments with OMFIT

One modeling framework for integrated tasks (OMFIT) is a comprehensive integrated modeling framework which has been developed to enable physics codes to interact in complicated workflows, and support scientists at all stages of the modeling cycle. The OMFIT development follows a unique bottom-up approach, where the framework design and capabilities organically evolve to support progressive integration of the components that are required to accomplish physics goals of increasing complexity. OMFIT provides a workflow for easily generating full kinetic equilibrium reconstructions that are constrained by magnetic and motional Stark effect measurements, and kinetic profile information that includes fast-ion pressure modeled by a transport code. It was found that magnetic measurements can be used to quantify the amount of anomalous fast-ion diffusion that is present in DIII-D discharges, and provide an estimate that is consistent with what would be needed for transport simulations to match the measured neutron rates. OMFIT was used to streamline edge-stability analyses, and evaluate the effect of resonant magnetic perturbation (RMP) on the pedestal stability, which have been found to be consistent with the experimental observations. The development of a five-dimensional numerical fluid model for estimating the effects of the interaction between magnetohydrodynamic (MHD) and microturbulence, and its systematic verification against analytic models was also supported by the framework. OMFIT was used for optimizing an innovative high-harmonic fast wave system proposed for DIII-D. For a parallel refractive index , the conditions for strong electron-Landau damping were found to be independent of launched and poloidal angle. OMFIT has been the platform of choice for developing a neural-network based approach to efficiently perform a non-linear multivariate regression of local transport fluxes as a function of local dimensionless parameters. Transport predictions for thousands of DIII-D discharges showed excellent agreement with the power balance calculations across the whole plasma radius and over a broad range of operating regimes. Concerning predictive transport simulations, the framework made possible the design and automation of a workflow that enables self-consistent predictions of kinetic profiles and the plasma equilibrium. It is found that the feedback between the transport fluxes and plasma equilibrium can significantly affect the kinetic profiles predictions. Such a rich set of results provide tangible evidence of how bottom-up approaches can potentially provide a fast track to integrated modeling solutions that are functional, cost-effective, and in sync with the research effort of the community.

Improved understanding of physics processes in pedestal structure, leading to improved predictive capability for ITER

Joint experiment/theory/modelling research has led to increased confidence in predictions of the pedestal height in ITER. This work was performed as part of a US Department of Energy Joint Research Target in FY11 to identify physics processes that control the H-mode pedestal structure. The study included experiments on C-Mod, DIII-D and NSTX as well as interpretation of experimental data with theory-based modelling codes. This work provides increased confidence in the ability of models for peeling–ballooning stability, bootstrap current, pedestal width and pedestal height scaling to make correct predictions, with some areas needing further work also being identified. A model for pedestal pressure height has made good predictions in existing machines for a range in pressure of a factor of 20. This provides a solid basis for predicting the maximum pedestal pressure height in ITER, which is found to be an extrapolation of a factor of 3 beyond the existing data set. Models were studied for a number of processes that are proposed to play a role in the pedestal ne and Te profiles. These processes include neoclassical transport, paleoclassical transport, electron temperature gradient turbulence and neutral fuelling. All of these processes may be important, with the importance being dependent on the plasma regime. Studies with several electromagnetic gyrokinetic codes show that the gradients in and on top of the pedestal can drive a number of instabilities.

Sustained suppression of type-I edge-localized modes with dominantly n = 2 magnetic fields in DIII-D

Type-I edge-localized modes (ELMs) have been suppressed in DIII-D (Luxon et al 2003 Nucl. Fusion 43 1813) H-mode discharges with a H98Y2 confinement factor near 1.0 using magnetic perturbations (MPs) with dominant toroidal mode number n = 2. This expands access to the ELM-suppressed regime, which was previously attainable in DIII-D only with n = 3 fields. ELM suppression is obtained with two rows of internal coils for 1.8 s with normalized beta of 1.9 and average triangularity of 0.53, corresponding to a scaled version of ITER scenario 2 at an ITER relevant electron collisionality of 0.2. The applied field reduces the pedestal pressure and edge current via the density without degrading the edge thermal transport barrier. ELITE calculations find that the resulting profiles are stable to intermediate-n peeling–ballooning modes. ELM suppression is found within different ranges of q95 depending on the coil configuration used to generate the MP. The edge safety factors associated with suppression do not correspond to those that maximize the pitch-resonant components of the applied vacuum field. Instead, ELM suppression is correlated with an increase in the amplification of kink-resonant components of the calculated ideal MHD plasma response field.

The limits and challenges of error field correction for ITERa)

Significant progress has been made in interpreting the effects of non-axisymmetric “error” fields on a plasma through ideal MHD stability and a dominant “least stable” ideal mode through which the fields couple to the tearing resonant surface. However, in contrast to expectations from such theories, experiments have found limited success in correcting error fields, with single correction coil arrays giving benefits of between 0% and ∼50% correction (in terms of improvement to a low density locked mode limit), dependent on the structure of the error and correcting fields. With additional coils up to ∼70% is possible. It was unclear whether this represented an intrinsic stability or control limit, or higher order toroidal or poloidal harmonic effects. Thus, studies on the DIII-D tokamak explored correction of a proxy error field, using two differently structured coil arrays. This enabled the principles of error correction to be tested at high amplitudes and operational densities, with known pure n=1 fields. Results showed substantial residual effects from the corrected n=1 field, with improvements of only ∼50% in the low density locked mode limit. This suggests that n=1 error fields must couple to more than one surface in the plasma, and this is conjectured to be through more than one ideal mode, thereby requiring precise correction. For ITER, updated predictions of field error have been obtained and compared with revised scalings for tearing mode thresholds, indicating 50% or better error field correction will be needed. This will likely require more than one well coupled correction coil array and sets a challenge for theory to model the behavior, in order to clarify the plasma response and braking mechanisms, and so the effectiveness of ITERs correction coils and the possible need for support from its edge localized mode control coils.

Modulation of prompt fast-ion loss by applied n=2 fields in the DIII-D tokamak

Energy and pitch angle resolved measurements of escaping neutral beam ions (E ≈ 80 keV) have been made during DIII-D L-mode discharges with applied, slowly rotating, n = 2 magnetic perturbations. Data from separate scintillator detectors (FILDs) near and well below the plasma midplane show fast-ion losses correlated with the internal coil (I-coil) fields. The dominant fast-ion loss signals are observed to decay within one poloidal transit time after beam turn-off indicating they are primarily prompt loss orbits. Also, during application of the rotating I-coil fields, outboard midplane edge density and bremsstrahlung emission profiles exhibit a radial displacement of up to δR ≈ 1 cm. Beam deposition and full orbit modeling of these losses using M3D-C1 calculations of the perturbed kinetic profiles and fields reproduce many features of the measured losses. In particular, the predicted phase of the modulated loss signal with respect to the I-coil currents is in close agreement with FILD measurements as is the relative amplitudes of the modulated losses for the co and counter-current beam used in the experiment. These simulations show modifications to the beam ion birth profile and subsequent prompt loss due to changes in the edge density; however, the dominant factor causing modulation of the losses to the fast-ion loss detectors is the perturbed magnetic field (δB/B ≈ 10−3 in the plasma). Calculations indicate total prompt loss to the DIII-D wall can increase with application of the n = 2 perturbation by up to 7% for co-current injected beams and 3% for counter-current injected beams depending on phase of the perturbation relative to the injected beam.

Experimental tests of linear and nonlinear three-dimensional equilibrium models in DIII-D

DIII-D experiments using new detailed magnetic diagnostics show that linear, ideal magnetohydrodynamics (MHD) theory quantitatively describes the magnetic structure (as measured externally) of three-dimensional (3D) equilibria resulting from applied fields with toroidal mode number n = 1, while a nonlinear solution to ideal MHD force balance, using the VMEC code, requires the inclusion of n ≥ 1 to achieve similar agreement. These tests are carried out near ITER baseline parameters, providing a validated basis on which to exploit 3D fields for plasma control development. Scans of the applied poloidal spectrum and edge safety factor confirm that low-pressure, n = 1 non-axisymmetric tokamak equilibria are determined by a single, dominant, stable eigenmode. However, at higher beta, near the ideal kink mode stability limit in the absence of a conducting wall, the qualitative features of the 3D structure are observed to vary in a way that is not captured by ideal MHD.

Fast-ion losses induced by ELMs and externally applied magnetic perturbations in the ASDEX Upgrade tokamak

Phase-space time-resolved measurements of fast-ion losses induced by edge localized modes (ELMs) and ELM mitigation coils have been obtained in the ASDEX Upgrade tokamak by means of multiple fast-ion loss detectors (FILDs). Filament-like bursts of fast-ion losses are measured during ELMs by several FILDs at different toroidal and poloidal positions. Externally applied magnetic perturbations (MPs) have little effect on plasma profiles, including fast-ions, in high collisionality plasmas with mitigated ELMs. A strong impact on plasma density, rotation and fast-ions is observed, however, in low density/collisionality and q95 plasmas with externally applied MPs. During the mitigation/suppression of type-I ELMs by externally applied MPs, the large fast-ion bursts observed during ELMs are replaced by a steady loss of fast-ions with a broad-band frequency and an amplitude of up to an order of magnitude higher than the neutral beam injection (NBI) prompt loss signal without MPs. Multiple FILD measurements at different positions, indicate that the fast-ion losses due to static 3D fields are localized on certain parts of the first wall rather than being toroidally/poloidally homogeneously distributed. Measured fast-ion losses show a broad energy and pitch-angle range and are typically on banana orbits that explore the entire pedestal/scrape-off-layer (SOL). Infra-red measurements are used to estimate the heat load associated with the MP-induced fast-ion losses. The heat load on the FILD detector head and surrounding wall can be up to six times higher with MPs than without 3D fields. When 3D fields are applied and density pump-out is observed, an enhancement of the fast-ion content in the plasma is typically measured by fast-ion D-alpha (FIDA) spectroscopy. The lower density during the MP phase also leads to a deeper beam deposition with an inward radial displacement of ?2?cm in the maximum of the beam emission. Orbit simulations are used to test different models for 3D field equilibrium reconstruction including vacuum representation, the free boundary NEMEC code and the two-fluid M3D-C1 code which account for the plasma response. Guiding center simulations predict the maximum level of losses, ?2.6%, with NEMEC 3D equilibrium. Full orbit simulations overestimate the level of losses in 3D vacuum fields with ?15% of lost NBI ions.