Search NASA⌕ Search

SEARCH · Search NASA

Results for “collisionality”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Compressional Alfvén eigenmodes excited by runaway electrons

Compressional Alfvén eigenmodes (CAEs) driven by energetic ions have been observed in magnetic fusion experiments. In this paper, we show that the modes can also be driven by runaway electrons formed in post-disruption plasma, which may explain kinetic instabilities observed in DIII-D disruption experiments with massive gas injection. The spatial structure is calculated, as are the frequencies which are in agreement with experimental observations. Using a runaway electron distribution function obtained from a kinetic simulation, the mode growth rates are calculated and found to exceed the collisional damping rate when the runaway electron density exceeds a threshold value. The excitation of CAEs poses a new possible approach to mitigate seed runaway electrons during the current quench and surpassing the avalanche.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bolometer tomography on Wendelstein 7-X for study of radiation asymmetry

The algorithm for bolometer tomography at Wendelstein 7-X (W7-X) has been recently improved using a novel regularization functional, based on relative gradient smoothing (RGS) of the sought radiation profile. It has been validated using radiation patterns provided by 3D modeling under real plasma conditions as phantoms and then applied to bolometer measurements performed during the first divertor operation phase of W7-X. The following results are presented: 1) edge-localized 2D radiation patterns with clearly resolved magnetic island radiation structures, 2) an up-down asymmetry in the impurity radiation that is not captured by the 3D edge plasma transport modeling, 3) reversal of the asymmetry with reversed magnetic field direction. Further analysis reveals a poloidal variation of the emissivity in the outer confined plasma region with a field-direction dependent asymmetry, also supported by the soft X-ray measurements. This asymmetry is considered to be related to asymmetric impurity distributions, driven by the pronounced ion-impurity friction force at the plasma edge where the collisionality of the W7-X plasma is sufficiently high to develop impurity asymmetry as predicted by neoclassical theory of parallel impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effects of zonal flows on ion temperature gradient instability in the scrape-off layer of a field-reversed configuration

Gyrokinetic simulations of long wavelength ion temperature gradient (ITG) turbulence in the scrape-off layer (SOL) of a field-reversed configuration (FRC) find that zonal flows are nonlinearly generated and are the dominant mechanism for the nonlinear saturation of the ITG instability. After the ITG saturation, zonal flows remain undamped and gradually suppress the turbulent transport to a very low level. In the simulations with collisions, collisional damping gradually reduces zonal flow amplitude to a lower level, which allows finite ITG turbulence intensity and ion heat transport in the SOL. The steady state turbulence intensity and ion heat transport are found to be proportional to the collision frequency. Furthermore, this favorable scaling suggests that minimizing collisions (e.g. increasing temperature, reducing impurity content, etc) and preserving toroidal symmetry could improve plasma confinement in the FRC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Identifying the microtearing modes in the pedestal of DIII-D H-modes using gyrokinetic simulations

Recent evidence points toward the microtearing mode (MTM) as an important fluctuation in the H-mode pedestal for anomalous electron heat transport. A study of the instabilities in the pedestal region carried out using gyrokinetic simulations to model an ELMy H-mode DIII-D discharge (USN configuration, 1.4 MA plasma current, and 3 MW heating power) is presented. The simulations produce MTMs, identified by predominantly electromagnetic heat flux, small particle flux, and a substantial degree of tearing parity. The magnetic spectrogram from Mirnov coils exhibits three distinct frequency bands—two narrow bands at lower frequency (~35–55 kHz and ~70–105 kHz) and a broader band at higher frequency (~300–500 kHz). Global linear GENE simulations produce MTMs that are centered at the peak of the ω* profile and correspond closely with the bands in the spectrogram. The three distinctive frequency bands can be understood from the basic physical mechanisms underlying the instabilities. For example (i) instability of certain toroidal mode numbers (n) is controlled by the alignment of their rational surfaces with the peak in the ω* profile, and (ii) MTM instabilities in the lower n bands are the conventional collisional slab MTM, whereas the higher n band depends on curvature drive. While many features of the modes can be captured with the local approximation, a global treatment is necessary to quantitatively reproduce the detailed band gaps of the low-n fluctuations. Notably, the transport signatures of the MTM are consistent with careful edge modeling by SOLPS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun.228 69–82). Here the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χ c and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λ q . From transport simulations, we obtain the critical thermal diffusivity χ c = 0.5 m 2 s –1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current I p = 5 MA, and χ c = 0.3 m 2 s –1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current I p = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity ${\chi }_{\text{c}}\propto {A}^{1/2}/(Z{\left(1+Z\right)}^{\frac{1}{2}}{B}_{\text{p}}^{2})$ using ITER scenarios with fixed safety factor q 95 , major radius R, aspect ratio R/a, and the separatrix temperature T sep , and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q 95 , R, R/a, and T sep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χ c , possibly yielding a transition from turbulence- to drift-dominant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

I-mode pedestal relaxation events in the Alcator C-Mod and ASDEX Upgrade tokamaks

Abstract In some conditions, I-mode plasmas can feature pedestal relaxation events (PREs) that transiently enhance the energy reaching the divertor target plates. To shed light on their appearance, characteristics and energy reaching the divertor targets, a comparative study between two tokamaks — Alcator C-Mod and ASDEX Upgrade — is carried out. It is found that PREs appear only in a subset of I-mode discharges, mainly when the plasma is close to the H-mode transition. Also, a growing oscillating precursor before the PRE onset is observed in the region close to the separatrix in both devices, and a discussion on a possible triggering mechanism is outlined. The PRE relative energy loss from the confined region is found to increase with decreasing pedestal top collisionality ν ped * . Similarly, also the relative electron temperature drop at the pedestal top, which is related to the conductive energy loss, rises with decreasing ν ped * . Based on these relations, the PRE relative energy loss in future devices such as DEMO and ARC is estimated. Finally, the divertor peak energy fluence due to the PRE is measured on each device. Those values are then compared to the model introduced in Eich et al (2017 Nucl. Mater. Energy 12 84–90) for type-I edge localized modes. The model is shown to provide an upper boundary for PRE energy fluence data, while a lower boundary is found by dividing the model by three. These two boundaries are used to make projections of the PRE divertor energy fluence to DEMO and ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental study on the role of the target electron temperature as a key parameter linking recycling to plasma performance in JET-ILW*

Changes in global and edge plasma parameters (H98(y,2), dimensionless collisionality ν*, core density peaking, separatrix density ne,sep) with variations in the D2 fueling rate and divertor configuration are unified into a single trend when mapped to $\langle$Te,ot$\rangle$, the spatially averaged spectroscopically derived outer target electron temperature. Dedicated JET with the ITER-like wall (JET-ILW) experiments in combination with an extended JET-ILW database of unseeded low-triangularity H-mode plasmas spanning a wide range of D2 fueling rates, Ip, Bt and heating power have demonstrated the importance of $\langle$Te,ot$\rangle$ as a key physics parameter linking the recycling particle source and detachment with plasma performance. The remarkably robust H98(y,2) trend with $\langle$Te,ot$\rangle$ is connected to a strong inverse correlation between $\langle$Te,ot$\rangle$, ne,sep and ν*, thus directly linking changes in the divertor recycling moderated by $\langle$Te,ot$\rangle$ with the previously established relationship between ν*, core density peaking and core pressure resulting in a degradation in core plasma performance with decreasing $\langle$Te,ot$\rangle$ (increasing ν*). Additionally, a strong inverse correlation between the separatrix to pedestal density ratio, ne,sep/ne,ped, and $\langle$Te,ot$\rangle$ is also established, with the rise in ne,sep/ne,ped saturating at $\langle$Te,ot$\rangle$ > 10 eV. A strong reduction in H98(y,2) is observed as $\langle$Te,ot$\rangle$ is driven from 30 to 10 eV via additional D2 gas fueling, while the divertor remains attached. Consequently, the pronounced performance degradation in attached divertor conditions has implications for impurity seeding radiative divertor scenarios, in which H98(y,2) is already low (~0.7) before impurities are injected into the plasma since moderate gas fueling rates are required to promote high divertor neutral pressure. A favorable pedestal pressure, pe,ped, dependence on Ip has also been observed, with an overall increase in pe,ped at Ip = 3.4 MA as $\langle$Te,ot$\rangle$ is driven down from attached to high-recycling divertor conditions. In contrast, pe,ped is reduced with decreasing $\langle$Te,ot$\rangle$ in the lower Ip branches. Further work is needed to (i) clarify the potential role of edge opacity on the observed favorable pedestal pressure Ip scaling; as well as to (ii) project the global and edge plasma performance trends with $\langle$Te,ot$\rangle$ to reactor-scale devices to improve predictive capability of the coupling between recycling and confined plasma fueling in what are foreseen to be more opaque edge plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of deuterium and carbon radiation transport in MAST-U tokamak advanced divertors

Modest effects of deuterium and carbon radiation opacity in the Super-X and snowflake divertor plasmas are predicted for MAST Upgrade tokamak with core plasma input power 2.5–5 MW and plasma current 1 MA. The radiation transport modeling is based on the SOLPS-EIRENE and UEDGE code divertor plasma predictions. Two radiation transport models are used: one is based on a full radiation transport equation implemented in the radiation transport and collisional-radiative code CRETIN (without feedback on the background plasma), and another is an internal self-consistent UEDGE model with ionization, recombination, and heating rates corrected for Ly α line trapping based on the escape probability model implemented in CRETIN. In MAST-U, the Super-X and snowflake divertor plasmas are predicted to reach detached regimes at lower upstream densities than the standard divertor, and the conclusion still holds with radiation transport effects included. At neutral densities m -3 , modest Ly α deuterium line trapping with optical depths 10–15 is predicted in the Super-X divertor. Divertor plasmas are optically thin to other Lyman and Balmer lines, as well as to strong C III and C IV lines that are responsible for most of divertor radiated power. Insignificant changes (within a few percent) to divertor deuterium ionization and recombination rates are found. Radiation fluxes on outer divertor target are modified within a factor of 2–3 when the radiation transport is accounted for, and a similar variation is found due to the line shape models that define the absorption and emission line profiles in the radiation transport modeling. The predicted Lyman and Balmer spectral intensities are significantly modified due to radiation trapping. A measurement of divertor radiation transport effects is discussed using the Ly β /Ba α line ratio. In the snowflake divertor configuration, divertor plasmas are found to be optically thin to Lyman series lines within a large range of parameter variations that include upstream density, divertor transport coefficients, and magnetic configurations. Finally, modest radiation transport effects are only found in a few cases with strongest divertor transport and magnetic configurations closest to the ideal snowflake configuration, however, the plasma background models that were used are yet to be validated with an experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hot-tail electrons’ impact on assimilation and injection penetration of D2 shattered pellet injections

The assimilation and the penetration of injected materials play significant roles in the mitigation efficiency of shattered pellet injection (SPI) as a disruption mitigation system (DMS), hence an accurate description of the fragment ablation rate is desirable. Current mainstream 3D MHD codes modelling SPIs mostly assume instantaneous thermalization between the previously hot ambient electrons and the newly released cold electrons. This could result in underestimation of the ablation rate if the hot electron thermalization time is comparable or even longer than the fragment flying time across flux tubes, putting the validity of these 3D SPI simulations into question. To resolve this doubt, we, in this study, investigate the thermalization dynamics and the overall hot-electron impact. The finite-time collisional thermalization of hot-tail electrons in a rapidly cooling plasma, as well as the so-called ‘self-limiting’ effect are considered. The former effect tends to deplete the colder population within a hot-tail species, while the latter is found to preferentially deplete the higher energy population. The combined result is found to cause an almost self-similar decay of the hot electron distribution function, while its shape does not deviate much from that of Maxwellian distribution and the mean energy does not change much during the thermalization process. Based on this observation, axisymmetric JOREK D2 SPI simulations were carried out with additional hot-tail contribution to evaluate their overall impact onto the injection assimilation and penetration. It is found that the hot-tail effect indeed causes enhanced assimilation and shallower penetration, although the overall effect depends on the exact injection configuration, with the slow injection showing negligible hot-tail effect while the fast single non-shattered pellet case shows drastic hot-tail ablation enhancement. For ITER-like SPI parameters, there is no significant deviation in the total assimilation, but some deviation in the injection penetration is observed for the fast injection velocity cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observation and analysis of lower-hybrid-current-drive density limit in EAST

Lower hybrid current drive (LHCD) experiments with line-averaged density up to ~5.1 × 10 19 m –3 were performed in EAST L-mode plasmas. When the line-averaged density rises above a critical value, the hard x-ray (HXR) emission falls to the noise level, indicating that the LHCD density limit is encountered. The experimental results show that the LHCD density limit can be increased with higher wave source frequency (ƒ 0 ) and higher magnetic field (B t ). Although a higher LHCD density limit is obtained by a higher magnetic field for both 2.45 GHz and 4.6 GHz waves, the results show a stronger dependence on the magnetic field for the 4.6 GHz case. Analysis suggests that, for normal operation with a relatively low magnetic field (1.6 T $\leqslant$ B t $\leqslant$ 2.5 T) on EAST, the dominant mechanisms responsible for the LHCD density limit are different between the 2.45 GHz and 4.6 GHz waves. The wave accessibility plays a more significant role during 4.6 GHz LHCD experiments, while parasitic losses due to parametric decay instability (PDIs) dominate the accessibility issue in the 2.45 GHz case. Collisional loss in the scrape-off layer (SOL) may explain the 4.6 GHz result when combined with the accessibility limit at high density and low temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The constraint of plasma power balance on runaway avoidance

Abstract In a post-thermal-quench plasma, mitigated or unmitigated, the plasma power balance is mostly between collisional or Ohmic heating and plasma radiative cooling. In a plasma of atomic mixture { n α } with α labeling the atomic species, the power balance sets the plasma temperature, ion charge state distribution { n α i } with i the charge number, and through the electron temperature T e and ion charge state distribution { n α i } , the parallel electric field E ∥ . Since the threshold electric field for runaway avalanche growth E a v is also set by the atomic mixture, ion charge state distribution and its derived quantity, the electron density n e , the plasma power balance between Ohmic heating and radiative cooling imposes a stringent constraint on the plasma regime for avoiding and minimizing runaways when a fusion-grade tokamak plasma is rapidly terminated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

How fluctuation intensity flux drives SOL expansion

Abstract Predictions of heat load widths λ q based on particle orbits alone are very pessimistic. This paper shows that pedestal peeling-ballooning (P-B) magnetohydrodynamic (MHD) turbulence broadens the stable scrape-off layer (SOL) by the transport, or spreading, of fluctuation energy from the pedestal. λ q is seen to increase with Γ ε , the fluctuation energy density flux. We elucidate the fundamental physics of the spreading process. Γ ε increases with pressure fluctuation correlation length. P-B turbulence is seen to be especially effective at spreading, on account of its large effective mixing length. Spreading is shown to be a multiscale process, which is enhanced by the synergy of large and small-scale modes. Pressure fluctuation skewness correlates well with the spreading flux–with the zero crossing of skewness and Γ ε spatially coincident–suggesting the role of coherent fluctuation structures and the presence of intermittency in λ q broadening. λ q ∼ B p − 1 scaling persists for the broadened SOL. We show that the spreading flux increases for increasing pedestal pressure gradient ∇ P 0 and for decreasing pedestal collisionality υ ped ∗ . This trend is due to the dominance of peeling modes for large ∇ P 0 and low υ ped ∗ . Ultimately, we see that a state of weak MHD turbulence, as for small ELMs, is very attractive for heat load management. Our findings have transformative implications for future fusion reactor designs and call for experimental investigations to validate the observed trends.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Density wavenumber spectrum measurements, synthetic diagnostic development, and tests of quasilinear turbulence modeling in the core of electron-heated DIII-D H-mode plasmas

Abstract Measurements of the turbulent density wavenumber spectrum, δ n ˆ e ( k ⊥ ) , using the Doppler Back-Scattering (DBS) diagnostic are reported from DIII-D H-mode plasmas with electron cyclotron heating as the only auxiliary heating method. These electron-heated plasmas have low collisionality, ν e ∗ < 1 , T e / T i > 1 , and zero injected torque—a regime expected to be relevant for future fusion devices. We probe density fluctuations in the core ( ρ ≈ 0.7) over a broad wavenumber range, 0.5 ⩽ k ⊥ ⩽ 16 cm −1 ( 0.1 ⩽ k ⊥ ρ s ⩽ 5 ), to characterize plasma instabilities and compare with theoretical predictions. We present a novel synthetic DBS diagnostic to relate the back-scattered power spectrum, P s ( k ⊥ ) —which is directly measured by DBS—to the underlying electron density fluctuation spectrum, δ n ˆ e ( k ⊥ ) . The synthetic DBS P s ( k ⊥ ) spectrum is calculated by combining the SCOTTY beam-tracing code with a model δ n ˆ e ( k ⊥ ) predicted either analytically or numerically. In this work we use the quasi-linear code Trapped Gyro-Landau Fluid (TGLF) to approximate the δ n ˆ e ( k ⊥ ) spectrum. We find that TGLF, using the experimental profiles, is capable of closely reproducing the DBS measurements. Both the DBS measurements and the TGLF-DBS synthetic diagnostic show a wavenumber spectrum with variable decay. The measurements show weak decay ( k −0.6 ) for k < 3.5 cm −1 , with k −2.6 at intermediate- k ( 3.5 ⩽ k ⩽ 8.5 cm −1 ), and rapid decay ( k −9.4 ) for k > 8.5 cm −1 . Scans of physics parameters using TGLF suggest that the normalized ∇ T e scale-length, R / L T e , is an important factor for distinguishing microturbulence regimes in these plasmas. A combination of DBS observations and TGLF simulations indicate that fluctuations remain peaked at ITG-scales (low k ) while R / L T e -driven TEM/ETG-type modes (intermediate/high k ) are marginally sub-dominant.

synthetic diagnostics↗
DOE OSTI2024

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗
DOE OSTI2024

NSTX-U research advancing the physics of spherical tokamaks

Abstract The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Abstract Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4–9 µ s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3–4 µ s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept as a last line of defense against runaway damage in ITER, SPARC, and other future burning-plasma tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗