Search

The ITER Research Plan envision operation around half of the nominal magnetic field (i.e. around B = 2.65 T) as a path to baseline operation. This work discusses constraints on the optimal range of magnetic field, which is bounded in the lower limit by the presence of the third-harmonic electron cyclotron resonance at half field, and on the upper limit by the loss of core heating and current drive. It will be shown that increasing the magnetic field by only 3%, i.e. to 2.75 T, eliminates the third harmonic parasitic absorption without compromising demonstration of access to H-mode, while operating at a magnetic field of 3.0 T - previously proposed for optimal use of the ion cyclotron system - would impair the use of the electron cyclotron system for core-heating and current drive. Operation at 2.65 T would still be possible if the polarization of the equatorial launcher is changed from X-mode to O-mode in the current flattop phase.

Pedestal collapse (i.e., the complete loss of the edge transport barrier (ETB)) in DIII-D H-mode plasmas occurs when resonant magnetic perturbations (RMPs) penetrate the steep gradient region at the plasma edge. Normally, RMP driven magnetic islands can occur at the top and bottom of the H-mode pedestal and these islands generate conditions consistent with edge-localized-mode (ELM) suppression and density pump-out, respectively, based on nonlinear two-fluid MHD simulations. In contrast, MHD simulations show that the steep pressure gradient region between the top and bottom of the DIII-D pedestal is generally immune to resonant field penetration due to large local E × B and diamagnetic flows. By this fortuitous circumstance, the edge-transport-barrier and H-mode confinement can be maintained while achieving ELM suppression. However, pedestal collapse can occur in DIII-D when the screening flows are inadequate to prevent field penetration in the steep gradient region of the pedestal. Non-linear two-fluid MHD simulations support the role of resonant field penetration in pedestal collapse for DIII-D H-mode plasmas with weak edge E × B and diamagnetic screening flows. ITER will likely have weaker edge screening flows than present experiments due to its much larger size, making it more susceptible to resonant field penetration in the steep gradient region of the pedestal. Analysis of model ITER equilibria demonstrates that resonant field penetration in the steep pressure gradient region is possible for RMP levels of the order required for ELM suppression. The effect of such penetration on the ITER pedestal will depend sensitively on the resulting degree of island overlap.

The local properties of the electron energy distribution function are investigated in the electron-cyclotron wave driven toroidal device TORPEX (Fasoli et al 2006 Phys. Plasmas 13 055902). The electron energy distribution function is measured by means of a time-resolved electrostatic energy analyser. The measurements indicate that electrons are accelerated to different energies at the electron-cyclotron and upper-hybrid resonances. Electrons from the upper-hybrid resonance have higher energies and can be described by a secondary maxwellian energy distribution function. The electron energy distribution function is found to be invariant with respect to experimental parameters such as the neutral gas pressure, the magnetic field configuration and the injected microwave power. The implications for a correct modelling of the particle source term for TORPEX plasmas are discussed.

The plasma particle source resulting from ionization by microwaves in the electron-cyclotron range of frequencies is investigated in the simple magnetized torus TORPEX ( Fasoli A et al 2006 Phys. Plasmas 13 055902), together with its dependence upon experimental parameters, such as the absorbed microwave power, the magnetic field configuration and the neutral gas pressure. The relative contributions to the ionization rate from thermal electrons and from suprathermal electrons accelerated at the plasma resonances are quantified through a global particle balance. An expression for the spatial profile of the particle source is derived from the experimental results and a simple numerical code.

Time-dependent simulations are used to evolve plasma discharges in combination with a modified Rutherford equation for calculation of neoclassical tearing mode (NTM) stability in response to electron cyclotron (EC) feedback control in ITER. The main application of this integrated approach is to support the development of control algorithms by analyzing the plasma response with physics-based models and to assess how uncertainties in the detection of the magnetic island and in the EC alignment affect the ability of the ITER EC system to fulfill its purpose. Simulations indicate that it is critical to detect the island as soon as possible, before its size exceeds the EC deposition width, and that maintaining alignment with the rational surface within half of the EC deposition width is needed for stabilization and suppression of the modes, especially in the case of modes with helicity (2, 1). A broadening of the deposition profile, for example due to wave scattering by turbulence fluctuations or not well aligned beams, could even be favorable in the case of the (2, 1)-NTM, by relaxing an over-focussing of the EC beam and improving the stabilization at the mode onset. Pre-emptive control reduces the power needed for suppression and stabilization in the ITER baseline discharge to a maximum of 5 MW, which should be reserved and available to the upper launcher during the entire flattop phase. Assuming continuous triggering of NTMs, with preemptive control ITER would be still able to demonstrate a fusion gain of Q = 10.

Non-activation phase operations in ITER in hydrogen (H) and helium (He) will be important for commissioning of tokamak systems, such as diagnostics, heating and current drive (HCD) systems, coils and plasma control systems, and for validation of techniques necessary for establishing operations in DT. The assessment of feasible HCD schemes at various toroidal fields (2.65–5.3 T) has revealed that the previously applied assumptions need to be refined for the ITER non-activation phase H/He operations. A study of the ranges of plasma density and profile shape using the JINTRAC suite of codes has indicated that the hydrogen pellet fuelling into He plasmas should be utilized taking the optimization of IC power absorption, neutral beam shine-through density limit and H-mode access into account. The EPED1 estimation of the edge pedestal parameters has been extended to various H operation conditions, and the combined EPED1 and SOLPS estimation has provided guidance for modelling the edge pedestal in H/He operations. The availability of ITER HCD schemes, ranges of achievable plasma density and profile shape, and estimation of the edge pedestal parameters for H/He plasmas have been integrated into various time-dependent tokamak discharge simulations. In this work, various H/He scenarios at a wide range of plasma current (7.5–15 MA) and field (2.65–5.3 T) have been developed for the ITER non-activation phase operation, and the sensitivity of the developed scenarios to the used assumptions has been investigated to provide guidance for further development.

Author(s): Menard, JE; Allain, JP; Battaglia, DJ; Bedoya, F; Bell, RE; Belova, E; Berkery, JW; Boyer, MD; Crocker, N; Diallo, A; Ebrahimi, F; Ferraro, N; Fredrickson, E; Frerichs, H; Gerhardt, S; Gorelenkov, N; Guttenfelder, W; Heidbrink, W; Kaita, R; Kaye, SM; Kriete, DM; Kubota, S; Leblanc, BP; Liu, D; Lunsford, R; Mueller, D; Myers, CE; Ono, M; Park, JK; Podesta, M; Raman, R; Reinke, M; Ren, Y; Sabbagh, SA; Schmitz, O; Scotti, F; Sechrest, Y; Skinner, CH; Smith, DR; Soukhanovskii, V; Stoltzfus-Dueck, T; Yuh, H; Wang, Z; Waters, I; Ahn, JW; Andre, R; Barchfeld, R; Beiersdorfer, P; Bertelli, N; Bhattacharjee, A; Brennan, D; Buttery, R; Capece, A; Canal, G; Canik, J; Chang, CS; Darrow, D; Delgado-Aparicio, L; Domier, C; Ethier, S; Evans, T; Ferron, J; Finkenthal, M; Fonck, R; Gan, K; Gates, D; Goumiri, I; Gray, T; Hosea, J; Humphreys, D; Jarboe, T; Jardin, S; Jaworski, MA; Koel, B; Kolemen, E; Ku, S; La Haye, RJ; Levinton, F; Luhmann, N; Maingi, R; Maqueda, R; McKee, G; Meier, E; Myra, J; Perkins, R | Abstract: The National Spherical Torus Experiment (NSTX) has undergone a major upgrade, and the NSTX Upgrade (NSTX-U) Project was completed in the summer of 2015. NSTX-U first plasma was subsequently achieved, diagnostic and control systems have been commissioned, the H-mode accessed, magnetic error fields identified and mitigated, and the first physics research campaign carried out. During ten run weeks of operation, NSTX-U surpassed NSTX record pulse-durations and toroidal fields (TF), and high-performance ∼1 MA H-mode plasmas comparable to the best of NSTX have been sustained near and slightly above the n = 1 no-wall stability limit and with H-mode confinement multiplier H98y,2 above 1. Transport and turbulence studies in L-mode plasmas have identified the coexistence of at least two ion-gyro-scale turbulent micro-instabilities near the same radial location but propagating in opposite (i.e. ion and electron diamagnetic) directions. These modes have the characteristics of ion-temperature gradient and micro-tearing modes, respectively, and the role of these modes in contributing to thermal transport is under active investigation. The new second more tangential neutral beam injection was observed to significantly modify the stability of two types of Alfven eigenmodes. Improvements in offline disruption forecasting were made in the areas of identification of rotating MHD modes and other macroscopic instabilities using the disruption event characterization and forecasting code. Lastly, the materials analysis and particle probe was utilized on NSTX-U for the first time and enabled assessments of the correlation between boronized wall conditions and plasma performance. These and other highlights from the first run campaign of NSTX-U are described.

This paper examines a method for real-time control of non-inductively sustained scenarios in NSTX-U by using TRANSP, a time-dependent integrated modeling code for prediction and interpretive analysis of tokamak experimental data, as a simulator. The actuators considered for control in this work are the six neutral beam sources and the plasma boundary shape. To understand the response of the plasma current, stored energy, and central safety factor to these actuators and to enable systematic design of control algorithms, simulations were run in which the actuators were modulated and a linearized dynamic response model was generated. A multi-variable model-based control scheme that accounts for the coupling and slow dynamics of the system while mitigating the effect of actuator limitations was designed and simulated. Simulations show that modest changes in the outer gap and heating power can improve the response time of the system, reject perturbations, and track target values of the controlled values.

High density plasmas (n0 ≈ 8 × 1020m-3) featuring steady improved core confinement have been obtained in FTU up to the maximum nominal toroidal field (8 T) by deep multiple pellet injection. These plasmas also feature high purity efficient electron-ion coupling and peaked density profiles sustained for several confinement times. Neutron yields in excess of 1 × 1013 n/s are measured, consistent with the reduction of the ion transport to neoclassical levels.

This work discusses sources of uncertainty in the validation of lower hybrid wave current drive simulations against experiments, by evolving self-consistently the magnetic equilibrium and the heating and current drive profiles, calculated with a combined toroidal ray tracing code and 3D Fokker–Planck solver. The simulations indicate a complex interplay of elements, where uncertainties in the input plasma parameters, in the models and in the transport solver combine and—in some cases—compensate each other. It is concluded that ray-tracing calculations should include a realistic representation of the density and temperature in the region between the confined plasma and the wall, which is especially important in regimes where the LH waves are weakly damped and undergo multiple reflections from the plasma boundary. Uncertainties introduced in the processing of diagnostic data as well as uncertainties introduced by model approximations are assessed. It is shown that, by comparing the evolution of the plasma parameters in self-consistent simulations with available data, inconsistencies can be identified and limitations in the models or in the experimental data assessed.

One of the goals of the National Spherical Torus Experiment Upgrade (NSTX-U) (Menard et al 2012 Nucl. Fusion 52 083015) is the demonstration of fully non-inductive start-up, current ramp-up and sustainment. This work discusses predictive simulations where the available heating and current drive systems are combined to maximize the non-inductive current and minimize the solenoidal contribution. Radio-frequency waves at harmonics higher than the ion cyclotron resonance (high-harmonic fast waves (HHFW)) and neutral beam injection are used to ramp the plasma current non-inductively starting from an initial Ohmic plasma. An interesting synergy is observed in the simulations between the HHFW and electron cyclotron (EC) wave heating. Furthermore, time-dependent simulations indicate that, depending on the phasing of the HHFW antenna, EC wave heating can significantly increase the effectiveness of the radio-frequency power, by heating the electrons and increasing the current drive efficiency, thus relaxing the requirements on the level of HHFW power that needs to be absorbed in the core plasma to drive the same amount of fast-wave current.

The MIT Plasma Science and Fusion Center and collaborators are proposing a high-performance Advanced Divertor and RF tokamak eXperiment (ADX)-a tokamak specifically designed to address critical gaps in the world fusion research programme on the pathway to next-step devices: fusion nuclear science facility (FNSF), fusion pilot plant (FPP) and/or demonstration power plant (DEMO). This high-field (>= 6.5 T, 1.5 MA), high power density facility (P/S similar to 1.5 MW m(-2)) will test innovative divertor ideas, including an 'X-point target divertor' concept, at the required performance parameters-reactor-level boundary plasma pressures, magnetic field strengths and parallel heat flux densities entering into the divertor region-while simultaneously producing high-performance core plasma conditions that are prototypical of a reactor: equilibrated and strongly coupled electrons and ions, regimes with low or no torque, and no fuelling from external heating and current drive systems. Equally important, the experimental platform will test innovative concepts for lower hybrid current drive and ion cyclotron range of frequency actuators with the unprecedented ability to deploy launch structures both on the low-magnetic-field side and the high-magneticfield side-the latter being a location where energetic plasma-material interactions can be controlled and favourable RF wave physics leads to efficient current drive, current profile control, heating and flow drive. This triple combination-advanced divertors, advanced RF actuators, reactor-prototypical core plasma conditions-will enable ADX to explore enhanced core confinement physics, such as made possible by reversed central shear, using only the types of external drive systems that are considered viable for a fusion power plant. Such an integrated demonstration of high-performance core-divertor operation with steady-state sustainment would pave the way towards an attractive pilot plant, as envisioned in the ARC concept (affordable, robust, compact) (Sorbom et al 2015 Fusion Eng. Des. submitted (arXiv: 1409.3540)) that makes use of high-temperature superconductor technology-a high-field (9.25 T) tokamak the size of the Joint European Torus that produces 270 MW of net electricity.

Steady state scenarios envisaged for ITER aim at optimizing the bootstrap current, while maintaining sufficient confinement and stability to provide the necessary fusion yield. Non-inductive scenarios will need to operate with internal transport barriers (ITBs) in order to reach adequate fusion gain at typical currents of 9 MA. However, the large pressure gradients associated with ITBs in regions of weak or negative magnetic shear can be conducive to ideal MHD instabilities, reducing the no-wall limit. The E × B flow shear from toroidal plasma rotation is expected to be low in ITER, with a major role in the ITB dynamics being played by magnetic geometry. Combinations of heating and current drive (H/CD) sources that sustain reversed magnetic shear profiles throughout the discharge are the focus of this work. Time-dependent transport simulations indicate that a combination of electron cyclotron (EC) and lower hybrid (LH) waves is a promising route towards steady state operation in ITER. The LH forms and sustains expanded barriers and the EC deposition at mid-radius freezes the bootstrap current profile stabilizing the barrier and leading to confinement levels 50% higher than typical H-mode energy confinement times. Using LH spectra with spectrum centred on parallel refractive index of 1.75–1.85, the performance of these plasma scenarios is close to the ITER target of 9 MA non-inductive current, global confinement gain H98 = 1.6 and fusion gain Q = 5.

Experiments on Alcator C-Mod from 2009?2012 have examined the rampup, flattop and rampdown phases of the proposed ITER 15?MA baseline scenario. Rampup studies show ICRF heating can significantly reduce the V-s requirement, and that an H-mode late in the ramp can reduce this further. ICRF modifications to li in L-mode are minimal, although the Te profile is peaked relative to ohmic in the plasma centre, and reduces the sawtooth onset times. Flattop plasmas targeting ITER baseline parameters have been sustained for 20?E or 8???13?CR, but only reach H98???0.6 at n/nGr?=?0.85, rising to 0.9 at n/nGr?=?0.65. Rampdown studies show H-modes can be routinely sustained with ICRF power injection, avoiding an OH coil over-current associated with the H?L transition. In addition, faster current rampdowns are preferred to avoid an over-current when an H?L transition ultimately does occur. In the H-mode rampdown the density is found to drop with Ip, preserving the n/nGr ratio, so long as ICRF power is injected.

Non-inductive steady-state scenarios on ITER will need to operate with internal transport barriers (ITBs) in order to reach adequate fusion gain at typical currents of 9 MA. The large pressure gradients at the location of the internal barrier are conducive to the development of ideal MHD instabilities that may limit the plasma performance and may lead to plasma disruptions. Fully non-inductive scenario simulations with five combinations of heating and current drive sources are presented in this work, with plasma currents in the range 7–10 MA. For each configuration the linear, ideal MHD stability is analysed for variations of the Greenwald fraction and of the pressure peaking factor around the operating point, aiming at defining an operational space for stable, steady-state operations at optimized performance. It is shown that plasmas with lower hybrid heating and current drive maintain the minimum safety factor above 1.5, which is desirable in steady-state operations to avoid neoclassical tearing modes. Operating with moderate ITBs at 2/3 of the minor radius, these plasmas have a minimum safety factor above 2, are ideal MHD stable and reach Q ≳ 5 operating above the ideal no-wall limit.

Electrostatic turbulence, related structures and their effect on particle, heat and toroidal momentum transport are investigated in TORPEX simple magnetized plasmas using high-resolution diagnostics, control parameters, linear fluid models and nonlinear numerical simulations. The nature of the dominant instabilities is controlled by the value of the vertical magnetic field, Bv, relative to that of the toroidal field, BT. For Bv/BT > 3%, only ideal interchange instabilities are observed. A critical pressure gradient to drive the interchange instability is experimentally identified. Interchange modes give rise to blobs, radially propagating filaments of enhanced plasma pressure. Blob velocities and sizes are obtained from electrostatic probe measurements using pattern recognition methods. The observed values span a wide range and are described by a single analytical expression, from the small blob size regime in which the blob velocity is limited by cross-field ion polarization currents, to the large blob size regime in which the limitation to the blob velocity comes from parallel currents to the sheath. As a first attempt at controlling the blob dynamical properties, limiter configurations with varying angles between field lines and the conducting surface of the limiter are explored. Mach probe measurements clearly demonstrate a link between toroidal flows and blobs. To complement probe data, a fast framing camera and a movable gas puffing system are installed. Density and light fluctuations show similar signatures of interchange activity. Further developments of optical diagnostics, including an image intensifier and laser-induced fluorescence, are under way. The effect of interchange turbulence on fast ion phase space dynamics is studied using movable fast ion source and detector in scenarios for which the development from linear waves into blobs is fully characterized. A theory validation project is conducted in parallel with TORPEX experiments, based on quantitative comparisons of observables that are defined in the same way in the data and in the output of numerical codes, including 2D and 3D local and global simulations.