SEARCH · Engineering Papers

Results for “collisional plasma”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Marginal stability constraint on runaway electron distribution

High-frequency kinetic instabilities of the strongly anisotropic runaway electrons (RE) can enhance the pitch-angle scattering of the runaways significantly. This wave-induced scattering can easily prevail over runaway scattering on high-Z impurities. In a steady state, collisional damping balances the kinetic drive of the unstable waves, keeping the RE distribution function at marginal stability. The marginal stability constraint limits the achievable RE densities and the shape of the RE distribution function. In this study, we consider whistler and compressional Alfvén waves as the primary source of enhanced elastic scattering of the runaways. By balancing the anomalous Doppler resonance drive with the collisional wave damping, we find the RE distribution function in the ultra-relativistic range of the phase space. We also derive an expression for the spectral energy density of the waves. We show that the power needed to compensate for the wave dissipation is negligible compared to the work of the electric field. The latter is in balance with the synchrotron losses of the runaway electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Turbulence-Driven Edge-Localized-Mode-Free High-Confinement Mode with Divertor Detachment in a Metal-Wall Tokamak

We report the first demonstration of a minute-scale, edge-localized-mode-free high-confinement plasma regime compatible with divertor partial detachment and enhanced pedestal performance in a metal-wall tokamak, the Experimental Advanced Superconducting Tokamak. This regime is enabled by a newly identified mechanism: during divertor partial detachment, reduced ionization and enhanced pumping in a closed divertor lead to less cooling of the pedestal by recycling neutrals and seeding impurities, resulting in an increased pedestal temperature gradient, which excites high-frequency broadband turbulence. Gyrokinetic simulations identify the high-frequency broadband turbulence as a temperature-gradient-driven trapped electron mode (𝜂 𝑒 -TEM), which drives outward transport of particles and heat, thereby maintaining the edge-localized-mode-free state. This pedestal regime is particularly promising for the International Thermonuclear Experimental Reactor, where the anticipated lower density gradient, reduced E×B shear, and lower collisionality in the pedestal will further facilitate 𝜂 𝑒 -TEM excitation. The achieved integrated scenario with a detached divertor and turbulence-dominated pedestal thus offers a compelling solution for managing heat loads and metal impurity sources for long-pulse high-performance operation in future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1970

Spatial cyclotron damping

To examine spatial electron cyclotron damping in a uniform Vlasov plasma, it is noted that the plasma response to a steady-state transverse excitation consists of several terms (dielectric-pole, free-streaming, and branch-cut), but that the cyclotron-damped pole term is the dominant term for z l = c/w sub ce provided (w sub pe/w sub ce) squared (c/a) is much greater than 1. If the latter inequality does not hold, then the free-streaming and branch-cut terms persist well past z = c/w sub ce as w sub 1 approaches w sub ce, making experimental measurement of cyclotron damping essentially impossible. Considering only (w sub pe/w sub ce) squared (c/a) is much greater than 1, it is shown how collisional effects should be estimated and how a finite-width excitation usually has little effect on the cyclotron-damped part of the response. Criteria is established concerning collisional damping, measurable damping length sizes, and allowed uncertainty in the magnetic field Beta. Results of numerical calculations, showing the regions in the appropriate parameter spaces that meet these criteria, are presented. From these results, one can determine the feasibility of, or propose parameter values for, an experiment designed to measure spatial cyclotron damping. It is concluded that the electron temperature T sub e should be at least 1 ev., and preferably 10 ev. or higher, for a successful experiment.

Olson, C. L.↗

DOE OSTI2020

Deep learning surrogate model for kinetic Landau-fluid closure with collision

In this work, the kinetic Landau-fluid (LF) closure with collision and periodic boundary condition is used in the development of the deep learning (DL) surrogate model. A classical neural network, namely, feedforward neural network or sometimes termed multilayer perceptron, is constructed and trained to learn the kinetic LF closure in the static limit and arbitrary mean free path in configuration space. The preliminary relation between best hyperparameters and critical parameters for data generation is found. Compared with the numerical approach (non-Fourier method) of the LF closure, the deep learning surrogate model shows an order of magnitude of improvement in terms of accuracy. Perhaps most importantly, the surrogate model closure has been integrated for the first time with fluid simulations. Our DL-enabled fluid simulations, for the first time, give the correct Landau damping rate for a wide range of wave vectors, while the Hammett–Perkins closure cannot produce the correct damping rate. We correctly connect the collisionless Hammett–Perkins closure and collisional Braginskii closure to reproduce the intrinsic nonlocal feature of the heat flux with DL techniques. We address the most concerning error accumulation problem and find that simulations with the deep learning surrogate model are as good as, if not better than, simulations with the analytic closure in terms of long-term numerical stability in the linear Landau damping test.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Investigation of divertor detachment induced through neon seeding and density ramp on HL-3

A new self-consistent 1D scrape-off layer model has been recently developed in BOUT++ framework, named SD1D, which includes equations for various particle species (e.g. main plasma, neutrals and impurities) and couples open databases like ADAS and AMJUEL. It is able to quickly and effectively simulate divertor detachment experiments. In this work, a typical detachment experiment (shot #6270) on HL-3 with neon seeding is simulated using the SD1D code. It is found that the target electron temperature and the target ion saturation current in the simulations are consistent with experimental results measured by Langmuir probes on the target plate. The variation of D α radiation intensity in the divertor is qualitatively similar to the measured D α signal. Following the experimental validations, different upstream densities are set in the simulations to study the impurity distribution under different plasma density conditions. It is found that increasing upstream density can be helpful for the control of the neon radiation front (closer to the target). In this work we also compare two detachment regimes in simulations. Based on the same initial experimental parameters (shot #6270) on HL-3, a scan of upstream density and a scan of neon seeding rate are carried out respectively. It is found that the role of atomic and molecular processes is different in the two detachment regimes. The current density roll-over is ascribed to a drop in the divertor ion source, and the variation of D α radiation intensity via different excitation channels is associated with the relevant collisional reaction sources.

BOUT++↗

DOE OSTI2022

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Transition characteristics and electron kinetics in microhollow cathode discharges

We demonstrate the transition characteristics and electron kinetics of microdischarges in a microgap with a cathode having microstructures using two-dimensional particle-in-cell/Monte Carlo collision (2d3v) simulations. It is found that the microdischarge is inside the hollow cavity at lower pressures, forming microhollow cathode discharges (MHCDs), while the dominant discharge moves outside the cavity at higher pressures, locating above the cathode rectangular protrusion tip. The spatial distributions of the microdischarge parameters (e.g., electric potential and electron density) are presented by capturing the transition characteristics. The electron kinetics of the microdischarges are examined based on the moment analysis of the electron Boltzmann equation. The collisional and non-collisional components of the electron power absorption are decomposed, which were seldomly investigated for microdischarges previously. The results provide fundamental insights into MHCD formation with a structured electrode at varying pressure conditions, which could be beneficial for the design and fabrication of microplasma devices in practical applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

A conservative discontinuous Galerkin algorithm for particle kinetics on smooth manifolds

A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.

Discontinuous Galerkin↗

DOE OSTI2024

Two-dimensional Thomson scattering measurements of misaligned electron density and temperature gradients and associated Biermann battery produced fields

We present optical Thomson scattering measurements of electron density and temperature in high Mach number laser-driven blast waves in homogeneous gases. Taylor–Sedov blast waves are launched in nitrogen (N2) or helium (He) at pressures between 0.4 mTorr and 10 Torr by ablating a solid plastic target with a high energy laser pulse (10 J, 1012 W/cm2). Experiments are performed at high repetition rate (1 Hz), which allows one-dimensional and two-dimensional Thomson scattering measurements over an area of several cm2 by automatically translating the scattering volume between shots. Electron temperature and density in the blast wave fronts were seen to increase with increasing background gas pressure. Measured electron density and temperature gradients were used to calculate ∂B/∂t∝∇Te×∇ne. The experimentally measured ∂B/∂t showed agreement with the magnetic field probe (B-dot) measurements, revealing that magnetic fields are generated in the observed blast waves via the Biermann battery effect. The results are compared to numerical three-dimensional collisional magnetohydrodynamic simulations performed with FLASH, and are discussed in the context of spontaneous magnetic field generation via the Biermann battery effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Atomic Oxygen (ATOX) simulation of Teflon FEP and Kapton H surfaces using a high intensity, low energy, mass selected, ion beam facility

A high intensity (greater than 10(exp 15) ions/sq cm) low energy (down to 5 eV) mass selected ion beam (MSIB) facility was used to study the effects of ATOX on two polymers commonly used for space applications (Kapton H and Teflon FEP). The polymers were exposed to O(+) and Ne(+) fluences on 10(exp 15) - 10(exp 19) ions/sq cm, using 30eV ions. A variety of analytical methods were used to analyze the eroded surfaces including: (1) atomic force microscopy (AFM) for morphology measurements; (2) total mass loss measurements using a microbalance; (3) surface chemical composition using x-ray photoelectron spectroscopy (XPS), and (4) residual gas analysis (RGA) of the released gases during bombardment. The relative significance of the collisional and chemical degradation processes was evaluated by comparing the effects of Ne(+) and O(+) bombardment. For 30 eV ions it was found that the Kapton is eroded via chemical mechanisms while Teflon FEP is eroded via collisional mechanisms. AFM analysis was found very powerful in revealing the evolution of the damage from its initial atomic scale (roughness of approx. 1 nm) to its final microscopic scale (roughness greater than 1 micron). Both the surface morphology and the average roughness of the bombarded surfaces (averaged over 1 micron x 1 micron images by the system's computer) were determined for each sample. For 30 eV a non linear increase of the Kapton roughness with the O(+) fluence was discovered (a slow increase rate for fluences phi less than 5 x 10(exp 17) O(+)/sq cm, and a rapid increase rate for phi greater than 5 x 10(exp 17) O(+)/sq cm). Comparative studies on the same materials exposed to RF and DC oxygen plasmas indicate that the specific details of the erosion depend on the simulation facility emphasizing the advantages of the ion beam facility.

Vered, R.↗

DOE OSTI2023

Neoclassical transport of impurities in tokamaks with non-axisymmetric perturbations

The effect of resonant magnetic perturbations (RMPs) on neoclassical transport of impurities is calculated in DIII-D and NSTX. Here, the neoclassical fluxes are evaluated using the NEO code with nonlinear one-fluid nonaxisymmetric equilibrium calculated using M3D-C1. Neoclassical fluxes of impurities show significant changes with RMPs if the impurity resides in low-collisionality regime, but are weakly affected by RMPs in the Pfirsch–Schlüter (P–S) regime. Charge number (Z) of the impurity affects the collisionality of impurity species, which then determines the effect of 3D fields on neoclassical impurity transport. This suggests that RMPs can be possibly used for exhaust of low-Z impurities in these facilities, but have little effect for impurities with charge number greater than 10 or so. Additionally, it is shown that it is the change in the convective fluxes which is a main contributor in modifying the neoclassical impurity fluxes in the presence of RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2002

A Chandra X-ray Study of Cygnus A: The Nucleus - 2

We report Chandra Advanced CCD Imaging Spectrometer and quasi-simultaneous Rossi X-Ray Timing Explorer (RXTE) observations of the nearby, powerful radio galaxy Cygnus A, with the present paper focusing on the properties of the active nucleus. In the Chandra observation, the hard (less than a few keV) X-ray emission is spatially unresolved with a size is approximately 1" (1.5 kpc, H(sub 0) = 50 km/s/Mpc) and coincides with the radio and near-infrared nuclei. In contrast, the soft (less than 2 keV) emission exhibits a bipolar nebulosity that aligns with the optical bipolar continuum and emission-line structures and approximately with the radio jet. In particular, the soft X-ray emission corresponds very well with the [O III] (lambda)5007 and H(alpha) + [N II] lambda(lambda)6548, 6583 nebulosity imaged with Hubble Space Telescope. At the location of the nucleus, there is only weak soft X-ray emission, an effect that may be intrinsic or result from a dust lane that crosses the nucleus perpendicular to the source axis. The spectra of the various X-ray components have been obtained by simultaneous fits to the six detectors. The compact nucleus is detected to 100 keV and is well described by a heavily absorbed power-law spectrum with Gamma(sub h) = 1.52(sup + 0.12, sub -0.12) (similar to other 0.12 narrow-line radio galaxies) and equivalent hydrogen column N(sub H)(nuc) = 2.0(sup +0.1, sub -0.1) x 10(exp 23)/sq cm. This 0.2 column is compatible with the dust obscuration to the near-infrared source for a normal gas-to-dust ratio. The soft (less than 2 keV) emission from the nucleus may be described by a power-law spectrum with the same index (i.e., Gamma(sub l) = Gamma(sub h), although direct fits suggest a slightly larger value for Gamma(sub l). Narrow emission lines from highly ionized neon and silicon, as well as a "neutral" Fe K(alpha) line, are detected in the nucleus and its vicinity (r approximately less than 2 kpc). The equivalent width (EW) of the Fe K(alpha) line (182(sup +40, sub -54) eV) is in good agreement with theoretical predictions for the EW versus N(sub H)(nuc) relationship in various geometries. An Fe K edge is also seen. The RXTE observations indicate a temperature of kT = 6.9(sup +0., sub -1.0) keV for the cluster gas (discussed in Paper III of this series) and cluster emission lines of Fe K(alpha) and Fe K(beta) and/or Ni K(alpha). We consider the possibility that the extended soft X-ray emission is electron-scattered nuclear radiation. Given that 1% of the unabsorbed 2 - 10 keV nuclear radiation would have to be scattered, the necessary gas column [N(sub H)(Scattering) approx. = 3.5 x 10(exp 22)/sq cm] would absorb the X-rays rather than scatter them if the gas is cold. Thus, the scattering plasma must be highly ionized. If this ionization is achieved through photoionization by the nucleus, the ionization parameter zeta greater than 1 ergs cm/s and the electron density n(sub e) approx. = 6 cc given the observed distance of the soft X-ray emission from the nucleus. The electron column density inferred from the X-ray observations is much too low to account for the extended optical scattered light, strongly suggesting that the polarized optical light is scattered by dust. The presence of highly ionized Ne lines in the soft X-ray spectrum requires 20 ergs cm/s approximately less than zeta approximately less than 300 ergs cm/s these lines may originate closer to the nucleus than the extended soft continuum or in a lower density gas. A collisionally ionized thermal model of the extended soft X-rays cannot be ruled out but is unattractive in view of the low metal abundance required (Z = 0.03 Z(mass)). The hard X-ray to far-infrared ratio for the nucleus of Cygnus A is similar to that seen in Seyfert 1 and unobscured radio galaxies. By means of the correlation between hard X-ray luminosity and nuclear optical absolute magnitude for these classes of object, we estimate M(sub B) = -22.4 for Cygnus A, near the .borderline between Seyfert galaxies and QSOs.

Young, Andrew J.↗

DOE OSTI2021

Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

Physics↗

DOE OSTI2026

Spectrally accurate, reverse-mode differentiable bounce-averaging algorithm and its applications

We present a fast, spectrally (exponentially) accurate, automatically differentiable bounce-averaging algorithm that is used to simplify kinetic models. Using this algorithm, implemented in the DESC stellarator optimisation suite, we can perform efficient optimisation of many objectives to improve stellarator performance, such as the effective ripple 𝜖 eff metric for the neoclassical transport coefficient in the low collisionality regime and proxies for energetic particle confinement. For the first time, we optimise a finite-beta stellarator to directly reduce neoclassical ripple transport using reverse-mode differentiation. This ensures the computational cost of differentiation is independent of the number of controllable parameters.

fusion plasma↗

DOE OSTI2026

First DIII-D-West hybrid scenario similarity experiments for iter-relevant long-pulse operation

For the first time, similarity experiments between DIII-D and WEST were performed in the ITER "hybrid-like" regime during dedicated campaigns in April and May 2025. The matched parameters include elongation, triangularity, ion ∇B drift direction toward the X-point, qprofile, and core normalized physics quantities in terms of normalized pressure, normalized gyroradius, electron collisionality, ratio of ion to electron temperature, T i /T e . Core transport physics is explored with different aspect ratio (R/a) values (typically 3 at DIII-D and 5 on WEST). DIII-D explored high-beta conditions (electromagnetic effect) with low torque injection (~0 ± 0.5 N•m) using high heating power (up to 6 MW NBI and 2 MW ECRH powers), while scanning the heating mix (ion vs electron), beta, T i /T e , core radiation via controlled tungsten injection using the Laser Blow-Off system. WEST extended operation toward long-duration pulses using its actively cooled tungsten divertor, achieving dominated electron heating regimes with reduced tungsten contamination. Boron impurity injection were scanned on WEST to control edge conditions and core performance. It is found that core confinement improves-manifested by higher electron temperature, total energy content, neutron rate, and ion temperatureunder conditions of low separatrix density, consistent with previous observations [Bourdelle et al., Nucl. Fusion 63 (2023) 056021]. Conditions for Hmode access and for ion heating in electron-dominated regimes in both WEST and DIII-D will be discussed and compared. The ratio of the thermal energy confinement time (τ E ) to the volume-averaged electron-ion collisional heat exchange time (τ e-i ) is a key parameter to enhance ion heating and potentially facilitate H-mode access in electron-heated regimes. These first-of-a-kind coordinated DIII-D and WEST experiments provide a unique multi-machine dataset to validate predictive models and to optimize ITER hybrid-scenario performance under diverse core and edge conditions.

DIII-D↗

NASA NTRS2017

Kinetic Theory and Fast Wind Observations of the Electron Strahl

We develop a model for the strahl population in the solar wind - a narrow, low-density and high-energy electron beam centred on the magnetic field direction. Our model is based on the solution of the electron drift-kinetic equation at heliospheric distances where the plasma density, temperature and the magnetic field strength decline as power laws of the distance along a magnetic flux tube. Our solution for the strahl depends on a number of parameters that, in the absence of the analytic solution for the full electron velocity distribution function (eVDF), cannot be derived from the theory. We however demonstrate that these parameters can be efficiently found from matching our solution with observations of the eVDF made by the Wind satellite's SWE strahl detector. The model is successful at predicting the angular width (FWHM) of the strahl for the Wind data at 1 au, in particular by predicting how this width scales with particle energy and background density. We find that the strahl distribution is largely determined by the local temperature Knudsen number γ ∼ |T dT/dx|/n, which parametrizes solar wind collisionality. We compute averaged strahl distributions for typical Knudsen numbers observed in the solar wind, and fit our model to these data. The model can be matched quite closely to the eVDFs at 1 au; however, it then overestimates the strahl amplitude at larger heliocentric distances. This indicates that our model may be improved through the inclusion of additional physics, possibly through the introduction of 'anomalous diffusion' of the strahl electrons.

Horaites, Konstantinos↗

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↗