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

Title: Investigation of the plasma shaping effects on the H-mode pedestal structure using coupled kinetic neoclassical/MHD stability simulations

Abstract

The effects of plasma shaping on the H-mode pedestal structure are investigated. High fidelity kinetic simulations of the neoclassical pedestal dynamics are combined with the magnetohydrodynamic (MHD) stability conditions for triggering edge localized mode (ELM) instabilities that limit the pedestal width and height in H-mode plasmas. We use the neoclassical kinetic XGC0 code [Chang et al., Phys. Plasmas 11, 2649 (2004)] to carry out a scan over plasma elongation and triangularity. As plasma profiles evolve, the MHD stability limits of these profiles are analyzed with the ideal MHD ELITE code [Snyder et al., Phys. Plasmas 9, 2037 (2002)]. In simulations with the XGC0 code, which includes coupled ion-electron dynamics, yield predictions for both ion and electron pedestal profiles. The differences in the predicted H-mode pedestal width and height for the DIII-D discharges with different elongation and triangularities are discussed. For the discharges with higher elongation, it is found that the gradients of the plasma profiles in the H-mode pedestal reach semi-steady states. In these simulations, the pedestal slowly continues to evolve to higher pedestal pressures and bootstrap currents until the peeling-ballooning stability conditions are satisfied. The discharges with lower elongation do not reach the semi-steady state, and ELM crashes aremore » triggered at earlier times. The plasma elongation is found to have a stronger stabilizing effect than the plasma triangularity. For the discharges with lower elongation and lower triangularity, the ELM frequency is large, and the H-mode pedestal evolves rapidly. It is found that the temperature of neutrals in the scrape-off-layer (SOL) region can affect the dynamics of the H-mode pedestal buildup. But the final pedestal profiles are nearly independent of the neutral temperature. The elongation and triangularity affect the pedestal widths of plasma density and electron temperature profiles differently. This provides a new mechanism of controlling the pedestal bootstrap current and the pedestal stability.« less

Authors:
 [1];  [2];  [2];  [3];  [4]; ORCiD logo [5]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Lehigh Univ., Bethlehem, PA (United States). Dept. of Physics
  3. National Fusion Research Inst., Daejeon (South Korea)
  4. General Atomics, San Diego, CA (United States)
  5. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Fusion Energy Sciences (FES) (SC-24)
OSTI Identifier:
1367977
Alternate Identifier(s):
OSTI ID: 1361916
Report Number(s):
LLNL-JRNL-718637
Journal ID: ISSN 1070-664X; PHPAEN
Grant/Contract Number:
AC52-07NA27344; SC0012174; SC0008605; FG02-92ER54141; SC0013977; AC02-05CH1131; AC02-05CH11231
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Physics of Plasmas
Additional Journal Information:
Journal Volume: 24; Journal Issue: 6; Journal ID: ISSN 1070-664X
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION

Citation Formats

Pankin, A. Y., Rafiq, T., Kritz, A. H., Park, G. Y., Snyder, P. B., and Chang, C. S.. Investigation of the plasma shaping effects on the H-mode pedestal structure using coupled kinetic neoclassical/MHD stability simulations. United States: N. p., 2017. Web. doi:10.1063/1.4984256.
Pankin, A. Y., Rafiq, T., Kritz, A. H., Park, G. Y., Snyder, P. B., & Chang, C. S.. Investigation of the plasma shaping effects on the H-mode pedestal structure using coupled kinetic neoclassical/MHD stability simulations. United States. doi:10.1063/1.4984256.
Pankin, A. Y., Rafiq, T., Kritz, A. H., Park, G. Y., Snyder, P. B., and Chang, C. S.. Thu . "Investigation of the plasma shaping effects on the H-mode pedestal structure using coupled kinetic neoclassical/MHD stability simulations". United States. doi:10.1063/1.4984256.
@article{osti_1367977,
title = {Investigation of the plasma shaping effects on the H-mode pedestal structure using coupled kinetic neoclassical/MHD stability simulations},
author = {Pankin, A. Y. and Rafiq, T. and Kritz, A. H. and Park, G. Y. and Snyder, P. B. and Chang, C. S.},
abstractNote = {The effects of plasma shaping on the H-mode pedestal structure are investigated. High fidelity kinetic simulations of the neoclassical pedestal dynamics are combined with the magnetohydrodynamic (MHD) stability conditions for triggering edge localized mode (ELM) instabilities that limit the pedestal width and height in H-mode plasmas. We use the neoclassical kinetic XGC0 code [Chang et al., Phys. Plasmas 11, 2649 (2004)] to carry out a scan over plasma elongation and triangularity. As plasma profiles evolve, the MHD stability limits of these profiles are analyzed with the ideal MHD ELITE code [Snyder et al., Phys. Plasmas 9, 2037 (2002)]. In simulations with the XGC0 code, which includes coupled ion-electron dynamics, yield predictions for both ion and electron pedestal profiles. The differences in the predicted H-mode pedestal width and height for the DIII-D discharges with different elongation and triangularities are discussed. For the discharges with higher elongation, it is found that the gradients of the plasma profiles in the H-mode pedestal reach semi-steady states. In these simulations, the pedestal slowly continues to evolve to higher pedestal pressures and bootstrap currents until the peeling-ballooning stability conditions are satisfied. The discharges with lower elongation do not reach the semi-steady state, and ELM crashes are triggered at earlier times. The plasma elongation is found to have a stronger stabilizing effect than the plasma triangularity. For the discharges with lower elongation and lower triangularity, the ELM frequency is large, and the H-mode pedestal evolves rapidly. It is found that the temperature of neutrals in the scrape-off-layer (SOL) region can affect the dynamics of the H-mode pedestal buildup. But the final pedestal profiles are nearly independent of the neutral temperature. The elongation and triangularity affect the pedestal widths of plasma density and electron temperature profiles differently. This provides a new mechanism of controlling the pedestal bootstrap current and the pedestal stability.},
doi = {10.1063/1.4984256},
journal = {Physics of Plasmas},
number = 6,
volume = 24,
place = {United States},
year = {Thu Jun 08 00:00:00 EDT 2017},
month = {Thu Jun 08 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on June 8, 2018
Publisher's Version of Record

Save / Share:
  • Multi-species kinetic neoclassical transport through the QH-mode pedestal and scrape-off layer on DIII-D is calculated using XGC0, a 5D full-f particle-in-cell drift-kinetic solver with self-consistent neutral recycling and sheath potentials. Quantitative agreement between the flux-driven simulation and the experimental electron density, impurity density, and orthogonal measurements of impurity temperature and flow profiles is achieved by adding random-walk particle diffusion to the guiding-center drift motion. The radial electric field (E{sub r}) that maintains ambipolar transport across flux surfaces and to the wall is computed self-consistently on closed and open magnetic field lines and is in excellent agreement with experiment. The E{submore » r} inside the separatrix is the unique solution that balances the outward flux of thermal tail deuterium ions against the outward neoclassical electron flux and inward pinch of impurity and colder deuterium ions. Particle transport in the pedestal is primarily due to anomalous transport, while the ion heat and momentum transport are primarily due to the neoclassical transport. The full-f treatment quantifies the non-Maxwellian energy distributions that describe a number of experimental observations in low-collisionallity pedestals on DIII-D, including intrinsic co-I{sub p} parallel flows in the pedestal, ion temperature anisotropy, and large impurity temperatures in the scrape-off layer.« less
  • Multi-species kinetic neoclassical transport through the QH-mode pedestal and scrapeoff layer on DIII-D is calculated using XGC0, a 5D full-f particle-in-cell drift-kinetic solver with self-consistent neutral recycling and sheath potentials. We achieved quantitative agreement between the fluxdriven simulation and the experimental electron density, impurity density and orthogonal measurements of impurity temperature and flow profiles by adding random-walk particle diffusion to the guiding-center drift motion. Furthermore, we computed the radial electric field (Er) that maintains ambipolar transport across flux surfaces and to the wall self-consistently on closed and open magnetic field lines, and is in excellent agreement with experiment. The Ermore » inside the separatrix is the unique solution that balances the outward flux of thermal tail deuterium ions against the outward neoclassical electron flux and inward pinch of impurity and colder deuterium ions. Particle transport in the pedestal is primarily due to anomalous transport, while the ion heat and momentum transport is primarily due to the neoclassical transport. The full-f treatment quantifies the non-Maxwellian energy distributions that describe a number of experimental observations in low-collisionallity pedestals on DIII-D, including intrinsic co-Ip parallel flows in the pedestal, ion temperature anisotropy and large impurity temperatures in the scrape-off layer.« less
  • The role of neoclassical, anomalous and neutral transport to the overall H-mode pedestal and scrape-off layer (SOL) structure in an ELM-free QH-mode discharge on DIII-D is explored using XGC0, a 5D full-f multi-species particle-in-cell drift-kinetic solver with self-consistent neutral recycling and sheath potentials. The work in this paper builds on previous work aimed at achieving quantitative agreement between the flux-driven simulation and the experimental electron density, impurity density and orthogonal measurements of impurity temperature and flow profiles. Improved quantitative agreement is achieved by performing the calculations with a more realistic electron mass, larger neutral density and including finite-Larmor-radius corrections self-consistentlymore » in the drift-kinetic motion of the particles. Consequently, the simulations provide stronger evidence that the radial electric field (E r) in the pedestal is primarily established by the required balance between the loss of high-energy tail main ions against a pinch of colder main ions and impurities. The kinetic loss of a small population of ions carrying a large proportion of energy and momentum leads to a separation of the particle and energy transport rates and introduces a source of intrinsic edge torque. Ion orbit loss and finite orbit width effects drive the energy distributions away from Maxwellian, and describe the anisotropy, poloidal asymmetry and local minimum near the separatrix observed in the T i profile.« less