Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Space plasma physics”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Propagation of ion cyclotron emission in the DIII-D tokamak

Experiments in L- and H-mode plasmas performed on the DIII-D tokamak explored ion cyclotron emission (ICE) propagation via the recently upgraded ICE diagnostic. The distance between the plasma and the outer wall was scanned to alter the evanescent region, which was seen to impact the detection of runaway electron-driven whistler modes in a comparable frequency range to typical ICE harmonics (f≳10fci) [Heidbrink et al., Plasma Phys. Controlled Fusion 61, 014007 (2018)]. In L-mode plasmas, core ICE did not exhibit clear trends as the plasma-wall distance was altered. Instead, inadvertent changes to the fast ion density correlated with different ICE fine structures; the most spectral bands (spaced ∼150–200 kHz apart) were observed at the highest fast ion density, and a just single band when nfast decreased by ∼ 40%. The full-wave Petra-M code simulated core ICE propagation through the plasma and to ICE probes, agreeing with experimental findings that these are likely compressional fast waves, which propagate to probe locations with no dependence on plasma-wall gap. Edge ICE in H-mode plasmas was similarly unaffected by plasma-wall distance. Unlike in core ICE cases, there were negligible changes to the edge fast ion distribution and edge bulk plasma profiles. However, changes in ELM frequency seemed most closely tied to ICE behavior. Ultimately, core and edge ICE did not exhibit clear dependencies on the plasma-wall distance, and ICE harmonics were observed with even the largest plasma-wall gaps. This independence bodes well for similar diagnostics in future tokamaks, which might be placed further from high radiation areas without adversely impacting measurement capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Three-dimensional plasmoid-mediated reconnection and turbulence in Hall magnetohydrodynamics

Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may transition into the collisionless regime as the fragmentation of the current sheet progresses toward kinetic scales. Hall magnetohydrodynamics (MHD) models are widely regarded as a simplified yet effective representation of the transition from collisional to collisionless reconnection. However, plasmoid instability in 2D Hall MHD simulations often leads to a single-X-line reconnection configuration, which significantly differs from fully kinetic particle-in-cell simulation results. This study shows that single-X-line reconnection is less likely to occur in 3D compared to 2D. Moreover, depending on the Lundquist number and the ratio between the system size and the kinetic scale, Hall MHD can also realize 3D self-generated turbulent reconnection. We analyze the features of the self-generated turbulent state, including the energy power spectra and the scale dependence of turbulent eddy anisotropy.

Physics↗
DOE OSTI2023

Multimodal and multiscale strengthening mechanisms in Al-Ni-Zr-Ti-Mn alloy processed by laser powder bed fusion additive manufacturing

The unique thermokinetics of laser-powder bed fusion additive manufacturing (L-PBFAM) has been exploited for development of a novel high-strength Al-Ni-Ti-Zr-Mn alloy. The addition of 0.5 wt% Mn leads to extraordinary improvement in ultimate tensile strength (502 MPa) and work hardening due to the activation of two Mn-induced strengthening mechanisms. First, by a bimodal particle strengthening effect due to Al 31 Mn 6 Ni 12 nano-quasi-crystal particles rejected in inter-dendritic spaces and fibrous Al 3 Ni eutectic dendritic channels, which predominately contributes to the strength improvement, and second by solid solution strengthening from remaining Mn entrapped in Al. These mechanisms supplement the particle strengthening effect imparted by coherent and incoherent Al 3 (Ti,Zr) co-precipitates present at melt pool boundaries, dislocation strengthening due to solidification induced strain, and Hall-Petch and backstress strengthening effect due to heterogenous grain size distribution occurring at various length scales. In conclusion, the solidification dynamics and hierarchical heat distribution that are associated with L-PBFAM resulted in complex spatial variations in these strengthening phenomena and were investigated via a high-throughput multiscale structure–property correlation that involved thermokinetic simulation, X-ray diffraction, high-resolution nanoindentation mapping, and site-specific transmission electron microscopy of the alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hamiltonian reduction of Vlasov–Maxwell to a dark slow manifold

In this work, we show that non-relativistic scaling of the collisionless Vlasov–Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov–Maxwell phase space. Vlasov–Maxwell dynamics restricted to the slow manifold recovers the Vlasov–Poisson and Vlasov–Darwin models as low-order approximations, and provides higher-order corrections to the Vlasov–Darwin model more generally. The slow manifold may be interpreted to all orders in perturbation theory as a collection of formal Vlasov–Maxwell solutions that do not excite light waves, and are therefore ‘dark’. We provide a heuristic lower bound for the time interval over which Vlasov–Maxwell solutions initialized optimally near the slow manifold remain dark. We also show how the dynamics on the slow manifold naturally inherits a Hamiltonian structure from the underlying system. After expressing this structure in a simple form, we use it to identify a manifestly Hamiltonian correction to the Vlasov–Darwin model. The derivation of higher-order terms is reduced to computing the corrections of the system Hamiltonian restricted to the slow manifold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling active control of resistive wall mode with power saturation and sensor noise on HL-2M

The resistive wall mode (RWM) control on the HL-2M tokamak is simulated with the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681), aiming at quantifying control current and voltage requirements when more realistic issues are taken into account, i.e. the control power saturation and the sensor signal noise. The fluid model predicts a narrow stability region for the n = 1 RWM without magnetic feedback, in the 2D parameter space of the plasma pressure versus the toroidal flow speed. Magnetic feedback can fully stabilize the RWM on HL-2M. Without considering the voltage limitation and the sensor signal noise, it is found that plasma flow helps active control of the mode, by reducing the required critical feedback gain for both flux-to-current and flux-to-voltage control schemes. In the absence of the sensor signal noise, the lowest control voltage saturation level, below which the RWM control is lost, is found to roughly satisfy a linear relation to the plasma flow frequency, indicating that subsonic plasma flow is effective in relaxing the control power requirement for the RWM feedback stabilization. The presence of the sensor signal noise substantially modifies the feedback results. A statistical study finds that the sensor signal noise, with the standard deviation of 0.1 G on HL-2M, roughly doubles the required control voltage for successful mode control. The synergistic stabilization effect due to plasma flow is somewhat weakened by the presence of the sensor signal noise. At a given rotation, the tolerable voltage limit generally increases with increasing feedback gain due to the sensor signal noise.

Physics↗
DOE OSTI2025

Intrinsic bulk viscosity of the one-component plasma

The intrinsic bulk viscosity of the one-component plasma (OCP) is computed and analyzed using equilibrium molecular dynamics simulations and the Green-Kubo formalism. It is found that bulk viscosity exhibits a maximum at Γ ≈ 1, corresponding to the condition that the average kinetic energy of particles equals the potential energy at the average interparticle spacing. The weakly coupled and strongly coupled limits are analyzed and used to construct a model that captures the full range of coupling strengths simulated: Γ ≈ 10 –2 –10 2 . Simulations are also run of the Yukawa one-component plasma (YOCP) in order to understand the impact of electron screening. It is found that electron screening leads to a smaller bulk viscosity due to a reduction in the excess heat capacity of the system. Bulk viscosity is shown to be at least an order of magnitude smaller than shear viscosity in both the OCP and YOCP. The generalized frequency-dependent bulk viscosity coefficient is also analyzed. This is found to exhibit a peak near twice the plasma frequency in strongly coupled conditions, which is associated with the oscillatory decay observed in the bulk viscosity autocorrelation function. As a result, the generalized shear and bulk viscosity coefficients are found to have a similar magnitude for ω ≳ 2⁢ω p at strongly coupled conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Growth of concomitant laser-driven collisionless and resistive electron filamentation instabilities over large spatiotemporal scales

Collective processes in plasmas often induce microinstabilities that play an important role in many space or laboratory plasma environments. Particularly of note is the Weibel-type current filamentation instability, which is believed to drive the creation of collisionless shocks in weakly magnetized astrophysical plasmas. In this paper, this instability class is studied through interactions of ultraintense and short laser pulses with solid foils, leading to localized generation of megaelectronvolt electrons. Proton radiographic measurements of both low- and high-resistivity targets show two distinct, superimposed electromagnetic field patterns arising from the interpenetration of the megaelectronvolt electrons and the background plasma. Particle-in-cell simulations and theoretical estimates suggest that the collisionless Weibel instability building up in the dilute expanding plasmas formed at the target surfaces causes the observed azimuthally symmetric electromagnetic filaments. For a sufficiently high resistivity of the target foil, an additional resistive instability is triggered in the bulk target, giving rise to radially elongated filaments. The data reveal the growth of both filamentation instabilities over large temporal (tens of picoseconds) and spatial (hundreds of micrometres) scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Multispecies Ion Acceleration in 3D Magnetic Reconnection with Hybrid-Kinetic Simulations

Magnetic reconnection drives multispecies particle acceleration broadly in space and astrophysics. Here we perform the first 3D hybrid simulations (fluid electrons, kinetic ions) that contain sufficient scale separation to produce nonthermal heavy-ion acceleration, with fragmented flux ropes critical for accelerating all species. We demonstrate the acceleration of all ion species (up to Fe) into power-law spectra with similar indices, by a common Fermi acceleration mechanism. The upstream ion velocities influence the first Fermi reflection for injection. The subsequent onsets of Fermi acceleration are delayed for ions with lower charge-mass ratios (Q/M), until growing flux ropes magnetize them. This leads to a species-dependent maximum energy/nucleon ∝(Q/M) α . These findings are consistent with in situ observations in reconnection regions, suggesting Fermi acceleration as the dominant multispecies ion acceleration mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Achieving a robust grassy-ELM operation regime in CFETR

We have identified a robust grassy-edge localized mode (ELM) operation regime for future tokamak reactors. The regime exists within a pedestal top electron collisionality (ν*) window at high global poloidal beta (β p ). The existence of an upper ν* limit for grassy-ELMs is consistent with results previously reported in experiments (Oyama et al 2010 Nucl. Fusion 50 064014), while the existence of a lower ν* limit has not been reported previously. Using EPED and BOUT++, a theoretical model that quantitatively explains the physics of the grassy-ELMs within the window, which distinguishes them from the small mixed-ELMs at lower ν*, is presented for the first time. A peeling-ballooning stability boundary is obtained by scanning the operating density space. The change in density corresponds to a change in ν* that affects the pedestal bootstrap current. High β p leads to a strong Shafranov shift, which affects the flux surface averaged pressure drive. The two effects combine to create a peeling-dominated window in intermediate ν* buffered by ballooning-dominated regimes. Only the peeling-dominated regime shows a cyclic behavior in the perturbed pressure during the nonlinear simulation of an ELM crash, reminiscent of grassy-ELM dynamics. Similarly, the energy released across the separatrix is demonstrated to be significantly smaller. The quick recovery of the ELM crash is explainable by the rapid rise of a low n kink-peeling instability when the pedestal current I ped exceeds a threshold at high β p . It minimizes the excursion beyond marginal stability and is absent in the ballooning-dominated regime. Comparison with recent experiments over a range of β p and ν* strongly supports the physical picture proposed by the modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

First fluctuation measurements using an Imaging Neutral Particle Analyzer on DIII-D

A recent upgrade to the Imaging Neutral Particle Analyzer (INPA) on DIII-D has allowed for the first fluctuation measurements using an INPA to be taken during a neoclassical tearing mode (NTM). The INPA signal tracked the mode over 150 ms as the mode frequency dropped to zero, capturing both the NTM with poloidal and toroidal mode numbers $m/n=2/1$ and a $3/2$ mode at double the frequency. Analysis shows that the signals originate from charge exchange events near the edge of the plasma, and relative fluctuation amplitudes are greater than 25% for the duration of the NTM. Filtered signals show frequency beating patterns that are phase-space dependent. Simulated signal is dominated by prompt transport from the neutral beams to the INPA sightline, while the contribution from the slowing down distribution is significantly lower. Simulated measurements in the range of pitches that the diagnostic is sensitive to (0.5≤|v∥/v|≤ 0.75) show the signal is dominated by trapped orbits that pass through magnetic islands near the edge of the plasma. Calculations of expected fluctuation levels show that only direct interaction with the NTM can provide the strong relative fluctuation levels seen experimentally. The prompt nature of the orbits and thin radial layer found to contribute to synthetic signals suggest INPA passive data may be used to measure the perturbations of confined orbits on a single pass through a plasma instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Radial coordinate maps, radial vectors, and binormal vectors for 5/6, 5/5 and 5/4 edge island domains in W7-X

Abstract The edge island domain in Wendelstein 7-X consists of divertor islands whose position coincide with the location of rational values of the rotational transform - ι ≈ ( 5 / 6 , 5 / 5 , 5 / 4 ) and surround the main confinement volume (the ‘main plasma’), providing the magnetic topology for an island divertor. The 5/5 edge consists of five individual islands that are unconnected. In contrast, a single island connects onto itself after 6 or 4 toroidal transits in the 5/6 or 5/4 edge, respectively. Many interesting phenomena are related to these islands and diagnostic analyses require a mapping from ‘laboratory’ or real space coordinates to the island coordinate system (ICS). A procedure is described here to calculate several scalar and vector quantities for closed island structures which can be utilized in fast interpolation schemes for inverse maps. For the 5/5 edge, a fixed-boundary vacuum (zero beta) magneto-hydrodynamic solution of the 5/5 island is found with the Variational Moments Equilibrium Code, VMEC. The solution is compatible with already existing routines which determine the radial and binormal vectors of VMEC solutions at arbitrary laboratory coordinates. VMEC does not support solutions for the 5/4 and 5/6 islands, but the radial and binormal vectors are available in a local 2D Fourier ICS.

Physics↗
DOE OSTI2023

Understanding the effects of fabrication process on BaZr 0.9 Y 0.1 O 3– δ grain-boundary chemistry using atom probe tomography

Proton-conducting oxides, such as BaZr 0.9 Y 0.1 O 3–δ (BZY10), are receiving significant scientific attention for application in a variety of electrochemical applications. It has been shown that the synthesis route has a significant effect on the ionic behavior of these materials, especially the grain boundaries (GB). In this study, laser-pulsed atom probe tomography (APT) is used to characterize the GB chemistry of BZY10 samples prepared by four different methods: (1) spark plasma sintering (SPS), (2) conventional sintering followed by high-temperature annealing (HT), (3) conventional sintering with NiO used as a sintering aid (SSR-Ni), and (4) solid-state reactive sintering of oxide precursor powders with NiO as a reactive sintering aid (SSRS-Ni). Oxygen depletion (corresponding to oxygen vacancy accumulation) was observed at every GB in this study, consistent with the positive space-charge GB region commonly assigned to these materials. In contrast to the consistent trends associated with oxygen depletion, cation constituents showed varying segregation/depletion behaviors. The samples that used NiO revealed Ni preferentially segregating to the GB. The most common impurities in BZY10 are Al, Fe, Mg, Si, and Sr, all of which generally accumulate at GBs. In order to quantify these chemical trends, the APT method was optimized for this material supplemented with understanding from density functional theory (DFT). Finally, this APT study reveals the complexity of BZY10 GB chemistry and shows that even adjacent GBs in the same material can sometimes show surprisingly different GB chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Effects of artificial collisions, filtering, and nonlocal closure approaches on Hermite-based Vlasov–Poisson simulations

Kinetic simulations of collisionless plasmas are computationally challenging due to phase-space mixing and filamentation, resulting in fine-scale velocity structures. This study compares three methods developed to reduce artifacts related to limited velocity resolution in Hermite-based Vlasov–Poisson simulations: artificial collisions, filtering, and nonlocal closure approaches. We evaluate each method's performance in approximating the linear kinetic response function and suppressing recurrence in linear and nonlinear regimes. Numerical simulations of Landau damping demonstrate that artificial collisions, particularly higher orders of the Lenard-Bernstein collisional operator, most effectively recover the correct damping rate across a range of wavenumbers. Moreover, Hou-Li filtering and nonlocal closures underdamp high wavenumber modes in linear simulations, and the Lenard-Bernstein collisional operator overdamps low wavenumber modes in both linear and nonlinear simulations. This study demonstrates that hypercollisions offer a robust approach to kinetic simulations, accurately capturing collisionless dynamics with limited velocity resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Operando XPS for Plasma Process Monitoring: A Case Study on the Hydrogenation of Copper Oxide Confined under h-BN

Here, we demonstrate that ambient pressure x-ray photoelectron spectroscopy (APXPS) can be used for in situ studies of dynamic changes in surface chemistry in a plasma environment. This opens a new and vast application space for XPS and greatly complements modern spectroscopy techniques to probe plasma-solid/liquid interactions relevant to process monitoring in the semiconductor industry, bio-medical plasma applications and plasma remediation technologies. Hexagonal boron nitride (h-BN) grown on Cu was used in this study as a well-defined model system for plasma process monitoring and because of its unique chemical, optical and electrical properties that make it a prospective material for advanced electronics. To better understand the stability and surface chemistry of h-BN during plasma assisted processing, we track in real time the plasma-induced chemical state changes of B, N and the underlying Cu substrate using APXPS equipped with an AC discharge plasma source operating at 13 Pa. Residual gas analysis (RGA) mass-spectra were concurrently collected during plasma-XPS to track reaction products formed during plasma exposure. A clear reduction of Cu x O is seen, while an h-BN layer remains intact, suggesting hydrogen radical (H • ) species can attack the exposed and h-BN covered Cu oxide patches and partially reduce the underlying substrate without significantly damaging the overlaying h-BN, which is of practical importance for development of h-BN encapsulated devices and interfaces. In addition to demonstration of plasma-XPS capabilities we discuss the observed challenges (e.g., parasitic plasma-chamber walls reactions and charging effects) and propose potential solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Phase-space sensitivity (weight functions) of 3 MeV proton diagnostics

Abstract The weight function W ( X ) is a diagnostic sensitivity to phase-space variables X that relates the measured signal C to the distribution function F ( X ) through the equation C = ∫ W ( X ) F ( X ) d X . In the present work, an algorithm to calculate W for a diagnostic that measures 3 MeV protons produced in d(d,p)t fusion reactions between a fast ion and a thermal deuteron is developed. The emitted protons escape the tokamak on curved orbits and are detected. These curved orbits constitute effective diagnostic ‘sightlines’. The presented algorithm accounts for the complications associated with these curved sightlines. An initial calculation of time-reversed proton orbits computes effective solid angles and sightlines for the relevant range of incident proton velocity vectors. These precomputed orbits are then used within the framework of FIDASIM (2020 Plasma Phys. Cont. Fusion 62 105008) to calculate the reactivity averaged over the thermal distribution of the ‘target’ deuterons and the probability that a fast ion of specified energy and pitch has a gyroangle that is consistent with the kinematic equations along each of the sightlines. Comparisons with analytic formulas and with independent calculations for the Mega Amp Spherical Tokamak 3 MeV proton diagnostic verify the algorithm.

Physics↗
DOE OSTI2025

Temporally-consistent koopman autoencoders for forecasting dynamical systems

Absence of sufficiently high-quality data often poses a key challenge in data-driven modeling of high-dimensional spatio-temporal dynamical systems. Koopman Autoencoders (KAEs) harness the expressivity of deep neural networks (DNNs), the dimension reduction capabilities of autoencoders, and the spectral properties of the Koopman operator to learn a reduced-order feature space with simpler, linear dynamics. However, the effectiveness of KAEs is hindered by limited and noisy training datasets, leading to poor generalizability. To address this, we introduce the Temporally-Consistent Koopman Autoencoder (tcKAE), designed to generate accurate long-term predictions even with limited and noisy training data. This is achieved through a consistency regularization term that enforces prediction coherence across different time steps, thus enhancing the robustness and generalizability of tcKAE over existing models. We provide analytical justification for this approach based on Koopman spectral theory and empirically demonstrate tcKAE’s superior performance over state-of-the-art KAE models across a variety of test cases, including simple pendulum oscillations, kinetic plasma, and fluid flow data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Low Loss, Fast Ramp HTS Solenoid Prototype for Compact Spherical Tokamaks

To explore AC loss characteristics of commercial Bi-2212 conductors made in the USA, a small ohmic heating (OH) prototype coil is being developed for compact tokamaks using high performance Bi-2212 wires newly manufactured by Bruker EST. Here, the goal is to build and test Bi-2212 prototype coils for fast ramp central solenoid (CS) operation of spherical tokamaks (ST). For a ST pilot plant, one of the key challenges is the tight space available for CS in a center column where in-board legs of the toroidal field (TF) and OH coils are assembled into a cryostat with a few millimeters clearance. A high current density Rutherford cable consisting of 17 legacy Bi-2212 wires, a contribution of the US Magnet Development Program (MDP), was fabricated at the Lawrence Berkeley National Laboratory (LBNL) for the development of fast ramp CS magnets in the ST Advanced Reactor (STAR) designed by PPPL. Characterization of the new and legacy HTS wires is presented, along with the coil design wound directly with new wire, and a Rutherford cable-wound solenoid using the legacy wire. A comparison of AC losses is made based on expected performance for a fast ramp ST operation. The HTS cable and cabled coil fabrication issues such as the heat treatment, coil mandrel design, assembly and 2212 quench characteristics are discussed, along with a test plan for CS design validation.

Zhai, Yuhu [Princeton Plasma Physics Laboratory (P↗
DOE OSTI2025

Turbulent Heating in Collisionless Low-beta Plasmas: Imbalance, Landau Damping, and Electron–Ion Energy Partition

An understanding of how turbulent energy is partitioned between ions and electrons in weakly collisional plasmas is crucial for modeling many astrophysical systems. Using theory and simulations of a four-dimensional reduced model of low-beta gyrokinetics (the “Kinetic Reduced Electron Heating Model”), we investigate the dependence of collisionless heating processes on plasma beta and imbalance (normalized cross-helicity). These parameters are important because they control the helicity barrier, the formation of which divides the parameter space into two distinct regimes with remarkably different properties. In the first, at lower beta and/or imbalance, the absence of a helicity barrier allows the cascade of injected power to proceed to small (perpendicular) scales, but its slow cascade rate makes it susceptible to significant electron Landau damping, in some cases leading to a marked steepening of the magnetic spectra on scales above the ion Larmor radius. In the second, at higher beta and/or imbalance, the helicity barrier halts the cascade, confining electron Landau damping to scales above the steep “transition-range” spectral break, resulting in dominant ion heating. We formulate quantitative models of these processes that compare well to simulations in each regime, and combine them with results of previous studies to construct a simple formula for the electron–ion heating ratio as a function of beta and imbalance. This model predicts a “winner takes all” picture of low-beta plasma heating, where a small change in the fluctuations' properties at large scales (the imbalance) can cause a sudden switch between electron and ion heating.

Interplanetary turbulence↗