Learning from Each Other: Cross-Cutting Diagnostic Development Activities Between Magnetic and Inertial Confinement Fusion
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Tinguely, R. A..
Explore the source record for details and available documents.
While much about Alfvén eigenmode (AE) stability has been explored in previous and current tokamaks, open questions remain for future burning plasma experiments, especially regarding exact stability threshold conditions and related isotope effects; the latter, of course, requiring good knowledge of the plasma ion composition. In the JET tokamak, eight in-vessel antennas actively excite stable AEs, from which their frequencies, toroidal mode numbers, and net damping rates are assessed. The effective ion mass can also be inferred using measurements of the plasma density and magnetic geometry. Thousands of AE stability measurements have been collected by the Alfvén Eigenmode Active Diagnostic in hundreds of JET plasmas during the recent Hydrogen, Deuterium, Tritium, DT, and Helium-4 campaigns. In this novel AE stability database, spanning all four main ion species, damping is observed to decrease with increasing Hydrogenic mass, but increase for Helium, a trend consistent with radiative damping as the dominant damping mechanism. These data are important for confident predictions of AE stability in both non-nuclear (H/He) and nuclear (D/T) operations in future devices. In particular, if radiative damping plays a significant role in overall stability, some AEs could be more easily destabilized in D/T plasmas than their H/He reference pulses, even before considering fast ion and alpha particle drive. Active MHD spectroscopy is also employed on select HD, HT, and DT plasmas to infer the effective ion mass, thereby closing the loop on isotope analysis and demonstrating a complementary method to typical diagnosis of the isotope ratio.
Dedicated experiments were performed in JET DTE2 plasmas for obtaining an α-particle bump-on-tail (BOT) distribution aiming at exciting Alfvén eigenmodes (AEs). Neutral beam injection-only heating with modulated power was used so that fusion-born α-particles were the only ions present in the MeV energy range in these DT plasmas. The beam power modulation on a time scale shorter than the α-particle slowing down time was chosen for modulating the α-particle source and thus sustaining a BOT in the α-particle distribution. High-frequency modes in the toroidicity-induced Alfven eigenmode (TAE) frequency range and multiple short-lived modes in a wider frequency range have been detected in these DT discharges with interferometry, soft x-ray cameras, and reflectometry. The modes observed were localised close to the magnetic axis and were not seen in the Mirnov coils. Analysis with the TRANSP and Fokker-Planck FIDIT codes confirms that α-particle distributions with BOT in energy were achieved during some time intervals in these discharges though no clear correlation was found between the times of the high-frequency mode excitation and the BOT time intervals. The combined magneto-hydrodynamic (MHD) and kinetic modelling studies show that the high-frequency mode in the TAE frequency range is best fitted with a TAE of toroidal mode number n= 9. This mode is driven mostly by the on-axis beam ions while the smaller drive due to the pressure gradient of α-particles allows overcoming the marginal stability and exciting the mode. The observed multiple short-lived modes in a wider frequency range are identified as the on-axis kinetic AEs predicted in Rosenbluth and Rutherford.
Asymmetric halo currents (HCs) can exert large net forces on the vacuum vessel and other components during disruptions on tokamaks. The displacements caused by these forces can then be amplified if these asymmetric forces rotate at frequencies resonant with the vessel. This paper reports on the investigation of a recently proposed scaling law for the disruption HC rotation frequency [Saperstein et al., “Halo current rotation scaling in post-disruption plasmas,” Nucl. Fusion 62, 026044 (2022)] that combines measurements on Alcator C-Mod with those on HBT-EP. We find that a new non-circular version of the scaling law [ ⟨ f rot ⟩ m / ⟨ m ⟩ ∝ 1 B T ( S / π )] takes into consideration the dependence of frot on the poloidal structure of the MHD instability (m) driving the asymmetry and describes the disruption-averaged rotation frequency on C-Mod. Disruption rotation is also found to be insensitive to the vertical position and impurity content of the plasma at the onset of the disruption. However, a stagnation in the time evolution of frot is occasionally observed. In conclusion, observations are consistent with the dominance of poloidal rotation during the disruption, which is motivated by the poloidal drift nature of the scaling law.
In Izzo et al (2022 Nucl. Fusion 62 096029), state-of-the-art modeling of thermal and current quench (CQ) magnetohydrodynamics (MHD) coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1MA compared to ~8.7MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.4 and with kinetic energies ~0.2–15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 m 2 s -1 , is needed in the core to fully mitigate these REs. Properly accounting for (a) the CQ electric field's effect on RE transport in islands and (b) the contribution of significant RE currents to disruption MHD may help achieve this.
Here, using the SPIRAL Monte Carlo, full particle-orbit simulation code, we investigate the effects of neoclassical tearing modes (NTMs) and toroidal field (TF) ripple on alpha power losses during steady-state operation of the SPARC primary reference discharge. Model perturbations for TF ripple and the $m/n = 2/1$ and $3/2$ NTMs with exaggerated widths selected based on an H-mode plasma approaching thermal quench are added to a simulated SPARC magnetic equilibrium through which marker particles are tracked. The 3/2 and 2/1 NTMs are located at $\rho_{\mathrm{pol}} \sim 0.76$ and $\rho_\mathrm{pol} \sim 0.86$ respectively, well positioned to increase alpha particle transport into and within an outer lossy region of the plasma beyond $\rho_\mathrm{pol} \sim 0.8$ where over 95% of lost alpha particles are born. Total alpha power losses are shown to increase modestly from 1.73% lost at a minimum to 2.34% lost at a maximum, and alpha particle surface power densities form localized hotspots on the first-wall near the low-field side midplane due to NTMs and TF ripple. We establish a conservative upper limit for first-wall alpha surface power densities on a toroidally symmetric wall for typical, flattop operation and motivate the consideration of NTMs in the design of three dimensional limiter surfaces for SPARC.
Here, we report the novel experimental observation of both unstable and stable toroidicity-induced Alfvén eigenmodes (TAEs) measured simultaneously in a JET tokamak plasma. The three-ion-heating scheme (D-DNBI-3He) is employed to accelerate deuterons to MeV energies, thereby destabilizing TAEs with toroidal mode numbers n = 3–5, each decreasing in mode amplitude. At the same time, the Alfvén eigenmode active diagnostic resonantly excites a stable n = 6 TAE with total normalized damping rate –γ/ω 0 ≈ 1%–4%. Hybrid kinetic-MHD modeling with codes NOVA-K and MEGA both find eigenmodes with similar frequencies, mode structures, and radial locations as in experiment. NOVA-K demonstrates good agreement with the n = 3, 4, and 6 TAEs, matching the damping rate of the n = 6 mode within uncertainties and identifying radiative damping as the dominant contribution. Improved agreement is found with MEGA for all modes: the unstable n = 3–5 and stable n = 2, 6 modes, with the latter two stabilized by higher intrinsic damping and lower fast ion drive, respectively. While some discrepancies remain to be resolved, this unique validation effort gives us confidence in TAE stability predictions for future fusion devices.
The operation of a 3D coil—passively driven by the current quench (CQ) loop voltage—for the deconfinement of runaway electrons (REs) is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear magnetohydrodynamic (MHD) modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for REs based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely et al (2021 Nucl. Fusion 61 124003), considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the n = 1 coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the thermal quench MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1–2 MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the CQ is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil.
The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0 , net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0 , γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.
Runaway electrons (REs) created during tokamak disruptions pose a threat to the reliable operation of future larger machines. Experiments using shattered pellet injection (SPI) have been carried out at the JET tokamak to investigate ways to prevent their generation or suppress them if avoidance is not sufficient. Avoidance is possible if the SPI contains a sufficiently low fraction of high-Z material, or if it is fired early in advance of a disruption prone to runaway generation. These results are consistent with previous similar findings obtained with Massive Gas Injection. Suppression of an already accelerated beam is not efficient using High-Z material, but deuterium leads to harmless terminations without heat loads. This effect is due to the combination of a large magnetohydrodynamic instability scattering REs on a large area and the absence of runaway regeneration during the subsequent current collapse thanks to the flushing of high-Z impurities from the runaway companion plasma. This effect also works in situations where the runaway beam moves upwards and undergoes scraping-off on the wall.
The interaction of Alfven Eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25-250 kHz and toroidal mode numbers |n| < 20. During the 2019-2020 deuterium campaign, almost 7500 resonances and their frequencies f0, net damping rates \gamma < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient's effect on accessibility and ELM-related noise's impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized Ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ICRH and NBI) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f0, \gamma, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.
Here, the potential formation of multi-mega-ampere beams of relativistic 'runaway' electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetric in-vessel coil will passively prevent RE beam formation during disruptions in the SPARC tokamak, a compact, high-field, high-current device capable of achieving a fusion gain Q > 2 in deuterium–tritium plasmas.
Recent JET deuterium experiments with an advanced tokamak scenario using an internal transport barrier (ITB) have been performed to clearly observe destabilized toroidicity-induced Alfvén eigenmodes (TAEs) by fast ions; interestingly, these also exhibit unstable electromagnetic (EM) perturbations in the sub-TAE frequency range. We identify such EM perturbations to be beta-induced ion temperature gradient (BTG) eigenmodes and not beta-induced Alfvén eigenmodes nor beta-induced Alfvén acoustic eigenmodes, which are often unstable in such high-beta plasmas with high-power neutral beam injection (NBI). The BTG modes are the most unstable modes due to the high thermal ion temperature gradient related to the ITB, high thermal ion temperature compared to thermal electron temperature (high T i /T e ), and a high ion beta regime. BTG mode experimental characteristics match the analytical theory, i.e., location in the vicinity of a rational magnetic surface with a low magnetic shear, mode frequency scaling with the ion diamagnetic frequency (ω$^*_i$), and a coupling among Alfvén and drift waves. We also perform linear gyrokinetic simulations with validated plasma profiles and equilibrium, and find a mode kinetically driven by thermal ions with similar characteristics as the experimental BTG modes.
Next generation High Performance (HP) tokamaks risk damage from unmitigated disruptions at high current and power. Achieving reliable disruption prediction for a device's HP operation based on its Low Performance (LP) data is a key to its success. In this letter, through explorative data analysis and dedicated numerical experiments on multiple existing tokamaks, we demonstrate how the operational regimes of tokamaks can affect the power of a trained disruption predictor. First, our results suggest data-driven disruption predictors trained on abundant LP discharges work poorly on the HP regime of the same tokamak, which is a consequence of the distinct distributions of the tightly correlated signals related to disruptions in these two regimes. Second, we find that matching operational parameters among tokamaks strongly improves cross-machine accuracy which implies our model learns from the underlying scalings of dimensionless physics parameters like q 95 , β p and confirms the importance of these parameters in disruption physics and cross machine domain matching from the data-driven perspective. Finally, our results show in the absence of HP data from the target devices, the best predictivity of the HP regime for the target machine can be achieved by combining LP data from the target with HP data from other machines. Furthermore, these results provide a possible disruption predictor development strategy for next generation tokamaks, such as ITER and SPARC, and highlight the importance of developing baseline scenario discharges of future tokamaks on existing machines to collect more relevant disruptive data.
This paper presents results of extensive analysis of mode excitation observed during the operation of the Alfven Eigenmode Active Diagnostic (AEAD) in the JET tokamak during the 2019-2020 deuterium campaign. Six of eight toroidally spaced antennas, each with independent power and phasing, were successful in actively exciting stable MHD modes in 479 plasmas. In total, 4768 magnetic resonances were detected with up to fourteen fast magnetic probes. In this work, we present the calculations of resonant frequencies f0, damping rates \gamma < 0, and toroidal mode numbers n, spanning the parameter range f0 = 30 - 250 kHz, -\gamma = 0 - 13 kHz, and |n| < 30. In general, good agreement is seen between the resonant and the calculated toroidal Alfven Eigenmode frequencies, and between the toroidal mode numbers applied by the AEAD and estimated of the excited resonances. We note several trends in the database: the probability of resonance detection decreases with plasma current and external heating power; the normalized damping rate increases with edge safety factor but decreases with external heating. These results provide key information to prepare future experimental campaigns and to better understand the physics of excitation and damping of Alfven Eigenmodes in the presence of alpha particles during the upcoming DT campaign, thereby extrapolating with confidence to future tokamaks.
This paper presents a dedicated study of plasma-antenna (PA) coupling with the Alfven Eigenmode Active Diagnostic (AEAD) in JET. Stable AEs and their resonant frequencies f, damping rates gamma < 0, and toroidal mode numbers n are measured for various PA separations and limiter versus X-point magnetic configurations. Two stable AEs are observed to be resonantly excited at distinct low and high frequencies in limiter plasmas. The values of f and n do not vary with PA separation. However, |gamma| increases with PA separation for the low-f, but not high-f, mode, yet this may be due to slightly different edge conditions. The high-f AE is detected throughout the transition from limiter to X-point configuration, though its damping rate increases; the low-f mode, on the other hand, becomes unidentifiable. The linear resistive MHD code CASTOR is used to simulate the frequency scan of an AEAD-like external antenna. For the limiter pulses, the high-f mode is determined to be an n = 0 GAE, while the low-f mode is likely an n = 2 TAE. During the transition from limiter to X-point configuration, CASTOR indicates that n = 1 and 2 EAEs are excited in the edge gap. These results extend previous experimental studies in JET and Alcator C-Mod; validate the computational work performed by Dvornova et al 2020 Phys. Plasmas 27 012507; and provide guidance for the optimization of PA coupling in upcoming JET energetic particle experiments, for which the AEAD will aim to identify the contribution of alpha particles to AE drive during the DT campaign.
In this paper, we present a new deep learning disruption prediction algorithm based on important findings from explorative data analysis which effectively allows knowledge transfer from existing devices to new ones, thereby predicting disruptions using very limited disruptive data from the new devices. The explorative data analysis conducted via unsupervised clustering techniques confirms that time-sequence data are much better separators of disruptive and non-disruptive behavior than the instantaneous plasma state data with further advantageous implications for a sequence-based predictor. Based on such important findings, we have designed a new algorithm for multi-machine disruption prediction that achieves high predictive accuracy on the C-Mod (AUC=0.801), DIII-D (AUC=0.947) and EAST (AUC=0.973). tokamaks with limited hyperparameter tuning. Through numerical experiments, we show that boosted accuracy (AUC=0.959) is achieved on EAST predictions by including in the training only 20 disruptive discharges, thousands of non-disruptive discharges from EAST, and combining this with more than a thousand discharges from DIII-D and C-Mod. The improvement of predictive ability obtained by combining disruptive data from other devices is found to be true for all permutations of the three devices. Furthermore, by comparing the predictive performance of each individual numerical experiment, we find that non-disruptive data are machine-specific while disruptive data from multiple devices contain device-independent knowledge that can be used to inform predictions for disruptions occurring on a new device.
Here, this paper presents a dedicated study of plasma-antenna (PA) coupling with the Alfvén Eigenmode Active Diagnostic (AEAD) in JET. Stable AEs and their resonant frequencies f, damping rates γ < 0, and toroidal mode numbers n are measured for various PA separations and limiter versus X-point magnetic configurations. Two stable AEs are observed to be resonantly excited at distinct low and high frequencies in limiter plasmas. The values of f and n do not vary with PA separation. However, |γ| increases with PA separation for the low-f, but not high-f, mode, yet this may be due to slightly different edge conditions. The high-f AE is detected throughout the transition from limiter to X-point configuration, though its damping rate increases; the low-f mode, on the other hand, becomes unidentifiable. The linear, resistive MHD code CASTOR is used to simulate the frequency scan of an AEAD-like external antenna. For the limiter pulses, the high-f mode is determined to be an n = 0 GAE, while the low-f mode is likely an n = 2 TAE. During the transition from limiter to X-point configuration, CASTOR indicates that n = 1 and 2 EAEs are excited in the edge gap. These results extend previous experimental studies in JET and Alcator C-Mod; validate the computational work performed by Dvornova et al; and provide guidance for the optimization of PA coupling in upcoming JET energetic particle experiments, for which the AEAD will aim to identify the contribution of alpha particles to AE drive during the DT campaign.