Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ionizing radiation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

KATRIN background due to surface radioimpurities

The goal of the KArlsruhe TRItrium Neutrino (KATRIN) experiment is the determination of the effective electron antineutrino mass with a sensitivity of 0.2 eV/c 2 at 90 % C.L. 1 . This goal can only be achieved with a very low background level in the order of 10 mcps 2 in the detector region of interest. A possible background source are α-decays on the inner surface of the KATRIN Main Spectrometer. Rydberg atoms, produced in sputtering processes accompanying the α-decays, are not influenced by electric or magnetic fields and freely propagate inside the vacuum of the Main Spectrometer. Here, they can be ionized by thermal radiation and the released electrons directly contribute to the KATRIN background. Two α-sources, 223 Ra and 228 Th, were installed at the Main Spectrometer with the purpose of temporarily increasing the background in order to study α-decay induced background processes. In this paper, we present a possible background generation mechanism and measurements performed with these two radioactive sources. Our results show a clear correlation between α-activity on the inner spectrometer surface and background from the volume of the spectrometer. Two key characteristics of the Main Spectrometer background – the dependency on the inner electrode offset potential, and the radial distribution – could be reproduced with this artificially induced background. These findings indicate a high contribution of α-decay induced events to the residual KATRIN background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Design of a new charge exchange recombination spectroscopy diagnostic for impurity transport experiments at Wendelstein 7-X

In this study, we present the design of a new charge exchange recombination spectroscopy (CXRS) system for Wendelstein 7-X (W7-X), which aims at measuring line radiation from highly ionized iron impurities after laser blow off injections over transport timescales. New fiber bundles with a core diameter of 600 µm have been added to one of the existing optical systems at W7-X used for CXRS. The fibers direct collected light to five newly developed single channel high frame rate F/2.8 spectrometers. The new custom-made spectrometers are built using entirely off-the-shelf and 3D printed components, and employ a round-to-linear fiber bundle, instead of an entrance slit, to maximize light throughput. The detector is an ultra-low readout noise EM-CCD camera capable of frame rates of up to 10 000 fps for about 500 consecutive frames when operated using a single readout channel. Such high frame rates will enable the study of inward convection of injected impurities. An initial sensitivity study is performed using a newly developed Markov chain Monte Carlo approach based on the pySTRAHL impurity transport code. This study indicates that, with the addition of the new CX diagnostic, impurity convection velocities can be inferred.

Instruments & Instrumentation↗
DOE OSTI2024

High-n Rydberg transition spectroscopy for heavy impurity transport studies in W7-X (invited)

Here, we present a novel spectroscopy approach to investigate impurity transport by analyzing line-radiation following high-n Rydberg transitions. While high-n Rydberg states of impurity ions are unlikely to be populated via impact excitation, they can be accessed by charge exchange (CX) reactions along the neutral beams in high-temperature plasmas. Hence, localized radiation of highly ionized impurities, free of passive contributions, can be observed at multiple wavelengths in the visible range. For the analysis and modeling of the observed Rydberg transitions, a technique for calculating effective emission coefficients is presented that can well reproduce the energy dependence seen in datasets available on the OPEN-ADAS database. By using the rate coefficients and comparing modeling results with the new high-n Rydberg CX measurements, impurity transport coefficients are determined with well-documented 2σ confidence intervals for the first time. This demonstrates that high-n Rydberg spectroscopy provides important constraints on the determination of impurity transport coefficients. By additionally considering Bolometer measurements, which provide constraints on the overall impurity emissivity and, therefore, impurity densities, error bars can be reduced even further.

Instruments & Instrumentation↗
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 OSTI2022

A Simple Condition for Sustained Super-Eddington Black Hole Growth

One of the most pressing questions in cosmology is how the black holes (BHs) powering quasars at high redshift grow to supermassive scales within a billion years of the Big Bang. Here we show that sustained super-Eddington accretion can be achieved for BHs with Eddington fractions f Edd ≳ 2/ϵ, where ϵ is the efficiency with which radiation is generated in the accretion process. In this regime, the radiation carries too little momentum to halt the accretion flow and the infalling gas traps the radiation. The BH growth then proceeds unimpeded until the gas supply is exhausted, in contrast to accretion at lower rates, which is limited by the radiation generated in the accretion process. The large gas supply available in massive high-redshift quasar host galaxies may be readily accreted onto seed BHs via this supply-limited mode of accretion, providing an explanation for how such supermassive BHs are assembled in the early universe. This sustained super-Eddington growth may also explain the short lifetimes inferred for the H ii regions surrounding high-redshift quasars, if the bulk of the BH growth occurs without the associated radiation escaping to ionize the intergalactic medium. It furthermore implies that a population of obscured rapidly growing BHs may be difficult to detect, perhaps explaining why so few quasars with Eddington fractions higher than a few have been observed. Finally, this simple condition for sustained super-Eddington growth can easily be implemented in cosmological simulations that can be used to assess in which environments it occurs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Revealing the Galaxy–Halo Connection through Machine Learning

Abstract Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. The complexity of these physics can prove difficult to disentangle and obfuscate how mass-dependent trends in the galaxy population originate. Here, we train a machine-learning method called Explainable Boosting Machines (EBMs) to infer how the stellar mass and star formation rate of nearly 6 million galaxies simulated by the Cosmic Reionization on Computers project depend on the physical properties of halo mass, the peak circular velocity of the galaxy during its formation history v peak , cosmic environment, and redshift. The resulting EBM models reveal the relative importance of these properties in setting galaxy stellar mass and star formation rate, with v peak providing the most dominant contribution. Environmental properties provide substantial improvements for modeling the stellar mass and star formation rate in only ≲10% of the simulated galaxies. We also provide alternative formulations of EBM models that enable low-resolution simulations, which cannot track the interior structure of dark matter halos, to predict the stellar mass and star formation rate of galaxies computed by high-resolution simulations with detailed baryonic physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

H α emission in local galaxies: star formation, time variability, and the diffuse ionized gas

ABSTRACT The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$10{{\ \rm per\ cent}}$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $\sim 40{{\ \rm per\ cent}}$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx 28{{\ \rm per\ cent}}$ and $f_{\rm He}\approx 9{{\ \rm per\ cent}}$, respectively. Together with an escape fraction of $f_{\rm esc}\approx 6{{\ \rm per\ cent}}$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx 57{{\ \rm per\ cent}}$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Advanced CMOS Reliability Update: Sub 20nm FinFET Assessment

Advances in FinFET design and fabrication enable manufacturing of denser, more compact integrated circuits (ICs) with substantially reduced leakage while shortening the channel-lengths. The same stress-induced leakage and breakdown degradation mechanisms that affect planar transistors also impact FinFET devices. Reliability concerns such as Bias Temperature Instability (BTI), Time Dependent Dielectric Breakdown (TDDB), and Hot Carrier Injection (HCI) become very important with changes to transistor geometry and fin sidewall crystal orientation. Recent testing has shown that FinFETs respond differently to radiation (radiation effects such as total ionizing dose) when compared to planar transistors. These reliability and radiation effects issues become very important when changing transistor geometry and scaling FinFETs towards smaller feature sizes (22-nm, 16-nm, 14- nm, 10-nm, and smaller critical dimensions). The comparable 2019 state of the art transistor densities in current high-volume manufacturing silicon-based foundries is 7-nm (ISMC, Samsung) and 10-nm (Intel) [www.anandtech.com,fuse.wikichip.org]. Released products include supporting components for the cellphone and commercial microprocessor markets respectively. Extensive development in the foundry industry is driving to a 5-nm technology node in late 2020.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Photonuclear Physics in SCALE

Radiation shielding modeling involves a wide range of applications, such as power reactors, critical experiments, portable sources, linear accelerators, and fusion systems. Particle energy ranges vary from thermal energies to fast energy spectra. MAVRIC, as the radiation shielding sequence of the SCALE code system [1] developed and maintained by Oak Ridge National Laboratory, is capable of modeling both neutron and photon interactions in matter and is widely used by the US Nuclear Regulatory Commission, Department of Energy, the National Nuclear Security Administration, and radiation shielding practitioners all around the world. Although most applications don’t involve highly energetic photons that can cause photonuclear reactions and create coupled chains of high energy, ionized electrons and bremsstrahlung radiation within a material medium, when present, these reactions can have a significant impact upon the downstream particle fluxes from incident photons above 20 MeV. For example, linear accelerators used for medical isotope production are examples of applications in which high-energetic photons are generated and used for driving or controlling other components of a complex facility. Radiation shielding analyses of such a facility would require accounting for any and all significant sources of radiation. Unfortunately, MAVRIC included no data and methods to model either photonuclear reactions or the bremsstrahlung photons until this work, which enabled accounting for both. Photonuclear reaction cross sections were imported from the Evaluated Nuclear Data File libraries, and corresponding methods were implemented in the SCALE code system; see Section 2 for how the nuclear data are imported and used. Bremsstrahlung photons were modeled by assuming thick-target approximation whereby all generated electrons and positrons due to photon interactions are assumed to be absorbed within the same material in which they were born, yielding a quick particle transport simulation method by using fast lookup tables for generating subsequent gammas; see Section 4 for details. A suite of verification models was used to test implemented methods and associated data. Tests yielded satisfactory results for the new features. Enhancements will be released as part of the next major SCALE release (v7.0), which is anticipated to be deployed in 2025, with beta releases available in the interim.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Effect of optical field ionisation on the generation of wake fields by femtosecond laser pulses in an inhomogeneous plasma

We consider the interaction of a high-intensity short laser pulse with an argon gas during optical field ionisation. Modelling in a three-dimensional cylindrical symmetric geometry is performed at various moderately relativistic intensities and positions of the focal plane of the laser beam in the case of both optical gas ionisation and a pre-ionised plasma. The influence of the processes occurring during optical field ionisation of the gas on the generation of wake waves is investigated, and the conditions are found under which intense wake fields are produced in a plasma formed from an inhomogeneous argon gas jet. The possibility of using a gas with a large number of electrons on the outer shell (argon) to excite intense wake waves and accelerate electrons is demonstrated. Despite significant ionisation refraction of the laser pulse on a radially inhomogeneous density profile of plasma electrons formed during optical field ionisation, the region of the parameters of the laser pulse and the gas target is determined at which ionisation refraction leads to the excitation of an intense wake wave. It is found that in a pre-ionised plasma a wake wave is not generated even at the same laser radiation intensity. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Nonlinear Thomson scattering of a tightly focused relativistically intense laser pulse by an ensemble of particles

We report a study of the process of nonlinear Thomson scattering (NTS) of a tightly focused relativistically intense laser pulse by an ensemble of electrons simulating the result of ionisation of an ultrathin nanofoil. The description of all six components of a laser pulse is based on the Stratton – Chu integrals describing radiation focused by an off-axis parabolic mirror into a spot down to the diffraction limit, which is of practical interest for modern experiments with extreme-intensity lasers, including diagnostics of their parameters. The interaction of a laser pulse with electrons is simulated by the test-particle method, and the resulting trajectories are used to calculate the characteristics of secondary electromagnetic radiation. The effect of the laser-pulse focal spot diameter in the range λ – 6λ on the angular distributions of the radiated energy and NTS spectra, as well as of a laser-beam peak intensity of 10{sup 21} – 10{sup 22} W cm{sup −2} on the spectral width of the secondary radiation is studied in relation to the experimentally implemented focusing scheme using the example of a laser pulse with an initial homogeneous spatial distribution and a duration of ∼25 fs. (interaction of laser radiation with matter. laser plasma)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Generation of a quasi-static magnetic field by a circularly polarised laser pulse due to tunnelling gas ionisation

We present a theoretical model of the quasi-static magnetic field generation in a laser channel, which is formed behind the front of a short laser pulse that ionises a gas. The generation of a magnetic field is caused by the appearance of the electron pressure anisotropy during tunnelling ionisation of atoms. In the considered case of subrelativistic laser light intensities, the generated magnetic field can reach ∼1 MG with an energy transformation ratio of about 1 %, which paves the way for identifying the proposed mechanism when use is made of a wide class of ultrashort pulse lasers. (interaction of laser radiation with matter. laser plasma)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Strategies for solving the excited-state self-consistent-field problem for highly excited and multiply ionized states

The dynamics of molecules exposed to intense x-ray radiation involve a large number of multiply ionized and highly excited electronic configurations. To model these dynamics a reliable and efficient electronic structure model is imperative. Employing the Hartree-Fock-Slater electronic structure model in combination with the maximum overlap method, we quantify the associated convergence failures when calculating electronic states of carbon monoxide with multiple vacancies in the core and valence levels. Here, we characterize these cases and describe strategies to overcome the convergence problems. The described techniques not only eliminate all convergence issues for CO but also result in a significant reduction of convergence failures for simulations of the x-ray-induced multiple ionization dynamics of the phenol molecule.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Integrated Circuit Design for Radiation-Hardened Charge-Sensitive Amplifier Survived up to 2 Mrad

According to the continuous development of metal-oxide semiconductor (MOS) fabrication technology, transistors have naturally become more radiation-tolerant through steadily decreasing gate-oxide thickness, increasing the tunneling probability between gate-oxide and channel. Unfortunately, despite this radiation-hardened property of developed transistors, the field of nuclear power plants (NPPs) requires even higher radiation hardness levels. Particularly, total ionizing dose (TID) of approximately 1 Mrad could be required for readout circuitry under severe accident conditions with 100 Mrad around a reactor in-core required. In harsh radiating environments such as NPPs, sensors such as micro-pocket-fission detectors (MPFD) would be a promising technology to be operated for detecting neutrons in reactor cores. For those sensors, readout circuits should be fundamentally placed close to sensing devices for minimizing signal interferences and white noise. Therefore, radiation hardening ability is necessary for the circuits under high radiation environments. This paper presents various integrated circuit designs for a radiation hardened charge-sensitive amplifier (CSA) by using SiGe 130 nm and Si 180 nm fabrication processes with different channel widths and transistor types of complementary metal-oxide-semiconductor (CMOS) and bipolar CMOS (BiCMOS). These circuits were tested under γ–ray environment with Cobalt-60 of high level activity: 490 kCi. The experiment results indicate amplitude degradation of 2.85%–34.3%, fall time increase of 201–1730 ns, as well as a signal-to-noise ratio (SNR) of 0.07–11.6 dB decrease with irradiation dose increase. These results can provide design guidelines for radiation hardening operational amplifiers in terms of transistor sizes and structures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Evaluation of Histone Deacetylase Inhibitors as Radiosensitizers for Proton and Light Ion Radiotherapy

Significant opportunities remain for pharmacologically enhancing the clinical effectiveness of proton and carbon ion-based radiotherapies to achieve both tumor cell radiosensitization and normal tissue radioprotection. We investigated whether pretreatment with the hydroxamate-based histone deacetylase inhibitors (HDACi) SAHA (vorinostat), M344, and PTACH impacts radiation-induced DNA double-strand break (DSB) induction and repair, cell killing, and transformation (acquisition of anchorage-independent growth in soft agar) in human normal and tumor cell lines following gamma ray and light ion irradiation. Treatment of normal NFF28 primary fibroblasts and U2OS osteosarcoma, A549 lung carcinoma, and U87MG glioma cells with 5–10 µM HDACi concentrations 18 h prior to cesium-137 gamma irradiation resulted in radiosensitization measured by clonogenic survival assays and increased levels of colocalized gamma-H2AX/53BP1 foci induction. We similarly tested these HDACi following irradiation with 200 MeV protons, 290 MeV/n carbon ions, and 350 MeV/n oxygen ions delivered in the Bragg plateau region. Unlike uniform gamma ray radiosensitization, effects of HDACi pretreatment were unexpectedly cell type and ion species-dependent with C-12 and O-16 ion irradiations showing enhanced G0/G1-phase fibroblast survival (radioprotection) and in some cases reduced or absent tumor cell radiosensitization. DSB-associated foci levels were similar for proton-irradiated DMSO control and SAHA-treated fibroblast cultures, while lower levels of induced foci were observed in SAHA-pretreated C-12 ion-irradiated fibroblasts. Fibroblast transformation frequencies measured for all radiation types were generally LET-dependent and lowest following proton irradiation; however, both gamma and proton exposures showed hyperlinear transformation induction at low doses (≤25 cGy). HDACi pretreatments led to overall lower transformation frequencies at low doses for all radiation types except O-16 ions but generally led to higher transformation frequencies at higher doses (>50 cGy). The results of these in vitro studies cast doubt on the clinical efficacy of using HDACi as radiosensitizers for light ion-based hadron radiotherapy given the mixed results on their radiosensitization effectiveness and related possibility of increased second cancer induction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Development, construction and tests of the Mu2e electromagnetic calorimeter mechanical structures

The “muon-to-electron conversion” (Mu2e) experiment at Fermilab will search for the charged lepton flavour violating neutrino-less coherent conversion of a muon into an electron in the field of an aluminum nucleus. The observation of this process would be the unambiguous evidence of the existence of physics beyond the standard model. Mu2e detectors comprise a straw-tracker, an electromagnetic calorimeter and an external veto for cosmic rays. In particular, the calorimeter provides excellent electron identification, a fast calorimetric online trigger, and complementary information to aid pattern recognition and track reconstruction. The detector has been designed as a state-of-the-art crystal calorimeter and employs 1348 pure Cesium Iodide (CsI) crystals readout by UV-extended silicon photosensors and fast front-end and digitization electronics. A design consisting of two identical annular matrices (named “disks”) positioned at the relative distance of 70 cm downstream the aluminum target along the muon beamline satisfies the Mu2e physics requirements. The hostile Mu2e operational conditions, in terms of radiation levels (total expected ionizing dose of 12 krad and a neutron fluence of 5 × 10$^{10}$ n/cm$^{2}$ @ 1 MeV$_{eq}$ (Si)/y), magnetic field intensity (1 T) and vacuum level (10$^{-4}$ Torr) have posed tight constraints on scintillating materials, sensors, electronics and on the design of the detector mechanical structures and material choice. The support structure of each 674 crystal matrix is composed of an aluminum hollow ring and parts made of open-cell vacuum-compatible carbon fiber. The photosensors and front-end electronics for the readout of each crystal are inserted in a machined copper holder and make a unique mechanical unit. The resulting 674 mechanical units are supported by a machined plate of vacuum-compatible plastic material. The plate also integrates the cooling system made of a network of copper lines flowing a low temperature radiation-hard fluid and placed in thermal contact with the copper holders to constitute a low resistance thermal bridge. The data acquisition electronics are hosted in aluminum custom crates positioned on the external lateral surface of the disks. The crates also integrate the electronics cooling system as lines running in parallel to the front-end system. In this paper we report on the calorimeter mechanical structure design, the mechanical and thermal simulations that have determined the design technological choices, and the status of component production, quality assurance tests and plans for assembly at Fermilab.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Ultra-short pulse generation from mid-IR to THz range using plasma wakes and relativistic ionization fronts

In this work, we discuss numerical and experimental results on frequency downshifting and upshifting of a 10 μ m infrared (IR) laser to cover the entire wavelength (frequency) range from λ = 1 to 150 μ m ( ν = 300–2 THz) using two different plasma techniques. The first plasma technique utilizes frequency downshifting of the drive laser pulse in a nonlinear plasma wake. Based on this technique, we have proposed and demonstrated that in a tailored plasma structure, multi-millijoule energy, single-cycle, long-wavelength IR (3–20 μ m) pulses can be generated by using an 810 nm Ti:sapphire drive laser. Here, we extend this idea to the THz frequency regime. We show that sub-joule, terawatts, single-cycle terahertz (2–12 THz or 150–25 μ m) pulses can be generated by replacing the drive laser with a picosecond 10 μ m CO 2 laser and a different shaped plasma structure. The second plasma technique employs frequency upshifting by colliding a CO 2 laser with a rather sharp relativistic ionization front created by ionization of a gas in less than half cycle (17 fs) of the CO 2 laser. Even though the electrons in the ionization front carry no energy, the frequency of the CO 2 laser can be upshifted due to the relativistic Doppler effect as the CO 2 laser pulse enters the front. The wavelength can be tuned from 1 to 10 μ m by simply changing the electron density of the front. While the upshifted light with 5 < λ ( μ m ) < 10 propagates in the forward direction, that with 1 < λ ( μ m ) < 5 is back-reflected. These two plasma techniques seem extremely promising for covering the entire molecular fingerprint region.

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↗