Abstract
In this paper, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical (NC) ion thermal transport simulated by NEO, account for 30%-60% of the total experimental heat flux. In addition, the NC electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne/L Te. Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for NC particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of NC particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. The implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.
Original language | English |
---|---|
Article number | 056005 |
Journal | Nuclear Fusion |
Volume | 61 |
Issue number | 5 |
DOIs | |
State | Published - May 2021 |
Funding
We would like to thank D.R. Hatch, P.B. Snyder and W. Dorland for useful discussions and S. Smith (GA) for assistance with some computational issues. We also thank an attentive referee for pointing out an incorrect interpretation regarding ETG saturation. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. This work was also supported by the U.S.Department of Energy underDE-AC02-09CH11466,DEFC02-04ER54698 and DE-FC02-06ER54873. Notice: This manuscript is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, and has been authored by Princeton University under Contract Number DE-AC02-09CH11466 with the U.S. Department of Energy.
Funders | Funder number |
---|---|
U.S.Department of Energy | DE-FC02-06ER54873, DEFC02-04ER54698 |
U.S. Department of Energy | DE-AC02-05CH11231 |
Office of Science | |
Fusion Energy Sciences | |
Princeton University | DE-AC02-09CH11466 |
Keywords
- ETG pedestal transport model
- gyrokinetic simulations
- pedestal transport
- validation