Search

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. The main aim of this work was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ∼20 cm in inner divertor and ∼40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ∼5 MW m−2 to ∼1 MW m−2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

An integrated modeling workflow capable of finding the steady-state plasma solution with self-consistent core transport, pedestal structure, current profile, and plasma equilibrium physics has been developed and tested against a DIII-D discharge. Key features of the achieved core-pedestal coupled workflow are its ability to account for the transport of impurities in the plasma self-consistently, as well as its use of machine learning accelerated models for the pedestal structure and for the turbulent transport physics. Notably, the coupled workflow is implemented within the One Modeling Framework for Integrated Tasks (OMFIT) framework, and makes use of the ITER integrated modeling and analysis suite data structure for exchanging data among the physics codes that are involved in the simulations. Such technical advance has been facilitated by the development of a new numerical library named ordered multidimensional arrays structure.

A new small angle slot (SAS) divertor concept has been developed to enhance neutral cooling across the divertor target by coupling a closed slot structure with appropriate target shaping. Initial tests on DIII-D find a strong interplay between such anticipated 'SAS' effects and cross-field drifts, favouring operation with the ion B × ∇B drift away from the X-point, as currently employed for advanced tokamaks. This offers the following key improvements relative to DIII-D's open lower divertor or partially-closed upper divertor: (i) SAS allows for transition to low temperature moderately detached divertor conditions with T e ≲ 10 eV at very low main plasma densities, lower than are usually attainable at all in DIII-D high confinement (H-mode) plasmas as used in these tests; (ii) Pedestal performance and core confinement are significantly improved with SAS. The final confinement collapse associated with the onset of X-point MARFE (multifaceted asymmetric radiation from the edge) following deep detachment occurs at significantly higher pedestal densities, thus widening the window of H-mode operation compatible with a dissipative divertor. For operation with the ion B × ∇B drift toward the X-point, the divertor plasma transitions to a bifurcative detached state at much higher densities, similar to other divertor configurations in DIII-D. These results highlight the strong interplay between divertor closure and drifts, and point to an interesting divertor optimization path to explore that offers potential for future fusion reactors.

The 3-D field induced relativistic runaway electron (RE) loss has been simulated for DIII-D and COMPASS plasmas, utilizing the MARS-F code incorporated with the recently developed and updated RE orbit module (REORBIT). Modeling shows effectively 100% loss of a post-disruption, high-current runaway beam in DIII-D due to the 1 kG level of magnetic field perturbation produced by a fast growing n = 1 resistive kink instability. This complete RE loss is shown to be independent of the particle energy or the initial location of particles in the configuration space. Applied resonant magnetic perturbation (RMP) fields from in-vessel coils are not effective for RE beam mitigation in DIII-D but do produce finite (>10%) RE loss in COMPASS post-disruption plasmas, consistent with experimental observations in the above two devices. The major reasons for this difference in RE control by RMP between these two devices are (i) the coil proximity to the RE beam and (ii) the effective coil current scaling vs the machine size and the toroidal magnetic field. In the modeling, the lost REs due to 3-D fields deposit onto the limiting surfaces of the devices. Distributions of the lost REs to the limiting surface show a poloidally peaked profile near the high-field-side in both DIII-D and COMPASS, covering about 100 ° poloidal angle. A higher perturbation field level and/or higher particle energy also result in REs being lost to the low-field-side of the limiting surface of these two devices, increasing the effective wetted area.

On DIII-D, the high βp scenario has an internal transport barrier (ITB), βN~βp~3,q95~10, and very high normalized confinement H98,y2~1.6. Recently, plasmas starting with these conditions have been dynamically driven to q95~6 and βp~3, where we find the ITB and high performance persist for multiple energy confinement times. These conditions are projected to meet the ITER steady-state goal of Q= 5. The ITB is maintained at lower βp with a strong reverse shear, consistent with predictions that negative central shear can lower the βp threshold for the ITB. There are two observed confinement states in the highβp scenario: H-mode confinement state with a high edge pedestal, and an enhanced confinement state with a low pedestal and an ITB. It has been observed in a scan of external resonant magnetic perturbation amplitude that when there are no large type-I ELMs, there is no transition to enhanced confinement. This is consistent with the proposed mechanism for ITB formation being a type-I ELM. Quasilinear gyro-Landau fluid predictive modeling of ITER suggests that only a modest reverse shear is required to achieve the ITB formation necessary for Q=5 when electromagnetic physics including the kinetic ballooning mode (KBM) is incorporated.

Fusion whole device modeling simulations require comprehensive models that are simultaneously physically accurate, fast, robust, and predictive. In this paper we describe the development of two neural-network (NN) based models as a means to perform a snon-linear multivariate regression of theory-based models for the core turbulent transport fluxes, and the pedestal structure. Specifically, we find that a NN-based approach can be used to consistently reproduce the results of the TGLF and EPED1 theory-based models over a broad range of plasma regimes, and with a computational speedup of several orders of magnitudes. These models are then integrated into a predictive workflow that allows prediction with self-consistent core-pedestal coupling of the kinetic profiles within the last closed flux surface of the plasma. The NN paradigm is capable of breaking the speed-accuracy trade-off that is expected of traditional numerical physics models, and can provide the missing link towards self-consistent coupled core-pedestal whole device modeling simulations that are physically accurate and yet take only seconds to run.

The runaway electron (RE) distributions driven by a large toroidal electric field induced by the drop in the temperature profile due to disruption or pellets are comprehensively simulated by the 3D Fokker–Planck (FP) solver CQL3D (Harvey and McCoy 1992 Proc. of IAEA TCM), recently coupled to the Ampere–Faraday (AF) equations. The evolution of the toroidal current in a plasma occurs on a resistive time scale, τ res = 4πa 2/(c 2 η), which is typically of the order of seconds in present tokamaks. Here, a and η are respectively plasma radius or radial extent of a current density perturbation, and Ohmic resistivity. From the Faraday EM equation, the toroidal electric field is proportional to the time derivative of the poloidal magnetic field, which, from the Ampere equation, is proportional to the toroidal current. Thus, the toroidal electric field rapidly increases due to an abrupt temperature drop decrease in conductivity, to prevent change in the toroidal current faster than τ res. This is a example of Lenz’s law. For example, in simulations with KPRAD (Whyte et al 2003 J. Nucl. Mater. 313–6 1239) of neon pellet injection into a DIII-D shot, T e drops from 2 keV to 10 eV in 0.1 ms and Z eff increases 1–4, giving that the electric field increases 3500× to 0.8 V cm−1. As described in Harvey et al (2000 PoP 7 4590), this places much of the tail electron distribution beyond the Dreicer runaway velocity, giving so-called ‘hot-tail runaways’ which for a time are the dominant source of runaways, more so than the knockon source. In this prior calculation, performed for a single flux surface, the toroidal current density is held constant, on the basis that τ res is large. Most of the initial current can be converted to runaway current, which is then dangerous, particularly for ITER. A more comprehensive A–F model recently implemented in CQL3D, taking into account the time-development of the full-plasma-width toroidal electric field on time-scales of order τ res applies an iterative technique for the toroidal field previously developed for a different application (Kupfer et al 1996 PoP 3 3644), maintaining the implicit-in-time evolution of CQL3D. The degree of runaway current formation is reduced in AF augmented CQL3D, but the basic mechanism of ‘hot-tail runaways’ remains a dominant contribution to the REs at early times after the T e drop in these simulations. On the other hand, a NIMROD (Sovinec et al 2004 J. Comput. Phys. 195 355) simulation of shattered-pellet shutdown of DIII-D plasma (Kim 2018 APS/DPP Meeting), gives a slower thermal quench; when the plasma profiles and electric field are coupled one-way to CQL3D, the ‘hot-tail’ REs are much less, and growth of RE is dominated by the knockon process.

An explanation is provided for the disruptive instability in diverted tokamaks when the safety factor$q$at the 95 % poloidal flux surface,$q_{95}$, is driven below 2.0. The instability is a resistive kink counterpart to the current-driven ideal mode that traditionally explained the corresponding disruption in limited cross-sections (Shafranov,Sov. Phys. Tech. Phys., vol. 15, 1970, p. 175) when$q_{edge}$, the safety factor at the outermost closed flux surface, lies just below a rational value$m/n$. Experimentally, external kink modes are observed in limiter configurations as the current in a tokamak is ramped up and$q_{edge}$decreases through successive rational surfaces. For$q_{edge}, the instability is always encountered and is highly disruptive. However, diverted plasmas, in which$q_{edge}$is formally infinite in the magnetohydrodynamic (MHD) model, have presented a longstanding difficulty since the theory would predict stability, yet, the disruptive limit occurs in practice when$q_{95}$, reaches 2. It is shown from numerical calculations that a resistive kink mode is linearly destabilized by the rapidly increasing resistivity at the plasma edge when$q_{95}, but$q_{edge}\gg 2$. The resistive kink behaves much like the ideal kink with predominantly kink or interchange parity and no real sign of a tearing component. However, the growth rates scale with a fractional power of the resistivity near the$q=2$surface. The results have a direct bearing on the conventional edge cutoff procedures used in most ideal MHD codes, as well as implications for ITER and for future reactor options.

Major developments on the Mega Amp Spherical Tokamak (MAST) have enabled important advances in support of ITER and the physics basis of a spherical tokamak (ST) based component test facility (CTF), as well as providing new insight into underlying tokamak physics. For example, L–H transition studies benefit from high spatial and temporal resolution measurements of pedestal profile evolution (temperature, density and radial electric field) and in support of pedestal stability studies the edge current density profile has been inferred from motional Stark effect measurements. The influence of the q-profile and E × B flow shear on transport has been studied in MAST and equilibrium flow shear has been included in gyro-kinetic codes, improving comparisons with the experimental data. H-modes exhibit a weaker q and stronger collisionality dependence of heat diffusivity than implied by IPB98(y,2) scaling, which may have important implications for the design of an ST-based CTF. ELM mitigation, an important issue for ITER, has been demonstrated by applying resonant magnetic perturbations (RMPs) using both internal and external coils, but full stabilization of type-I ELMs has not been observed. Modelling shows the importance of including the plasma response to the RMP fields. MAST plasmas with q > 1 and weak central magnetic shear regularly exhibit a long-lived saturated ideal internal mode. Measured plasma braking in the presence of this mode compares well with neo-classical toroidal viscosity theory. In support of basic physics understanding, high resolution Thomson scattering measurements are providing new insight into sawtooth crash dynamics and neo-classical tearing mode critical island widths. Retarding field analyser measurements show elevated ion temperatures in the scrape-off layer of L-mode plasmas and, in the presence of type-I ELMs, ions with energy greater than 500 eV are detected 20 cm outside the separatrix. Disruption mitigation by massive gas injection has reduced divertor heat loads by up to 70%.

We report tests of an alternate technique for constraining MHD equilibrium analysis in tokamak plasmas using internal magnetic field measurements based on |B| measurements from the motional Stark splitting of Dα spectral lines emitted by a neutral heating beam (MSE-LS). We compare results using MSE-LS with those of the standard equilibrium analysis technique based on line polarization of the Dα emission (MSE-LP). An alternative to MSE-LP is needed in future devices such as ITER where MSE-LP will be difficult due to a plasma-induced coating of the first optical element. The tests utilized data from 10 DIII-D shots with 7 MSE-LS and 14 MSE-LP views covering a range of radii along the outer midplane of the plasma. Seven MSE-LS measurements can contribute significantly to the equilibrium reconstruction of pressure and q profiles using both synthetic and experimental DIII-D MSE-LS data. For example, 7 MSE-LS plus seven MSE-LP measurements give a fit quality that is as good as the same cases with 14 MSE-LP measurements. Analyzing synthetic data for 14 MSE-LS measurements shows significant improvement in fitting quality over the case with 7 MSE-LS locations.We report tests of an alternate technique for constraining MHD equilibrium analysis in tokamak plasmas using internal magnetic field measurements based on |B| measurements from the motional Stark splitting of Dα spectral lines emitted by a neutral heating beam (MSE-LS). We compare results using MSE-LS with those of the standard equilibrium analysis technique based on line polarization of the Dα emission (MSE-LP). An alternative to MSE-LP is needed in future devices such as ITER where MSE-LP will be difficult due to a plasma-induced coating of the first optical element. The tests utilized data from 10 DIII-D shots with 7 MSE-LS and 14 MSE-LP views covering a range of radii along the outer midplane of the plasma. Seven MSE-LS measurements can contribute significantly to the equilibrium reconstruction of pressure and q profiles using both synthetic and experimental DIII-D MSE-LS data. For example, 7 MSE-LS plus seven MSE-LP measurements give a fit quality that is as good as the same cases with 14 MSE-LP meas...

Experiments have demonstrated improved energy confinement due to the spontaneous formation of an internal transport barrier in high bootstrap fraction discharges. Gyrokinetic analysis, and quasilinear predictive modeling, demonstrates that the observed transport barrier is caused by the suppression of turbulence primarily from the large Shafranov shift. It is shown that the Shafranov shift can produce a bifurcation to improved confinement in regions of positive magnetic shear or a continuous reduction in transport for weak or negative magnetic shear. Operation at high safety factor lowers the pressure gradient threshold for the Shafranov shift-driven barrier formation. Two self-organized states of the internal and edge transport barrier are observed. It is shown that these two states are controlled by the interaction of the bootstrap current with magnetic shear, and the kinetic ballooning mode instability boundary. Election scale energy transport is predicted to be dominant in the inner 60% of the profile. Evidence is presented that energetic particle-driven instabilities could be playing a role in the thermal energy transport in this region.

An integrated-modeling workflow has been developed for the purpose of performing predict-first analysis of transient-stability experiments. Starting from an existing equilibrium reconstruction from a past experiment, the workflow couples together the EFIT Grad-Shafranov solver [L. Lao et al., Fusion Sci. Technol. 48, 968 (2005)], the EPED model for the pedestal structure [P. B. Snyder et al., Phys. Plasmas 16, 056118 (2009)], and the NEO drift-kinetic-equation solver [E. A. Belli and J. Candy, Plasma Phys. Controlled Fusion 54, 015015 (2012)] (for bootstrap current calculations) in order to generate equilibria with self-consistent pedestal structures as the plasma shape and various scalar parameters (e.g., normalized β, pedestal density, and edge safety factor [q95]) are changed. These equilibria are then analyzed using automated M3D-C1 extended-magnetohydrodynamic modeling [S. C. Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)] to compute the plasma response to three-dimensional magnetic perturbations. This workflow was created in conjunction with a DIII-D experiment examining the effect of triangularity on the 3D plasma response. Several versions of the workflow were developed, and the initial ones were used to help guide experimental planning (e.g., determining the plasma current necessary to maintain the constant edge safety factor in various shapes). Subsequent validation with the experimental results was then used to revise the workflow, ultimately resulting in the complete model presented here. We show that quantitative agreement was achieved between the M3D-C1 plasma response calculated for equilibria generated by the final workflow and equilibria reconstructed from experimental data. A comparison of results from earlier workflows is used to show the importance of properly matching certain experimental parameters in the generated equilibria, including the normalized β, pedestal density, and q95. On the other hand, the details of the pedestal current did not significantly impact the plasma response in these equilibria. A comparison to the experimentally measured plasma response shows mixed agreement, indicating that while the equilibria are predicted well, additional analysis tools may be needed. Finally, we note the implications that these results have for the success of future predict-first studies, particularly the need for scans of uncertain parameters and for close collaboration between experimentalists and theorists.

Measurements show fast ion losses correlated with applied three-dimensional (3D) fields in a variety of plasmas ranging from L-mode to resonant magnetic perturbation (RMP) edge localized mode (ELM) suppressed H-mode discharges. In DIII-D L-mode discharges with a slowly rotating n = 2 magnetic perturbation, scintillator detector loss signals synchronized with the applied fields are observed to decay within one poloidal transit time after beam turnoff indicating they arise predominantly from prompt loss orbits. Full orbit following using M3D-C1 calculations of the perturbed fields and kinetic profiles reproduce many features of the measured losses and points to the importance of the applied 3D field phase with respect to the beam injection location in determining the overall impact on prompt beam ion loss. Modeling of these results includes a self-consistent calculation of the 3D perturbed beam ion birth profiles and scrape-off-layer ionization, a factor found to be essential to reproducing the experimental measurements. Extension of the simulations to full slowing down timescales, including fueling and the effects of drag and pitch angle scattering, show the applied n = 3 RMPs in ELM suppressed H-mode plasmas can induce a significant loss of energetic particles from the core. With the applied n = 3 fields, up to 8.4% of the injected beam power is predicted to be lost, compared to 2.7% with axisymmetric fields only. These fast ions, originating from minor radii ρ > 0.7, are predicted to be primarily passing particles lost to the divertor region, consistent with wide field-of-view infrared periscope measurements of wall heating in n = 3 RMP ELM suppressed plasmas. Edge fast ion Dα (FIDA) measurements also confirm a large change in edge fast ion profile due to the n = 3 fields, where the effect was isolated by using short 50ms RMP-off periods during which ELM suppression was maintained yet the fast ion profile was allowed to recover. The role of resonances between fast ion drift motion and the applied 3D fields in the context of selectively targeting regions of fast ion phase space is also discussed.

The National Spherical Torus Experiment (NSTX) is currently being upgraded to operate at twice the toroidal field and plasma current (up to 1 T and 2 MA), with a second, more tangentially aimed neutral beam (NB) for current and rotation control, allowing for pulse lengths up to 5 s. Recent NSTX physics analyses have addressed topics that will allow NSTX-Upgrade to achieve the research goals critical to a Fusion Nuclear Science Facility. These include producing stable, 100% non-inductive operation in high-performance plasmas, assessing plasma–material interface (PMI) solutions to handle the high heat loads expected in the next-step devices and exploring the unique spherical torus (ST) parameter regimes to advance predictive capability. Non-inductive operation and current profile control in NSTX-U will be facilitated by co-axial helicity injection (CHI) as well as radio frequency (RF) and NB heating. CHI studies using NIMROD indicate that the reconnection process is consistent with the 2D Sweet–Parker theory. Full-wave AORSA simulations show that RF power losses in the scrape-off layer (SOL) increase significantly for both NSTX and NSTX-U when the launched waves propagate in the SOL. Toroidal Alfven eigenmode avalanches and higher frequency Alfven eigenmodes can affect NB-driven current through energy loss and redistribution of fast ions. The inclusion of rotation and kinetic resonances, which depend on collisionality, is necessary for predicting experimental stability thresholds of fast growing ideal wall and resistive wall modes. Neutral beams and neoclassical toroidal viscosity generated from applied 3D fields can be used as actuators to produce rotation profiles optimized for global stability. DEGAS-2 has been used to study the dependence of gas penetration on SOL temperatures and densities for the MGI system being implemented on the Upgrade for disruption mitigation. PMI studies have focused on the effect of ELMs and 3D fields on plasma detachment and heat flux handling. Simulations indicate that snowflake and impurity seeded radiative divertors are candidates for heat flux mitigation in NSTX-U. Studies of lithium evaporation on graphite surfaces indicate that lithium increases oxygen surface concentrations on graphite, and deuterium–oxygen affinity, which increases deuterium pumping and reduces recycling. In situ and test-stand experiments of lithiated graphite and molybdenum indicate temperature-enhanced sputtering, although that test-stand studies also show the potential for heat flux reduction through lithium vapour shielding. Non-linear gyro kinetic simulations have indicated that ion transport can be enhanced by a shear-flow instability, and that non-local effects are necessary to explain the observed rapid changes in plasma turbulence. Predictive simulations have shown agreement between a microtearing-based reduced transport model and the measured electron temperatures in a microtearing unstable regime. Two Alfven eigenmode-driven fast ion transport models have been developed and successfully benchmarked against NSTX data. Upgrade construction is moving on schedule with initial physics research operation of NSTX-U planned for mid-2015.

The major objective of the National Spherical Torus Experiment (NSTX) is to understand basic toroidal confinement physics at low aspect ratio and high βT in order to advance the spherical torus (ST) concept. In order to do this, NSTX utilizes up to 7.5 MW of neutral beam injection, up to 6 MW of high harmonic fast waves (HHFWs), and it operates with plasma currents up to 1.5 MA and elongations of up to 2.6 at a toroidal field up to 0.45 T. New facility, and diagnostic and modelling capabilities developed over the past two years have enabled the NSTX research team to make significant progress towards establishing this physics basis for future ST devices. Improvements in plasma control have led to more routine operation at high elongation and high βT (up to ~40%) lasting for many energy confinement times. βT can be limited by either internal or external modes. The installation of an active error field (EF) correction coil pair has expanded the operating regime at low density and has allowed for initial resonant EF amplification experiments. The determination of the confinement and transport properties of NSTX plasmas has benefitted greatly from the implementation of higher spatial resolution kinetic diagnostics. The parametric variation of confinement is similar to that at conventional aspect ratio but with values enhanced relative to those determined from conventional aspect ratio scalings and with a BT dependence. The transport is highly dependent on details of both the flow and magnetic shear. Core turbulence was measured for the first time in an ST through correlation reflectometry. Non-inductive start-up has been explored using PF-only and transient co-axial helicity injection techniques, resulting in up to 140 kA of toroidal current generated by the latter technique. Calculated bootstrap and beam-driven currents have sustained up to 60% of the flat-top plasma current in NBI discharges. Studies of HHFW absorption have indicated parametric decay of the wave and associated edge thermal ion heating. Energetic particle modes, most notably toroidal Alfven eigenmodes and fishbone-like modes result in fast particle losses, and these instabilities may affect fast ion confinement on devices such as ITER. Finally, a variety of techniques has been developed for fuelling and power and particle control.

An overview of the research capabilities and the future plans on the MA-class National Spherical Torus Experiment (NSTX) at Princeton is presented. NSTX research is exploring the scientific benefits of modifying the field line structure from that in more conventional aspect ratio devices, such as the tokamak. The relevant scientific issues pursued on NSTX include energy confinement, MHD stability at high beta, non-inductive sustainment, solenoid-free start-up, and power and particle handling. In support of the NSTX research goal, research tools are being developed by the NSTX team. In the context of the fusion energy development path being formulated in the US, an ST-based Component Test Facility (CTF) and, ultimately a high beta Demo device based on the ST, are being considered. For these, it is essential to develop high performance (high beta and high confinement), steady-state (non-inductively driven) ST operational scenarios and an efficient solenoid-free start-up concept. We will also briefly describe the Next-Step-ST (NSST) device being designed to address these issues in fusion-relevant plasma conditions.

This work investigates the potential roles played by the scrape-off-layer current (SOLC) in MHD activity of tokamak plasmas, including effects on stability. SOLCs are found to be an integral part of the MHD activity, with a propensity to flow in a non-axisymmetric pattern and with magnitude potentially large enough to play a role in the MHD stability. Candidate mechanisms that can drive these SOLCs are identified: (a) non-axisymmetric thermoelectric potential, (b) electromotive force from MHD activity, and (c) flux swing, both toroidal and poloidal, of the plasma column. Other potential magnetic consequences of the SOLC are identified: (a) its error field can introduce complications in feedback control schemes for stabilizing MHD activity, and (b) its non-axisymmetric field can be falsely identified as an axisymmetric field by the tokamak control logic and in equilibrium reconstruction. The radial profile of an SOLC observed during a quiescent discharge period is determined, and found to possess polarity reversals as a function of radial distance.

Steady-state operation of the advanced tokamak reactor relies on maintaining plasma stability with respect to the resistive wall mode ~RWM!. Active magnetic feedback and plasma rotation are the two methods proposed and demonstrated for this purpose. A comprehensive modeling effort including both magnetic feedback and plasma rotation is needed for understanding the physical mechanisms of the stabilization and to project to future devices. For plasma with low rotation, a complete solution for the feedback issue is obtained by assuming the plasma obeys ideal magnetohydrodynamics ~MHDs! and utilizing a normal mode approach ~NMA! @M. S. Chu et al., Nucl. Fusion 43, 441 ~2003!#. It is found that poloidal sensors are more effective than radial sensors and coils inside of the vacuum vessel more effective than outside. For plasmas with non-negligible rotation, a comprehensive linear nonideal MHD code, the MARS-F has been found to be suitable. MARS-F @Y. Q. Liu et al., Phys. Plasmas 7, 3681 ~2000!# has been benchmarked in the ideal MHD limit against the NMA. The effect of rotation stabilization of the plasma depends on the plasma dissipation model. Broad qualitative features of the experiment are reproduced. Rotation reduces the feedback gain required for RWM stabilization. Reduction is significant when rotation is near the critical rotation speed needed for stabilization. The International Thermonuclear Experimental Reactor ~ITER! @R. Aymar et al., Plasma Phys. Controlled Fusion 44, 519 ~2002!# ~scenario IV for advanced tokamak operation! may be feedback stabilized with babove the no wall limit and up to an increment of ;50% towards the ideal limit. Rotation further improves the stability.

Recent progress in the development of electron cyclotron heating⁄current drive systems in regard to power and pulse duration have allowed the extension of two normalized parameters (toroidal electric field Eφ and wave power density p). In the extended regime, the conventional theory requires consideration of a distorted electron distribution function. The first validation of the ECCD theory in the extended regime having large Eφ and p is presented. Linear calculation has a tendency to overestimate the EC driven current, as normalized parameters for Eφ and p increase. While the EC driven current IEC obtained by a linearized Fokker-Planck calculation (1.1 MA) did not agree with the measured EC driven current (0.74 ± 0.06MA), non-linear calculation of the Fokker-Planck equation considering the effect of Eφ (0.76 MA) shows close agreement with the experimental result. Calculations show that the decreasing effect of Eφ on IEC was stronger than the increasing effect of p on IEC in the experiment.

A major research goal of the national spherical torus experiment is establishing long-pulse, high beta, high confinement operation and its physics basis. This research has been enabled by facility capabilities developed during 2001 and 2002, including neutral beam (up to 7?MW) and high harmonic fast wave (HHFW) heating (up to 6?MW), toroidal fields up to 6?kG, plasma currents up to 1.5?MA, flexible shape control, and wall preparation techniques. These capabilities have enabled the generation of plasmas with of up to 35%. Normalized beta values often exceed the no-wall limit, and studies suggest that passive wall mode stabilization enables this for H mode plasmas with broad pressure profiles. The viability of long, high bootstrap current fraction operations has been established for ELMing H mode plasmas with toroidal beta values in excess of 15% and sustained for several current relaxation times. Improvements in wall conditioning and fuelling are likely contributing to a reduction in H mode power thresholds. Electron thermal conduction is the dominant thermal loss channel in auxiliary heated plasmas examined thus far. HHFW effectively heats electrons, and its acceleration of fast beam ions has been observed. Evidence for HHFW current drive is obtained by comparision of the loop voltage evolution in plasmas with matched density and temperature profiles but varying phases of launched HHFW waves. Studies of emissions from electron Bernstein waves indicate a density scale length dependence of their transmission across the upper hybrid resonance near the plasma edge that is consistent with theoretical predictions. A peak heat flux to the divertor targets of 10?MW?m?2 has been measured in the H mode, with large asymmetries being observed in the power deposition between the inner and outer strike points. Non-inductive plasma startup studies have focused on coaxial helicity injection. With this technique, toroidal currents up to 400?kA have been driven, and studies to assess flux closure and coupling to other current drive techniques have begun.

Research on the spherical torus (or spherical tokamak) (ST) is being pursued to explore the scientific benefits of modifying the field line structure from that in more moderate aspect ratio devices, such as the conventional tokamak. The ST experiments are being conducted in various US research facilities including the MA-class National Spherical Torus Experiment (NSTX) at Princeton, and three medium sized ST research facilities: PEGASUS at University of Wisconsin, HIT-II at University of Washington, and CDX-U at Princeton. In the context of the fusion energy development path being formulated in the US, an ST-based Component Test Facility (CTF) and, ultimately a Demo device, are being discussed. For these, it is essential to develop high performance, steady-state operational scenarios. The relevant scientific issues are energy confinement, MHD stability at high beta (?), non-inductive sustainment, Ohmic-solenoid-free start-up, and power and particle handling. In the confinement area, the NSTX experiments have shown that the confinement can be up to 50% better than the ITER-98-pby2 H-mode scaling, consistent with the requirements for an ST-based CTF and Demo. In NSTX, CTF-relevant average toroidal beta values ?T of up to 35% with a near unity central ?T have been obtained. NSTX will be exploring advanced regimes where ?T up to 40% can be sustained through active stabilization of resistive wall modes. To date, the most successful technique for non-inductive sustainment in NSTX is the high beta poloidal regime, where discharges with a high non-inductive fraction (~60% bootstrap current+NBI current drive) were sustained over the resistive skin time. Research on radio-frequency (RF) based heating and current drive utilizing high harmonic fast wave and electron Bernstein wave is also pursued on NSTX, PEGASUS, and CDX-U. For non-inductive start-up, the coaxial helicity injection, developed in HIT/HIT-II, has been adopted on NSTX to test the method up to Ip ~ 500?kA. In parallel, start-up using a RF current drive and only external poloidal field coils are being developed on NSTX. The area of power and particle handling is expected to be challenging because of the higher power density expected in the ST relative to that in conventional aspect-ratio tokamaks. Due to its promise for power and particle handling, liquid lithium is being studied in CDX-U as a potential plasma-facing surface for a fusion reactor.

Recently a number of major, unanticipated effects have been reported in tokamak edge research raising the question of whether we understand the controlling physics of the edge. This report is on the first part – here focused on the outer divertor – of a systematic study of the simplest possible edge plasma – no ELMs, no detachment, etc. – for a set of 10 repeat, highly diagnosed, single-null, divertor discharges in DIII-D. For almost the entire, extensive data set so far evaluated, the matches of experiment and model are so close as to imply that the controlling processes at the outer divertor for these simple plasma conditions have probably been correctly identified and quantitatively characterized in the model. The principal anomaly flagged so far relates to measurements of Te near the target, potentially pointing to a deficiency in our understanding of sheath physics in the tokamak environment.

The National Spherical Torus Experiment (NSTX) at Princeton University, Princeton, NJ, is designed to assess the potential of the low-aspect-ratio spherical torus concept for magnetic plasma confinement. The plasma has been heated by up to 7 MW of neutral beam injection (NBI) at an injection energy of 100 keV and up to 6 MW of high harmonic fast wave (HHFW) at 30 MHz. NSTX has achieved /spl beta//sub T/ of 32%. A variety of MHD phenomena have been observed to limit /spl beta/. NSTX has now begun addressing /spl tau//sub E/ scaling, /spl beta/ limits, and current drive issues. During the NBI heating experiments, a broad T/sub i/ profile with T/sub i/ up to 2 keV, T/sub i/>T/sub e/ and a large toroidal rotation were observed. Transport analysis suggests that the impurity ions have diffusivities approaching neoclassical. For L-Mode plasmas, /spl tau//sub E/ is up to two times the ITER97L L-Mode scaling and exceeds the ITER98pby2 H-Mode scaling in some cases. Transitions to H-Mode have been observed which result in an approximate doubling of /spl tau//sub E/ after the transition in some conditions. During HHFW heating, T/sub e/>T/sub i/ and T/sub e/ up to 3.5 keV were observed. Current drive has been studied using both coaxial helicity injection with up to 390 kA of toroidal current and HHFW. HHFW has produced H-modes with significant bootstrap current fraction at low I/sub p/, high q, and high /spl beta//sub p/.

Electron cyclotron current drive (ECCD) experiments on the DIII-D tokamak are solidifying the physics basis for localized, off-axis current drive, the goal being to validate a predictive model for ECCD. The ECCD profiles are determined from the magnetic field pitch angles measured by motional Stark effect polarimetry. The measured ECCD switches from the co- to the counter-direction as the toroidal injection angle is varied with a profile width that is in accordance with ray tracing calculations. Tests of electron trapping in low beta plasmas show that the ECCD efficiency decreases rapidly as the deposition is moved off-axis and towards the outboard side of the plasma, but the detrimental effects of electron trapping on the current drive are greatly reduced in high beta plasmas. Overall, the measured ECCD is in good agreement with theoretical calculations using a quasilinear Fokker–Planck code over a wide range of injection angles and plasma parameters.

Some analyses that survey experimental data are carried out at a sparse sample rate between pulses during tokamak operation and/or completed as a batch job overnight because the complete analysis on a single fast workstation cannot fit in the narrow time window between two pulses. Scientists therefore miss the opportunity to use these results to guide experiments quickly. With a dedicated Beowulf type cluster at a cost less than that of a workstation, these analyses can be accomplished between pulses and the analyzed data made available for the research team during the tokamak operation. A Linux PC cluster comprised of 12 processors was installed at DIII-D National Fusion Facility in CY00 and expanded to 24 processors in CY01 to automatically perform between-pulse magnetic equilibrium reconstructions using the efit code written in fortran , cer analyses using cerquick code written in interactive data language ( idl) and full profile fitting analyses (ne, Te, Ti, Vr, Zeff) using idl code zipfit . This paper reports the current status of the system, the details of the between-pulse profile fitting analyses, and discusses some problems and concerns raised during the implementation and expansion of the system.

Ideal MHD stability limits of anticipated plasma configurations for NSTX [Ono M. et al 2000 Nucl. Fusion 40 557] and the dependence on the parameters defining the MHD equilibrium are evaluated. The study provides a quantitative computational evaluation of the stability limit variations induced by changes to the equilibrium of NSTX high β plasmas. The analysis is based on a reference free-boundary equilibrium with β = 41.5%, monotonic safety factor q profile (qa = 12.1, q0 = 2.8) and broad pressure profile p (peaking factor Fp≡p(0)/p = 1.7). On this reference target local variation of the plasma boundary, and the safety factor q and pressure profiles p are imposed. Localized inflection of the outboard plasma boundary, produced by near field effects from poloidal shaping field coils, weaken the stability due to the destabilization of high n ballooning modes. Variation of the q profile at different radial locations can also degrade stability. Both experimental profiles from existing tokamaks and spherical torus machines and profiles generated from transport modelling of anticipated neutral beam heated plasmas are used. Degraded stability is found at increasing pressure peaking factor due to the destabilization of n = 1 kink/ballooning modes. Direct access to the second region of stability is found in certain configurations and, for the entire set of variations considered, the lower calculated β limit values are still in the range of 20.0% without considering the stabilizing effect of the passive conducting structures.

Substantial progress has been made towards both understanding and control of internal transport barriers (ITBs) on DIII-D, resulting in the discovery of a new sustained high performance operating mode termed the quiescent double barrier (QDB) regime. The QDB regime combines core transport barriers with a quiescent ELM-free H mode edge (termed QH mode), giving rise to separate (double) core and edge transport barriers. The core and edge barriers are mutually compatible and do not merge, resulting in broad core profiles with an edge pedestal. The QH mode edge is characterized by ELM-free behaviour with continuous multiharmonic MHD activity in the pedestal region and has provided density and radiated power control for longer than 3.5 s (25τE) with divertor pumping. QDB plasmas are long pulse high performance candidates, having maintained a βN H89 product of 7 for five energy confinement times (Ti ≤ 16 keV, βN ≤ 2.9, H89 ≤ 2.4 τE ≤ 150 ms, DD neutron rate Sn ≤ 4 × 1015 s-1). The QDB regime has only been obtained in counter-NBI discharges (injection antiparallel to the plasma current) with divertor pumping. Other results include successful expansion of the ITB radius using (separately) both impurity injection and counter-NBI, and the formation of ITBs in the electron thermal channel using both ECH and strong negative central shear (NCS) at high power. These results are interpreted within a theoretical framework in which turbulence suppression is the key to ITB formation and control, and a decrease in core turbulence is observed in all cases of ITB formation.

Experimental and modeling investigations on the DIII-D and EAST tokamaks show the attractive transport and stability properties of fully noninductive, high poloidal-beta (β P ) plasmas, and their suitability for steady-state operating scenarios in ITER and CFETR. A key feature of the high-β P regime is the large-radius (ρ > 0.6) internal transport barrier (ITB), often observed in all channels (ne, Te, Ti, rotation), and responsible for both excellent energy confinement quality and excellent stability properties. Experiments on DIII-D have shown that, with a large-radius ITB, very high β N and β P values (both ≥ 4) can be reached by taking advantage of the stabilizing effect of a nearby conducting wall. Synergistically, higher plasma pressure provides turbulence suppression by Shafranov shift, leading to ITB sustainment independent of the plasma rotation. Experiments on EAST have been used to assess the long pulse potential of the high-β P regime. Using RF-only heating and current drive, EAST achieved minute-long fully noninductive steady state H-mode operation with strike points on an ITER-like tungsten divertor. Improved confinement (relative to standard H-mode) and steady state ITB features are observed with a monotonic q-profile with q min ~ 1.5. Separately, experiments have shown that increasing the density in plasmas driven by lower hybrid wave broadens the q-profile, a technique that could enable a large radius ITB. These experimental results have been used to validate MHD, current drive, and turbulent transport models, and to project the high-β P regime to a burning plasma. These projections suggest the Shafranov shift alone will not suffice to provide improved confinement (over standard H-mode) without rotation and rotation shear. However, increasing the negative magnetic shear (higher q on axis) provides a similar turbulence suppression mechanism to Shafranov shift, and can help devices such as ITER and CFETR achieve their steady-state fusion goals.

Transport modeling of a proposed ITER steady-state scenario based on DIII-D high poloidal-beta () discharges finds that ITB formation can occur with either sufficient rotation or a negative central shear q-profile. The high scenario is characterized by a large bootstrap current fraction (80%) which reduces the demands on the external current drive, and a large radius internal transport barrier which is associated with excellent normalized confinement. Modeling predictions of the electron transport in the high scenario improve as approaches levels similar to typical existing models of ITER steady-state and the ion transport is turbulence dominated. Typical temperature and density profiles from the non-inductive high scenario on DIII-D are scaled according to 0D modeling predictions of the requirements for achieving a steady-state fusion gain in ITER with 'day one' heating and current drive capabilities. Then, TGLF turbulence modeling is carried out under systematic variations of the toroidal rotation and the core q-profile. A high bootstrap fraction, high scenario is found to be near an ITB formation threshold, and either strong negative central magnetic shear or rotation in a high bootstrap fraction are found to successfully provide the turbulence suppression required to achieve .

An efficient technique for building plasma equilibria with preset or constrained q profiles by solving the nonlinear Grad–Shafranov equation iteratively is proposed. The technique uses a set of dynamic virtual poloidal field probes in the plasma core region. Numerical examples are shown for monotonic q profiles and negative central magnetic shear equilibria which are obtained with this technique. As a proof of principle, the technique is applied on the EAST tokamak, to constrain the plasma equilibrium reconstruction with the experimental q = 1 surface information from the ECE diagnostic in addition to the external magnetic measurements.

SOLPS-EIRENE edge code analysis shows that a gas-tight slot divertor geometry with a small-angle (glancing-incidence) target, named the small angle slot (SAS) divertor, can achieve cold, dissipative/detached divertor conditions at relatively low values of plasma density at the outside midplane separatrix. SAS exhibits the following key features: (1) strong enhancement of the buildup of neutral density in a localized region near the plasma strike point on the divertor target; (2) spreading of the cooling front across the divertor target with the slot gradually flaring out from the strike point, thus effectively reducing both heat flux and erosion on the entire divertor target surface. Such a divertor may potentially provide a power and particle handling solution for long pulse advanced tokamaks.

A new diagnostic has been developed on DIII-D that determines the impurity poloidal rotation from the poloidal asymmetry in the toroidal angular rotation velocity. This asymmetry is measured with recently added tangential charge exchange viewchords on the high-field side of the tokamak midplane. Measurements are made on co- and counter-current neutral beams, allowing the charge exchange cross section effect to be measured and eliminating the need for atomic physics calculations. The diagnostic implementation on DIII-D restricts the measurement range to the core (r/a < 0.6) where, relative to measurements made with the vertical charge exchange system, the spatial resolution is improved. Significant physics results have been obtained with this new diagnostic; for example, poloidal rotation measurements that significantly exceed neoclassical predictions.

Energy and pitch angle resolved measurements of escaping neutral beam ions (E approximate to 80 keV) have been made during DIII-D L-mode discharges with applied, slowly rotating, n = 2 magnetic perturbations. Data from separate scintillator detectors (FILDs) near and well below the plasma midplane show fast-ion losses correlated with the internal coil (I-coil) fields. The dominant fast-ion loss signals are observed to decay within one poloidal transit time after beam turn-off indicating they are primarily prompt loss orbits. Also, during application of the rotating I-coil fields, outboard midplane edge density and bremsstrahlung emission profiles exhibit a radial displacement of up to delta R approximate to 1 cm. Beam deposition and full orbit modeling of these losses using M3D-C1 calculations of the perturbed kinetic profiles and fields reproduce many features of the measured losses. In particular, the predicted phase of the modulated loss signal with respect to the I-coil currents is in close agreement with FILD measurements as is the relative amplitudes of the modulated losses for the co and counter-current beam used in the experiment. These simulations show modifications to the beam ion birth profile and subsequent prompt loss due to changes in the edge density; however, the dominant factor causing modulation of the losses to the fast-ion loss detectors is the perturbed magnetic field (delta B/B approximate to 10(-3) in the plasma). Calculations indicate total prompt loss to the DIII-D wall can increase with application of the n = 2 perturbation by up to 7% for co-current injected beams and 3% for counter-current injected beams depending on phase of the perturbation relative to the injected beam. US Department of Energy DE-FC02-04ER54698, SC-G903402, DEAC02- 09CH11466, DE-FG02-04ER54761, DE-FG02- 05ER54808

Edge conditions in DIII-D are being quantified in order to provide insight into the physics of the H?mode regime. Several studies show that electron temperature is not the key parameter that controls the L-H transition. Gradients of edge temperature and pressure are much more promising candidates for elements of such parameters. They systematically increase during the L phases of discharges which make a transition to H?mode, and these increases are typically larger than the increases in the underlying quantities. The quality of H?mode confinement is strongly correlated with the height of the H?mode pedestal for the pressure. The gradient of the pressure is limited by MHD modes, in particular by ideal kink ballooning modes with finite mode number n. For a wide variety of discharges, the width of the barrier for electron pressure is well described by a relationship that is proportional to (?pedp)1/2. A new regime of confinement, called the quiescent H?mode, which provides steady state operation with no ELMs, low radiated power and normal H?mode confinement, has been discovered. A coherent edge MHD mode provides adequate particle transport to control the plasma density while permitting the pressure pedestal to remain almost identical to that observed in ELMing discharges.

Significant progress in obtaining high performance discharges lasting many energy confinement times in the DIII-D tokamak has been realized in recent experimental campaigns. Normalized performance ∼10 has been sustained for more than 5>E with qmin > 1.5. (The normalized performance is measured by the product βNH89, indicating the proximity to the conventional β limits and energy confinement quality, respectively.) These H mode discharges have an ELMing edge and β < 5%. The limit to increasing β is a resistive wall mode, rather than the tearing modes as previously observed. Confinement remains good despite qmin > 1. The global parameters were chosen to optimize the potential for fully non-inductive current sustainment at high performance, which is a key program goal for the DIII-D facility. Measurement of the current density and loop voltage profiles indicate that ≈75% of the current in the present discharges is sustained non-inductively. The remaining ohmic current is localized near the half-radius. The electron cyclotron heating system is being upgraded to replace this remaining current with ECCD. Density and β control, which are essential for operating advanced tokamak discharges, were demonstrated in ELMing H mode discharges with βNH89 ≈ 7 for up to 6.3 s or ≈34τE. These discharges appear to have stationary current profiles with qmin ≈ 1.05 in agreement with the current profile relaxation time ≈1.8 s.

The main aim of the National Spherical Torus Experiment (NSTX) is to establish the fusion physics principles of the spherical torus (ST) concept. The NSTX device began plasma operations in February 1999 and the plasma current Ip was successfully brought up to the design value of 1 MA on 14 December 1999. The planned plasma shaping parameters, elongation κ = 1.6-2.2 and triangularity δ = 0.2-0.4, were achieved in inner wall limited, and single null and double null diverted configurations. The coaxial helicity injection (CHI) and high harmonic fast wave (HHFW) experiments were also initiated. CHI current of 27 kA produced up to 260 kA toroidal current without using an ohmic solenoid. With the injection of 2.3 MW of HHFW power, using 12 antennas connected to six transmitters, electrons were heated from a central temperature of 400 eV to 900 eV at a central density of 3.5 × 1013 cm-3, increasing the plasma energy to 59 kJ and the toroidal β, βT, to 10%. The NBI system commenced operation in September 2000. The initial results with two ion sources (PNBI = 2.8 MW) show good heating, producing a total plasma stored energy of 90 kJ corresponding to βT ≈ 18% at a plasma current of 1.1 MA.

Coaxial helicity injection (CHI) on the National Spherical Torus Experiment (NSTX) has produced 240 kA of toroidal current without the use of the central solenoid. Values of the current multiplication ratio (CHI produced toroidal current/injector current) up to 10 were obtained, in agreement with predictions. The discharges, which lasted for up to 200 ms, limited only by the programmed waveform, are more than an order of magnitude longer in duration than any CHI discharges previously produced in a spheromak or a spherical torus.

Scrape-off layer plasma simulation modeling has been carried out to assess the effect of tightly closing the lower divertor in DIII-D, which at present is almost fully open, on the achievement of cold dissipative/detached divertor conditions. To isolate the impact of other factors on the divertor plasma solution and to make direct comparisons, most of the parameters including the meshes were kept as similar as possible. Only the neutral baffling was modified to compare a fully open divertor with a tightly closed one. The modeling shows that the tightly closed divertor greatly improves trapping of recycling neutrals, thereby increasing radiative and charge exchange losses in the divertor and reducing the electron temperature T et and deposited power density q dep at the target plate. Furthermore, the closed structure enables the divertor plasma to enter into highly dissipative and detached divertor conditions at a significantly lower upstream density. The effects of divertor closure on the neutral density and pressure, and their correlation with the divertor plasma conditions are also demonstrated. The effect of molecular D2−ion D+ elastic collisions and neutral–neutral collisions on the divertor plasma solution are assessed.

Recent progress in the development of long pulse, high performance discharges on the DIII-D tokamak is reviewed. This is highlighted by a discharge achieving simultaneously a ?N H of 9, bootstrap current fraction of 0.5 and non-inductive current fraction of 0.75, sustained for 16 energy confinement times. The physics challenge has changed in the long pulse regime. Non-ideal MHD modes limit the stability, fast ion driven modes may play a role in fast ion transport, which limits the stored energy, and plasma edge behaviour can affect the global performance. New control tools are being developed to address these issues.

The characteristics of the H-mode are studied in discharges with varying triangularity and squareness. The pressure at the top of the H-mode pedestal increases strongly with triangularity, primarily due to an increase in the margin by which the edge pressure gradient exceeds the ideal ballooning mode first stability limit. Two models are considered for how the edge may exceed the ballooning mode limit. In one model, access to the ballooning mode second stable regime allows the edge pressure gradient and associated bootstrap current to continue to increase until an edge localized, low toroidal mode number, ideal kink mode is destabilized. In the second model, the finite width of the H-mode transport barrier and diamagnetic effects raise the pressure gradient limit above the ballooning mode limit. We observe a weak inverse dependence of the width of the H-mode transport barrier, Δ, on triangularity relative to the previously obtained scaling Δ∝(βPEP)1/2. The energy loss for Type I ELMs increases with triangularity in proportion to the pedestal energy increase. At low density, the energy confinement of high-triangularity discharges is larger than discharges with low triangularity, as a result of an increase in the energy in the H-mode pedestal. At high density, both the change in pedestal pressure and the response of the density profile are found to play a role in setting the energy confinement. The highest energy confinement at high density was obtained in low-triangularity discharges.

Two approaches to achieving long timescale stabilization of the ideal kink mode with a real, finite conductivity wall are considered: plasma rotation and active feedback control. DIII-D experiments have demonstrated stabilization of the resistive wall mode (RWM) by sustaining ? greater than the no-wall limit for up to 200?ms, much longer than the wall penetration time of a few milliseconds. These plasmas are typically terminated by an m = 3, n = 1 mode as the plasma rotation slows below a few kilohertz. Recent temperature profile data show an ideal MHD mode structure, as expected for the RWM at ? above the no-wall limit. The critical rotation rate for stabilization is in qualitative agreement with recent theories for dissipative stabilization in the absence of magnetic islands. However, drag by small amplitude RWMs or damping of stable RWMs may contribute to an observed slowing of rotation at high ?, rendering rotational stabilization more difficult. An initial open loop active control experiment, using non-axisymmetric external coils and a new array of saddle loop detectors, has yielded encouraging results indicating a delayed onset of the RWM.

Improving confinement and β limits simultaneously in long pulse ELMy H mode discharges is investigated. The product βN H98y serves as a useful figure of merit for performance, where βN ≡ β/(I/aB) and H98y is the ratio of the thermal confinement time to the most recent ELMy H mode confinement scaling established by the ITER confinement database working group. In discharges with q0 ≈ 1 (no sawteeth) and discharges with qmin > 1.5 and negative central magnetic shear, βN ≈ 2.9 and H98y ≈ 1.4 are sustained for up to 2 s. Although peaked profiles are observed, steep internal transport barriers are not present. Further increases in βN in these discharges are limited by neoclassical tearing modes (NTMs) in the positive shear region. In another recently developed regime, βN ≈ 3.8 and H98y ≈ 2 have been sustained during large infrequent ELMs in non-sawtoothing discharges with q0 ≈ 1. This level of performance is similar to that obtained in ELM free regimes such as VH mode. The limitation on βN and pulse length in these discharges is also the onset of NTMs.