Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fusion Reactor”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

W-band system-on-chip electron cyclotron emission imaging system on DIII-D

Monolithic, millimeter wave “system-on-chip” (SoC) technology has been employed in heterodyne receiver integrated circuit radiometers in a newly developed Electron Cyclotron Emission Imaging (ECEI) system on the DIII-D tokamak for 2D electron temperature profile and fluctuation evolution diagnostics. Here, a prototype module operating in E-band (72-80 GHz) was first employed in a 2 x 10 element array that demonstrated significant improvements over the previous quasi-optical Schottky diode mixer arrays during the 2018 operational campaign of the DIII-D tokamak. For compatibility with International Thermonuclear Experimental Reactor (ITER) relevant scenarios on DIII-D, the SoC ECEI system was upgraded with 20 horn-waveguide receiver modules. Each individual module contains a University of California at Davis designed W-band (75 -110 GHz) receiver die that integrates a broadband low noise amplifier (LNA), a double balanced down-converting mixer, and a x4 multiplier on the local oscillator (LO) chain. A x2 multiplier and two IF amplifiers are packaged and selected to further boost the signal strength as well as downconvert the signal frequency. The upgraded W-band array exhibits > 30 dB additional gain and 20x improvement in noise temperature comparing with the previous Schottky diode radio frequency (RF) mixer input systems; an internal 8 times multiplier chain is used to bring down the LO frequency below 12 GHz, thereby obviating the need for a large aperture for quasi-optical LO coupling and replacing it with coaxial connectors. The horn-waveguide shielding housing avoids out-of-band noise interference on each individual module. The upgraded ECEI system plays an important role for absolute electron temperature evolution and fluctuation measurements for edge and core region transport physics studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Benchmarking core turbulence and transport predictions for an inductive compact tokamak reactor plasma

Motivated by the need for accurate, timely, and efficient calculations of plasma transport, predictions of plasma turbulence properties made using different TGLF saturation rules are benchmarked against corresponding predictions from linear and nonlinear gyrokinetic CGYRO simulations. This benchmarking is carried out using parameters taken from an inductive burning plasma scenario in a hypothetical compact high-field (R maj = 4 m, B T = 8 T) tokamak, lying in a much different regime of parameter space than either the TGLF calibration regime or current-day experiments. The core turbulent transport in this scenario is predicted to be dominated by ion temperature gradient (ITG) turbulence. In general, the ITG critical gradients predicted by various TGLF saturation rules are quite close to the CGYRO predictions. Both codes predict similar linear ITG growth rates and frequency spectra, as well as their scaling with R/L T i = −Rd ln(T i )/dr. However, TGLF systematically predicts unstable trapped-electron modes (TEMs) above k y ρ s ≃ 0.5 not seen by CGYRO for the same parameters, due to TGLF predicting a lower threshold in R/L T e than CGYRO for TEM onset. It is shown that for this scenario, nonlinear CGYRO simulations predict stiffer ITG turbulence than the TGLF SAT0 and SAT1 saturation rules, with energy fluxes close in magnitude and scaling with R/L T i to what is predicted by the SAT2 saturation rule. Self-consistent core profiles calculated using nonlinear CGYRO flux predictions and the PORTALS transport solver are shown to agree fairly well with corresponding predictions made using the TGLF SAT2 model, including a similar level of density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design considerations for high entropy alloys in advanced nuclear applications

The demanding operating environments of advanced nuclear reactors require the development of new nuclear materials that can withstand their increased physical, chemical, thermal, and radiation-related challenges. High-entropy alloys (HEAs) have shown often-impressive mechanical, thermomechanical, and corrosion-resistant properties, and offer a massive, unexplored compositional space which allows for the targeted development of application-specific materials. Furthermore, although still in a nascent stage, research has shown that HEAs may exhibit unique irradiation tolerance, including reduced defect production and resistance to irradiation-induced swelling and hardening. The mechanisms behind this increased tolerance are not yet well-understood, although the HEA-specific attributes of a complex energy landscape, reduced thermal conductivity, and shift in defect migration energies and pathways provide promising explanations. Here this work assesses the current and future challenges faced by structural nuclear materials, identifying the specific applications in which HEAs may provide a competitive advantage compared to industry-standard materials with the aid of Ashby material selection maps. Considerations are provided for the design of future nuclear HEAs, including calculations of nuclear-relevant properties to assist in the initial down-selection of elements depending on application requirements (e.g., low neutron capture for in-core applications), narrowing the existing compositional space of HEAs to a manageable scope.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Beryllium-driven structural evolution at the divertor surface

Erosion of the beryllium first wall material in tokamak reactors has been shown to result in transport and deposition on the tungsten divertor. Experimental studies of beryllium implantation in tungsten indicate that mixed W–Be intermetallic deposits can form, which have lower melting temperatures than tungsten and can trap tritium at higher rates. To better understand the formation and growth rate of these intermetallics, we performed cumulative molecular dynamics (MD) simulations of both high and low energy beryllium deposition in tungsten. In both cases, a W–Be mixed material layer (MML) emerged at the surface within several nanoseconds, either through energetic implantation or a thermally-activated exchange mechanism, respectively. While some ordering of the material into intermetallics occurred, fully ordered structures did not emerge from the deposition simulations. Targeted MD simulations of the MML to further study the rate of Be diffusion and intermetallic growth rates indicate that for both cases, the gradual re-structuring of the material into an ordered intermetallic layer is beyond accessible MD time scales(≤1 μs). However, the rapid formation of the MML within nanoseconds indicates that beryllium deposition can influence other plasma species interactions at the surface and begin to alter the tungsten material properties. Therefore, beryllium deposition on the divertor surface, even in small amounts, is likely to cause significant changes in plasma-surface interactions and will need to be considered in future studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Full-wave simulations on helicon and parasitic excitation of slow waves near the edge plasma

Helicon waves are thought to be promising in various tokamaks, such as DIII-D, because they can penetrate reactor-grade high-density cores and drive the off-axis current with higher efficiency. In the frequency regime ~476 MHz, both slow electrostatic and fast electromagnetic helicon waves can coexist in DIII-D. If the antenna parasitically excites the slow mode, these waves can propagate along the magnetic field line into the scrape-off layer (SOL). Although the importance of the misalignment of the Faraday screen and the electron density in the SOL on the excitation and propagation of slow modes is well known, the conditions for minimizing slow mode excitation have yet to be optimized. Using the Petra-M simulation code in the 2D domain, we analyze the effects of the misalignment of the antenna in the poloidal direction, the misalignment of the Faraday screen in the toroidal direction, and the density in front of the antenna on slow mode generation. Our results suggest that the misalignment of the Faraday screen is a critical factor in reducing the slow mode and that the misalignment angle should be below ~5° to minimize the slow wave excitation. When the electron density is higher than 3.5 x 10 18 m –3 in the SOL, the generation of the slow mode from the antenna is minimized and unaffected by the misalignment of the Faraday screen.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Magnetic probe diagnostics of nonlinear standing waves and bulk ohmic electron power absorption in capacitive discharges

It is recognized that standing wave effects appearing in large-area, very-high-frequency capacitively coupled plasma (CCP) reactors cause center-high plasma non-uniformity. Using a high-frequency magnetic probe, we present a direct experimental diagnostic of the nonlinear standing waves and bulk ohmic electron power absorption dynamics in low pressure CCP discharges for different driving frequencies of 13.56, 30, and 60 MHz. The design, principle, calibration, and validation of the probe are described in detail. Spatial structures of the harmonics of the magnetic field, determined by the magnetic probe, were used to calculate the distributions of the harmonic current and the corresponding ohmic electron power deposition, providing insights into the behavior of nonlinear harmonics. At a low driving frequency, i.e. 13.56 MHz, no remarkable nonlinear standing waves were identified and the bulk ohmic electron power absorption was observed to be negligible. Here, the harmonic magnetic field/current was found to increase dramatically with the driving frequency, due to decreased sheath reactance and more remarkable nonlinear standing waves at a higher driving frequency, leading to the enhancements of the ohmic heating and the plasma density in the bulk, specifically at the electrode center. At a high driving frequency, i.e. 60 MHz, the high-order harmonic current density and the corresponding ohmic electron power absorption exhibited a similar node structure, with the main peak on axis, and one or more minor peaks between the electrode center and the edge, contributing to the center-high profile of the plasma density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Power matching to pulsed inductively coupled plasmas

Matching of power delivery to nonlinear loads in plasma processing is a continuing challenge. Plasma reactors used in microelectronics fabrication are increasingly multi-frequency and/or pulsed, producing a non-linear and, in many cases, non-steady state electrical termination that can complicate efficient power coupling to the plasma. This is particularly the case for pulsed inductively coupled plasmas where the impedance of the plasma can significantly change during the start-up-transient and undergo an E–H (capacitive-to-inductive) transition. In this paper, we discuss the results from a computational investigation of the dynamics of power matching to pulsed inductively coupled plasmas (Ar/Cl 2 mixtures of tens of mTorr pressure) using fixed component impedance matching networks and their consequences on plasma properties. In this investigation, we used set-point matching where the components of the matching network provide a best-case impedance match (relative to the characteristic impedance of the power supply) at a chosen time during the pulsed cycle. Matching impedance early during the pulse enables power to feed the E-mode, thereby emphasizing capacitive coupling and large excursions in the plasma potential. This early power coupling enables a more rapid ramp-up in plasma density while being mismatched during the H-mode later in the pulse. The early match also produces more energetic ion bombardment of surfaces. Matching late in the pulse diminishes power dissipated in the E-mode at the cost of also reducing the rate of increase in plasma density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Improvements to the High-Field-Side Transient CHI System on QUEST

Transient coaxial helicity injection (CHI) is a promising candidate for solenoid-free plasma current start-up in a low-aspect-ratio tokamak in support of developing fully non-inductive scenarios. The aim of the transient CHI research on QUEST is to develop a reactor-relevant CHI design. On QUEST, a CHI discharge is initiated by driving current along magnetic field lines that connect an electrically insulated electrode plate (which is referred to as a bias electrode) to a vessel component at the bottom region of the spherical tokamak. In the first application of CHI on QUEST, the electrically insulated electrode plate was biased with respect to the outer vessel in a configuration referred to as low-field-side (LFS) injection. To maintain a narrow injector flux footprint width throughout the discharge, QUEST is now developing a high-field-side (HFS) injection configuration, in which the electrically insulated electrode plate is biased with respect to the inner vessel components. Through the implementation of a CHI-dedicated gas injection system, studies in the HFS injection configuration have now demonstrated good magnetic flux evolution into the vacuum vessel. Toroidal currents up to 43 kA have been achieved, and the generated current has increased in proportion to the magnetic flux connecting the electrodes. Finally, these results which show agreement with the CHI-scaling suggest that much higher levels of toroidal current can be generated on QUEST in an optimized CHI system in which the magnetic flux connecting the CHI electrodes is further increased.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Investigation of Drift Effects in UEDGE Simulations of NSTX-U Edge Plasma With Lithium Divertors

Lithium is a low-Z material, and lithium-based plasma-facing components (PFCs) are planned for the National Spherical Torus Experiment Upgrade (NSTX-U) to explore potential benefits for divertor power exhaust and core plasma management. NSTX-U is a medium-sized spherical tokamak with up to 12 MW of auxiliary heating, capable of generating reactor-relevant plasma conditions. This work presents boundary plasma simulations for NSTX-U with lithium PFCs using the UEDGE code, incorporating full magnetic and 𝐄 ×𝐁 drift physics. The simulations show that drifts strongly influence heat and particle transport: they enhance convective transport, broaden the scrape-off layer heat-flux width 𝜆 𝑞 , and reduce the anomalous heat diffusivity 𝜒 required to reproduce predicted SOL heat-flux width. 𝐄 ×𝐁 drifts provide poloidal transport, while ∇𝐵 (which includes both gradB and curvature) drifts provide radial heat and particle transport. Lithium transport is also affected by drifts, with lithium ions migrating from the outer divertor to the inner divertor through the private flux region (PFR) following the 𝐄 ×𝐁 drifts flow, lowering upstream impurity lithium densities. UEDGE is self-consistently coupled with the Wall-Li model to study plasma lithium PFC interactions depending on the local lithium sourcing based on local plasma conditions and lithium surface temperature. In these simulations, lithium evaporation shows a vapor-shielding effect that reduces divertor heat flux and increases radiative losses once surface temperatures exceed 450°C. This research work provides a first step toward self-consistent modeling of lithium PFCs in NSTX-U, demonstrating the impact of drift-driven plasma transport in SOL and divertor regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry

Plasma facing components (PFCs) in the next generation of tokamak devices will operate in challenging environments, with heat loads predicted to exceed 10 MWm -2 . The magnitude of these heat loads is set by the width of the channel, the ‘scrape-off layer’ (SOL), into which heat is exhausted, and can be characterised by an e-folding length scale for the decay of heat flux across the channel. It is expected this channel will narrow as tokamaks move towards reactor relevant conditions. Understanding the processes involved in setting the SOL heat flux width is imperative to be able to predict the heat loads PFCs must handle in future devices. Measurements of the SOL width are performed on the high-field spherical tokamak, ST40, using a newly commissioned infrared thermography system. With its high on-axis toroidal magnetic field (≥1.5 T) ST40 is uniquely positioned to investigate the influence of toroidal field on the heat flux width in spherical tokamaks, whilst also extending measurements of the SOL width in spherical tokamaks to increased poloidal field (≥0.3 T). Due to the divertor on ST40 having a low degree of axisymmetry, it is necessary for a set of radial measurements of the heat flux to be taken across the divertor, made possible using an automated toolchain that fully incorporates its 3D geometry. These radial profiles are combined with the magnetic geometry of the plasma to infer the width of the SOL, with both Eich and double exponential profiles of heat flux observed. A reduction in the heat flux is observed toroidally across part of the divertor, along with increased heat loads observed locally around the edges of the tiles. Future work in characterising the impact of tile misalignment and uncertainties in the reconstructed divertor magnetic geometry is required in order to further understand the observed heat flux patterns, as are additional investigations into the role potentially being played by an inhomogeneous sheath electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

DECAF Cross-device characterization of tokamak disruptions indicated by abnormalities in plasma vertical position and current

Abnormal (deviating from target) variations in the plasma vertical position Z and current I p (such as vertical displacements, transient I p 'spikes' and quenches) constitute common elements of a disruption, a phenomenon that is to be mitigated, or ultimately avoided in future reactor-relevant tokamaks. While those abnormalities are generally recognized cross-shot and cross-device, details in terms of appearance (or not) and order of those abnormalities in disruption event chains is bound with the plasma state at the time of the chain initiation. Detection of those abnormalities is thus indicative not only of the onset of the plasma collapse itself, but also of the disruption driving cause that is promoted at the particular plasma state. Here, occurrence of disruptions, explored via detection of a Ip quench, and analysis of disruption event chains constituted by Ip and Z abnormalities, is reported for in total 7 full years of operation of 3 devices (KSTAR, MAST-U and NSTX-U) using the DECAFTM code expanded tools and capabilities. It is shown that the disruption occurrence depends not only on details of the plasma state, but also on (device-dependent) technical elements of the shot exit scenario. A year-to-year change in main disruption causes and a reduction of the disruptivity rate, bound with device and operation upgrades, are reported. Particular trigger instances of disruption event chains (and the full chains, when applicable) are shown to occupy different parts of the operation space diagrams, in accordance with prior expectation. Plasma elongation is identified as an important factor influencing details of the chains and its role will be further explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

High performance power handling in the absence of an H-mode edge in negative triangularity DIII-D plasmas

Abstract Experiments performed during strongly-shaped high-power diverted negative triangularity (NT) experiments in DIII-D achieved detached divertor conditions and a transient-free edge, showcasing the potential for application of NT to a core-edge integrated reactor-like scenario and providing the first characterization of the parametric dependencies for detachment onset. Detached divertor conditions will be required in future devices to mitigate divertor heat fluxes. Access to dissipative divertor conditions was investigated via an increase in upstream density. Detachment onset at the outer strike point was achieved with H-mode level confinement H 98 − y 2 ∼ 1 and reactor-relevant normalized pressures β N ∼ 2 . Confinement degradation was observed with deeper detachment, associated with the loss of an electron temperature pedestal. Differences in geometry, radial transport, impact of cross field drifts are discussed to explain differences in access to detachment in NT discharges. Higher normalized densities, with respect to equivalent discharges in positive triangularity, were necessary to achieve detachment, partially explained by the shorter parallel connection length to the targets. The effect of cross-field particle drifts ( E × B , B × ∇ B ) on access to detachment was demonstrated by the lower upstream density needed to access detachment with ion B × ∇ B drift directed outside of the active divertor (Greenwald fraction f Gw ∼ 0.9–1.0 vs f Gw ∼ 1.3). The upstream density at detachment onset was observed to increase linearly with plasma current with ion B × ∇ B drift into the divertor, consistent with the observed narrowing of the scrape-off layer heat flux width λ q . Edge fluid simulations capture separatrix densities needed to achieve detachment in NT plasma and their dependence on drift direction. The ability to reproduce detachment dynamics in NT plasma increases the confidence in future design studies for NT divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

From MAST-U to STEP: Power exhaust control challenges and opportunities

This work explores the challenges and opportunities for power exhaust control in STEP, informed by dedicated MAST-U experiments. The STEP system primarily relies on hydrogenic and argon gas injection into the divertors, compensating for transients originating in the plasma core. Fast transients (>10 Hz) that evolve too quickly for actuators to compensate must be absorbed passively; the enhanced transient buffering of the foreseen long-legged divertor provides a clear advantage over conventional geometries, as demonstrated in the MAST-U Super-X divertor. Experiments further indicate that STEP’s tight divertor baffling enables near-independent control of the upper and lower divertors. The most challenging transients are expected from power-sharing fluctuations, as MAST-U experiments observed extremely fast dynamics. Fluctuations from core pellet fuelling remain relatively benign due to the small pellet size according to simulations. Exhaust control in a reactor like STEP requires an integrated approach befitting its machine-critical nature. A predictive control element, integrated with core plasma control, allows pre-emptive preparation of the divertor for incoming transients. An observer-based diagnostic approach is foreseen to monitor the divertor in the challenging reactor environment, supported by an extended diagnostic set in the non-nuclear phase to validate the required dynamic models.

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