Prospects for Disruption Handling in a Tokamak-Based Fusion Reactor
Not provided.
Nuclear Science & Technology↗
SEARCH · Engineering Papers
Results for “tokamak disruption”
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.
Observation of non-thermal electron formation during the thermal quench of shattered pellet injection shutdowns in DIII-D
Formation of non-thermal (hot) electrons is studied during the thermal quench of shattered pellet injection (SPI) discharge shutdowns in the DIII-D tokamak. Both pure neon and mixed neon/deuterium shutdowns are studied. High-energy radiation (SXR and EUV) brightness levels show indications of non-thermal electron formation during the TQ. The non-thermal densities are of order 10 18 m –3 and non-thermal temperatures are estimated to be of order 10 keV. The non-thermal electrons appear to be mostly lost during the subsequent current quench, as no signal of relativistic runaway electrons is typically seen at the end of the current quench. In conclusion, the non-thermals appear to have an effect on the deposition profile of the SPI pellets, especially at the end of the TQ, due to increased pellet ablation.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The constraint of plasma power balance on runaway avoidance
Abstract In a post-thermal-quench plasma, mitigated or unmitigated, the plasma power balance is mostly between collisional or Ohmic heating and plasma radiative cooling. In a plasma of atomic mixture { n α } with α labeling the atomic species, the power balance sets the plasma temperature, ion charge state distribution { n α i } with i the charge number, and through the electron temperature T e and ion charge state distribution { n α i } , the parallel electric field E ∥ . Since the threshold electric field for runaway avalanche growth E a v is also set by the atomic mixture, ion charge state distribution and its derived quantity, the electron density n e , the plasma power balance between Ohmic heating and radiative cooling imposes a stringent constraint on the plasma regime for avoiding and minimizing runaways when a fusion-grade tokamak plasma is rapidly terminated.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Runaway electron current reconstitution after a nonaxisymmetric magnetohydrodynamic flush
Benign termination of mega-ampere (MA) level runaway current has been convincingly demonstrated in recent JET and DIII-D experiments, establishing it as a leading candidate for runaway mitigation on ITER. This comes in the form of a runaway flush by parallel streaming loss along stochastic magnetic field lines formed by global magnetohydrodynamic instabilities, which are found to correlate with a low-Z injection that purges the high-Z impurities from a post-thermal-quench plasma. Here, in this study, we show the competing physics that govern the postflush reconstitution of the runaway current in an ITER-like reactor where significantly higher current is expected. The trapped “runaways” are found to dominate the seeding for runaway reconstitution, and the incomplete purge of high-Z impurities helps drain the seed but produces a more efficient avalanche, two of which compete to produce a 2–3 MA step in current drop before runaway reconstitution of the plasma current.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Coupled 2D MHD and runaway electron fluid simulations of SPARC disruptions
Runaway electrons (REs) generated during disruption events in tokamaks can carry mega-Ampère level currents, potentially causing damage to plasma-facing components. Understanding RE evolution during disruption events is important for evaluating strategies to mitigate RE damage. Using two-dimensional toroidally symmetric magnetohydrodynamic (MHD) simulations in M3D-C1, which incorporates a fluid RE model evolved self-consistently with the bulk MHD fluid, we examine the seeding and avalanching of REs during disruptions in the SPARC tokamak – a compact, high-field, high-current device designed to achieve a fusion gain Q > 2 in deuterium–tritium plasmas. The M3D-C1 simulations of unmitigated disruptions demonstrate RE plateau formation and peaking of the final current density, which agree well with the results of lower-fidelity reduced RE fluid models. This work provides the first systematic comparison and benchmarking of different primary sources, including activated tritium beta decay and Compton scattering, in SPARC disruption simulations with self-consistent MHD and RE coupling.
Compton scattering↗
Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks
Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Density regulation with disruption avoidance in next-generation tokamaks using a safe reinforcement learning-based controller
Not Available
Paruchuri, Sai Tej (ORCID:0000000323723888)↗
Model predictive current profile control in tokamaks by exploiting spatially moving electron cyclotron current drives
The plasma control systems in next-generation tokamaks like ITER will balance competing control objectives to achieve the desired level of performance in advanced scenarios while preventing magnetohydrodynamic instabilities and disruptions. During normal tokamak operation, the points of incidence of the electromagnetic waves generated by the electron cyclotron heating and current drives (EC H&CDs) are usually fixed in space. However, the points of incidence can be modified in real-time by changing the angles of the mirrors that reflect the EC H&CD waves. Altering the points of incidence, in turn, varies the ability of the plasma control system to regulate a plasma property. For instance, changing the EC H&CD wave incidence location may place the power demands necessary to achieve a particular plasma target within saturation limits. Therefore, using the EC H&CD deposition location, which is related to the EC H&CD mirror angle, as a supplementary controllable variable may facilitate access to a given target scenario. However, active scenario-control algorithms have not been designed so far to fully exploit this capability in real time. In this work, a model predictive controller that can handle actuation locations as control inputs is developed. In particular, the controller is designed to regulate both the auxiliary powers and the EC H&CD deposition locations in a pre-defined optimal sense to achieve the control objective of attaining and sustaining a target current profile. Further, the proposed controller is tested for a DIII-D tokamak scenario in nonlinear simulations using the Control Oriented Transport SIMulator (COTSIM).
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Non-disruptive error field measurement in DIII-D low safety factor plasmas and projection to ITER
Previous experiments in DIII-D (Paz-Soldan et al 2022 Nucl. Fusion 62 126007) introduced a method to identify intrinsic error fields (EFs) in tokamaks with minimal disruption risk by promptly healing driven magnetic islands during the conventional 'compass scan'. This paper presents recent experimental and numerical advancements in extending this approach to low q 95 plasmas, and projects its applicability to ITER. Non-disruptive EF measurement is achieved at q 95 = 4.5 and 3.9 without any initial EF correction (EFC) by reducing the time between the occurrence of the locked mode (LM) and control action to 10 ms and increasing the density 50%–100%. However, 50% correction of the intrinsic EF is required to achieve island healing at q 95 = 3.2 with 10 ms delay for the control action. Nonlinear two-fluid modeling with the TM1 code reproduces the DIII-D experimental observations, indicating that promptly turning off the 3D coil current reduces both magnetic island width and electromagnetic force, while raising the density increases plasma viscosity, facilitating magnetic island healing. The simulations show that for scenarios with q 95 = 3.2, lowering the control action time to 5 ms will lead to island healing without EFC. TM1 simulations are extended to future ITER scenarios with 5 MA and 7.5 MA plasma currents, predicting the dependence of required density rise on action time and EF amplitude. These simulations indicate that, benefiting from the much longer resistive time, island healing can be successfully achieved in ITER when taking control action 100–500 ms after a LM occurrence.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Non-disruptive error field identification based on magnetic island healing
Here a technique to identify intrinsic error fields (EFs) in tokamaks with minimized risk of disruption is demonstrated on the DIII-D tokamak. The method extends the conventional driven magnetic island ‘compass scan’ approach by modifying asynchronous control waveforms to enable prompt healing of the island instability. Healing of the island is achieved by reducing the imposed non-axisymmetric coil current and raising the density (here via gas fueling). The method is also shown to support multiple island threshold measurements per pulse, thus reducing the number of dedicated pulses necessary to conduct an EF identification. Non-linear modeling with the TM1 code reproduces the experimental results and approximately recovers the critical density required for island healing. Island healing is explained in the non-linear modeling by an increase in the viscous coupling between the static island and the nearby flowing plasma, thus healing the island as it accelerates into the plasma frame. Due to both simplicity and risk minimization, this technique is suitable for plasma-based EF identification in the early commissioning stages of future disruption-averse tokamaks such as ITER and SPARC.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
Experimental inference of flux tunneling between magnetic island chains in tokamaks
We report confinement bifurcations due to flux tunneling in magnetic islands in toroidal plasmas. This occurs in 2/1 islands, triggered by decoupling of non-overlapping adjacent 5/2 islands. While subject to electron cyclotron current drive (ECCD), the temperature at the 2/1 O-point is flat before but peaks after decoupling. Here, this is shown to be due to a bifurcation from stochastic to nested magnetic structure, caused by the removal of flux tunneling through intersecting manifolds of overlapping heteroclinic tangles. 2/1 island stability correlates with coupling events in an ITER baseline scenario plasma where ECCD can not stabilize the coupled islands, showing the critical impact of flux tunneling on disruption avoidance in tokamaks.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Tech-X contributions to the Center for Tokamak Transients (CTTS) Scientific Discovery through Advanced Computing (SciDAC) project (Final Report)
The Center for Tokamak Transients Simulations (CTTS) consists of computational plasma physicists, computer scientists, and applied mathematicians who are together developing and applying simulation software to improve our understanding of and ability to predict and control transient events in tokamak discharges, particularly those that lead to violent disruptions. We have major simulation activities concerning disruptions that are caused by violating ideal MHD stability criteria, vertical displacement events, resistive wall modes, and tearing modes. We are performing comprehensive modeling of disruption mitigation techniques, particularly by shattered pellet injection (SPI). The Center is built around two mature extended MHD codes: NIMROD and M3D-C1. These codes are being extended in a number of ways as required by the demands of disruption modeling, and are being optimized to work well for larger problem size and modern computer architectures as required to model tokamak plasmas with realistic parameters over the required timescales. The codes will be coupled to the FronTier code for state-of-the-art modeling of SPI. Validation with experimental data on DIII-D, NSTX, and JET is being facilitated by rigorous uncertainty quantification (UQ). The high-level objectives of the center are to: (1) Develop better understanding and improved predictive capability of when crossing a linear stability boundary will likely lead to a disruption, (2) Provide a tool for calculating vessel forces, both axisymmetric and non-axisymmetric, for worst case VDE and other disruptions in tokamaks, (3) Develop increased understanding and improved predictive capability for locked mode disruptions and how best to avoid them, (4) Develop, verify, and validate 3D models for disruption mitigation by shattered pellets that can be used to design and optimize a system for future tokamaks, and (5) Increase the efficiency and scalability of the M3D-C1 and NIMROD codes so that they can model disruptions in ITER scale devices using the most powerful DOE HPC hardware, and (6) Demonstrate code-coupling via FronTier coupling with NIMROD and M3D-C1 and by interfacing these codes within a WDM framework.
97 MATHEMATICS AND COMPUTING↗
Effects of electron viscosity on resonant layer responses to non-axisymmetric magnetic perturbations
The resonant field penetration to magnetic islands is the central MHD mechanism of non-axisymmetric plasma responses in a tokamak such as disruptive locking or favorable ELM stabilization. The resonant field penetration can be induced by any non-ideal processes as manifested in the delicate balance under the generalized Ohm's law. Here, we show that the viscous effects by electrons are not ignorable in the field penetration unlike previous presumption, even if the electron viscosity is as small as the square root of its mass compared to the ions. It is clear that its effects become only bigger if the electron viscosity becomes anomalously large. The work strictly follows the three-field model in the linear regime targeting the prediction of the onset of the field penetration and successfully extending it with electron viscosity and identifying new regimes. The results also indicate that the error field thresholds become more strongly dependent on plasma density than ones predicted in the linear regimes without the electron viscosity, which is consistent with experimental observations and thus a significant implication.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Issues in Formation of Cryogenic Pellets for Fusion Applications
Cryogenic pellets are used for injection into fusion plasmas to add fuel to build up density and replace the ions lost from fusion reactions and imperfect confinement in the plasma. These pellets are formed at cryogenic temperatures with pure hydrogenic isotopes or mixtures of the isotopes. Technology to make these pellets and inject them into plasmas has been under development for many years, and various methods using freezing or desublimation have been shown to produce high-quality solid pellets suitable for injection. The throughput needed and possible impurity content from the necessary recirculation of fusion exhaust gases are two of the key issues to overcome for fusion pellet fueling systems in long-pulse burning plasmas. Here, we describe the technical challenges associated with these issues and the capability of pellet formation extruders to overcome them.Cryogenic pellets of deuterium, neon, and argon are also used in fusion tokamak devices for disruption mitigation in the form of large pellets that can be injected on demand to quickly dissipate the plasma thermal energy through radiation and add significant density in order to prevent runaway electron formation. Here, the issue is not throughput as with the fueling pellets but rather is the time it takes to form pellets of the size needed and the ability to dislodge them immediately on demand when needed to mitigate a disruption. The method used to make these pellets by desublimation is described, and examples related to how pellet size and input gas parameters affect the formation time are provided.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Pressure Response Optimization of an Eddy Current-Driven Flyer Plate Valve for the ITER Shattered Pellet Injection System
One technique for mitigating disruptions in a tokamak is shattered pellet injection (SPI). SPI is a process in which a large solid pellet consisting of deuterium, neon, or argon is desublimated in a pipe gun barsrel and launched downstream. Pellets are shattered just before entering the plasma by an impact with an angled tube. Injection of these materials into the plasma radiates stored thermal energy, limits current decay rates, suppresses the generation of runaway electrons, and dissipates runaway electrons if necessary. A critical element of the SPI system is a fast-acting valve that releases high-pressure gas to dislodge and accelerate pellets directly, or indirectly via a mechanical punch. A prototype valve sized for the ITER SPI system has been designed and fabricated. A pulsed high-voltage power supply energizes the valve’s internal magnetic coil, which induces eddy currents in the adjacent flyer plate resulting in a repulsive force between the flyer plate and the coil. The flyer plate action lifts a valve seat, allowing high-pressure gas to flow from the valve plenum to the downstream (breech) location of the pellet or mechanical punch. All of the valve’s internal components are designed to operate in ITER-level static background magnetic fields. Here, a study was conducted to optimize the downstream pressure response for a range of valve sizes and operating pressures. In particular, the study analyzes the breech pressure response associated with varying plenum pressures as well as varying breech volumes. A computational fluid dynamics simulation was built in STAR-CCM+ and validated against data from laboratory experiments. The resulting simulation outputs, in the form of downstream responses for a variety of initial plenum pressures and breech volumes, will be used as a complement to experimental data to ensure the pressure pulse is suitable for pellet survivability. These data, combined with studies on pellet shear strength and shock response, will be applied to optimization of overall operating parameters of the ITER SPI system.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The magnetohydrodynamics of current sheets
Examples of current sheets are summarized and their formation is described. A universal phenomenon in cosmic plasmas is the creation of sheets off intense current near X-type neutral points (where the magnetic field vanishes). These sheets are important as sites where the magnetic-field energy is converted efficiently into heat and bulk kinetic energy and where particles can be accelerated to high energies. Examples include disruptions in laboratory tokamaks, substorms in the earth's magnetosphere, and flares on the sun. The basic behavior of a one-dimensional sheet is presented, together with an account of the linear tearing-mode instability that can cause the field lines in such a sheet to reconnect. Such reconnection may develop in different ways: it may arise from a spontaneous instability or it may be driven, either from outside by motions or locally by a resistivity enhancement. Various processes are described that may occur during the nonlinear development of tearing, along with the many numerical and laboratory experiments that are aiding our understanding of this intriguing cosmical process.
Priest, E. R.↗