An H-mode edge pedestal plasma transport benchmarking exercise was undertaken for a single DIII-D pedestal. Transport modelling codes used include 1.5D interpretive (ONETWO, GTEDGE), 1.5D predictive (ASTRA) and 2D ones (SOLPS, UEDGE). The particular DIII-D discharge considered is 98889, which has a typical low density pedestal. Profiles for the edge plasma are obtained from Thomson and charge-exchange recombination data averaged over the last 20% of the average 33.53 ms repetition time between type I edge localized modes. The modelled density of recycled neutrals is largest in the divertor X-point region and causes the edge plasma source rate to vary by a factor ∼10 2 on the separatrix. Modelled poloidal variations in the densities and temperatures on flux surfaces are small on all flux surfaces up to within about 2.6 mm (ρ N > 0.99) of the mid-plane separatrix. For the assumed Fick's-diffusion-type laws, the radial heat and density fluxes vary poloidally by factors of 2-3 in the pedestal region; they are largest on the outboard mid-plane where flux surfaces are compressed and local radial gradients are largest. Convective heat flows are found to be small fractions of the electron (10%) and ion (25%) heat flows in this pedestal. Appropriately averaging the transport fluxes yields interpretive 1.5D effective diffusivities that are smallest near the mid-point of the pedestal. Their 'transport barrier' minima are about 0.3 (electron heat), 0.15 (ion heat) and 0.035 (density) m 2 s −1. Electron heat transport is found to be best characterized by electron-temperaturegradient-induced transport at the pedestal top and paleoclassical transport throughout the pedestal. The effective ion heat diffusivity in the pedestal has a different profile from the neoclassical prediction and may be smaller than it. The very small effective density diffusivity may be the result of an inward pinch flow nearly balancing a diffusive outward radial density flux. The inward ion pinch velocity and density diffusion coefficient are determined by a new interpretive analysis technique that uses information from the force balance (momentum conservation) equations; the paleoclassical transport model provides a plausible explanation of these new results. Finally, the measurements and additional modelling needed to facilitate better pedestal plasma transport modelling are discussed.

Abstract. Recent DIII-D pedestal studies provide improved characterization of pedestal scaling for comparison with models. A new pedestal model accurately predicts the maximum achieved pedestal width and height in TYPE-I ELMing discharges over a large range of DIII-D operational space, including ITER demonstration discharges. The model is a combination of the peeling-ballooning theory for the MHD stability limits on the pedestal with a simple empirical pedestal width scaling, which is proportional to the square root of the pedestal poloidal beta. Width scalings based on the ion toroidal or poloidal gyroradius are much poorer descriptions of DIII-D data. A mass scaling experiment in H and D provides additional support for a poloidal beta scaling and is not consistent with an ion poloidal gyroradius scaling. Studies of pedestal evolution during the inter-ELM cycle provide evidence that the pedestal width and height increase during pedestal buildup. These studies show a correlation of pedestal width with both the square root of the pedestal poloidal beta and the square root of the pedestal ion temperature during the pedestal buildup.

An analytic model, derived from coupled continuity equations for the electron and neutral deuterium densities, is consistent with many features of edge electron density profiles in the DIII-D tokamak [J. L. Luxon et al., Plasma Physics and Controlled Fusion Research, 1986, Vol. I (International Atomic Energy Agency, Vienna, 1987), p. 159]. For an assumed constant particle diffusion coefficient, the model shows that particle transport and neutral fueling produce electron and neutral density profiles that have the same characteristic scale lengths at the plasma edge. For systematic variations of density in high-mode (H-mode) discharges, the model predicts that the width of the electron density transport barrier decreases and the maximum gradient increases, as observed in the experiments. The widths computed from the model agree quantitatively with the experimental widths for conditions in which the model is valid. These results support models of transport barrier formation in which the H-mode particle barrier is driven by the edge particle flux and the width of the barrier is approximately the neutral penetration length.

No abstract

To study the effect of edge neutrals on the low-to-high confinement transition threshold, a broad range of plasma discharges has been analyzed. These discharges vary by gas puffing and pumping rates, position of the X point, and line-averaged density. It is shown that the determination of the neutral density (or neutral pressure) in the scrape-off layer (SOL) can give a misleading indication of the neutral population inside the separatrix. An increase of neutral density in the SOL creates an increase of plasma density that, in turn, increases the opacity to the neutrals and results in reduced neutral penetration. At a constant magnetic field, the transition power divided by the density appears to be a function of a single parameter measuring the neutrals effect. From this analysis, this parameter cannot be uniquely identified. For instance, it may be the radial decay length of the neutral profile or the charge-exchange damping rate at about r/a≈0.95. A similar correlation exists between these neutral parameters and local plasma parameters such as electron and ion temperature. This indicates that a missing parameter linked to the neutrals exists in the power threshold scaling laws.