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

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗
DOE OSTI2023

The application of parallel kinetic simulations to laser and electron transport through plasmas (Final technical report)

This is a final report for the grant entitled, “The application of parallel kinetic simulations to laser and electron transport through plasmas”. The objectives of this grant were to significantly advance the fundamental understanding of the nonlinear optics of plasmas and electron transport in high-energy-density laboratory plasmas (HEDLP), including conditions of relevance to Inertial Fusion Energy (IFE). The ultimate goal was to use the understanding to determine how to fully control laser plasma interactions. The primary research tools were our own kinetic particle-in-cell software, OSIRIS, that includes kinetic physics and can run effectively on leadership class computing facilities. Therefore, one objective was to ensure that OSIRIS in continually improved so that it was more accurate and could effectively utilize state-of-the-art computing facilities. Another objective was to attract and train young researchers into the field of high energy density plasma physics. To meet the research objectives, the funds from this proposal were used to conduct research on stimulated Raman scattering (SRS) and enhance our PIC software OSIRIS. It was found that small normalized magnetic fields can in some cases mitigate SRS and that speckles can mutually interact through SRS. It was also found that it is possible for instabilities drive near quarter critical (e.g., the high frequency hybrid instability-HFHI) can generate light waves that propagate back down a density gradient where they can rescatter into the HFHI at 1/16 of the original quarter critical density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗
DOE OSTI2024

Kinetic corrections to heat-flow and Nernst advection for laser heated plasmas

Reduced models for approximating the impact of kinetic electron behavior on the transport of thermal energy and magnetic field are investigated. The thermal flux limiter has improved agreement with Vlasov–Fokker–Planck data when a harmonic form is used that adjusts the electron mean free path to account for electron–electron collisions; these results apply to both unmagnetized and magnetized plasmas. Once a magnetic field is incorporated, the mean free path should also be modified using the electron gyroradius. A flux limiter on Nernst advection of magnetic fields is also required; a form that limits Nernst by the same fraction as the thermal heat-flow best reproduces kinetic simulations. A flux limiter form for the cross terms (Righi–Leduc and cross-gradient-Nernst) is also suggested. Hohlraum simulations relevant to fusion experiments on the National Ignition Facility are found to be sensitive to all of these details.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

1D kinetic study of pinch formation in a dense plasma focus: Transition from collisional to collisionless regimes

Here, the pinch-formation stage of a deuterium dense plasma focus, and associated “shock-flash” neutron yield, is studied using 1D kinetic simulations considering a plasma column with initial pressure P, initial radius R, and the compression to be driven by a constant current I. The relative behavior of the compression is shown to be similar for fixed ratios of the characteristic ion mean free path to the radius of the plasma column at stagnation, $λ_{st}/R_{st}$. This dimensionless parameter is shown to scale like $I^4/(P^3R^5)$. The compression ratio, $R/R_{st}$, is found to be a minimum when $λ_{st}/R_{st} ≈$ 1 and is the largest in the collisionless limit where $λ_{st} \gg R_{st}$. This behavior is in contrast to the analogous planar pinch where R/Rst decreases from one constant for λst/Rst$\ll$1 to a smaller constant for $λ_{st}/R_{st} \gg$ 1. The yield in the collisionless regime is shown to fall between the two well-known $I^4$ scaling laws. Furthermore, this regime exhibits qualities that potentially make it appealing for radiography applications, such as increased localization in time and space of the neutron formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Interplanetary gas dynamics.

Observations of the interplanetary medium are considered together with the structure of the solar atmosphere and the kinetic properties of the interplanetary plasma. Methods of gas kinetics and continuum flows are examined, taking into account microscopic and macroscopic representations, the kinetic flow equations, the formal solution of the Vlasov equation, and the continuum flow equations. The collective particle behavior of that interplanetary gas is discussed along with the hydrodynamic coronal expansion, a wave-pump problem, the free expansion phenomenon, and a generalized free expansion problem.

Liu, V. C.↗
NASA NTRS2011

The Nonlinear Magnetosphere: Expressions in MHD and in Kinetic Models

Like most plasma systems, the magnetosphere of the Earth is governed by nonlinear dynamic evolution equations. The impact of nonlinearities ranges from large scales, where overall dynamics features are exhibiting nonlinear behavior, to small scale, kinetic, processes, where nonlinear behavior governs, among others, energy conversion and dissipation. In this talk we present a select set of examples of such behavior, with a specific emphasis on how nonlinear effects manifest themselves in MHD and in kinetic models of magnetospheric plasma dynamics.

Hesse, Michael↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Kinetic effects of self-organization in low-temperature magnetized plasmas (Final Report for DE-SC0019045: 9/1/2018-8/31/2019)

The research goal of this Early Career Research Program project is to advance the science of low temperature magnetized plasmas using first-principles kinetic, e.g., solving the Boltzmann equation, models. A better understanding of the plasma transport properties, such as electrical resistivity and thermal conductivity must be obtained to better control and design magnetized plasma sources. A key physical phenomenon in low-temperature magnetized plasmas is the self-organization of plasma flows, e.g., coherent plasma structures and transition between different discharge modes. Understanding the physics of coherent plasma structures has been identified as one of the four important frontier topics in a recent workshop by the Office of Fusion Energy Science. These self-emerging plasma patterns results from the nonlinear coupling between plasma constituents (ions, electrons, and neutral atoms), electromagnetic fields, and materials. The physics of low-temperature magnetized plasmas is particularly complex because (i) the collisionless instabilities affect collisional phenomena (ionization, transport, etc.) and vice versa and (ii) the presence of materials (sheath, electron injection, etc.) influences the bulk plasma properties in low-temperature plasma (LTP) devices. Computational models of these plasma flows remain challenging because the plasma density can vary several orders of magnitudes, a wide range of temporal and spatial scales must be resolved, the phenomena are inherently three-dimensional, and various physical processes, such as wall interaction, instabilities, inelastic collisions, etc., must be simultaneously and self-consistently taken into account. Understanding the fundamental transport mechanisms in cross-field configurations enables additional controllability of plasma properties and electron energy distribution functions (EEDFs).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effective viscosity, resistivity, and Reynolds number in weakly collisional plasma turbulence

ABSTRACT We examine dissipation and energy conversion in weakly collisional plasma turbulence, employing in situ observations from the Magnetospheric Multiscale mission and kinetic particle-in-cell simulations of proton–electron plasma. A previous result indicated the presence of viscous-like and resistive-like scaling of average energy conversion rates – analogous to scalings characteristic of collisional systems. This allows for extraction of collisional-like coefficients of effective viscosity and resistivity, and thus also determination of effective Reynolds numbers based on these coefficients. The effective Reynolds number, as a measure of the available bandwidth for turbulence to populate various scales, links turbulence macroscale properties with kinetic plasma properties in a novel way.

Astronomy & Astrophysics↗
DOE OSTI2023

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2016

Investigating Global Ion and Neutral Atom Populations with IBEX and Voyager

The main objective of this project was to investigate pickup ion (PUI) production in the solar wind and heliosheath (the region between the termination shock and the heliopause) and compute the distributed energetic neutral atom fluxes throughout the helioshpere. The simulations were constrained by comparing the model output against observations from Ulysses, New Horizons, Voyager 1 and 2, and IBEX space probes. As evidenced by the number of peer reviewed journal publications resulting from the project (13 plus three submitted) and their citation rate (156 citations over three years), the project has made a lasting contribution to the field. The outcome is a significant improvement of our understanding of the pickup ion production and distribution in the distant heliosphere. The team has accomplished the entire set of tasks A-H set forth in the proposal. Namely, the transport modeling framework has been augmented with two populations of pickup ions (PUIs), the boundary conditions for the plasma and interstellar neutral hydrogen were verified against Ulysses and New Horizons PUI and an optimal set of velocity diffusion parameters established. The multi-component fluxes of PUIs were computed and isotropic velocity distributions generated for each cell in the computer simulation that covered the heliosphere from 1.5 AU to the heliopause. The distributions were carefully compared with in situ measurements at 3 AU (Ulysses), 12 AU (New Horizons), and 80-90 AU (Voyager 1 and 2) as well as those inferred from ENA fluxes measured by Cassini and IBEX (Wu et al., 2016). Some examples of modeldata comparison are shown in Figure 1. We have used coupled MHD-plasma and kinetic-neutral code to investigate the likely range of plasma and magnetic field parameters in the local interstellar medium (LISM), based on the assumption that the shape of the IBEX ribbon could be used to determine the orientation of the interstellar magnetic field. While the magnetic field is believed to be oriented toward the center of the ribbon, constraining its strength requires comparing the model-predicted angular diameter and circularity of the ribbon with the observations. The study, published in Heerikhuisen et al. (2014), found that the most likely range for the LISM magnetic field strength is between 0.2 and 0.3 nT, which is less than previously thought. Figure 2 shows the IBEX data (left) and compares it to the simulation with a 0.2 nT interstellar magnetic field (center) and a 0.4 nT (right).

Florinski, Vladimir↗
DOE OSTI2024

Anticathode effect on electron kinetics in electron beam generated E × B plasma

Electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E x B) fields are promising for low damage processing of materials with applications to microelectronics and quantum information systems. In cylindrical e-beam E x B plasmas, radial confinement of electrons and ions is achieved by an axial magnetic field and radial electric field, respectively. To control the axial confinement of electrons, such e-beam generated plasma sources may incorporate a conducting boundary known as an anticathode, which is placed on the axially opposite side of the plasma from the cathode. In this work, it is shown that varying the anticathode voltage bias can control the degree to which the anticathode collects or repels incident electrons, allowing control of warm electron (electron energies in 10–30 eV range) and beam electron population confinement. It is suggested that the effect of the anticathode bias on the formation of these distinct electron populations is also associated with the transition between weak turbulence and strong Langmuir turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Linear simulation of kinetic electromagnetic instabilities in a tokamak plasma with weak magnetic shear

Gyrokinetic simulation and eigenvalue calculation of electromagnetic instabilities are carried out for an experimentally observed low-n mode in weak magnetic shear discharge. With different magnetic shear values, the Gyrokinetic Toroidal Code simulation of the ion temperature gradient mode is consistent with the eigenvalue code (HD7) calculation. Due to the sensitivity of the kinetic ballooning mode (KBM) to global equilibrium, the simulation of the KBM deviates from the eigenvalue results, for the ballooning representation used in HD7 satisfies its spatial scale separations. Under a flat safety factor profile, the KBM is more unstable and its mode structure tends to move with the peak of the ion temperature drive. Further simulation of the KBM in an HL-2A-like equilibrium shows that the β excitation threshold of the mode is lower than 0.2% and the dominant toroidal mode number is n = 4, which is consistent with the measured experimental spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1971

Theory of the unmagnetized plasma.

The Vlasov mathematical model of a plasma, which has come to be thought more useful than any other in describing the dynamical behavior of the majority of plasmas of interest, is first examined. Macroscopic variables and moment equations; linear electrostatics solutions; plasma oscillations, ion acoustic waves, and linear instabilities are treated, as well as external fields, 'test' charges, and nonlinear Vlasov phenomena. Plasmas are statistically described, and attention is given to the kinetic theory of the stable, uniform plasma and the Balescu-Lenard equation; two-time ensemble averages and fluctuation spectra in stable plasmas; the kinetic theory of the unstable plasma; and ensembles of Vlasov plasmas. Some illustrative experiments are described. Four appendixes deal with the electrostatic approximation and transverse waves; solution of the linearized Vlasov equation in a magnetic field; estimates of correlation functions from thermal equilibrium; and equivalence of spatially uniform BBGKY and Klimontovich correlations.

Montgomery, D. C.↗
DOE OSTI2024

Unraveling Nonequilibrium Generation of Atomic Nitrogen and Hydrogen in Plasma-Aided Ammonia Synthesis

In this study, nonequilibrium generation of atomic nitrogen and hydrogen governing ammonia production both in the gas phase and on the catalyst surface is critical to plasma-aided ammonia synthesis. Here, this work studies the nonequilibrium generation of atomic nitrogen and hydrogen with a focus on the kinetic role of vibrational energy transfer of hydrogen molecules in plasma-aided ammonia synthesis. By combining two-photon absorption laser-induced fluorescence measurements and plasma kinetic modeling, we found that plasma not only generates ammonia but also produces critical H, N, and NH radicals via both electron impact and vibrational energy transfer excitations. The vibrational energy transfer from the excited hydrogen H 2 (v = 1) to higher vibrational levels H 2 (v = 2–3) via the V–V exchange (H 2 (v)–H 2 (v)) and V–V' exchange (N 2 (v)–H 2 (v)) significantly enhances the H and NH production and then promotes the coupling between N and NH for ammonia synthesis both in the gas phase and on the catalyst surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗