Kinetic simulations of piston-driven collisionless shock formation in magnetized laboratory plasmas
Explore the source record for details and available documents.
D. B. Schaeffer↗
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.
A kinetic study of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas
We present particle-in-cell simulations with Monte Carlo collisions of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas. We study the energy balance in the one-dimensional expansion of a hot-spot by simulating Coulomb collisions, fusion reactions, and bremsstrahlung emission with a Monte Carlo model and inverse bremsstrahlung absorption using a new PIC model. This allows us to self-consistently capture the alpha particle heating and radiative losses in the expanding hot-spot and surrounding cold fuel. After verifying our model in a code-to-code comparison with both kinetic and fluid codes for the case of a deuterium–tritium hot-spot, we simulate the expansion of a proton–boron hot-spot initialized at 200 keV and 1,000 g/cm 3 . Our model predicts that energy radiated by the hot-spot is recaptured by the surrounding high-density opaque fuel reducing the expansion work done by the propagating burn wave. As a result, we find the net fusion energy produced over the course of $20$~ps is twice the initial hot-spot energy independent of whether radiation physics is included.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma
A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.
neoclassical transport↗
Solar Wind Access to Lunar Polar Craters: Feedback Between Surface Charging and Plasma Expansion
Determining the plasma environment within permanently shadowed lunar craters is critical to understanding local processes such as surface charging, electrostatic dust transport, volatile sequestration, and space weathering. In order to investigate the nature of this plasma environment, the first two-dimensional kinetic simulations of solar wind expansion into a lunar crater with a self-consistent plasma-surface interaction have been undertaken. The present results reveal how the plasma expansion into a crater couples with the electrically-charged lunar surface to produce a quasi-steady wake structure. In particular, there is a negative feedback between surface charging and ambipolar wake potential that allows an equilibrium to be achieved, with secondary electron emission strongly moderating the process. A range of secondary electron yields is explored, and two distinct limits are highlighted in which either surface charging or ambipoiar expansion is responsible for determining the overall wake structure.
Zimmerman, M. I.↗
Excitation of fishbone-like mode in tokamaks due to bounce resonances of trapped thermal ions
When the drift kinetic effect of thermal ions is taken into account, a high-frequency fishbone-like mode (FLM) is found to be driven unstable by trapped thermal ions (TTIs) in tokamak plasmas, according to self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid modeling utilizing the MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)] as well as an analytic theory. It is found that, similar to energetic particles, TTIs can also stabilize the internal kink mode, whereas the FLM is excited when the effective beta of TTIs exceeds a threshold value. The real frequency of the FLM is comparable to the bounce frequency of TTIs. The mode structure of the FLM can be significantly different from the conventional step-like function for the associated plasma radial displacement. This drift kinetic induced modification of the mode structure near the q = 1 surface is captured by non-perturbative MHD-kinetic hybrid computations with MARS-K. Furthermore, the FLM can only be triggered by TTIs at sufficiently high thermal temperatures. Both the FLM and the internal kink can be stabilized by sufficiently fast plasma toroidal rotation and parallel sound wave damping. These two conditions of high thermal temperature and (fast) flow stabilization, though making it challenging to observe the TTI-driven FLM in present day experiments, are favorable for the mode excitation in future reactor scale devices.
36 MATERIALS SCIENCE↗
Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas
Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.
Plasma kinetic theory↗
Effects of wave damping and finite perpendicular scale on three-dimensional Alfvén wave parametric decay in low-beta plasmas
Shear Alfvén wave parametric decay instability (PDI) provides a potential path toward significant wave dissipation and plasma heating. However, fundamental questions regarding how PDI is excited in a realistic three-dimensional (3D) open system and how the finite perpendicular wave scale—as found in both laboratory and space plasmas—affects the excitation remain poorly understood. Here, we present the first 3D, open-boundary, hybrid kinetic-fluid simulations of kinetic Alfvén wave PDI in low-beta plasmas. Key findings are that the PDI excitation is strongly limited by the wave damping present, including electron–ion collisional damping (represented by a constant resistivity) and geometrical attenuation associated with the finite-scale Alfvén wave, and ion Landau damping of the child acoustic wave. The perpendicular wave scale alone, however, plays no discernible role: waves of different perpendicular scales exhibit similar instability excitation as long as the magnitude of the parallel ponderomotive force remains unchanged. These findings are corroborated by theoretical analysis and estimates. This new understanding of 3D kinetic Alfvén wave PDI physics is essential for laboratory study of the basic plasma process and may also aid future evaluation of the relevance/role of PDI in low-beta space plasma.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
An acceleration mechanism for cometary plasma tails
Cometary plasma tail formation by the interaction between the solar wind plasma flow and the plasma at the head of the coma is discussed using the unipolar electric generation theory. The plasma in the 'plasma tail' is almost directly accelerated from the cometary ionopause along the sun-nucleus line where the tail current flows. For steady state solar wind conditions, the cometary 'plasma tail' velocity distribution is obtained self-consistently. The solution of a kinetic equation gives the velocity of the cometary plasma as a function of the cometary tail position. The characteristic length is 3 x 10 to the 6th km when the plasma density near the nucleus of the comet is 10 to the 6th/cu cm and the component of the interplanetary magnetic field perpendicular to the solar wind flow is 3 nT. The tailward cometary plasma is finally accelerated to the speed of the solar wind. The theory is compared with the observational plasma velocities in the tails of comet Bennett (1970II) and Comet Halley (1985).
Minami, S.↗
Plasma edge and scrape-off layer turbulence in gyrokinetic simulations of negative triangularity plasmas
Gyrokinetic simulations in the long-wavelength or drift-kinetic limit are carried out of DIII-D inner-wall-limited (IWL) plasmas to investigate the effect of triangularity on edge and scrape-off layer (SOL) turbulence. The effect of neutral interactions and triangularity on plasma blobs is explored due to the impact blobs can have in setting the SOL width or introducing impurities through interactions with plasma-facing components. Seeded blob simulations with neutrals in shaped SOL scenarios demonstrate that increasing elongation, triangularity, or Shafranov shift decreases radial blob velocities, but neutral interactions have a minor effect. Fully turbulent simulations of DIII-D IWL plasmas include both open- and closed-field-line regions. The negative triangularity (NT) simulation has lower average core T e , lower normalized T e fluctuations, and lower fluxes, but a greater number of coherent structures (blobs) identified with increased size and velocity, on average. Density and electron temperature profiles are within a factor of 2 of experimental values. In conclusion, the increased trapped electron particle fraction in NT simulations is consistent with previous studies.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Upper-Hybrid Waves Driven by Meandering Electrons Around Magnetic Reconnection X Line
Magnetic reconnection is a fundamental process in collisionless space plasma environment, and plasma waves relevant to the kinetic interactions can have a significant impact on the multiscale behavior of reconnection. Here, we present Magnetospheric Multiscale (MMS) observations during an encounter of an X line of symmetric magnetic reconnection in the magnetotail. The X line is characterized by reversals of ion and electron jets and electromagnetic fields, agyrotropic electron velocity distribution functions (VDFs), and an electron-scale current sheet. MMS observe large-amplitude nonlinear upper-hybrid (UH) waves on both sides of the neutral line, and the wave amplitudes have highly localized distribution along the normal direction. Here, the inbound meandering electrons drive the UH waves, releasing the free energy stored from the reconnection electric field along the meandering trajectories. The interaction between the meandering electrons and the UH waves may modify the balance of the reconnection electric field around the X line.
58 GEOSCIENCES↗
RAM-SCB Runs on Request at CCMC
The Ring current–Atmosphere interactions Model (RAM) with Self-Consistent magnetic (B) field (SCB) combines a large-scale kinetic model of ring current plasma with a three-dimensional (3-D) force-balanced model of the terrestrial magnetic field. RAM-SCB simulates the evolution of major ion species (H + , O + , and He + by default) and electrons as a function of azimuth, radial distance, energy, and pitch angle. Simulation outputs include the Dst index, and pressure and differential flux of the modeled species, thus providing benefit as a science code and to inform surface charging hazard within the model domain. Version 2.2 of this open-source simulation code has now—as of 8 August 2023—been made available for Runs-On-Request via the Community Coordinated Modeling Center.
79 ASTRONOMY AND ASTROPHYSICS↗
Advances in laser-based bremsstrahlung x-ray sources. I. Optimizing laser-accelerated electrons
In this work, we have performed a suite of kinetic simulations of relativistic laser–plasma interaction under settings relevant to recent and planned experiments on a variety of laser systems. The goal of the study is to illuminate the physics of laser–target coupling and to provide guidance for how to optimize these sources for applications. It is shown that the production of relativistic electrons is maximized when conditions of relativistic induced transparency (RIT) in dense plasmas can be achieved over a large interaction volume at the time of arrival of most intense part of the laser pulse. RIT is shown to enhance both the numbers of relativistic electrons and the energies of the electrons, leading to an increased x-ray dose. A variety of approaches to enhancing laser–target coupling are considered. These include optimizing the effects of low-density pre-plasma (arising either from finite laser pedestal or from the use of foam coatings) and of modifying the laser focusing geometry to reduce effects of filamentation and self-focusing. Evidence of a novel approach to achieving stable laser propagation over distances of tens of micrometers in a plasma gradient is also presented. These conditions coincide with plasma and laser conditions explored in recent experiments on the Omega EP laser system and compare favorably with an analytic criterion for stable laser propagation in relativistically underdense plasma obtained from a nonlinear Wentzel–Kramers–Brillouin analysis.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Energy deposition in parallel-plate plasma accelerators
Parallel plate plasma accelerator energy deposition, considering kinetic and thermal modes based on flow velocity, temperature and Mach number measurements
Di Capua, M. S.↗
Correlations between plasma variables and the deposition process of Si films from chlorosilanes in low pressure RF plasma of argon and hydrogen
The dissociation of chlorosilanes to silicon and its deposition on a solid substrate in a RF plasma of mixtures of argon and hydrogen were investigated as a function of the macrovariables of the plasma. The dissociation mechanism of chlorosilanes and HCl as well as the formation of Si in the plasma state were studied by sampling the plasma with a quadrupole mass spectrometer. Macrovariables such as pressure, net RF power input and locations in the plasma reactor strongly influence the kinetics of dissociation. The deposition process of microcrystalline silicon films and its chlorine contamination were correlated to the dissociation mechanism of chlorosilanes and HCl.
Avni, R.↗
Helium escape from the Earth's atmosphere - The charge exchange mechanism revisited
We have studied the escape of neutral helium from the terrestrial atmosphere through exothermic charge exchange reactions between He(+) ions and the major atmospheric constituents N2, O2 and O. Elastic collisions with the neutral background particles were treated quantitatively using a recently developed kinetic theory approach. An interhemispheric plasma transport model was employed to provide a global distribution of He(+) ions as a function of altitude, latitude and local solar time and for different levels of solar ionization. Combining these ion densities with neutral densities from an MSIS model and best estimates for the reaction rate coefficients of the charge exchange reactions, we computed the global distribution of the neutral He escape flux. The escape rates show large diurnal and latitudinal variations, while the global average does not vary by more than a factor of three over a solar cycle. We find that this escape mechanism is potentially important for the overall balance of helium in the Earth's atmosphere. However, more accurate values for the reaction rate coefficients of the charge exchange reactions are required to make a definitive assessment of its importance.
Lie-Svendsen, O.↗
Heat flux and viscosity of ions in the collisionless solar wind
Between 1 and 2 solar radii, the Coulomb-collision mean free path for thermal ions exceeds the scale height of the solar atmosphere. The expanding solar plasma becomes collisionless and the kinetics of the solar wind are no longer dominated by thermalizing collisions. The usual Braginskii-type expressions for solar wind ion heat flux and viscosity are no longer valid. However, another microscale still exists in the solar wind, dictated by the gyro-radius of ions in the turbulent embedded solar wind magnetic field. Wave-particle interactions will act to isotropize (but not thermalize) particle distributions, and the relevant microscale for this process is the ion gyro-radius. The ion distribution can be modelled as undergoing isotropizing 'collisions,' with the relevant mean free path scaling with gyro-radius. Here, the author presents the heat flux and viscosity expected for solar wind protons which are relaxing to isotropy on a microscale that scales with gyro-radius. The collisionless viscosity and heat flux have a functional dependence different than their collisional analogs. The collisional expressions for ion viscosity and heat flux drastically overestimate the efficiency of diffusive energy and momentum transport actually operative in the solar wind.
Williams, L. L.↗
Electron-Scale Measurements of Magnetic Reconnection in Space
Magnetic reconnection is a fundamental physical process in plasmas whereby stored magnetic energy is converted into heat and kinetic energy of charged particles. Reconnection occurs in many astrophysical plasma environments and in laboratory plasmas. Using measurements with very high time resolution, NASA's Magnetospheric Multiscale (MMS) mission has found direct evidence for electron demagnetization and acceleration at sites along the sunward boundary of Earth's magnetosphere where the interplanetary magnetic field reconnects with the terrestrial magnetic field. We have (i) observed the conversion of magnetic energy to particle energy; (ii) measured the electric field and current, which together cause the dissipation of magnetic energy; and (iii) identified the electron population that carries the current as a result of demagnetization and acceleration within the reconnection diffusion/dissipation region.
Burch, J. L.↗
Direct Energy Conversion for Nuclear Propulsion at Low Specific Mass
The project will continue the FY13 JSC IR&D (October-2012 to September-2013) effort in Travelling Wave Direct Energy Conversion (TWDEC) in order to demonstrate its potential as the core of a high potential, game-changing, in-space propulsion technology. The TWDEC concept converts particle beam energy into radio frequency (RF) alternating current electrical power, such as can be used to heat the propellant in a plasma thruster. In a more advanced concept (explored in the Phase 1 NIAC project), the TWDEC could also be utilized to condition the particle beam such that it may transfer directed kinetic energy to a target propellant plasma for the purpose of increasing thrust and optimizing the specific impulse. The overall scope of the FY13 first-year effort was to build on both the 2012 Phase 1 NIAC research and the analysis and test results produced by Japanese researchers over the past twenty years to assess the potential for spacecraft propulsion applications. The primary objective of the FY13 effort was to create particle-in-cell computer simulations of a TWDEC. Other objectives included construction of a breadboard TWDEC test article, preliminary test calibration of the simulations, and construction of first order power system models to feed into mission architecture analyses with COPERNICUS tools. Due to funding cuts resulting from the FY13 sequestration, only the computer simulations and assembly of the breadboard test article were completed. The simulations, however, are of unprecedented flexibility and precision and were presented at the 2013 AIAA Joint Propulsion Conference. Also, the assembled test article will provide an ion current density two orders of magnitude above that available in previous Japanese experiments, thus enabling the first direct measurements of power generation from a TWDEC for FY14. The proposed FY14 effort will use the test article for experimental validation of the computer simulations and thus complete to a greater fidelity the mission analysis products originally conceived for FY13.
Scott, John H.↗