Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space 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.

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

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

Schoeffler, K. M. (ORCID:0000000276447381)↗
DOE OSTI2024

Understanding extraction limits of plasma cathodes with experiment and simulation

The project is focused on computations for enhanced ionization near the exit orifice of a plasma cathode and how the cathode extraction electrode geometry and stray magnetic field structure influence available extractable current. The computational effort will take place at the Princeton Collaborative Research Laboratory (PCRF). The specific modeling platform to be utilized from PCRF is the electrostatic direct implicit particle in cell (EDIPIC) code. EDIPIC will be used to study ionization processes in the electron extraction region of the cathode. The control of electron energy is a key aspect to minimizing power losses in the plasma. Exiting electrons ionize gas leaving the orifice and in the process produce the dense plasma from which large amounts of current can be extracted. Project activities include studying processes that determine the energy distribution of the exiting electrons thus providing insight into how to optimize the source. The model will be validated and developed using experimental measurements as warranted. The familiarity and expertise of the team at PPPL with this type of plasma and the challenges posed greatly enhances the likelihood of gaining insight into the extraction sheath under the duration of the proposed activity. The effort is expected to shed a great deal of insight into electron extraction processes and pathways for optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Saturation of Fishbone Instability by Self-Generated Zonal Flows in Tokamak Plasmas

Gyrokinetic simulations of the fishbone instability in DIII-D tokamak plasmas find that self-generated zonal flows can dominate the nonlinear saturation by preventing coherent structures from persisting or drifting in the energetic particle phase space when the mode frequency down-chirps. Results from the simulation with zonal flows agree quantitatively, for the first time, with experimental measurements of the fishbone saturation amplitude and energetic particle transport. Moreover, the fishbone-induced zonal flows are likely responsible for the formation of an internal transport barrier that was observed after fishbone bursts in this DIII-D experiment. Finally, gyrokinetic simulations of a related ITER baseline scenario show that the fishbone induces insignificant energetic particle redistribution and may enable high performance scenarios in ITER burning plasma experiments. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗
DOE OSTI2022

Oxidation in laser-generated metal plumes

The temporal evolution of atoms and molecules in a laser-produced plasma was investigated using optical emission spectroscopy for several metal targets (i.e., Al, Ti, Fe, Zr, Nb, and Ta). Plasmas from metal targets were generated by focusing 1064 nm, 6 ns pulses from an Nd:YAG laser. Gas-phase oxidation/plasma chemistry was initiated by adding O2 (partial pressures up to ≈20%) to an N2 environment where the total background pressure was kept at a constant 1 atmosphere. Temporally resolved emission spectral features were used to track the gas-phase oxidation. The dynamics of atomic and molecular species were monitored using space-resolved time-of-flight emission spectroscopy. Our results highlight that the partial pressure of O2 strongly influences spectral features and molecular formation in laser-produced plasmas. Atoms and molecules co-exist in plasmas, although with different temporal histories depending on the target material due to differences in thermo- and plasma chemical reactions occurring in the plume.

Kautz, Elizabeth J. (ORCID:0000000263389223)↗
DOE OSTI2023

Modeling of the cold electron plasma density for radiation belt physics

This review focusses strictly on existing plasma density models, including ionospheric source models, empirical density models, physics-based and machine-learning density models. This review is framed in the context of radiation belt physics and space weather codes. The review is limited to the most commonly used models or to models recently developed and promising. A great variety of conditions is considered such as the magnetic local time variation, geomagnetic conditions, ionospheric source regions, radial and latitudinal dependence, and collisional vs. collisionless conditions. These models can serve to complement satellite observations of the electron plasma density when data are lacking, are for most of them commonly used in radiation belt physics simulations and can improve our understanding of the plasmasphere dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Understanding the limitations and potential of micro tensile testing of tungsten and needs for crystal plasticity modelling

Tungsten is the leading plasma facing material candidate due to its exceptional properties. Understanding the response to neutron irradiation is crucial for the lifetime evaluation of tungsten. The limited space in nuclear reactors and the high levels of radioactivity of the specimens after irradiation are significant barriers to accurately measuring the mechanical properties after neutron irradiation. Testing of micro tensile specimens is one approach to reduce the total amount of irradiated material needed for a set of mechanical testing experiments. However, micro samples are not necessarily measuring the bulk material properties, as size effects produce higher measured mechanical properties than is observed in engineering size counterparts. Determining the minimal sample size for reliable bulk property measurement, combined with modeling efforts, is essential. We conducted room temperature tensile tests on tungsten micro tensile specimens fabricated with focused ion beam, plasma focused ion beam and femto-second laser ablation system to dimension of 2×2×7?µm³, 7×7×18?µm³, and 80×100×233?µm³ (width×thickness×gauge-length), respectively. Elevated tensile testing was performed up to 475°C on 5×5×18?µm³ specimens fabricated with focused ion beam. The smallest specimens exhibited a high degree of ductility and strength, whereas the largest specimens demonstrated behavior akin to bulk tungsten. To investigate the effects introduced by the micro specimen fabrication processes and to obtain the necessary bulk dislocation density for the crystal plasticity model, we employed X-ray diffraction depth profiling. This measurement was performed using different X-ray sources with varying penetration depths. Finally, the potential and limitations of micro mechanical tests for tungsten will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

Generation of waves in magnetized plasma [PowerPoint]

The presentation goes through assumptions made for the study, sections on plasma dispersion, whistler modes for RBR, excitation of plane waves with antenna, dipole vs loop antenna, verification of results, excitation of whistler waves with electron beams, assumptions for the model, distribution of current in space in general case, Fourier transform of current distribution, resonant excitation of waves, and ending with a discussion of generic debunching.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A self-consistent Hamiltonian model of the ponderomotive force and its structure preserving discretization

In the presence of an inhomogeneous oscillatory electric field, charged particles experience a net force, averaged over the oscillatory timescale, known as the ponderomotive force. We derive a one-dimensional Hamiltonian model which self-consistently couples the electromagnetic field to a plasma which experiences the ponderomotive force. We derive a family of structure preserving discretizations of the model of varying order in space and time using conforming and broken finite element exterior calculus spectral element methods. In all variants of our discretization framework, the method is found to conserve the Casimir invariants of the continuous model to machine precision and the energy to the order of the splitting method used.

Physics↗
DOE OSTI2021

Optical shock-enhanced self-photon acceleration

Photon accelerators can spectrally broaden laser pulses with high efficiency in moving electron density gradients driven in a rapidly ionizing plasma. When driven by a conventional laser pulse, the group velocity walk-off experienced by the accelerated photons and deterioration of the gradient from diffraction and plasma-refraction limit the extent of spectral broadening. In this work, we show that a laser pulse with a shaped space-time and transverse intensity profile overcomes these limitations by creating a guiding density pro le at a tunable velocity. Self-photon acceleration in this profile leads to dramatic spectral broadening and intensity steepening, forming an optical shock that further enhances the rate of spectral broadening. In this new regime, multi-octave spectra extending from 400nm to 60nm wavelengths, which support near-transform limited < 400 as pulses, are generated over < 100 μm of interaction length.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2026

Verification of global gyrokinetic simulation of low frequency mode excited by thermal plasma in spherical tokamak

A low frequency Alfven eigenmode in the ST40 spherical tokamak (ST) is found to be resonantly excited by thermal ions from global gyrokinetic GTC simulations. The analyses of phase-space resonances and wave-particle energy exchanges show that the Beta-induced Alfvén-Acoustic Eigenmode (BAAE) is primarily driven unstable by thermal ions, but is stabilized by energetic particles (EPs). When all kinetic effects are artificially suppressed, GTC simulations in the ideal magnetohydrodynamic (MHD) limit confirm the existence of the MHD eigenmode with a frequency in the BAAEs gap, as previously found by the ideal MHD code NOVA. The kinetic effects of thermal ions induce a frequency downshift compared to the stable BAAE in the MHD simulations. These results support the verification for global gyrokinetic simulations of the BAAE in the ST. However, the simulations do not explain the experimental observations of the BAAE in the ST40, possibly owing to the limitation of the collisionless gyrokinetic Monte Carlo model for the EPs.

Alfvén eigenmode↗
DOE OSTI2024

Power handling in a highly-radiative negative triangularity pilot plant

Abstract This work explores detailed power handling solutions for a class of high-field, highly-radiative negative triangularity (NT) reactors based around the MANTA concept (Rutherford et al 2024 Plasma Phys. Control. Fusion ). The divertor design is kept as simple as possible, opting for a standard divertor with standard leg length. FreeGS is used to create an equilibrium for the boundary region, prioritizing a short outer leg length of only ∼50 cm (∼40% of the minor radius). The UEDGE code package is used for the boundary plasma solution to track plasma temperatures and fluxes to the divertor targets. It is found that for P SOL = 25 MW and n sep = 0.96 × 10 20 m −3 , conditions consistent with initial core transport modeling, little additional power mitigation is necessary. For a fixed impurity fraction of just 0.13% Ne in the plasma, the peak heat flux density at the more heavily loaded outer targets falls to 7.8 MW m −2 , while the electron temperature T e remains just under 5 eV. Scans around the parameter space reveal that even at densities lower than in the primary operating scenario, P SOL can be increased up to 50 MW, so long as a slightly higher fraction of extrinsic radiator is used. With less than 1% neon (Ne) impurity content, the divertor still experiences less than 10 MW m −2 at the outer target. Design of the plasma-facing components includes a close-fitting vacuum vessel with a tungsten inner surface as well as FLiBe-carrying cooling channels fashioned into the VV wall directly behind the divertor targets. For the seeded heat flux profile, Ansys Fluent heat transfer simulations estimate that the outer target temperature remains at just below 1550 ∘ C. Initial scoping of advanced divertor designs shows that for an X-divertor, detachment of the outer target becomes much simpler, and plasma fluxes to the targets drop considerably with only 0.01% Ne content.

Miller, M. A. (ORCID:0000000265406533)↗
DOE OSTI2023

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

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 OSTI2021

Production and Persistence of Extreme Two-temperature Plasmas in Radiative Relativistic Turbulence

Turbulence is a predominant process for energizing electrons and ions in collisionless astrophysical plasmas, and thus is responsible for shaping their radiative signatures (luminosity, spectra, and variability). To better understand the kinetic properties of a collisionless radiative plasma subject to externally driven turbulence, we investigate particle-in-cell simulations of relativistic plasma turbulence with external inverse Compton cooling acting on the electrons. We find that ions continuously heat up while electrons gradually cool down (due to the net effect of radiation), and hence the ion-to-electron temperature ratio T i /T e grows in time. Furthermore, we show that Ti/Te is limited only by the size and duration of the simulations (reaching ${T}_{i}/{T}_{e}\sim {10}^{3}$), indicating that there are no efficient collisionless mechanisms of electron–ion thermal coupling. This result has implications for models of radiatively inefficient accretion flows, such as observed in the Galactic center and in M87, for which so-called two-temperature plasmas with ${T}_{i}/{T}_{e}\gg 1$ have been invoked to explain their low luminosity. Additionally, we find that electrons acquire a quasi-thermal distribution (dictated by the competition of turbulent particle energization and radiative cooling), while ions undergo efficient nonthermal acceleration (acquiring a harder distribution than in equivalent nonradiative simulations). There is a modest nonthermal population of high-energy electrons that are beamed intermittently in space, time, and direction; these beamed electrons may explain rapid flares in certain high-energy astrophysical systems (e.g., in the Galactic center). These numerical results demonstrate that extreme two-temperature plasmas can be produced and maintained by relativistic radiative turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗