Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “massive gas injection”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparison of different noble gas injections by massive gas injection on plasma disruption mitigation on Experimental Advanced Superconducting Tokamak

Massive gas injection (MGI) is a traditional plasma disruption mitigation method. This method directly injected massive gas into the pre-disruption plasma and had been developed on the Experimental Advanced Superconducting Tokamak (EAST). Different noble gas injection experiments, including He, Ne, and Ar, were performed to compare the mitigation effect of plasma disruption by evaluating the key parameters such as flight time, pre-thermal quench (pre-TQ), and current quench (CQ). The flight time was shorter for low atomic number ( Z ) gas, and the decrease in flight time by increasing the amount of gas was insignificant. However, both pre-TQ and CQ durations decreased considerably with the increase in gas injection amount. The effect of atomic mass on pre-TQ and CQ durations showed the opposite trend. The observed trend could help in controlling CQ duration in a reasonable area. Moreover, the analysis of radiation distribution with different impurity injections indicated that low Z impurity could reduce the asymmetry of radiation, which is valuable in mitigating plasma disruption. These results provided essential data support for plasma disruption mitigation on EAST and future fusion devices.

Physics↗
DOE OSTI2021

MHD simulations of cold bubble formation from 2/1 tearing mode during massive gas injection in a tokamak

Massive gas injection (MGI) experiments have been carried out in many tokamaks to study disruption dynamics and mitigation schemes. Two events often observed in those experiments are the excitation of the m = 2, n = 1 magnetohydrodynamic mode, and the formation of cold bubble structure in the temperature distribution before the thermal quench (TQ). Here m is the poloidal mode number, n the toroidal mode number. The physics mechanisms underlying those phenomena, however, have not been entirely clear. In this work, our recent NIMROD simulations of the MGI process in a tokamak have reproduced the main features of both events, which has allowed us to examine and establish the causal relation between them. In these simulations, the 3/1 and 2/1 islands are found to form successively after the arrival of impurity ion cold front at the corresponding q = 3 and q = 2 rational surfaces. At the interface between impurity and plasma, a local thin current sheet forms due to an enhanced local pressure gradient and moves inward following the gas cold front, this may contribute to the formation of a dominant 2/1 mode. Following the growth of the 2/1 tearing mode, the impurity penetration into the core region inside the q = 2 surface gives rise to the formation of the cold bubble temperature structure and initiates the final TQ. Here, a subdominant 1/1 mode developed earlier near the q = 1 surface alone does not cause such a cold bubble formation, however, the exact manner of the preceding impurity penetration depends on the nature of the 1/1 mode: kink-tearing or quasi-interchange.

2/1 tearing mode↗
DOE OSTI2026

Radiation asymmetry during shattered pellet and massive gas injection in DIII-D

Infrared thermography of the first wall in DIII-D is used to show the toroidal peaking of plasma radiation during mitigated disruptions with shattered pellet injection (SPI) and with massive gas injection (MGI). During MGI, the radiation peak location is shown to be due to the n = 1 magnetohydrodynamic (MHD) activity and continuously variable in toroidal phase based on externally applied error fields, consistent with previous experimental results. Furthermore, the measured toroidal peaking factor (TPF) is in agreement with the previous estimates based on radiometry, with a total TPF of 1.07 ± 0.05 when integrated over the entire duration of the disruption, and a value of 1.3 ± 0.1 during the thermal quench (TQ). For SPI, the location of the radiation peak is found to be determined predominantly by the particle source location, and thus can be varied on DIII-D by injecting from either of two toroidally separated injectors. 3D non-linear MHD simulations with the NIMROD code support this picture. The localized SPI particle source leads to higher peaking than for MGI, with a total TPF of 1.3 ± 0.1 over the full disruption and a TQ TPF of 1.9 +0.5/−0.3, a value also consistent with the NIMROD simulations. This TQ value is close to or potentially exceeding the allowable limit for the previously planned Be first wall in ITER, pointing to the importance of radiation asymmetries depending on the choice of first wall material.

Physics - Plasma physics↗
DOE OSTI2024

Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC

Recent developments to the M3D-C1 code enable higher fidelity modeling of disruptions, and can be applied in the design verification of reactor-scale tokamaks. Among these new capabilities is a method to mesh conducting vessel structures such as coils and passive plates, packing of the toroidal mesh around gas injectors, as well as anisotropic resistivity inside the vessel structures. We present extended-magnetohydrodynamic (MHD) simulations of disruption mitigation via massive gas injection (MGI) in SPARC. The goal of this study is to inform the disruption mitigation layout of SPARC and aid in the design of an effective gas injector configuration. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures around the plasma. A localized mixture of deuterium and neon with a small toroidal and poloidal width is injected in up to six locations. We demonstrate that M3D-C1 can model a rapid shutdown via MGI using narrow and more realistic gas plumes than in previous simulations. As a result of the q = 1 surface in the SPARC baseline case a sawtooth is observed early in the simulations. Despite the sawtooth and the onset of edge MHD instabilities, the impurity distribution remains localized around the injector locations, but enables a radiative shutdown of the plasma. We find that using the maximum of six gas injectors results in a lower peaking factor and leads to a more even distribution of radiation toroidally than using two injectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Study of argon expulsion from the post-disruption runaway electron plateau following low-Z massive gas injection in DIII-D

Herein, a 1D radial diffusion model is developed to study the observed rapid expulsion of argon from the runaway electron plateau in the DIII-D tokamak following secondary massive low-Z (D 2 or He) gas injection. The expulsion of argon is found to be caused by further cooling of the background plasma due to the added neutrals, accompanied by recombination of argon ions and greatly increased outward radial transport rate of argon (now dominantly in neutral form) out of the runaway electron beam. After argon expulsion, power loss out of the runaway electron plateau is found to be dominated by neutral transport to the wall (rather than by radiation); this result resolves the power balance discrepancy highlighted in previous work on argon expulsion out of the runaway electron plateau.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Enhanced plasma current spike formation due to onset of 1/1 kink-tearing reconnection during a massive gas injection process

The formation of the plasma current spike at the end of the thermal quench (TQ) phase is studied systematically, which is found to strongly correlate with the onset of the 1/1 kink-tearing reconnection only in the simulation results presented here. Here, the magnetohydrodynamic activity on the q = 1 surface plays a critical role in the spike formation and the disruption process. In particular, when the safety factor in the magnetic axis q 0 exceeds 1, the plasma major disruption transits into successive minor disruptions, and the start of the TQ phase is delayed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

ARC disruption physics and strategy

Commonwealth Fusion Systems (CFS) plans to operate a tokamak power plant called ARC in the early 2030s. Tokamak plasmas have stability limits that, if crossed, lead to a rapid termination of the plasma, referred to as a disruption. Disruptions pose a melt risk to the first wall resulting from thermal and non-thermal particle heat fluxes, and an electromagnetic loading risk on all metal components within the equilibrium coils. A comprehensive set of models is used herein to provide an assessment of both mitigated and unmitigated ARC disruption loads. A preliminary massive gas injection system is baselined and a runaway electron mitigation coil option is proposed to close possible gaps in the baseline. It is predicted that all ARC disruption loads are within a factor of 2 of the disruption loads in SPARC, a tokamak presently under construction by CFS, and therefore SPARC provides an opportunity to calibrate models, test solutions and inform the design of ARC. The goal for ARC is disruption-free operation, however, the pragmatic design target is to withstand one mitigated disruption per day, and to restart the plasma following mitigation in tens of seconds without interrupting the power output. Unmitigated disruptions must be rare, and experience with unmitigated disruption impacts in SPARC will better define what rare means. The implications of this strategy for plasma disruptivity and disruption prediction are discussed, and operating the ARC scenario on SPARC is expected to refine the ARC final design and operational plan.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

First results of Ne shattered pellet injection for mitigating plasma disruption with full metal wall in EAST tokamak

Disruption mitigation poses a significant and unresolved challenge for ITER and future fusion reactor devices. To mitigate the effect of plasma disruption, a Shattered Pellet Injection (SPI) system has been developed and the first rapid shutdown experiments have been successfully performed in the Experimental Advanced Superconducting Tokamak with full metal wall. The experiments confirmed the system’s capability to inject 5 mm diameter neon (Ne) pellets at velocities ranging from 100 to 400 m s −1 , with a maximum injected quantity of 13.2 Pa m 3 . Compared to unmitigated disruptions, the total radiation power was significantly higher with the respective use of SPI and Massive Gas Injection (MGI). Additionally, the radiation distribution and divertor heat flux were compared between SPI and MGI methods. The results demonstrated that SPI exhibited shorter cooling time, stronger core radiation and more uniform poloidal radiation distribution compared to MGI, indicating deeper deposition using SPI. Furthermore, during plasma shutdown, the electron temperature and peak heat flux near outer divertor strike points were reduced by 40% and 50% respectively, with SPI, in comparison to MGI. These findings serve as a valuable reference for implementing SPI technology as the baseline approach for disruption mitigation in ITER.

EAST↗
DOE OSTI2020

Pellet sublimation and expansion under runaway electron flux

This work provides a qualitative description of the pellet response to the ambient runaway electrons. For ITER-relevant parameters, our estimates suggest that the cryogenic pellets will be sublimated instantly at the edge of the runaway beam. The subsequent rapid expansion of the sublimated material spreads the impurities over the poloidal cross-section of a tokamak on a millisecond time scale prior to the complete ionization of the expanding cloud. Here, the injected solid pellet turns into a rapidly expanding gas cloud before it reaches the core of the runaway beam. As a result, the pellet acts similar to the massive gas injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

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

Shattered Pellet Injection (SPI) was previously 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. On account 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

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 OSTI2024

Simulation of DIII-D disruption with argon pellet injection and runaway electron beam

Abstract The next generation of large tokamaks, including ITER, will be equipped with a disruption mitigation system (DMS) that can be activated if a disruption is deemed to be imminent. Introducing impurities by pellet (large or shattered) or massive gas injection has been shown to be an effective mitigation mechanism on many tokamaks. The goal of the mitigation is to lessen the thermal and electromagnetic loads from the disruption without generating enough high-energy (runaway) electrons to damage the device. Variations of this mitigation process with impurity injection are presently being tested on many experiments. We have modeled one such impurity injection experiment on DIII-D using the M3D-C1 nonlinear 3D extended MHD code (Jardin et al 2012 Comput. Sci. Discovery 6 014002), The model includes an argon large pellet injection and ablation model, impurity ionization, recombination, and radiation, and runaway electron formation and subsequent evolution, including both Dreicer and avalanche sources. We obtain reasonable agreement with the experimental results for the timescale of the thermal and current quench and for the magnitude of the runaway electron plateau formed during the mitigation. This is the first 3D full MHD simulation with pellets and REs to simulate the disruption process and it also provides a partial validation of the M3D-C1 DMS model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

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’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2022

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or –90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s –1 . The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

divertor↗