The behaviour of impurities in fusion plasmas is of crucial importance for achieving sustained fusion reactions, and understanding similarities and differences between Deuterium (D) and Tritium (T) plasmas is needed to assess potential changes from DD to DT in ITER and future reactors. The first dimensionless and dimensional isotope identity experiments between Deuterium (D) and Tritium (T) L-mode plasmas were conducted at the JET W/Be wall. In the first approach, the discharges with matched rho & lowast;, nu & lowast;, beta n, q, and Te/Ti were compared to emphasize direct isotope effects, while in the dimensional approach engineering parameters such as toroidal magnetic field BT, plasma current Ip, plasma electron density and NBI power PNBI were matched. The dimensionless isotope scaling showed an improvement in global confinement and local transport in T plasmas in comparison to the matched D one (Cordey et al 1999 Nucl. Fusion 39 301). Detailed impurity analyses using VUV, visible spectroscopy, SXR cameras, and bolometry revealed that T plasmas exhibited higher radiation and impurity content, particularly Ni and W, compared to D plasmas. Understanding the origin of the increased impurity content is addressed in this paper. The dimensionless experiments showed differences in impurity transport. The Be source behaviour varied: D plasmas had higher Be influx in the dimensionless approach due to lower electron density and enhanced sputtering (Saibene et al 1999 Nucl. Fusion 39 1133), while T plasmas showed a higher Be source in the dimensional experiments, highlighting isotope mass effects. W in the divertor region was not sputtered by hydrogen isotopes. W in the divertor region was not sputtered by hydrogen isotopes. In the dimensionless experiments, W sputtering was primarily influenced by Ni in T plasmas and by Be in D plasmas. However, in the dimensional approach, Be played a more significant role in W sputtering within T plasmas. MHD instabilities, including ST oscillations, were present in all cases other ones were correlated with NBI power levels; higher NBI power led to elevated levels of Be, Ni, and W impurities. The comprehensive comparison underscores the necessity of accounting for isotope mass effects in predictive modelling and optimization of plasma performance in fusion reactors.
Tokamak a configuration variable (TCV), recently celebrating 30 years of near-continual operation, continues in its missions to advance outstanding key physics and operational scenario issues for ITER and the design of future power plants such as DEMO. The main machine heating systems and operational changes are first described. Then follow five sections: plasma scenarios. ITER Base-Line (IBL) discharges, triangularity studies together with X3 heating and N2 seeding. Edge localised mode suppression, with a high radiation region near the X-point is reported with N-2 injection with and without divertor baffles in a snowflake configuration. Negative triangularity (NT) discharges attained record, albeit transient, beta(N) similar to 3 with lower turbulence, higher low-Z impurity transport, vertical stability and density limits and core transport better than the IBL. Positive triangularity L-Mode linear and saturated ohmic confinement confinement saturation, often-correlated with intrinsic toroidal rotation reversals, was probed for D, H and He working gases. H-mode confinement and pedestal studies were extended to low collisionality with electron cyclotron heating obtaining steady state electron iternal transport barrier with neutral beam heating (NBH), and NBH driven H-mode configurations with off-axis co-electron cyclotron current drive. Fast particle physics. The physics of disruptions, runaway electrons and fast ions (FIs) was developed using near-full current conversion at disruption with recombination thresholds characterised for impurity species (Ne, Ar, Kr). Different flushing gases (D2, H2) and pathways to trigger a benign disruption were explored. The 55 kV NBH II generated a rich Alfvenic spectrum modulating the FI fas ion loss detector signal. NT configurations showed less toroidal Alfven excitation activity preferentially affecting higher FI pitch angles. Scrape-off layer and edge physics. gas puff imaging systems characterised turbulent plasma ejection for several advanced divertor configurations, including NT. Combined diagnostic array divertor state analysis in detachment conditions was compared to modelling revealing an importance for molecular processes. Divertor physics. Internal gas baffles diversified to include shorter/longer structures on the high and/or low field side to probe compressive efficiency. Divertor studies concentrated upon mitigating target power, facilitating detachment and increasing the radiated power fraction employing alternative divertor geometries, optimised X-point radiator regimes and long-legged configurations. Smaller-than-expected improvements with total flux expansion were better modelled when including parallel flows. Peak outer target heat flux reduction was achieved (>50%) for high flux-expansion geometries, maintaining core performance (H-98 > 1). A reduction in target heat loads and facilitated detachment access at lower core densities is reported. Real-time control. TCV's real-time control upgrades employed MIMO gas injector control of stable, robust, partial detachment and plasma beta feedback control avoiding neoclassical tearing modes with plasma confinement changes. Machine-learning enhancements include trajectory tracking disruption proximity and avoidance as well as a first-of-its-kind reinforcement learning-based controller for the plasma equilibrium trained entirely on a free-boundary simulator. Finally, a short description of TCV's immediate future plans will be given.
Superconducting (SC) tokamak JT-60SA plays an essential role in fusion research and development by supporting and complementing the ITER project, providing directions to the DEMO design activity and fostering next generation scientists and engineers. Since the short circuit incident at the terminal joints of equilibrium field coil #1 during the integrated commissioning (IC) in March 2021, both EU and JA implementing agencies (IAs) have examined how to ensure safe operation of JT-60SA by mitigating the risk of possible discharge occurrence inside the cryostat. Based on the experience of the global Paschen tests, the IAs have established a strategy of risk mitigation measures, which is a combination of (i) reinforcement of insulation, (ii) avoiding unnecessary voltage application to the coil systems and (iii) immediate de-energization of the coils when deteriorated vacuum conditions are detected. Thanks to the considerable efforts of the Integrated Project Team members, the IC restarted in May 2023. After confirmation of the SC state of the coil systems (TF, EF and CS), the coil energization test and the plasma operation phase 1 (OP-1) started. The first plasma was successfully achieved on 23 October 2023 with a limited value of voltage and current applied to the coils. The plasma configuration control was also confirmed with low plasma current and low auxiliary heating power conditions. Based on the IO-F4E-QST collaboration, activities of JT-60SA have been shared with the IO and provided an important lesson for ITER assembly and commissioning, and will provide an outstanding contribution to fusion research at large. After OP-1, maintenance & enhancement phase 1 (M/E-1) starts from January 2024, in which in-vessel components are installed, and heating and diagnostic systems are extensively upgraded to allow a high power heating experiment planned in OP-2. In order to make the best use of JT-60SA, a newly organized JT-60SA experiment team will refine the research plan for the future high heating power operation phase.
In 2021 JET exploited its unique capabilities to operate with T and D-T fuel with an ITER-like Be/W wall (JET-ILW). This second major JET D-T campaign (DTE2), after DTE1 in 1997, represented the culmination of a series of JET enhancements-new fusion diagnostics, new T injection capabilities, refurbishment of the T plant, increased auxiliary heating, in-vessel calibration of 14 MeV neutron yield monitors-as well as significant advances in plasma theory and modelling in the fusion community. DTE2 was complemented by a sequence of isotope physics campaigns encompassing operation in pure tritium at high T-NBI power. Carefully conducted for safe operation with tritium, the new T and D-T experiments used 1 kg of T (vs 100 g in DTE1), yielding the most fusion reactor relevant D-T plasmas to date and expanding our understanding of isotopes and D-T mixture physics. Furthermore, since the JET T and DTE2 campaigns occurred almost 25 years after the last major D-T tokamak experiment, it was also a strategic goal of the European fusion programme to refresh operational experience of a nuclear tokamak to prepare staff for ITER operation. The key physics results of the JET T and DTE2 experiments, carried out within the EUROfusion JET1 work package, are reported in this paper. Progress in the technological exploitation of JET D-T operations, development and validation of nuclear codes, neutronic tools and techniques for ITER operations carried out by EUROfusion (started within the Horizon 2020 Framework Programme and continuing under the Horizon Europe FP) are reported in (Litaudon et al Nucl. Fusion accepted), while JET experience on T and D-T operations is presented in (King et al Nucl. Fusion submitted).
Within the 9th European Framework programme, since 2021 EUROfusion is operating five tokamaks under the auspices of a single Task Force called ‘Tokamak Exploitation’. The goal is to benefit from the complementary capabilities of each machine in a coordinated way and help in developing a scientific output scalable to future largre machines. The programme of this Task Force ensures that ASDEX Upgrade, MAST-U, TCV, WEST and JET (since 2022) work together to achieve the objectives of Missions 1 and 2 of the EUROfusion Roadmap: i) demonstrate plasma scenarios that increase the success margin of ITER and satisfy the requirements of DEMO and, ii) demonstrate an integrated approach that can handle the large power leaving ITER and DEMO plasmas. The Tokamak Exploitation task force has therefore organized experiments on these two missions with the goal to strengthen the physics and operational basis for the ITER baseline scenario and for exploiting the recent plasma exhaust enhancements in all four devices (PEX: Plasma EXhaust) for exploring the solution for handling heat and particle exhaust in ITER and develop the conceptual solutions for DEMO. The ITER Baseline scenario has been developed in a similar way in ASDEX Upgrade, TCV and JET. Key risks for ITER such as disruptions and run-aways have been also investigated in TCV, ASDEX Upgrade and JET. Experiments have explored successfully different divertor configurations (standard, super-X, snowflakes) in MAST-U and TCV and studied tungsten melting in WEST and ASDEX Upgrade. The input from the smaller devices to JET has also been proven successful to set-up novel control schemes on disruption avoidance and detachment.
JET’s frequency-modulated continuous wave (FMCW) reflectometers have been operating well with the current design since 2005, and density profiles have been automatically calculated intershot since then. However, the calculated profiles had long suffered from several shortcomings: poor agreement with other diagnostics, sometimes inappropriately moving radially by several centimeters, elevated levels of radial jitter, and persistent wriggles (strong unphysical oscillations). In this research, several techniques are applied to the reflectometry data analysis, and the shortcomings are significantly improved. Starting with improving the equilibrium reconstruction that estimates the background magnetic field, adding a ripple correction in the reconstructed magnetic field profile, and adding new inner-wall reflection positions estimated through ray-tracing, these changes not only improve the agreement of reconstructed profiles to other diagnostics but also solve density profile wriggles that were present during band transitions. Other smaller but also persistent wriggles were also suppressed by applying a localized correction to the measured beat frequency where persistent oscillations are present. Finally, the burst analysis method, as introduced by Varela et al. [Nucl. Fusion 46 S693 (2006)], has been implemented to extract the beat frequency from stacked spectrograms. Due to the strong suppression of spurious reflections, the radial jitter that sometimes would span several centimeters has been strongly reduced. The stacking of spectrograms has also been shown to be very useful for stacking recurring events, like small gas puff modulations, and extracting transport coefficients that would otherwise be below the noise level.
The DIII-D tokamak has elucidated crucial physics and developed projectable solutions for ITER and fusion power plants in the key areas of core performance, boundary heat and particle transport, and integrated scenario operation, with closing the core-edge integration knowledge gap being the overarching mission. New experimental validation of high-fidelity, multi-channel, non-linear gyrokinetic turbulent transport models for ITER provides strong confidence it will achieve Q 10 operation. Experiments identify options for easing H-mode access in hydrogen, and give new insight into the isotopic dependence of transport and confinement. Analysis of 2,1 islands in unoptimized low-torque IBS demonstration discharges suggests their onset time occurs randomly in the constant beta phase, most often triggered by non-linear 3-wave coupling, thus identifying an NTM seeding mechanism to avoid. Pure deuterium SPI for disruption mitigation is shown to provide favorable slow cooling, but poor core assimilation, suggesting paths for improved SPI on ITER. At the boundary, measured neutral density and ionization source fluxes are strongly poloidally asymmetric, implying a 2D treatment is needed to model pedestal fuelling. Detailed measurements of pedestal and SOL quantities and impurity charge state radiation in detached divertors has validated edge fluid modelling and new self-consistent 'pedestal-to-divertor' integrated modeling that can be used to optimize reactors. New feedback adaptive ELM control minimizes confinement reduction, and RMP ELM suppression with sustained high core performance was obtained for the first time with the outer strike point in a W-coated, compact and unpumped small-angle slot divertor. Advances have been made in integrated operational scenarios for ITER and power plants. Wide pedestal intrinsically ELM-free QH-modes are produced with more reactor-relevant conditions, Low torque IBS with W-equivalent radiators can exhibit predator-prey oscillations in T-e and radiation which need control. High-beta(P) scenarios with q(min) > 2, q(95)-7.9, beta(N) > 4, beta(T)-3.3% and H-98y2 > 1.5 are sustained with high density ((n) over bar = 7E19 m(-3), f(G)-1) for 6 tau(E), improving confidence in steady-state tokamak reactors. Diverted NT plasmas achieve high core performance with a non-ELMing edge, offering a possible highly attractive core-edge integration solution for reactors.
This study employs the established momentum transport analysis at ASDEX Upgrade [Zimmermann et al., Nucl. Fusion 63, 124003 (2023)] to investigate the parametric variations of the momentum transport coefficients in the core of H-mode plasmas. These experimental results are compared to a comprehensive database of gyrokinetic calculations. Generally, good agreement between predicted and measured diffusive and convective transport coefficients is found. The predicted and measured Prandtl numbers correlate most dominantly with the magnetically trapped particle fraction. The experimentally inferred pinch numbers strongly depend on the logarithmic density gradient and magnetic shear, consistent with the theoretical predictions of the Coriolis pinch. The intrinsic torque from residual stress in the inner core is small, scales with the local logarithmic density gradient, and the data indicate a possible sign reversal. In the outer periphery of the core, the intrinsic torque is always co-current-directed and scales with the pressure gradient. This is consistent with prior experimental findings and global, non-linear gyrokinetic predictions. It suggests that profile shearing effects generate the intrinsic torque in the inner core. Toward the outer core, most likely, effects from E×B-shearing become more influential. These results offer the first comprehensive picture of this transport channel in the core plasma and contribute to validating the corresponding theoretical understanding. The derived scaling laws are used to construct a reduced momentum transport model, which has been validated against an additional dataset. This demonstrates that the model captures the essential contributions to momentum transport in the core of H-mode plasmas.
Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98( p,y ) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’ q -profile) were established on AUG and characterised concerning their plasma performance.
In the JET DTE2 campaign a new method was successfully tested to detect the heating of bulk electrons by α-particles, using the dynamic response of the electron temperature T e to the modulation of ion cyclotron resonance heating (ICRH). A fundamental deuterium (D) ICRH scheme was applied to a tritium-rich hybrid plasma with D-neutral beam injection (NBI). The modulation of the ion temperature T i and of the ICRH accelerated deuterons leads to modulated α -heating with a large delay with respect to other modulated electron heating terms. A significant phase delay of ∼40° is measured between central T e and T i , which can only be explained by α -particle heating. Integrated modelling using different models for ICRH absorption and ICRH/NBI interaction reproduces the effect qualitatively. Best agreement with experiment is obtained with the European Transport Solver/Heating and Current Drive workflow.
The simplification of collision operators is necessary for quasilinear turbulence modelling used with integrated modelling frameworks, such as the gyrokinetic code QuaLiKiz. The treatment of collisions greatly impacts the accuracy of trapped electron mode (TEM) modelling, which is necessary to predict the electron heat flux and the balance between inward and outward particle fluxes. In particular, accurate particle flux predictions are necessary to successfully model density peaking in the tokamak core. We explored two ways of improving collisional TEM model reduction for tokamak core plasmas. First, we carried out linear GENE simulations to study the complex interplay between collisions and trapped electrons. We then used these simulations to define an effective trapped fraction to characterize the collisional TEM based on two key parameters, the local inverse-aspect ratio $\epsilon$ and the collisionality $\nu ^\ast$ . One aspect missing from analytical TEM research is that the collisional frequency and the bounce-transit frequency are both velocity dependent; this effective trapped fraction takes both into account. In doing so, we determined that two parameters are not enough to model the collisional TEM, as an additional third free parameter was necessary. We determined that this model, as currently formulated, is not suitable for integrated modelling purposes. Second, we directly improved QuaLiKiz's Krook operator, which relies on two free parameters. We determined that these parameters required adjustments against higher-fidelity collisional models. In order to improve density profile predictions when paired with integrated models, we refined the Krook operator by using GENE simulations as a higher-fidelity point of comparison. We then demonstrate strong improvement of density peaking predictions of QuaLiKiz within the integrated modelling framework JETTO.
The dimensionless isotope mass scaling experiment between pure Deuterium and pure Tritium plasmas with matched ρ ∗ , ν ∗ , β n , q and T e / T i has been achieved in JET L-mode with dominant electron heating (NBI+ohmic) conditions. 28% higher scaled energy confinement time B t τ E , t h / A is found in favour of the Tritium plasma. This can be cast in the form of the dimensionless energy confinement scaling law as Ω i τ E , t h ∼ A 0.48 ± 0.16 . This significant isotope mass scaling is consequently seen in the scaled one-fluid heat diffusion coefficient A χ e f f / B t which is around 50% lower in the Tritium plasma throughout the whole plasma radius. The isotope mass dependence in the particle transport channel is negligible, supported also by the perturbative particle transport analysis with gas puff modulation. The comparison of the edge particle fuelling or ionisation profiles from the EDGE2D-EIRENE simulations show that the absolute density differences that are necessary for the dimensionless match in the confined plasma dominate over any isotope mass dependencies of particle fuelling and ionization profiles at the plasma edge. Local GENE simulation results indicate a mild anti-gyroBohm effect at ρ t o r = 0.6 and thereby a small isotope mass dependence in favour of Tritium on heat transport and a negligible effect on particle transport. A significant fraction of the isotope scaling and reduced heat transport observed in the Tritium plasma is not captured in the GENE and ASTRA-TGLF-SAT2 simulations by simply changing the isotope mass for the same input profiles.
A new approach to infer the momentum transport in tokamak core plasmas via perturbation experiments is presented. For the first time, the analysis self-consistently includes all momentum transport components and their time dependencies, which are essential to separate the momentum fluxes and closely match the experiment. The quantitative agreement between the experimentally inferred transport coefficients and the gyrokinetic predictions provides an unprecedented validation. This work shows that the new methodology and gyrokinetic predictions can now be utilized on the route to physics-based prediction of momentum transport in future reactor plasmas.
Advanced momentum transport analysis is used to study matched hydrogen (H) and deuterium (D) plasmas in the core of ASDEX Upgrade. The aim is to validate gyrokinetic theory and assess a possible isotope dependence. The methodology extracts momentum diffusion, convection, and intrinsic torque as a function of time from experiments employing neutral beam injection (NBI) modulation. H and D plasma scenarios with comparable ion heat fluxes, NBI torque, electron densities, and several dimensionless parameters were designed to highlight any mass dependency. Linear gyrokinetic simulations predict that, for similar background gradients, the Prandtl and pinch numbers should be similar for H and D. This was confirmed by the experimental momentum transport analyses. The assessed intrinsic torques were found to be similar between H and D, co-current directed and located near the outermost region of the plasma core. The strength of the intrinsic torque is correlated with the amplitude of the plasma pressure gradient in the pedestal. Finally, a robust error analysis demonstrates the uniqueness of the parameters obtained together with their uncertainties. Neglecting the intrinsic torque, or its time dependence, systematically distorts the assessed momentum diffusion and convection. This is the first method to separate all three transport mechanisms from experimental data by retaining their time dependencies, that is found to match, quantitatively, the gyrokinetic predictions for Prandtl and pinch numbers, within experimental uncertainties.
Strike-point sweeping, a technique often used to spread heat loads on divertor targets, was employed in JET experiments for the first time to generate an edge-localized modulated particle source for investigating plasma fuelling and particle transport in the edge region. This approach was motivated by the possibility of achieving higher modulation frequencies than those available from traditional gas puff modulation at JET. Higher frequencies would enable the collection of more edge-localized information from the electron density response to the modulated particle source. Various sweeping frequencies, up to 18.5 Hz, were commissioned and utilized in the experiments. Both strong and weak electron density responses were observed in H-mode plasmas, depending on the strike-point configuration and the distance the strike-points moved during the sweep cycle. The electron density response exhibited complex and unconventional behaviour (compared to gas puff modulation), which presented challenges for interpretation. In this study, we analyse one experiment in detail using an optimization framework in which transport and particle source parameters are determined by fitting our forward model parameters to the experimental electron density measurements. We demonstrate that a consistent picture emerges and that our approach can provide new insights into these complex data. However, we note that while strike-point sweeping generates the desired modulated edge-localized particle source, it also modifies the properties of the edge transport barrier. Therefore, the strike-point sweeping methodology is a promising but challenging way to study edge particle transport and edge fuelling properties, requiring very precise measurements.
The prediction of plasma rotation is of high interest for fusion research due to the effects of the rotation upon magnetohydrodynamic (MHD) instabilities, impurities, and turbulent transport in general. In this work, an analysis method was studied and validated to reliably extract momentum transport coefficients from neutral beam injection (NBI) modulation experiments. To this end, a set of discharges was created with similar background profiles for the ion and electron temperatures, the heat fluxes, the electron density, and the plasma rotation that, therefore, should exhibit similar momentum transport coefficients. In these discharges, a range of temporal perturbations were imposed by modulating and varying the power deposition of the NBI, electron-cyclotron-resonance heating (ECRH), and ion-cyclotron-resonance heating (ICRH). The transport model including diffusion, convection, and residual stress was implemented within the ASTRA code. The Prandtl number P r = χ φ / χ i was assessed via the GKW code. A convective Coriolis pinch was fitted and the intrinsic torque from the residual stress was estimated. The obtained transport coefficients agree within error bars for sufficiently small imposed temperature perturbations, as would be expected, from the similar background profiles. This successful validation of the methodology opens the door to study the parametric dependence of the diffusive and convective momentum transport of the main ions of the plasma as well as the turbulent intrinsic torque in a future work.
DIII-D physics research addresses critical challenges for the operation of ITER and the next generation of fusion energy devices. This is done through a focus on innovations to provide solutions for high performance long pulse operation, coupled with fundamental plasma physics understanding and model validation, to drive scenario development by integrating high performance core and boundary plasmas. Substantial increases in off-axis current drive efficiency from an innovative top launch system for EC power, and in pressure broadening for Alfven eigenmode control from a co-/counter- I p steerable off-axis neutral beam, all improve the prospects for optimization of future long pulse/steady state high performance tokamak operation. Fundamental studies into the modes that drive the evolution of the pedestal pressure profile and electron vs ion heat flux validate predictive models of pedestal recovery after ELMs. Understanding the physics mechanisms of ELM control and density pumpout by 3D magnetic perturbation fields leads to confident predictions for ITER and future devices. Validated modeling of high- Z shattered pellet injection for disruption mitigation, runaway electron dissipation, and techniques for disruption prediction and avoidance including machine learning, give confidence in handling disruptivity for future devices. For the non-nuclear phase of ITER, two actuators are identified to lower the L–H threshold power in hydrogen plasmas. With this physics understanding and suite of capabilities, a high poloidal beta optimized-core scenario with an internal transport barrier that projects nearly to Q = 10 in ITER at ∼ 8 MA was coupled to a detached divertor, and a near super H-mode optimized-pedestal scenario with co- I p beam injection was coupled to a radiative divertor. The hybrid core scenario was achieved directly, without the need for anomalous current diffusion, using off-axis current drive actuators. Also, a controller to assess proximity to stability limits and regulate β N in the ITER baseline scenario, based on plasma response to probing 3D fields, was demonstrated. Finally, innovative tokamak operation using a negative triangularity shape showed many attractive features for future pilot plant operation.
An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L–H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts.
Density peaking has been studied between an ICRH and NBI identity plasma in JET. The comparison shows that 8 MW of NBI heating/fueling increases the density peaking by a factor of two, being R/L n = 0.45 for the ICRH pulse and R/L n = 0.93 for the NBI one averaged radially over ρ tor = 0.4, 0.8. The dimensionless profiles of q, ρ *, υ *, β n and T i/T e ≈ 1 were matched within 5% difference except in the central part of the plasma (ρ tor < 0.3). The difference in the curvature pinch (same q-profile) and thermo-pinch (T i = T e) between the ICRH and NBI discharges is virtually zero. Both the gyro-kinetic simulations and integrated modelling strongly support the experimental result where the NBI fuelling is the main contributor to the density peaking for this identity pair. It is to be noted here that the integrated modeling does not reproduce the measured electron density profiles, but approximately reproduces the difference in the density profiles between the ICRH and NBI discharge. Based on these modelling results and the analyses, the differences between the two pulses in impurities, fast ions (FIs), toroidal rotation and radiation do not cause any such changes in the background transport that would invalidate the experimental result where the NBI fuelling is the main contributor to the density peaking. This result of R/L n increasing by a factor of 2 per 8 MW of NBI power is valid for the ion temperature gradient dominated low power H-mode plasmas. However, some of the physics processes influencing particle transport, like rotation, turbulence and FI content scale with power, and therefore, the simple scaling on the role of the NBI fuelling in JET is not necessarily the same under higher power conditions or in larger devices.
The use of Panda-type polarization-maintaining (PM) fiber for the localized sensing of high temperatures was analyzed with simulations and experiments up to 900°C. Accuracy and repeatability of the results started to decline above 800°C. Fused silica optical fiber melts at 1700°C, which sets an ultimate limit for measurable temperatures. In practice, optical fiber birefringence restricts the maximum temperature to 1060°C where PM fiber loses its ability to maintain polarization. Three sensor fibers (4, 5 and 10 cm long) were spliced at 45° angles to input/output fibers and calibrated in an oven from room temperature to 850–900°C temperature range. Two superluminescent light-emitting diodes (SLEDs) were coupled together as a broadband light source. Birefringence-induced change of the polarization in the sensor fiber was measured with a polarization splitter and an optical spectrum analyzer (OSA) as a function of the wavelength. Temperature-dependent birefringence generates a sinusoidal reflection spectrum from the input polarization mode to the orthogonal output polarization mode. Temperature changes could be concluded from variations in these spectra. Finally, a small fusion device, NORTH, at DTU, Denmark was successfully used as a testbed to make sure that the sensors can handle transportation and the instrumentation required for vacuum operation and still produce sensible data from a harsh environment.