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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electrical conductivity of a warm neutron star crust in magnetic fields: Neutron-drip regime

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at nonzero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve protoneutron stars born in supernova explosions, binary neutron star mergers, and accreting neutron stars. The temperature-density range studied covers the transition from a semidegenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (i) the screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (ii) the correlations of the ionic component in a one-component plasma, and (iii) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropy introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust, which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviations of relaxation time and components of the conductivity tensor from the average are below ≤25% except close to crust-core transition, where nonspherical nuclear structures are expected. Finally, our results can be used in dissipative magnetohydrodynamics simulations of warm compact stars.

Physics↗
DOE OSTI2020

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗
DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetic reconnection and kinetic waves generated in the Earth's quasi-parallel bow shock

Magnetic reconnection in quasi-parallel shocks, relevant to the Earth’s bow shock, is studied by means of two-dimensional full particle-in-cell simulations. As the Alfvénic Mach number increases, the propagation direction of the waves excited in the transition region changes, and the shock becomes more turbulent with more reconnection sites. In the higher Mach number shock, abundant electron-only reconnection sites are generated with scales on the order of the ion skin depth or less. Non-reconnecting current sheets can also generate electron jets and energy dissipation can occur there as well. However, non-reconnecting current sheets with the magnetic field reversal typically show a smaller energy dissipation rate than reconnecting current sheets. In the shock transition region, two types of waves are responsible for driving reconnection: one has a wavelength around three ion skin depths (d i ), and the other has a wavelength less than 1 d i . Here, electron and ion distribution functions show that in regions where the former type of waves is excited, there are two ion beams and a single-peaked electron distribution. In contrast, in regions where the latter type of waves is excited, there are multiple electron and ion beams. The waves propagating obliquely to the magnetic field bend the magnetic field lines, and magnetic reconnection occurs where oppositely directed field lines come into contact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effective Resistivity in Relativistic Collisionless Reconnection

Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing nonthermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations, we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without a guide field. In X-points, where diverging flows result in a nondiagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain, we identify the mechanisms driving the nonideal electric field using a full Ohm law by means of a statistical analysis based on our PIC simulations. We show that the nonideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated nonuniform effective resistivity model that is negligible on global scales and becomes significant only locally in X-points. It captures the properties of collisionless reconnection with the aim of mimicking its essentials in nonideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of NBI-induced energetic ions on internal kink stability in HL-2M

A systematic investigation is carried out, studying the effect of the neutral beam injection induced energetic particles (EPs) on the n = 1 (n is the toroidal mode number) internal kink (IK) instability in the HL-2M tokamak, utilizing the MARS-F/K code [Liu et al., Phys. Plasmas 7, 3681 (2000) and 15, 112503 (2008)]. A high-beta sawteething HL-2M scenario, simulated by the TRANSP code [Breslau et al. Computer Software (2018)], is chosen for this study. Compared to the fluid model, non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid computations with MARS-K show a generally stabilization effect on the IK, due to drift kinetic resonances associated with EPs. Further, the bounce resonance of trapped EPs has minor influence on the mode stability. In the absence of the plasma equilibrium flow and with the assumed particle pitch distribution, the transit resonance of co-current (countercurrent) passing EPs destabilizes (stabilizes) the IK. With plasma flow, both co- and countercurrent passing EPs tend to stabilize the mode, but the effect is stronger with the countercurrent particles. These modeling results provide useful guidance for interpreting MHD instabilities in the future high-performance experiments in HL-2M.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Global flow structure and exact formal transseries of the Gubser flow in kinetic theory

In this work we introduce the generic conditions for the existence of a non-equilibrium attractor that is an invariant manifold determined by the long-wavelength modes of the physical system. We investigate the topological properties of the global flow structure of the Gubser flow for the Israel-Stewart theory and a kinetic model for the Boltzmann equation by employing Morse-Smale theory. We present a complete classification of the invariant submanifolds of the flow and determine all the possible flow lines connecting any pair of UV/IR fixed points. The formal transseries solutions to the Gubser dynamical system around the early-time (UV) and late-time (IR) fixed points are constructed and analyzed. It is proven that these solutions are purely perturbative (or power-law asymptotic) series with a finite radius of convergence. Based on these analyses, we find that Gubser-like expanding kinetic systems do not hydrodynamize owing to the failure of the hydrodynamization process which heavily relies on the classification of (non)hydrodynamic modes in the IR regime. This is in contrast to longitudinal boost-invariant plasmas where the asymptotic dynamics is described by a few terms of the hydrodynamic gradient expansion. We finally compare our results for both Bjorken and Gubser conformal kinetic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
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↗