Engineering topics

Sweeney, R.

Publications and source records attributed to Sweeney, R..

14 records

DOE OSTI2023

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil

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.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Physics of runaway electrons with shattered pellet injection at JET

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Shattered pellet injection experiments at JET in support of the ITER disruption mitigation system design

A series of experiments have been executed at JET to assess the efficacy of the newly installed shattered pellet injection (SPI) system in mitigating the effects of disruptions. Issues, important for the ITER disruption mitigation system, such as thermal load mitigation, avoidance of runaway electron (RE) formation, radiation asymmetries during thermal quench mitigation, electromagnetic load control and RE energy dissipation have been addressed over a large parameter range. The efficiency of the mitigation has been examined for the various SPI injection strategies. Here, the paper summarises the results from these JET SPI experiments and discusses their implications for the ITER disruption mitigation scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Disruption thermal load mitigation with shattered pellet injection on the Joint European Torus (JET)

Disruption mitigation remains a critical, unresolved challenge for ITER. To aid in addressing this challenge, a shattered pellet injection (SPI) system was installed on JET and experiments conducted at a range of thermal energy fractions and stored energies in excess of 7 MJ. The primary goals of these experiments were to investigate the efficacy of the SPI on JET and the ability of the plasma to assimilate multiple pellets. Single pellet injections produced a saturation in total radiated energy (Wrad) with increasing injected neon content, suggesting total radiation of stored thermal energy. Further increases in injected neon quantities resulted in reduced cooling times and current quench (CQ) durations, indicating higher impurity assimilation. No significant variation in CQ duration or W rad was observed when varying the deuterium content at fixed neon quantities. Additionally, higher assimilation, inferred by shorter CQ durations, was measured when a mechanical punch was used to launch the pellets and this was attributed to a lower pellet velocity leading to higher solid content in the pellet plume and larger fragments penetrating deeper into the plasma. Radiation asymmetries averaged over the cooling time were inferred from Emis3D and ranged from 1.6 to 1.9. Asymmetries averaged over the entire disruption sequence were found to increase at higher thermal energy fractions. The radiated energy fractions decreased with increasing thermal energy fractions but this trend was eliminated when toroidal asymmetries were accounted for with Emis3D. Pure deuterium pellets were able to produce cooling times of up to 75 ms with a gradual loss in thermal stored energy of up to 80%. Experiments with multiple pellet injection indicated W rad can be increased through pellet superposition and density can be increased with an additional D2 injection without a reduction in W rad . KPRAD modelling accurately reproduced the cooling times and the CQ duration at high thermal energies. Assimilation estimates from KPRAD indicated CQ rates scale strongly whilst W rad scales weakly and saturates with assimilated neon content. Comparable W rad can be achieved with lower assimilated neon quantities as longer cooling times are attained. Thus reduced neon content can be preferential in a thermal load mitigation scheme as it may reduce radiation asymmetries and prevent flash melting.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Modeling the complete prevention of disruption-generated runaway electron beam formation with a passive 3D coil in SPARC

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Shattered pellet penetration in low and high energy plasmas on DIII-D

Shattered pellet injection (SPI) has been adopted as the baseline disruption mitigation system for ITER, as the radiative payload penetration into DIII-D plasmas from SPI is superior to those using the massive gas injection (MGI) method. Because of the substantial differences in the energy content of ITER plasma and those in present experiments, reliable 3D MHD modeling, benchmarked against present experiments is needed to project to ITER plasmas. In support of these needs, the depth of SPI fragment penetration in DIII-D plasmas was investigated by injecting SPI into two discharges with vastly different energy content and pedestal height. 400 Torr-L pure Ne fragmented pellets at a velocity of about 200 m s−1 were injected into a 0.2 MJ L-mode discharge and a 2 MJ super H-mode discharge. Results show deep penetration of SPI fragments into low-energy plasmas in DIII-D. SPI fragment penetration is reduced as the plasma energy content increases, with some discharges exhibiting penetration that is confined to the outer regions of the plasma. The injected SPI fragments are also spread out over a distance of about 20 cm, which results in some fragments arriving near the end of or after the thermal quench is over.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Scenario adaptive disruption prediction study for next generation burning-plasma tokamaks

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Hybrid deep learning architecture for general disruption prediction across tokamaks

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Methods to determine the radiated power in SPI-mitigated disruptions in JET

Here we present techniques for evaluating the radiated power in JET disruptions. Disrupting plasmas are shown to have non-axisymmetric radiation profiles, motivating the re-evaluation of the standard techniques for calculating the total radiated power at JET using bolometry. Four single-channel bolometers at different toroidal locations are exploited to quantify the radiation asymmetry. Toroidal radiation peaking factors integrated over the entire disruption of up to 1.5 have been observed when varying the quantity of neon in pellets used in disruptions mitigated by shattered pellet injection. Using synthetic bolometer diagnostics developed with the Cherab spectroscopy modeling framework, we can estimate the systematic error on total power calculations for relevant radiation profiles and improve estimates of the total radiated power. We show that the component of the systematic error on the total power due to the poloidal radiation profile can be reduced from 70% to 10% with suitable assumptions about the structure of the poloidal profile.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

MHD stability and disruptions in the SPARC tokamak

SPARC is being designed to operate with a normalized beta of $β_N$ = 1.0, a normalized density of $n_G$ = 0.37 and a safety factor of $q_{95}$ ≈ 3.4, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $β_p$ = 0.19 at the safety factor $\textit{q}$ = 2 surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ~80 %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order 10 –2 that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Overview of the SPARC tokamak

The SPARC tokamak is a critical next step towards commercial fusion energy. SPARC is designed as a high-field ( $B_0 = 12.2$ T), compact ( $R_0 = 1.85$ m, $a = 0.57$ m), superconducting, D-T tokamak with the goal of producing fusion gain $Q>2$ from a magnetically confined fusion plasma for the first time. Currently under design, SPARC will continue the high-field path of the Alcator series of tokamaks, utilizing new magnets based on rare earth barium copper oxide high-temperature superconductors to achieve high performance in a compact device. The goal of $Q>2$ is achievable with conservative physics assumptions ( $H_{98,y2} = 0.7$ ) and, with the nominal assumption of $H_{98,y2} = 1$ , SPARC is projected to attain $Q \approx 11$ and $P_{\textrm {fusion}} \approx 140$ MW. SPARC will therefore constitute a unique platform for burning plasma physics research with high density ( $\langle n_{e} \rangle \approx 3 \times 10^{20}\ \textrm {m}^{-3}$ ), high temperature ( $\langle T_e \rangle \approx 7$ keV) and high power density ( $P_{\textrm {fusion}}/V_{\textrm {plasma}} \approx 7\ \textrm {MW}\,\textrm {m}^{-3}$ ) relevant to fusion power plants. SPARC's place in the path to commercial fusion energy, its parameters and the current status of SPARC design work are presented. This work also describes the basis for global performance projections and summarizes some of the physics analysis that is presented in greater detail in the companion articles of this collection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Projections of H-mode access and edge pedestal in the SPARC tokamak

In order to inform core performance projections and divertor design, the baseline SPARC tokamak plasma discharge is evaluated for its expected H-mode access, pedestal pressure and edge-localized mode (ELM) characteristics. A clear window for H-mode access is predicted for full field DT plasmas, with the available 25 MW of design auxiliary power. Additional alpha heating is likely needed for H-mode sustainment. Pressure pedestal predictions in the developed H-mode are surveyed using the EPED model. The projected SPARC pedestal would be limited dominantly by peeling modes and may achieve pressures in excess of 0.3 MPa at a density of approximately 3 × 10 20 m -3 . High pedestal pressure is partially enabled by strong equilibrium shaping, which has been increased as part of recent design iterations. Edge-localized modes (ELMs) with >1 MJ of energy are projected, and approaches for reducing the ELM size, and thus the peak energy fluence to divertor surfaces, are under consideration. The high pedestal predicted for SPARC provides ample margin to satisfy its high fusion gain (Q) mission, so that even if ELM mitigation techniques result in a 2x reduction of the pedestal pressure, Q > 2 is still predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Multi-device study of temporal characteristics of magnetohydrodynamic modes initiating disruptions

Disruptions in tokamaks are often preceded by magnetohydrodynamic (MHD) instabilities that can rotate or become locked to the wall. Measurements from magnetic diagnostics in the presence of MHD mode precursors to disruptions can yield potentially valuable input to the plasma control system, with a view to disruption avoidance, prediction and mitigation. This report presents an exploratory analysis of the growth of MHD modes and corresponding time scales on the basis of magnetic measurements in multiple tokamaks. To this end, a database was compiled using disruptive discharges from COMPASS, ASDEX Upgrade, DIII-D and JET, manually classified according to disruption root cause, and characterized by a great diversity of operational conditions and mode dynamics. The typical time during which a mode can be detected using saddle coils and the duration of the locked mode phase in the database both extend over several orders of magnitude, but generally the time scales increase with plasma size. Several additional factors are discussed that can influence these durations, including the disruption root cause. A scaling law for the locked phase duration was estimated, yielding predictions toward ITER of the order of hundreds of milliseconds or even seconds. In addition, a scaling law for the mode amplitude at the disruption onset, proposed earlier by de Vries et al. (2016), is applied to the database, and its predictive capabilities are assessed. Despite significant uncertainty on the predictions from both scaling laws, encouraging trends are observed of the fraction of disruptions that may be detected with sufficient warning time to allow mitigation or even avoidance, based solely on observations of MHD mode dynamics. When combined with similar analysis of measurements from diagnostics that are sensitive to other disruption precursors, our analysis methods and results may contribute to the reliability, robustness and generalization of disruption warning schemes for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗