Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA DECAY”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗
DOE OSTI2024

Characterization of the image plate multi-scan response to mono-energetic x-rays

Image plates (IPs), or phosphor storage screens, are a technology employed frequently in inertial confinement fusion (ICF) and high energy density plasma (HEDP) diagnostics because of their sensitivity to many types of radiation, including, x rays, protons, alphas, beta particles, and neutrons. Prior studies characterizing IPs are predicated on the signal level remaining below the scanner saturation threshold. Since the scanning process removes some signal from the IP via photostimulated luminescence, repeatedly scanning an IP can bring the signal level below the scanner saturation threshold. This process, in turn, raises concerns about the signal response of IPs after an arbitrary number of scans and whether such a process yields, for example, a constant ratio of signal between the nth and n + 1st scan. Here, the sensitivity of IPs is investigated when scanned multiple times. It is demonstrated that the ratio of signal decay is not a constant with the number of scans and that the signal decay depends on the x-ray energy. As such, repeatedly scanning an IP with a mixture of signal types (e.g., x ray, neutron, and protons) enables ICF and HEDP diagnostics employing IPs to better isolate a particular signal type.

Instruments & Instrumentation↗
DOE OSTI2026

Post-puff SOL broadening on MAST-U under high-recycling conditions: evidence consistent with cross-field transport changes

Transient broadening of the scrape-off layer (SOL) density profile can modify main-chamber first-wall particle fluxes and divertor loading, yet its control parameters remain debated between divertor-regime transitions, neutral dynamics and changes in cross-field transport. We investigate fueling-driven SOL density-profile evolution and post-fueling relaxation on MAST-Upgrade (MAST-U) in ohmic L-mode, using two otherwise similar double-null Conventional Divertor discharges (I p = 450 kA, B T = 0.33 T) with identical 50 ms low-field side gas puffs; the only intentional difference is the puff start time. Upstream Thomson scattering shows that both discharges develop a transient far-SOL density profile modification, expressed as an increased SOL-width metric λ n e and a far-SOL enhancement consistent with a shoulder-like signature. In the earlier-puff case, the SOL-width metric remains elevated after puff termination and the post-puff decay time scales are systematically longer across the analysed radii. Outer-target Langmuir probes indicate high-recycling conditions during the analysed window (few-eV T e with radially peaked j sat,∥ and q ∥ profiles, without signatures of deep detachment). The outer-divertor collisionality proxy Λ div is elevated in both cases and does not discriminate between the different post-puff persistence. The target-integrated ion flux proxy ∫ J sat dA evolves nearly identically in both discharges when aligned to the puff start. Taken together, these observations suggest that the late/post-puff upstream SOL evolution is not set by parallel exhaust to the outer targetalone and is more consistent with upstream cross-field redistribution, with a possible role for plasma–neutral coupling and fueling geometry.

MAST-U↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: II. Time evolution

Dedicated DIII-D experiments studied instabilities below the ion cyclotron frequency f ci and at cyclotron harmonics by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H and D at five different values of toroidal field. The resulting database of 2529 conditions is described in Part 1: Stability. This paper focuses on the temporal evolution. For some beam geometries, the ion cyclotron emission (ICE) is much stronger at the beginning of injection, suggesting the importance of ‘bump-on-tail’ energy drive. Sawtooth crashes cause the emission to cease for some beam types and emission bands, but the opposite effect occurs for others. When unstable, for virtually all conditions and frequency bands, the emission bursts in a cycle reminiscent of predator-prey oscillations. Evidence for nonlinear interaction between ICE harmonics is observed. Strongly unstable modes typically appear within 1 ms from the start of beam injection and often decay within 0.1 ms when the beam turns off, but both the initial rise and final decay of the mode power exhibit a wide variety of behavior.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2025

Relaxation times for disoriented isospin condensates in high energy heavy ion collisions

Fluctuations between charged and neutral kaons measured by the ALICE Collaboration in Pb-Pb collisions at the CERN Large Hadron Collider (LHC) exceed conventional explanations. Previously it was shown that if the scalar condensate is accompanied by an electrically neutral isospin-1 field then the combination can produce large equilibrium fluctuations where ⟨$\overline{𝑢}$ ⁢𝑢⟩≠⟨ $\overline{𝑑}$𝑑⟩. Hadronizing strange and antistrange quarks might then strongly fluctuate between charged ($\overline{𝑢}$⁢𝑠 or $\overline{𝑠}$⁢𝑢) and neutral (𝑑⁢$\overline{𝑠}$⁢ or 𝑠⁢$\overline{𝑑}$) kaons. Here, we estimate the times for the condensates to achieve their equilibrium probability distributions within causal volumes in high energy heavy ion collisions. This is achieved by modeling the temperature dependence of the condensates, mesonic collective excitations, decay rates of the associated fields, and employing the Langevin and Fokker-Planck equations. Within this model, we find that the equilibration times are short compared with the expansion time, making disoriented isospin condensates a viable explanation for the anomalous fluctuations observed at the LHC.

Quark-gluon plasma↗
DOE OSTI2024

Inter-ELM pedestal turbulence dynamics dependence on q 95 and temperature gradient

A series of dedicated experiments from DIII-D tokamak provide spatially and temporally resolved measurements of electron density and temperature, and multiscale and multichannel fluctuations over a wide range of conditions. Measurements of long wavelength density fluctuations in the type-I ELMing H-mode pedestal routinely revealed a coexistence of multiple instabilities that exhibit dramatic different dynamic behaviors as q 95 and temperature gradient are varied, apparently responsible for limiting pedestal temperature profiles. Two distinct frequency bands of density fluctuations are modulated with ELM cycle with frequency above 200kHz propagating in the electron diamagnetic direction in the lab frame (electron mode) and below 200kHz propagating in the ion diamagnetic direction (ion mode). The electron mode amplitude peaks near the electron temperature gradient region and increases with q 95 which seems to be correlated with the increased χ e at higher q 95 , similar to the characteristics expected for Micro-tearing Mode (MTM). At higher q 95 , during the inter-ELM period, the ion mode decays at later phase of the ELM cycle. Consistently, the poloidal correlation length of the ion mode is also found to reduce which suggests the possible E×B flow shear suppression of ion mode at later phase of the ELM cycle as the Er well recovers. In contrast, the electron mode grows during the ELM cycle and reaches saturation at around 50-60% of ELM period. Linear gyrokinetic simulations find the MTMs as the most unstable mode in the pedestal electron temperature gradient region. The higher q 95 and lower magnetic shear destabilize MTMs. These observations provide key insights of the underlying physics of multifield properties and rich dataset of experimental 'fingerprints' that enable new tests of the theoretical pedestal models and lead to developing a predictive model for pedestal formation on ITER and future burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Constraints on dark photon dark matter from Lyman- α forest simulations and an ultrahigh signal-to-noise quasar spectrum

The ultralight dark photon is a well-motivated, hypothetical dark matter candidate. In a dilute plasma, they can resonantly convert into photons, and heat up the intergalactic medium between galaxies. In this work, we explore the dark photon dark matter parameter space by comparing synthetic Lyman- α forest data from cosmological hydrodynamical simulations to observational data from VLT/UVES of the quasar HE0940-1050 ( z em = 3.09 ). We use a novel flux normalization technique that targets underdense gas, reshaping the flux probability distribution. Not only do we place robust constraints on the kinetic mixing parameter of dark photon dark matter, but notably our findings suggest that this model can still reconcile simulated and observed Doppler parameter distributions of z ∼ 0 Lyman- α lines, as seen by HST/COS. This work opens new pathways for the use of the Lyman- α forest to explore new physics, and can be extended to other scenarios such as primordial black hole evaporation, dark matter decay, and annihilation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Particle signal considerations for isotope ratio analysis with single particle multi-collector inductively coupled plasma mass spectrometry

Particle analysis has benefitted from the advent of single particle inductively coupled plasma-mass spectrometry (spICP-MS) due to its robustness, sensitivity, and high-throughput nature. Previous methods of spICP-MS have typically utilized quadrupole or time-of-flight mass analyzers and therefore employ electron multiplier-based detectors (such as secondary electron multipliers or microchannel plates). However, to obtain precise measurements on elemental or isotopic ratios within individual particles, multi-collector ICP-MS (MC-ICP-MS) can be used. Here, we investigate Ce isotope ratios, specifically 142 Ce/ 140 Ce, by spMC-ICP-MS using an all-Faraday cup collector array. Using 1 μm (diameter) cerium dioxide particles, integration times of the Faraday cup detectors were varied from 50–500 ms. The signal from the cerium isotopes in the particles was used to determine isotope ratios, which closely matched the expected natural isotopic abundances. Due to the signal decay response from the Faraday cups, the signal from particles lasts much longer than the expected 1–2 ms (up to 100 s of ms). To explore this effect on isotope ratio analysis, multiple ratio analysis methods were used to determine how to obtain optimal precision and accuracy. Relative differences were around 2% for methods that calculated isotope ratios from summing the total signal of an individual particle before calculating the ratio (rather than using every data point individually). It was found that summing all data points per particle, or integrating under the signal peak, yielded both accurate and precise isotope ratios within the particle population. Particles were also sampled off a solid substrate via microextraction, and isotope ratios were determined with relative differences of 0.13% to 9%. This demonstrates the ability to use spMC-ICP-MS to obtain isotope ratios on particles, with little to no relative difference in comparison to the expected ratio, even when operating Faraday detectors at fast 50 ms integration times.

Szakas, Sarah E. [Oak Ridge National Laboratory (O↗
DOE OSTI2024

Evolution of highly multimodal Rayleigh–Taylor instabilities

Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Approach to startup inventory for viable commercial fusion power plant

With the increasing efforts to commercialize fusion power, private and government organizations are investing heavily in the development of technology to support a viable fusion power plant. Deuterium-Tritium (DT) fueled reactors are more prevalent than other proposed designs, requiring tritium processing and handling technology for safe operations and self-sufficiency. Further, each fusion power plant will need a specific-to-design startup inventory of tritium to begin operations. This startup inventory is required prior to breeding and is the minimum tritium inventory required to fill each processing component in the fuel cycle, to offset radioactive decay losses, and to avoid a zero-fuel situation for continuous operation. We present an approach to calculate the startup tritium inventory for a 500 MW th reactor, with considerations for reserve inventory for maintenance and commissioning. A baseline startup inventory was calculated to be approximately 327 gs. This value was obtained using modest assumptions about the technology and operating parameters of a fusion power plant. The required operating reserve inventory or the inventory necessary to keep a fusion power plant operational using only direct internal recycling for 24 h for the same plant design is approximately 642 gs. The approach and findings of this paper will enable fusion energy stakeholders to better utilize the existing scarce global tritium supply.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Establishing Pb-203 production from electrodeposited Tl targets at Brookhaven National Laboratory

Background: Promising developments in Pb-212 radiopharmaceutical therapies have increased demand for Pb-203 diagnostic agents. Building on previous work from various isotope production facilities, this study optimized Pb-203 production from electrodeposited Tl targets at Brookhaven National Laboratory (BNL). The additional supply of Pb-203 may help meet growing preclinical and clinical demands. Results: Two Tl targets were irradiated at the Brookhaven Linac Isotope Producer facility with 30 ± 1 MeV protons, measured using previously published cross section data. Distribution coefficients for Pb Resin in acetate media were investigated for both Na + and K + cations, where potassium acetate was ~ 4 times more effective at stripping Pb from the Pb Resin. The Tl electrodeposition was optimized to deposit 350 mg of Tl (~ 60 mg/cm 2 ) on Au backing in under 6 h. The proposed separation process was completed in < 1.5 h and achieved > 98% and 92 ± 3% recovery of Tl and Pb, respectively, with an overall Tl-Pb separation factor of 6 × 10 5 . The experimentally measured half-life of Pb-203 was 52.4 ± 0.7 h, agreeing with 51.93 ± 0.02 h reported by the National Nuclear Data Center. The radioisotopic purity of the Pb fraction at 24 h post end of bombardment (EOB) from a 24 h irradiation was 66% Pb-203, 28% Pb-201, and 6% Pb-200. Following chemical separation, the Pb-203 produced in this work (21 MBq Pb-203 EOB) achieved apparent molar activities of 10 ± 5 and 0.9 ± 0.5 GBq/µmol for [ 203 Pb]Pb-DOTAM and [ 203 Pb]Pb-DO3A, respectively, decay corrected to EOB. Data derived from this work suggests BNL can produce > 10’s GBq Pb-203 with > 99% radiochemical and radioisotopic purity from Tl-205 for worldwide distribution. Conclusions: The production and separation of Pb-203 from natural Tl target material was successfully demonstrated at BNL. Existing methods were adapted and optimized for the facilities at BNL. Results from this work will guide future large-scale Pb-203 production opportunities at BNL for clinical applications.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Bottom quark energy loss and hadronization with B + and ${\textrm{B}}_{\textrm{s}}^0$ nuclear modification factors using pp and PbPb collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV

The production cross sections of $ {\textrm{B}}_{\textrm{s}}^0 $ and B$^{+}$ mesons are reported in proton-proton (pp) collisions recorded by the CMS experiment at the CERN LHC with a center-of-mass energy of 5.02 TeV. The data sample corresponds to an integrated luminosity of 302 pb$^{−1}$. The cross sections are based on measurements of the $ {\textrm{B}}_{\textrm{s}}^0 $→ J/ψ(μ$^{+}$μ$^{−}$)ϕ(1020)(K$^{+}$K$^{−}$) and B$^{+}$→ J/ψ(μ$^{+}$μ$^{−}$)K$^{+}$ decay channels. Results are presented in the transverse momentum (p$_{T}$) range 7–50 GeV/c and the rapidity interval |y| < 2.4 for the B mesons. The measured p$_{T}$-differential cross sections of B$^{+}$ and $ {\textrm{B}}_{\textrm{s}}^0 $ in pp collisions are well described by fixed-order plus next-to-leading logarithm perturbative quantum chromodynamics calculations. Using previous PbPb collision measurements at the same nucleon-nucleon center-of-mass energy, the nuclear modification factors, R$_{AA}$, of the B mesons are determined. For p$_{T}$> 10 GeV/c, both mesons are found to be suppressed in PbPb collisions (with R$_{AA}$ values significantly below unity), with less suppression observed for the $ {\textrm{B}}_{\textrm{s}}^0 $ mesons. In this p$_{T}$ range, the R$_{AA}$ values for the B$^{+}$ mesons are consistent with those for inclusive charged hadrons and D$^{0}$ mesons. Below 10 GeV/c, both B$^{+}$ and $ {\textrm{B}}_{\textrm{s}}^0 $ are found to be less suppressed than either inclusive charged hadrons or D$^{0}$ mesons, with the $ {\textrm{B}}_{\textrm{s}}^0 $R$_{AA}$ value consistent with unity. The R$_{AA}$ values found for the B$^{+}$ and $ {\textrm{B}}_{\textrm{s}}^0 $ are compared to theoretical calculations, providing constraints on the mechanism of bottom quark energy loss and hadronization in the quark-gluon plasma, the hot and dense matter created in ultrarelativistic heavy ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Investigating the role of edge neutrals in exciting tearing mode activity and achieving flat temperature profiles in LTX-β

Abstract We present observations, numerical simulations, and analysis from experiments in the Lithium Tokamak Experiment-Beta (LTX- β ) in which the electron temperature profile ( T e ( r )) shifts from flat to peaked and a tearing mode is also destabilized when the average density ( n e ave ) exceeds ∼10 19 m −3 . Flat T e ( r ) is obtained routinely in LTX- β , with a lithium coated, low-recycling first wall, once the external fueling is stopped and density decays [Boyle et al 2023 Nucl. Fusion 63 056020]. In the present experiment, flat T e profiles can be sustained while maintaining constant n e ave below a line averaged density threshold ( n e ave th ) of ∼10 19 m −3 . Above n e ave th , T e ( r ) shifts from flat to peaked and a tearing mode is destabilized. Due to low recycling, the achieved n e ave can be controlled precisely by external fueling and hence, a certain threshold of the edge neutral inventory from the external fueling is experimentally manifested through n e ave th . The goal of the present work is to investigate the role of edge neutrals in determining T e ( r ) and MHD stability in the unique low-recycling regime of LTX- β . Our hypothesis is that the peaking of T e ( r ) beyond n e ave th is due ultimately to the edge cooling by the cold neutrals beyond a critical fueling flux. At lower fueling flux, flat T e ( r ) results in broader pressure profile and lower resistivity, which in turn stabilizes the tearing mode. This hypothesis is supported by edge neutral density estimation by DEGAS 2 code. Mode analysis by singular value decomposition confirms the tearing mode structure to be m / n = 2/1 ( m and n being the poloidal and toroidal mode numbers). Linear tearing stability analysis with M3D-C1 predicts that plasmas with n e ave > 10 19 are highly susceptible to a n = 1 tearing mode. ORBIT simulations, however, confirmed that the tearing modes do not contribute to the loss of fast ions from neutral beam injection. This study shows for the first time that the neutral inventory at the edge could be one of the deciding factors for the achievability of the unique operation regime of flat T e ( r ) and the excitation of tearing activity that could be disruptive for the plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Experimental observations of bifurcated power decay lengths in the near Scrape-Off Layer of ST40 High Field Spherical Tokamak

The scrape-off layer parallel heat flux decay lengths measured at ST40, a high field, low aspect ratio spherical tokamak, have been observed to bifurcate into two groups. The wide group follows established H-mode scalings (ranging between 2 to 8 mm) while the narrow group falls up to 10 times below these scalings (between 0.2 and 0.8 mm), being comparable to the ion total Larmor radius rather than the ion poloidal Larmor radius. The heat flux profiles of the latter group can only be described by a multi-exponential function, rather than the single exponential function convoluted with a Gaussian. The onset of the narrow scrape-off layer width is observed to be associated with suppressed magnetic fluctuations, suggesting reduced electromagnetic turbulence levels in the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development and validation of fully open-source R2S shutdown dose rate capabilities in OpenMC *

We present the first fully open-source capabilities for shutdown dose rate (SDR) calculations of fusion energy facilities based on the Rigorous 2-Step (R2S) methodology. These capabilities have been implemented in the OpenMC Monte Carlo particle transport code, building on its existing capabilities while also leveraging new features that have been added to the code to support SDR calculations, such as decay photon source generation. Each of the individual physics components in the R2S workflow—neutron transport, activation, decay photon source generation, and photon transport—have been verified through code-to-code comparisons with MCNP6.2 and FISPACT-II 4.0. These comparisons generally demonstrate excellent agreement between codes for each of the physics components. The full cell-based R2S workflow was validated by performing a simulation of the first experimental campaign from the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For short cooling times, the dose calculated by OpenMC agrees with the experimental measurements within the stated experimental uncertainties. For longer cooling times, an overprediction of the shutdown dose was observed relative to experiment, which is consistent with previous studies in the literature. Altogether, these features constitute a combination of capabilities in a single, open-source codebase to provide the fusion community with a readily-accessible option for SDR calculations and a platform for rapidly analyzing the performance of fusion technology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Diagnostic development and needs for laser driven MeV x-ray radiography

Laser-driven MeV x-ray radiography of dynamic, dense objects demands a small, high flux source of energetic x-rays to generate an image with sufficient quality. Understanding the multi-MeV x-ray spectrum underscores the ability to extrapolate from the current laser sources to new future lasers that might deploy this radiography modality. Here, in this study, we present a small study of the existing x-ray diagnostics and techniques. We also present work from National Ignition Facility-Advanced Radiographic Capability, where we deploy three diagnostics to measure the x-ray spectrum up to 30 MeV. Finally, we also discuss the needs and developments of two new diagnostics: a single crystal scintillator spectrometer and a fast decay activation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Comparison of reduced model predictions for divertor detachment onset and reattachment timescales in ASDEX Upgrade and JET experiments

Building on prior analysis of ASDEX Upgrade (AUG) experiments (Henderson et al 2023 Nucl. Fusion63 086024), this study compares simple analytical formula predictions for divertor detachment onset and reattachment timescales in JET experiments. Detachment onset primarily scales with divertor neutral pressure, impurity concentration, power directed to the targets, machine size, and integral perpendicular power decay length. JET experiments, focusing on seeding mixtures of Ne and Ar, align with the detachment onset predictions. Radiation efficiencies among the impurities show good agreement with the model predictions, contrasting with AUG observations which suggested higher efficiency for Ar and lower efficiency for Ne. The time taken to re-ionise the neutral volume in front of the outer target in fully detached divertor conditions was measured following both abrupt increases in injected neutral beam power and, separately, cutting of the impurity gas flow. Re-ionisation of the neutrals occurs within approximately 1 s on JET, which aligns with the simple model prediction derived from AUG data. While the AUG results are not new, their comparison with the JET results enhances understanding, reinforcing confidence in using simple models to predict future reactor scenarios.

JET↗
DOE OSTI2024

Corrigendum: Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation (2024 Nucl. Fusion 64 056019)

In the original article, we inadvertently omitted a reference which previously introduced the concept of tungsten injection for 'Disruption Mitigation in Tokamak Reactor via Reducing the Seed Electrons of Avalanche.' In that reference, the authors have proposed injection of a tungsten cylinder of 80 mm by 8 mm dimensions, by rail guns, to absorb the runaway seed population. This can be contrasted with another interesting idea of using small tungsten pellets coated with a low Z material for depleting the runaway seeds as an option for upgrade of the ITER Disruption Mitigation System. Interested readers are referred to those two papers on the tungsten injection for runaway seed removal/reduction, and the related issue of not shortening the current quench excessively. The application focus of was not on reducing the runaway seeds in the initial Ohmic-to-runaway current conversion phase, but on safe termination of a fully formed and likely decaying runaway beam that is about to scrape off against the first wall, which can be accelerated by vertical displacement events. This scheme thus serves as a last line of defense against potential wall damage. The choice here is locally released tungsten particulates, similar to the previous idea of using tungsten particulates as a dust shield for the divertor. The most significant finding of is that strong pitch angle scattering, in addition to energy attenuation and absorption, of high-energy runaways by tungsten particulates, provides a transport mechanism by which runaways damage of the first wall can be mitigated by reducing and spreading the runaway wall load. We regret the omission in the original article, and thank the authors of the cited article for bringing this matter to our attention. We also find the possibility of using the injected tungsten rod from to terminate the runaways, as opposed to the cloud of tungsten particulates in intriguing. In that context, it is of interest to mention another possibility of a retractable tungsten metal arm that can be swung out for deployment in runaway termination, which would remove the need for post-mitigation recovery of tungsten debris in injection schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗