Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic plasmas”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Collaborative Research: Effects of the magnetic field shear and flow shear on the kinetic physics and particle acceleration in relativistic magnetic reconnection (Final Report)

Magnetic reconnection is a fundamental plasma process that allows rapid changes of magnetic field topology and the conversion of magnetic energy into plasma kinetic energy. There has been a strong surge of interest in relativistic reconnection over the past few years in plasma astrophysics, but the rich physics of collisionless reconnection and its associated particle acceleration in the relativistic regime remain less studied compared to the non-relativistic counterparts. Magnetic reconnection in a realistic system is often accompanied with the magnetic field shear and shear flows. This setting is more general and likely the more common situation where most reconnection takes place. Within the past few decades, researchers had made significant progress in understanding how magnetic reconnection works under these conditions in the non-relativistic regime. However, how the effect of special relativity, among other effects, influences the dynamics of reconnection in the strongly magnetized astrophysical plasmas under these conditions is largely unknown. Our proposed research is to use fully kinetic simulations to understand kinetic physics and particle acceleration in relativistic magnetic reconnection in more realistic configurations and physical conditions. The primary goal of this project is to identify the fundamental effects caused by special relativity under the presence of field shear and flow shear.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Collisional stopping power of ions in warm dense matter

Here, a model for the collisional stopping of ions on free electrons in warm dense matter is developed and explored. It is based on plasma kinetic theory, but with modifications to address the warm dense matter regime. Specifically, it uses the Boltzmann-Uehling-Uhlenbeck kinetic equation to incorporate effects of Fermi degeneracy of electrons. The cross section is computed from quantum scattering of electrons and ions occurring via the potential of mean force derived from an average atom model, which incorporates effects of strong Coulomb correlations. Predictions from this model show comparable accuracy to results from time-dependent density functional theory calculations for deuterium near solid density and a temperature of several electronvolts, at a fraction of the computational cost. Further, the model captures the transition of a plasma from the classical limit to the degenerate limit, including qualitative behaviors of solid state theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Local reduced-order modeling for electrostatic plasmas by physics-informed solution manifold decomposition

Despite advancements in high-performance computing and modern numerical algorithms, computational cost remains prohibitive for multi-query kinetic plasma simulations. Here, in this work, we develop data-driven reduced-order models (ROMs) for collisionless electrostatic plasma dynamics, based on the kinetic Vlasov-Poisson equation. Our ROM approach projects the equation onto a linear subspace defined by the proper orthogonal decomposition (POD) modes. We introduce an efficient tensorial method to update the nonlinear term using a precomputed third-order tensor. We capture multiscale behavior with a minimal number of POD modes by decomposing the solution manifold into multiple time windows and creating temporally local ROMs. We consider two strategies for decomposition: one based on the physical time and the other based on the electric field energy. Applied to the 1D1V Vlasov–Poisson simulations, that is, prescribed E-field, Landau damping, and two-stream instability, we demonstrate that our ROMs accurately capture the total energy of the system both for parametric and time extrapolation cases. The temporally local ROMs are more efficient and accurate than the single ROM. In addition, in the two-stream instability case, we show that the energy-windowing reduced-order model (EW-ROM) is more efficient and accurate than the time-windowing reduced-order model (TW-ROM). With the tensorial approach, EW-ROM solves the equation approximately 90 times faster than Eulerian simulations while maintaining a maximum relative error of 7.5% for the training data and 11% for the testing data.

Electrostatic plasmas↗
DOE OSTI2021

Kinetic equilibrium reconstruction and the impact on stability analysis of KSTAR plasmas

High fidelity kinetic plasma equilibrium reconstructions are an essential requirement for accurate stability and disruption prediction analyses to support continuous operation of high beta tokamak plasmas. The present kinetic equilibrium reconstructions of plasmas in the KSTAR device include plasma density and temperature profiles from Thomson scattering and ion temperature from charge exchange spectroscopy diagnostics, and allowance for fast particle pressure. In addition, up to 25 channels of motional Stark effect diagnostic data are used to constrain the magnetic field pitch angle profile in the plasma to produce a reliable computation of the safety factor, q, profile. In this work, H-mode plasmas exhibit clear pedestal characteristics in the reconstructed pressure profile compared to internal transport barrier or L-mode plasmas. The plasma configuration and vertical position of inner strike points are validated by CCD and infrared camera images. Ideal and resistive magnetohydrodynamic (MHD) stability analyses using the DCON and resistive DCON codes utilize these kinetic equilibrium reconstructions to compare to experimental plasma stability. Equilibria with sufficiently low convergence error can provide reliable computation of ideal and resistive MHD stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Scale Separation Effects on Simulations of Plasma Turbulence

Abstract Understanding plasma turbulence requires a synthesis of experiments, observations, theory, and simulations. In the case of kinetic plasmas such as the solar wind, the lack of collisions renders the fluid closures such as viscosity meaningless and one needs to resort to higher-order fluid models or kinetic models. Typically, the computational expense in such models is managed by simulating artificial values of certain parameters such as the ratio of the Alfvén speed to the speed of light ( v A / c ) or the relative mass ratio of ions and electrons ( m i / m e ). Although, typically care is taken to use values as close as possible to realistic values within the computational constraints, these artificial values could potentially introduce unphysical effects. These unphysical effects could be significant at sub-ion scales, where kinetic effects are the most important. In this paper, we use the 10-moment fluid model in the Gkeyll framework to perform controlled numerical experiments, systematically varying the ion–electron mass ratio from a small value down to the realistic proton–electron mass ratio. We show that the unphysical mass ratio has a significant effect on the kinetic range dynamics as well as the heating of both plasma species. The dissipative process for both ions and electrons becomes more compressive in nature, although the ions remain nearly incompressible in all cases. The electrons move from being dominated by incompressive viscous-like heating/dissipation to very compressive heating/dissipation dominated by compressions/rarefactions. While the heating change is significant for the electrons, a mass ratio of m i / m e ∼ 250 captures the asymptotic behavior of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2012

Simulation Tool for Dielectric Barrier Discharge Plasma Actuators at Atmospheric and Sub-Atmospheric Pressures: SBIR Phase I Final Report

This report is the final report of a SBIR Phase I project. It is identical to the final report submitted, after some proprietary information of administrative nature has been removed. The development of a numerical simulation tool for dielectric barrier discharge (DBD) plasma actuator is reported. The objectives of the project were to analyze and predict DBD operation at wide range of ambient gas pressures. It overcomes the limitations of traditional DBD codes which are limited to low-speed applications and have weak prediction capabilities. The software tool allows DBD actuator analysis and prediction for subsonic to hypersonic flow regime. The simulation tool is based on the VORPAL code developed by Tech-X Corporation. VORPAL's capability of modeling DBD plasma actuator at low pressures (0.1 to 10 torr) using kinetic plasma modeling approach, and at moderate to atmospheric pressures (1 to 10 atm) using hydrodynamic plasma modeling approach, were demonstrated. In addition, results of experiments with pulsed+bias DBD configuration that were performed for validation purposes are reported.

Likhanskii, Alexandre↗
NASA NTRS1995

The microscopic state of the solar wind: Links between composition, velocity distributions, waves and turbulence

An overview is given of the microscopic state of the solar wind with emphasis on recent Ulysses high-latitude observations and previous Helios in-ecliptic observations. The possible links between composition, ionization state. velocity distribution functions of electrons, protons and heavy ions. kinetic plasma waves and MHD-scale turbulence are elaborated. Emphasis is placed on a connection of interplanetary kinetic-scale phenomena with their generating microscopic processes in the corona. The fast streams seem to consist of mesoscale pressure-balanced plasma filaments and magnetic flux tubes, reminiscent of the supergranular-size structures building the open corona, from which copious Alfven waves emanate. The wind from the magnetically structured and active corona shows developed compressive turbulence and considerable abundance and ionization state variations. Some modelling attempts to explain the observed element fractionation are briefly reviewed. The causes of the nonthermal particle features. such as proton-ion differential streaming, ion beams. temperature anisotropies, and skewed distributions associated with collisionless heat conduction, are ultimately to be searched in the fact, that the corona is never quiescent but fundamentally variable in space and time. Consequently, the radial evolution of the internal state of the wind resembles at all latitudes a complicated relaxation process, in the course of which the free (in comparison with LTE conditions) particle kinetic energy is converted into plasma waves and MHD turbulence on a wide range of scales. This leads to intermittent wave-particle interactions and unsteady anomalous transport, mixed with the weak effects of the rare Coulomb collisions. Spherical expansion and large-scale inhomogeneity forces the wind to attain microscopically a complex internal state of dynamic statistical equilibrium.

Marsch, E.↗
DOE OSTI2023

Ponderomotive recoil for electromagnetic waves

When waves damp or amplify on resonant particles in a plasma, nonresonant particles experience a recoil force that conserves the total momentum between particles and electromagnetic fields. This force is important to understand, as it can completely negate current drive and rotation drive mechanisms that are predicted on the basis of only resonant particles. Here, the existing electrostatic theory of this recoil force is extended to electromagnetic waves. While the result bears close similarity to historical fluid theories of laser–plasma interactions, it now incorporates both resonant and nonresonant particles, allowing momentum conservation to be self-consistently proven. Furthermore, the result is shown to be generally valid for kinetic plasmas, which is verified through single-particle hot-plasma simulations. Here, the new form of the force provides physical insight into the nature of the generalized Minkowski (plasmon) momentum of geometrical optics, which is shown to correspond to the momentum gained by the field and nonresonant particles as the wave is self-consistently ramped up from vanishing amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Plasma-assisted catalysis for ammonia synthesis in a dielectric barrier discharge reactor: key surface reaction steps and potential causes of low energy yield

Ammonia synthesis experiments were carried out in a coaxial dielectric barrier discharge (DBD) reactor packed with several different supports and metal catalysts. There was a marked increase in the reaction rate, over that obtained in an empty DBD plasma reactor, upon introduction of a packed bed of γ-Al 2 O 3 , Ru/γ-Al 2 O 3 , or SiO 2 particles. The difference in the reaction rates over γ-Al 2 O 3 and Ru/γ-Al 2 O 3 was minimal. Complementary zero-dimensional plasma kinetic model analysis was also performed using inputs from experimental data. This kinetic analysis allowed for gas phase reactions, Eley–Rideal (E–R) reactions, and direct adsorption of radical species on the γ-Al 2 O 3 surface. On the metal surface, dissociative adsorption of N 2 and H 2 , and Langmuir–Hinshelwood reactions were also included. This analysis revealed that, under the conditions of our experiments, ammonia synthesis proceeds principally by the formation of reactive radicals in the gas phase, which then adsorb and participate in E–R reactions on both the metal and support material surfaces. Furthermore, this finding illustrates a challenge for substantially increasing the energy yield for plasma-assisted ammonia synthesis in typical DBD reactors containing packed catalyst beads.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines

Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.

Singh, Nagendra↗
NASA NTRS1993

Magnetic reconnection with large separatrix angles

The magnetic reconnection process is studied here using incompressible MHD simulations with different inflow boundary conditions and different magnetic Reynolds numbers R(m). The angle between the magnetic separatrices is in steady state reconnection depends mainly on the normal magnetic field on the inflow boundary. In steady state nonuniform reconnection with large separatrix angles, field-aligned plasma jets appear slightly downstream of the magnetic separatrices. The field-aligned plasma jet are stronger when R(m) is larger. Each field-aligned plasma jet consists of two parts: a slow shock and a fast-mode compressional wave. The slow shock converts the magnetic energy into plasma kinetic energy by acceleration and heating. The fast-mode compressional wave decelerates the plasma to a smaller outflow speed and heats it further. Nearly all the magnetic energy flowing into the diffusion region is converted into other forms. The length and width of the diffusion region depend on the values of the reconnection rate, R(m), and the normal magnetic field on the inflow boundary.

Yan, M.↗
NASA NTRS2012

Plasma Wake Simulations and Object Charging in a Shadowed Lunar Crater During a Solar Storm

Within a permanently shadowed lunar crater the horizontal flow of solar wind is obstructed by upstream topography, forming a plasma wake that electrostatically diverts ions toward the crater floor and generates a surface potential that can reach kilovolts. In the present work kinetic plasma simulations are employed to investigate the morphology of a lunar crater wake during passage of a solar storm. Results are cast in terms of leading dimensionless ratios including the ion Mach number, ratio of crater depth to plasma Debye length, peak secondary electron yield, and electron temperature vs. electron impact energy at peak secondary yield. This small set of ratios allows generalization to a much wider range of scenarios. The kinetic simulation results are fed forward into an equivalent-circuit model of a roving astronaut. In very low-plasma-current environments triboelectric charging of the astronaut suit becomes effectively perpetual, representing a critical engineering concern for roving within shadowed lunar regions. Finally, simulated ion fluxes are used to explore sputtering and implantation processes within an idealized crater. It is suggested that the physics of plasma mini-wakes formed in the vicinity of permanently shadowed topography may play a critical role in modulating the enigmatic spatial distribution of volatiles at the lunar poles.

Zimmerman, Michael I.↗
NASA NTRS1973

Physics of planetary ionospheres

The fundamental physical and chemical processes in an idealized planetary ionosphere are considered as a general abstraction, with actual planetary ionospheres representing special cases. After describing the structure of the neutral atmospheres (the barosphere, the thermosphere, and the exosphere) and noting the principal ionizing radiations responsible for the formation of planetary ionospheres, a detailed study is made of the thermal structure of these ionospheres and of the chemical processes and plasma-transport processes occurring in them. The features of equilibrium and realistic models of planetary ionospheres are discussed, and an attempt is made to determine the extent of these ionospheres. Considering the ionosphere as a plasma, a plasma kinetic approach is developed for determining the effects of interactions between individual particles and waves in this plasma. The use of remote-sensing radio techniques and direct measurement or in situ techniques is discussed. Finally, the observed properties of the ionospheres of the Earth, Mars, Venus, and Jupiter are reviewed.

Bauer, S. J.↗
DOE OSTI2023

Electron Injection via Modified Diffusive Shock Acceleration in High-Mach-number Collisionless Shocks

The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Radiative properties of a plasma moving across a magnetic field. I - Theoretical analysis. II - Numerical results

The early-time evolution of plasmas moving across a background magnetic field is addressed with a 2D model in which a plasma cloud is assumed to have formed instantaneously with a velocity across a uniform background magnetic field and with a Gaussian density profile in the two dimensions perpendicular to the direction of motion. This model treats both the dynamics associated with the formation of a polarization field and the generation and propagation of electromagnetic waves. In general, the results indicate that, to zeroth order, the plasma cloud behaves like a large dipole antenna oriented in the direction of the polarization field which oscillates at frequencies defined by the normal mode of the system. Radiation damping is shown to play an important role in defining the plasma cloud evolution, causing a rapid decay of the polarizaiton field and a loss of plasma kinetic energy and momentum on time scales comprable to several ion gyroperiods. Scaling laws are derived for the plasma momentum and energy loss rates, and predictions for the braking time, the amplitude and spectrum of the radiation field, and the total radiated power are presented for conditions relevant to the recent Combined Release and Radiation Effects Satellite experiments.

Roussel-Dupre, Robert↗
NASA NTRS1979

A fluid description of plasma double-layers

The space-charge double-layer that forms between two plasmas with different densities and thermal energies was investigated using three progressively realistic models which are treated by fluid theory, and take into account four species of particles: electrons and ions reflected by the double-layer, and electrons and ions transmitted through it. The two plasmas are assumed to be cold, and the self-consistent potential, electric field and space-charge distributions within the double-layer are determined. The effects of thermal velocities are taken into account for the reflected particles, and the modifications to the cold plasma solutions are established. Further modifications due to thermal velocities of the transmitted particles are examined. The applicability of a one dimensional fluid description, rather than plasma kinetic theory, is discussed. Theoretical predictions are compared with double layer potentials and lengths deduced from laboratory and space plasma experiments.

Levine, J. S.↗