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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Polar Spacecraft Observations of the Turbulent Outer Cusp/Magnetopause Boundary Layer of Earth

The orbit of the Polar spacecraft has been ideally suited for studying the turbulent region of the cusp that is located near or just outside the magnetopause current sheet at 7-9 R(sub E). The wave data obtained in this region show that electromagnetic turbulence is dominant in the frequency range 1-10 Hz. The waves responsible for this turbulence usually propagate perpendicular to the local magnetic field and have an index of refraction that generally falls between the estimated cold plasma theoretical values of the electromagnetic lower hybrid and whistler modes and may be composed of both modes in concert with kinetic Alfven waves and/or fast magnetosonic waves. Fourier spectra of the higher frequency wave data also show the electromagnetic turbulence at frequencies up to and near the electron cyclotron frequency. This higher frequency electromagnetic turbulence is most likely associated with whistler mode waves. The lower hybrid drift and current gradient instabilities are suggested as possible mechanisms for producing the turbulence. The plasma and field environment of this turbulent region is examined and found to be extremely complex. Some of the wave activity is associated with processes occurring locally, such as changes in the DC magnetic field, while others are associated with solar wind and interplanetary magnetic field changes.

Pickett, J. S.↗
DOE OSTI2022

Theory, observations, and simulations of kinetic entropy in a magnetotail electron diffusion region

In this work, we examine velocity-space kinetic entropy, a spatially local measure of entropy for systems out of thermal equilibrium, during an encounter of an electron diffusion region at a magnetic reconnection site in Earth’s magnetotail by the Magnetospheric Multiscale (MMS) mission. We start by generalizing the theory of kinetic entropy to the case of non-uniform velocity space grids and transforming the equations into spheri- cal energy coordinates useful to experimental plasma detectors. The theory is then applied to MMS data and compared to particle-in-cell simulations of reconnection. We demonstrate that the entropy-based non-Maxwellianity measure from the MMS data is of sufficiently high precision to reliably identify non-Maxwellian distributions and therefore the measurements when kinetic effects are most significant. By comparing two different non-Maxwellian measures, we show that total entropy density suffers from “information loss” because it lacks a dependence on the velocity space grid, and so has lost information about how well a distribution function is resolved. Local velocity-space kinetic entropy density recovers this information. We quantify information loss and argue that the considerations needed to minimize it are crucial for instruments designed to measure distribution functions in situ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗
NASA NTRS1990

Cyclotron maser and plasma wave growth in magnetic loops

Cyclotron maser and plasma wave growth which results from electrons accelerated in magnetic loops are studied. The evolution of the accelerated electron distribution is determined by solving the kinetic equation including Coulomb collisions and magnetic convergence. It is found that for modest values of the column depth of the loop the growth rates of instabilities are significantly reduced and that the reduction is much larger for the cyclotron modes than for the plasma wave modes. The large decrease in the growth rate with column depth suggests that solar coronal densities must be much lower than commonly accepted in order for the cyclotron maser to operate. The density depletion has to be similar to that which occurs during auroral kilometric radiation events in the magnetosphere. The resulting distributions are much more complicated than the idealized distributions used in many theoretical studies, but the fastest growing mode can still simply be determined by the ratio of electron plasma to gyrofrequency, U=omega(sub p)/Omega(sub e). However, the dominant modes are different than for the idealized situations with growth of the z-mode largest for U approximately less than 0.5, and second harmonic x-mode (s=2) or fundamental o-mode (s=1) the dominant modes for 0.5 approximately less than U approximately less than 1. The electron distributions typically contain more than one inverted feature which could give rise to wave growth. It is shown that this can result in simultaneous amplification of more than one mode with each mode driven by a different feature and can be observed, for example, by differences in the rise times of the right and left circularly polarized components of the associated spike bursts.

Hamilton, Russell J.↗
DOE OSTI2025

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nonlinear gyrokinetic predictions of SPARC burning plasma profiles enabled by surrogate modeling

Multi-channel, nonlinear predictions of core temperature and density profiles are performed for the SPARC tokamak accounting for both kinetic neoclassical and fully nonlinear gyro-kinetic turbulent fluxes. A series of flux-tube, nonlinear, electromagnetic simulations using the CGYRO code with six gyrokinetic species are coupled to a nonlinear optimizer using Gaussian process regression techniques. The simultaneous evolution of energy sources, including alpha heat, radiation, and energy exchange, coupled with these high fidelity models and techniques, leads to a converged solution in electron temperature, ion temperature and electron density channels with a minimal number of expensive gyrokinetic simulations without compromising accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Acceleration of Oxygen Ions In Dipolarization Events: 2. PSBL Distributions

Abstract This paper represents the second part of an investigation of the acceleration of energetic oxygen ions from encounters with a dipolarization front (DF), based on test particle tracing in the fields of an MHD simulation. In this paper, we focus on distributions in the plasma sheet boundary layer (PSBL). O + beams close to the plasma sheet boundary are found to be less pronounced and/or delayed against the H + beams. The reason is that these particles are accelerated by nonadiabatic motion in the duskward electric field such that O + ions gain the same amount of energy, but only 1/4 of the speed of protons. This causes a delay and larger equatorward displacement by the E × B drift. In contrast, the O + beams somewhat deeper inside the plasma sheet, where previously multiple proton beams were found, are accelerated at an earthward propagating DF just like H + , forming a field‐aligned beam at a similar speed as the lowest‐energy H + beam. We found that the source location depends on the adiabaticity of the orbit. For larger adiabaticity, the beam ions originate initially from the outer plasma sheet, but later from the opposite PSBL or lobe, but for low adiabaticity, sources are well inside the plasma sheet. The energy gained from a single encounter of a DF is comparable to the kinetic energy associated with the front speed. Assuming maximum speeds of 500–1,000 km/s, this yields a mass dependent acceleration of about 1–5 keV for protons and 20–80 keV for oxygen ions, independent of their charge state.

Birn, J.↗
DOE OSTI2021

Features of ion generation by a picosecond laser in the range of 10 11 –10 13 W cm -2 power densities

Picosecond lasers (ps-lasers) have significant advantages for the generation of low charge state ions compared to nanosecond lasers because the influence of heat conductivity on a solid target is almost negligible in the case of ps-laser ablation for laser pulse durations less than 10 ps. However, there is no comprehensive data on ion yields for different elements and target irradiation conditions for laser power densities at the target surface around and below 10 13 W cm –2 , which is of interest to our study of such plasmas as a source of low charge state ions for various applications, particularly for external injection of those ions into an Electron Beam Ion Source (EBIS). We investigated ion generation from Al, Ti, Cu, Nb and Ta target elements by a ps-laser with power densities in the range of 10 11 –10 13 W cm –2 at the target surface. A ps-laser with 1.27 mJ maximum energy within an 8 ps pulse and repetition rate up to 400 Hz has been used to generate a laser-ablated plasma. Dependencies of ion current versus time, total charge of registered ions as well as ion kinetic energy distributions are characterized using a Faraday cup. Significant difference in ion current dynamics between first, second and following shots onto the same target spot was found for all five target elements. The total charge of ions registered by the Faraday cup increases linearly with increasing laser pulse energy and is almost independent of the target element and number of shots onto the same target spot for all five target elements studied. The results obtained give us a basis for specification and design of the source of low charge state ions for external injection into EBIS.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Study of Ion Generation by a ps-laser for External Injection into Electron Beam Ion Source

We investigated ion generation from Al, Ti, Cu, Nb and Ta target elements by a picosecond-laser in the range of 10 11 – 10 13 W/cm 2 power densities at the target surface. A ps-laser with 1.27 mJ maximum energy within a 8 ps pulse and repetition rate up to 400 Hz has been used to generate a laser-ablated plasma. Dependencies of ion current vs time, total charge of registered ions as well as ion kinetic energy distributions are characterized using a Faraday cup. Significant difference in ion current dynamics between first, second and following shots onto the same target spot was found for all five target elements. The total charge of ions registered by the Faraday cup increase linearly with increasing laser pulse energy and are almost independent on the target element and number of shots onto the same target spot for all five target elements studied. Target craters have been examined using a commercially available microscope (Nikon Eclipse LV150) and crater cross section areas were found to be in the range of 0.0076 – 0.0172 mm 2 within the first five shots onto the same target spot for different target elements. A linear dependence of crater cross section square on laser pulse energy has been observed for the Ta target. Ion energy distributions for all target elements have a plateau or a slight peak in the energy range of 10 - 100 eV and long low-intensity energy tail up to tens of keV.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Krash-ESL-Final report

During the period of the current Grant, the effort of the UCSD researchers was focused on: i) The application of the gyrokinetic code COGENT code to study the kinetic drift wave instability and its impact on blob dynamics, and ii) The modeling of edge plasma transport in collisional and semi-collisional regimes, and iii) Impurity transport in edge plasma and the incorporation of plasma-impurity collision operator (hereafter unlike particle collision operator) into COGENT. Over the years four UCSD graduate students (J. Angus, W. Lee, Y. Zhang, and A. Knyazev), funded under this Grant, contributed to the COGENT development/utilization and defended their PhD theses, which are based, partly, on their research related to the COGENT. Three of them (J. Angus, W. Lee, and Y. Zhang) have received “The Best MAE Graduate Student of the Year” awards. Whereas W. Lee after defending his PhD thesis returned back to his home country (South Korea), J. Angus, and Y. Zhang are continuing their research in plasma physics in the LLNL and LANL, whereas A. Knyazev currently is a postdoctoral fellow at Columbia University working on the theoretical aspects of advanced stellarators. The results obtained in our studies were disseminated in many papers published in first-class refereed scientific journals and presented on prestigious scientific meetings including the Annual APS DPP Meeting, Sherwood Fusion Theory Conference, IAEA Fusion Energy Conferences, International Conferences on Plasma-Surface Interactions in Controlled Fusion Devices, International Workshops on Edge Plasma Theory in Fusion Devices, etc. In what follows we highlight the major accomplishments of the UCSD Team in the main areas of our research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Self-Consistent Relativistic Electron Scattering using the Sherlock Scattering Model for X-ray Diagnostics

We present on a new, self-consistent, arbitrary-temperature Romberg integration scheme for modeling electron scattering in materials in a LANL Lagrangian Shock Hydro (LSH) code. Electron beam-target interactions are fundamental to a wide range of scientific and technological applications. When high-energy electron beams hit their target, they may scatter, deposit energy, or ionize the source. These processes govern the behavior and outcomes in nanotechnology manufacturing, electron microscopy, and modern X-ray diagnostics. Simulating these interactions is essential for interpreting experimental results, predicting material responses, and designing efficient tools and experiments. At Los Alamos, this is done using a LSH code, which is a multi-dimension, multi-material, massively parallel, multi-physics code used to simulate applications from asteroid impacts to electron beam interactions. By effectively and efficiently modeling the way that electrons scatter from the beam we can bolster these simulations and more accurately predict experimental outcomes. The model currently implemented in the LSH of interest is based on work by Papp and does not self-consistently preserve momentum in the slightly relativistic regime; here we adopt a model proposed by Braams and Karney and implement a Romberg integration scheme to compute the diffusion tensor. In this paper we will provide background on the Braams-Karney diffusion tensor as well as the Romberg integration scheme we employed to numerically solve for it. We will show that our integration scheme is accurate in solving for the set of scalar potentials used to re-express the diffusion tensor in differential form, and in solving for the diffusion coefficients in the larger LSH code. By using this diffusion tensor rather than the existing Papp one, and numerically integrating it with a Romberg method, we produce much more accurate, self-consistent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Data for "Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium"

This contains the dataset used in "Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium" published in Physics of Plasmas in Oct 2024. It consists of plaintext files as well as .bp files (which may be read by ADIOS open-source software) used for creating the plots that appear in the paper. These are derived from simulation outputs from XGC (X-point Gyrokinetic Code). The data mainly consists of growth rate and frequency measurements as well as poloidal cross-section data.

gyrokinetic↗
NASA NTRS1976

Diamagnetic boundary layers - A kinetic theory

A kinetic theory is presented for boundary layers associated with MHD tangential 'discontinuities' in a collisionless magnetized plasma, such as those observed in the solar wind. The theory consists of finding self-consistent solutions of Vlasov's equation and Maxwell's equation for stationary one-dimensional boundary layers separating two Maxwellian plasma states. Layers in which the current is carried by electrons are found to have a thickness of the order of a few electron gyroradii, but the drift speed of the current-carrying electrons is found to exceed the Alfven speed, and accordingly such layers are not stable. Several types of layers in which the current is carried by protons are discussed; in particular, cases are considered in which the magnetic-field intensity, direction, or both, changed across the layer. In every case, the thickness was of the order of a few proton gyroradii, and the field changed smoothly, although the characteristics depended somewhat on the boundary conditions. The drift speed was always less than the Alfven speed, consistent with stability of such structures. These results are consistent with observations of boundary layers in the solar wind near 1 AU.

Lemaire, J.↗
NASA NTRS1990

Numerical methods for TVD transport and coupled relaxing processes in gases and plasmas

The construction of second-order upwind schemes for nonequilibrium plasmas, for both one- and two-fluid formulations is demonstrated. Coupled relaxation processes, including ionization kinetics and radiative processes and their algorithms for nonequilibrium, multiple temperature conditions are described as well. The paper applies the numerical techniques on some simple test cases, points out critical problems and their solutions, and makes qualitative comparisons with known results, whenever possible.

Cambier, Jean-Luc↗
NASA NTRS1993

Remote plasma enhanced chemical vapor deposition of GaP with in situ generation of phosphine precursors

Thin homoepitaxial films of gallium phosphide (GaP) were grown by remote plasma enhanced chemical vapor deposition utilizing in situ generated phosphine precursors. The GaP forming reaction is kinetically controlled with an activation energy of 0.65 eV. The increase of the growth rate with increasing radio frequency (rf) power between 20 and 100 W is due to the combined effects of increasingly complete excitation and the spatial extension of the glow discharge toward the substrate, however, the saturation of the growth rate at even higher rf power indicates the saturation of the generation rate of phosphine precursors at this condition. Slight interdiffusion of P into Si and Si into GaP is indicated from GaP/Si heterostructures grown under similar conditions as the GaP homojunctions.

Choi, S. W.↗
NASA NTRS1996

Investigation of ion-beam machining methods for replicated x-ray optics

The final figuring step in the fabrication of an optical component involves imparting a specified contour onto the surface. This can be expensive and time consuming step. The recent development of ion beam figuring provides a method for performing the figuring process with advantages over standard mechanical methods. Ion figuring has proven effective in figuring large optical components. The process of ion beam figuring removes material by transferring kinetic energy from impinging neutral particles. The process utilizes a Kaufman type ion source, where a plasma is generated in a discharge chamber by controlled electric potentials. Charged grids extract and accelerate ions from the chamber. The accelerated ions form a directional beam. A neutralizer outside the accelerator grids supplies electrons to the positive ion beam. It is necessary to neutralize the beam to prevent charging workpieces and to avoid bending the beam with extraneous electro-magnetic fields. When the directed beam strikes the workpiece, material sputters in a predicable manner. The amount and distribution of material sputtered is a function of the energy of the beam, material of the component, distance from the workpiece, and angle of incidence of the beam. The figuring method described here assumes a constant beam removal, so that the process can be represented by a convolution operation. A fixed beam energy maintains a constant sputtering rate. This temporally and spatially stable beam is held perpendicular to the workpiece at a fixed distance. For non-constant removal, corrections would be required to model the process as a convolution operation. Specific figures (contours) are achieved by rastering the beam over the workpiece at varying velocities. A unique deconvolution is performed, using series-derivative solution developed for the system, to determine these velocities.

Drueding, Thomas W.↗