Development and validation of a critical gradient energetic particle driven Alfven eigenmode transport model for DIII-D tilted neutral beam experiments

Recent experiments with the DIII-D tilted neutral beam injection (NBI) varying the beam energetic particle (EP) source profiles have provided strong evidence that unstable Alfven eigenmodes (AE) drive stiff EP transport at a critical EP density gradient [Heidbrink et al 2013 Nucl. Fusion 53 093006]. Here the critical gradient is identified by the local AE growth rate being equal to the local ITG/TEM growth rate at the same low toroidal mode number. The growth rates are taken from the gyrokinetic code GYRO. Simulation show that the slowing down beam-like EP distribution has a slightly lower critical gradient than the Maxwellian. The ALPHA EP density transport code [Waltz and Bass 2014 Nucl. Fusion 54 104006], used to validate the model, combines the low-n stiff EP critical density gradient AE mid-core transport with the Angioni et al (2009 Nucl. Fusion 49 055013) energy independent high-n ITG/TEM density transport model controling the central core EP density profile. For the on-axis NBI heated DIII-D shot 146102, while the net loss to the edge is small, about half the birth fast ions are transported from the central core r/a < 0.5 and the central density is about half the slowing down density. These results are in good agreement with experimental fast ion pressure profiles inferred from MSE constrained EFIT equilibria

Adriana Lunardi: a vendedora de fósforos

Recent experiments with the DIII-D tilted neutral beam injection (NBI) varying the beam energetic particle (EP) source profiles have provided strong evidence that unstable Alfven eigenmodes (AE) drive stiff EP transport at a critical EP density gradient [Heidbrink et al 2013 Nucl. Fusion 53 093006]. Here the critical gradient is identified by the local AE growth rate being equal to the local ITG/TEM growth rate at the same low toroidal mode number. The growth rates are taken from the gyrokinetic code GYRO. Simulation show that the slowing down beam-like EP distribution has a slightly lower critical gradient than the Maxwellian. The ALPHA EP density transport code [Waltz and Bass 2014 Nucl. Fusion 54 104006], used to validate the model, combines the low-n stiff EP critical density gradient AE mid-core transport with the Angioni et al (2009 Nucl. Fusion 49 055013) energy independent high-n ITG/TEM density transport model controling the central core EP density profile. For the on-axis NBI heated DIII-D shot 146102, while the net loss to the edge is small, about half the birth fast ions are transported from the central core r/a < 0.5 and the central density is about half the slowing down density. These results are in good agreement with experimental fast ion pressure profiles inferred from MSE constrained EFIT equilibria

Progress toward fully noninductive, high beta conditions in DIII-D

The DIII-D Advanced Tokamak (AT) program in the DIII-D tokamak [J. L. Luxon, Plasma Physics and Controlled Fusion Research, 1986, Vol. I (International Atomic Energy Agency, Vienna, 1987), p. 159] is aimed at developing a scientific basis for steady-state, high-performance operation in future devices. This requires simultaneously achieving 100% noninductive operation with high self-driven bootstrap current fraction and toroidal beta. Recent progress in this area includes demonstration of 100% noninductive conditions with toroidal beta, ΒT =3.6%, normalized beta, ΒN =3.5, and confinement factor, H89 =2.4 with the plasma current driven completely by bootstrap, neutral beam current drive, and electron cyclotron current drive (ECCD). The equilibrium reconstructions indicate that the noninductive current profile is well aligned, with little inductively driven current remaining anywhere in the plasma. The current balance calculation improved with beam ion redistribution that was supported by recent fast ion diagnostic measurements. The duration of this state is limited by pressure profile evolution, leading to magnetohydrodynamic (MHD) instabilities after about 1 s or half of a current relaxation time (τCR). Stationary conditions are maintained in similar discharges (∼90% noninductive), limited only by the 2 s duration (1 τCR) of the present ECCD systems. By discussing parametric scans in a global parameter and profile databases, the need for low density and high beta are identified to achieve full noninductive operation and good current drive alignment. These experiments achieve the necessary fusion performance and bootstrap fraction to extrapolate to the fusion gain, Q=5 steady-state scenario in the International Thermonuclear Experimental Reactor (ITER) [R. Aymar, Fusion Energy Conference on Controlled Fusion and Plasma Physics, Sorrento, Italy (International Atomic Energy Agency, Vienna, 1987), paper IAEA-CN-77/OV-1]. The modeling tools that have been successfully employed to both plan and interpret the experiment are used to plan future DIII-D experiments with higher power and longer pulse ECCD and fast wave and co- and counterneutral beam injection in a pumped double-null configuration. The models predict our ability to control the current and pressure profiles to reach full noninductivity with increased beta, bootstrap fraction, and duration. The same modeling tools are applied to ITER, predicting favorable prospects for the success of the ITER steady-state scenario. © 2006 American Institute of Physics