skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information
  1. Anomalous transport in the H-mode pedestal of Alcator C-Mod discharges

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

    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 simulationsmore » 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.« less
  3. The impact of collisionality, FLR, and parallel closure effects on instabilities in the tokamak pedestal: Numerical studies with the NIMROD code

    The extended-MHD NIMROD code [C. R. Sovinec and J. R. King, J. Comput. Phys. 229, 5803 (2010)] is verified against the ideal-MHD ELITE code [H. R. Wilson et al., Phys. Plasmas 9, 1277 (2002)] on a diverted tokamak discharge. When the NIMROD model complexity is increased incrementally, resistive and first-order finite-Larmour radius effects are destabilizing and stabilizing, respectively. The full result is compared to local analytic calculations which are found to overpredict both the resistive destabilization and drift stabilization in comparison to the NIMROD computations.
  4. MHD modeling of a DIII-D low-torque QH-mode discharge and comparison to observations

    Extended-MHD modeling of DIII-D tokamak quiescent H-mode (QH-mode) discharges with nonlinear NIMROD simulations saturates into a turbulent state but does not saturate when the steady-state flow inferred from measurements is not included. This is consistent with the experimental observations of the quiescent regime on DIII-D. The simulation with flow develops into a saturated turbulent state where the n Φ = 1 and 2 toroidal modes become dominant through an inverse cascade. Each mode in the range of n Φ = 1–5 is dominant at a different time. Consistent with experimental observations during QH-mode, the simulated state leads to large particlemore » transport relative to the thermal transport. Analysis shows that the amplitude and phase of the density and temperature perturbations differ resulting in greater fluctuation-induced convective particle transport relative to the convective thermal transport. As a result, comparison to magnetic-coil measurements shows that rotation frequencies differ between the simulation and experiment, which indicates that more sophisticated extended-MHD two-fluid modeling is required.« less
  5. The impact of collisionality, FLR, and parallel closure effects on instabilities in the tokamak pedestal: Numerical studies with the NIMROD code

    The extended-MHD NIMROD code [C. R. Sovinec and J. R. King, J. Comput. Phys. 229, 5803 (2010)] is verified against the ideal-MHD ELITE code [H. R. Wilson et al., Phys. Plasmas 9, 1277 (2002)] on a diverted tokamak discharge. When the NIMROD model complexity is increased incrementally, resistive and first-order finite-Larmour radius effects are destabilizing and stabilizing, respectively. Lastly, the full result is compared to local analytic calculations which are found to overpredict both the resistive destabilization and drift stabilization in comparison to the NIMROD computations.
  6. Microtearing modes in tokamak discharges

    Microtearing modes (MTMs) have been identified as a source of significant electron thermal transport in tokamak discharges. In order to describe the evolution of these discharges, it is necessary to improve the prediction of electron thermal transport. This can be accomplished by utilizing a model for transport driven by MTMs in whole device predictive modeling codes. The objective of this paper is to develop the dispersion relation that governs the MTM driven transport. A unified fluid/kinetic approach is used in the development of a nonlinear dispersion relation for MTMs. The derivation includes the effects of electrostatic and magnetic fluctuations, arbitrarymore » electron-ion collisionality, electron temperature and density gradients, magnetic curvature, and the effects associated with the parallel propagation vector. An iterative nonlinear approach is used to calculate the distribution function employed in obtaining the nonlinear parallel current and the nonlinear dispersion relation. The third order nonlinear effects in magnetic fluctuations are included, and the influence of third order effects on a multi-wave system is considered. An envelope equation for the nonlinear microtearing modes in the collision dominant limit is introduced in order to obtain the saturation level. In the limit that the mode amplitude does not vary along the field line, slab geometry, and strong collisionality, the fluid dispersion relation for nonlinear microtearing modes is found to agree with the kinetic dispersion relation.« less
  7. Fusion power production in International Thermonuclear Experimental Reactor baseline H-mode scenarios

    Self-consistent simulations of 15 MA ITER H-mode DT scenarios, from ramp-up through flat-top, are carried out. Electron and ion temperatures, toroidal angular frequency, and currents are evolved, in simulations carried out using the predictive TRANSPort and integrated modeling code starting with initial profiles and equilibria obtained from tokamak simulation code studies. Studies are carried out examining the dependence and sensitivity of fusion power production on electron density, argon impurity concentration, choice of radio frequency heating, pedestal temperature without and with E × B flow shear effects included, and the degree of plasma rotation. The goal of these whole-device ITER simulationsmore » is to identify dependencies that might impact ITER fusion performance.« less
  8. Integrated modeling of temperature profiles in L-mode tokamak discharges

    Simulations of doublet III-D, the joint European tokamak, and the tokamak fusion test reactor L-mode tokamak plasmas are carried out using the PTRANSP predictive integrated modeling code. The simulation and experimental temperature profiles are compared. The time evolved temperature profiles are computed utilizing the Multi-Mode anomalous transport model version 7.1 (MMM7.1) which includes transport associated with drift-resistive-inertial ballooning modes (the DRIBM model [T. Rafiq et al., Phys. Plasmas 17, 082511 (2010)]). The tokamak discharges considered involved a broad range of conditions including scans over gyroradius, ITER like current ramp-up, with and without neon impurity injection, collisionality, and low and highmore » plasma current. The comparison of simulation and experimental temperature profiles for the discharges considered is shown for the radial range from the magnetic axis to the last closed flux surface. The regions where various modes in the Multi-Mode model contribute to transport are illustrated. In the simulations carried out using the MMM7.1 model it is found that: The drift-resistive-inertial ballooning modes contribute to the anomalous transport primarily near the edge of the plasma; transport associated with the ion temperature gradient and trapped electron modes contribute in the core region but decrease in the region of the plasma boundary; and neoclassical ion thermal transport contributes mainly near the center of the discharge.« less
  9. Kinetic modeling of divertor heat load fluxes in the Alcator C-Mod and DIII-D tokamaks

    The guiding-center kinetic neoclassical transport code, XGC0 [Chang et al., Phys. Plasmas 11, 2649 (2004)], is used to compute the heat fluxes and the heat-load width in the outer divertor plates of Alcator C-Mod and DIII-D tokamaks. The dependence of the width of heat-load fluxes on neoclassical effects, neutral collisions, and anomalous transport is investigated using the XGC0 code. The XGC0 code includes realistic X-point geometry, a neutral source model, the effects of collisions, and a diffusion model for anomalous transport. It is observed that the width of the XGC0 neoclassical heat-load is approximately inversely proportional to the total plasmamore » current I{sub p.} The scaling of the width of the divertor heat-load with plasma current is examined for an Alcator C-Mod discharge and four DIII-D discharges. The scaling of the divertor heat-load width with plasma current is found to be weaker in the Alcator C-Mod discharge compared to scaling found in the DIII-D discharges. The effect of neutral collisions on the 1/I{sub p} scaling of heat-load width is shown not to be significant. Although inclusion of poloidally uniform anomalous transport results in a deviation from the 1/I{sub p} scaling, the inclusion of the anomalous transport that is driven by ballooning-type instabilities results in recovering the neoclassical 1/I{sub p} scaling. The Bohm or gyro-Bohm scalings of anomalous transport do not strongly affect the dependence of the heat-load width on plasma current. The inclusion of anomalous transport, in general, results in widening the width of neoclassical divertor heat-load and enhances the neoclassical heat-load fluxes on the divertor plates. Understanding heat transport in the tokamak scrape-off layer plasmas is important for strengthening the basis for predicting divertor conditions in ITER.« less
...

Search for:
All Records
Creator / Author
"Pankin, A. Y."

Refine by:
Resource Type
Availability
Publication Date
Creator / Author
Research Organization