Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA DYNAMICS”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulating single-particle dynamics in magnetized plasmas: The RMF code

The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Runaway electron beam dynamics at low plasma density in DIII-D: energy distribution, current profile, and internal instability

Parameters of the post-disruption runaway electron (RE) beam in the collisionless background plasma achieved after deuterium injection are investigated in DIII-D. The spatially resolved RE energy distribution function is measured for the first time during the RE plateau stage by inverting hard X-ray bremsstrahlung spectra. It has maximum energy up to 20 MeV and a non-monotonous feature at 5{6 MeV observed only in the core of the beam supporting the possibility of kinetic instabilities. The RE current profile is reconstructed for the first time using the spatially resolved RE energy distribution. It is found to be more peaked than the pre-disruption plasma current, with higher internal inductance, suggesting preferential formation of REs in the core plasma or potentially a radially inward motion of the REs. The accessed relatively low-current (180 kA) RE beam is found stable, likely due to its elevated q profile. From this base stable equilibrium, instability is accessed by ramping the solenoid current. Under a large applied accelerating loop voltage an internal kink mode is observed. The kink mode leads to sawtooth-like relaxation of the RE current profile, but drives no RE loss. Under conditions of a large decelerating voltage, chirping Alfvénic instabilities driven by REs during the RE plateau stage are observed for the first time. The low-frequency (0.1{7 MHz) Alfvénic instabilities, supposedly compressional Alfvén eigenmodes, are found to be correlated with a RE loss signal. Finally, this provides novel plasma conditions to study these instabilities and investigate opportunities to utilize them for RE control. Bifurcated electron cyclotron emission (ECE) spectrum is measured during the collisionless RE plateau, with a break point at ≈100 GHz suggesting resonant absorption of the ECE at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Strong Coulomb coupling influences ion and neutral temperatures in atmospheric pressure plasmas

Molecular dynamics simulations are used to model ion and neutral temperature evolution in partially-ionized atmospheric pressure plasma at different ionization fractions. Results show that ion–ion interactions are strongly coupled at ionization fractions as low as 10 –5 and that the temperature evolution is influenced by effects associated with the strong coupling. Specifically, disorder-induced heating is found to rapidly heat ions on a timescale of the ion plasma period (~10 s ps) after an ionization pulse. This is followed by the collisional relaxation of ions and neutrals, which cools ions and heats neutrals on a longer (~ns) timescale. Slight heating then occurs over a much longer (~100 s ns) timescale due to ion-neutral three-body recombination. An analytic model of the temperature evolution is developed that agrees with the simulation results. Finally, a conclusion is that strong coupling effects are important in atmospheric pressure plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Accuracy of the time-averaged ponderomotive approximation for laser-plasma accelerator modeling

Reliable modeling of laser-plasma accelerators, where a short and intense laser pulse propagates in an underdense plasma over long distances, is a computationally challenging task. This is due to the great disparity among the scales involved in the modeling, ranging from the micrometer scale of the laser wavelength to, for instance, the meter scale of the laser-plasma interaction length for a multi-GeV-class laser-plasma accelerator. To reduce such imbalance, the time-averaged ponderomotive approximation may be used, where the plasma particle dynamics is analytically averaged over the laser frequency, and only spatiotemporal scales associated with the laser envelope are retained in the calculations, resulting in significant computational savings. Here, we characterize the accuracy and robustness of the time-averaged ponderomotive approximation for a range of laser parameters of interest for present and future laser-plasma accelerators, and we show that the error introduced by the averaging process is small in all relevant cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Whistler–Alfvén Turbulence in a Nonneutral Ultrarelativistic Pair Plasma

The large-scale dynamics of most conventional space and astrophysical plasmas are predominantly governed by Alfvén modes, which are low-frequency magnetohydrodynamic modes existing in magnetized media. At scales smaller than the ion gyroscale or frequencies exceeding the ion cyclotron frequency, the Alfvén modes transform into kinetic-Alfvén or whistler modes that significantly contribute to plasma dynamics. However, this scenario reverses in nonneutral pair plasmas, such as those found in the magnetospheres of pulsars and magnetars, around rotating black holes, and in their relativistic jets, as well as in certain laboratory plasmas. In these systems, the large-scale dynamics are governed by hybrid whistler–Alfvén modes, which transform into pure Alfvén modes at smaller scales. We derive the nonlinear equations that describe the dynamics of whistler–Alfvén modes in ultrarelativistic nonneutral magnetically dominated pair plasma and discuss the spectrum of turbulence governed by these equations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effects of low-frequency voltage on nonlinear standing wave excitation, plasma uniformity, and ion dynamics in dual-frequency asymmetric capacitive discharges

Abstract It is known that in very-high-frequency (VHF) capacitively coupled plasmas, the higher harmonics generated by nonlinear sheath motion can enhance the standing wave effect (SWE), which can lead to center-peaked plasma density profiles. In this work, an improved nonlinear electromagnetic model incorporating a transmission line model, an electron momentum balance model, a bulk plasma model, a collisionless nonlinear numerical sheath model, and an ion Monte-Carlo collision (MCC) model is developed to study the effects of low-frequency (LF) voltage V L on the nonlinear standing wave excitation, plasma uniformity, and ion energy and angular distribution functions (IEDFs and IADFs) in dual-frequency (DF) asymmetric capacitive argon discharges at relatively low pressure of 3 Pa. The plasma diffusion in the radial direction and ion dynamics within the LF oscillating sheath are self-consistently considered. The LF voltage V L at 2 MHz varies from 0 to 700 V while the HF voltage V H at 60 MHz is fixed at 100 V. Simulation results indicate that without the addition of an LF source (i.e. V L = 0 V), there are a considerable number of high-order harmonics with short wavelengths, leading to significant SWE and central peak in the radial plasma density profile. Nevertheless, the high-order harmonic excitations tend to be weakened and merely occur around the phase of the full LF sheath collapse due to a shorter characteristic damping time of the surface waves as V L increases. This, combined with increased surface wavelengths of both the driving frequency and the higher harmonics at a higher V L , leads to suppressed standing waves and improved plasma uniformity. Meanwhile, the simulations show that both the low and the high energy peaks of IEDF move towards higher energies, and the energy peak separation width ΔEbecomes wider with the increase of V L . The IEDF at the radial center of the powered electrode exhibits a broader ΔEthan that at the edge. For the IADF, an increased V L results in more ions incident on the electrode with a smaller deflection angle. Because of a thinner sheath and a higher sheath voltage at the electrode center, the peak value of IADF at the electrode center is greater than that at the edge.

Physics↗
DOE OSTI2024

Kinetic simulations comparing quasi-parallel and quasi-perpendicular piston-driven collisionless shock dynamics in magnetized laboratory plasmas

Magnetized collisionless shocks are common in astrophysical systems, and scaled versions can be created in laboratory experiments by utilizing laser-driven piston plasmas to create these shocks in a magnetized background plasma. A key parameter for these experiments is the angle θB between the shock propagation direction and the background magnetic field. We performed quasi-1D piston-driven shock simulations to explore shock formation, evolution, and key observables relevant to laboratory experiments for a range of shock angles between θB=90° to θB=30°. Our results show that the spatial and temporal scales of shock formation for all angles considered are similar when expressed in terms of the perpendicular component of the magnetic field. In a steady state, ion and electron temperatures become more isotropic, and the electron-to-ion temperature ratio is higher for smaller θB. At θB=30°, ion heating parallel to the magnetic field becomes dominant, associated with more ions being reflected at one discontinuity and subsequently trapped by the next discontinuity due to shock reformation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Dynamics of the water-plasma interface in various discharge modes of atmospheric-pressure plasmas

This paper presents an experimental study of the dynamics of the water-plasma interface during the interaction high-voltage (HV) atmospheric pressure discharge with the water surface. Under the influence of the voltage applied to the pin type electrode, discharges are formed at the air layer between the sharp-tip HV electrode and deionized water. Three discharge regimes were identified from synchronized electrical waveforms and imaging: (i) a weak linear regime at 3–10.6 kV, (ii) a branching streamer regime at the maximum applied voltage of 12.6 kV, and (iii) a continuous arc-like regime occurring during the transition to lower voltage and higher current (e.g., U ≈ 3.2 kV, I ≈ 4.07 mA). High-speed shadowgraph images showed a symmetric interfacial cavity whose depth increased nonlinearly with voltage from h 0 ≈ 0.1 mm at 3 kV to h 0 ≈ 2.7 mm at 10.6 kV, reaching a maximum depth of h 0 ≈ 5.9 mm at 12.6 kV, while the cavity disappeared in the continuous-channel regime and was replaced by outward-propagating wave-like motion. The deformation of the water surface during the discharge-water interaction is governed by the balance between electric field forces, surface tension, gravitational forces, and electrohydrodynamic forces, whose relative contributions vary with the applied voltage. This effect was explained by quantitatively estimating the magnitudes of the acting forces and establishing the force balance. Optical emission spectroscopy in the continuous-channel regime indicated air plasma signatures and yielded a gas temperature of approximately 400 K, while prolonged operation (≈ 20 min) increased the water temperature to ~ 70 °C, reduced the water-layer thickness from 6 mm to 3 mm, and decreased pH from 7 to 4. These regime-resolved, quantitative results clarify how the dominant interfacial forcing shifts with discharge mode and provide a mechanistic basis for controlling plasma-water interactions in atmospheric-pressure applications.

atmospheric pressure plasma↗
DOE OSTI2020

Measurement of Kinetic-Scale Current Filamentation Dynamics and Associated Magnetic Fields in Interpenetrating Plasmas

Here, we present the first local, quantitative measurements of ion current filamentation and magnetic field amplification in interpenetrating plasmas, characterizing the dynamics of the ion Weibel instability. The interaction of a pair of laser-generated, counterpropagating, collisionless, supersonic plasma flows is probed using optical Thomson scattering (TS). Analysis of the TS ion-feature revealed anticorrelated modulations in the density of the two ion streams at the spatial scale of the ion skin depth c/ω pi = 120 μm, and a correlated modulation in the plasma current. The inferred current profile implies a magnetic field amplitude ~30 ± 6 T , corresponding to ~1% of the flow kinetic energy, indicating that magnetic trapping is the dominant saturation mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Developments and Challenges in the Design of the ITER DRGA

We report the ITER Diagnostic Residual Gas Analyzer (DRGA) will measure the distribution of gas species, i.e., deuterium (D), tritium (T), and impurities, in the divertor exhaust stream and in the plasma periphery, with time resolution relevant to fusion plasma–wall particle dynamics. The uniqueness of the DRGA, over previous implementations of plasma dynamics residual gas analysis, is an integrated approach, combining mass and low-temperature plasma-activated optical spectroscopy, in a differentially pumped analysis station. A further unique feature of the ITER divertor-specific DRGA is an ~8-m separation of the analysis station from the sampled pumping duct, while still maintaining a ~1-s response time for hydrogen isotopic concentrations. ITER DRGA final design activities are strongly benefiting from testing of prototypical DRGA components and methods on present fusion devices, most currently on JET and W7-X. DRGA systems are implemented on both these devices and include sensors (and pumping methods) that are directly relevant to the ITER DRGA design. The recent JET-DTE2 campaign has provided the first experience on operating the combined ITER DRGA sensors with D-T plasmas. While enhancing system design for ITER, this experience on operating devices has also revealed additional engineering challenges, which further guide the continuing final design project. Meanwhile, the recent determination that the ITER DRGA, with slight optimization, will resolve the helium isotopes well enough to support an ITER pre-DT, He-3-based heating scheme, has greatly increased ITER Research Program interest in the DRGA and its implementation well ahead of the DT phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Expansion dynamics of femtosecond laser-induced plasmas: Influence of thermophysical plasma properties

This study investigates the expansion dynamics of femtosecond laser-induced plasmas, emphasizing the impact of plasma thermophysical properties and ambient gas composition. Through shadowgraphy experiments and multiphase computational fluid dynamics (CFD) simulations, the influence of parameters such as heat capacity, molecular weight, and thermal conductivity on plume morphology, shockwave evolution, and energy dissipation mechanisms is examined. A mixture multiphase model is implemented to capture the interaction between the plasma and the surrounding gas. Simulation results reveal that plasma expansion is strongly inertia-driven. Results show that differences in plasma properties and ambient conditions affect the shape and temperature distribution of the expanding plume. The early-stage dynamics are primarily dictated by pressure forces, whereas thermal and viscous effects play a growing role in the plume's behavior during later stages of expansion. The CFD findings show the necessity of accurate initial condition characterization, including crater geometry and plasma pressure and temperature, for reliable modeling of plasma evolution in laser ablation processes.

CFD modeling↗
DOE OSTI2020

General formulas for adiabatic invariants in nearly periodic Hamiltonian systems

While it is well known that every nearly periodic Hamiltonian system possesses an adiabatic invariant, extant methods for computing terms in the adiabatic invariant series are inefficient. The most popular method involves the heavy intermediate calculation of a non-unique near-identity coordinate transformation, even though the adiabatic invariant itself is a uniquely defined scalar. Here, a less well-known method, developed by S. Omohundro, avoids calculating intermediate sequences of coordinate transformations but is also inefficient as it involves its own sequence of complex intermediate calculations. In order to improve the efficiency of future calculations of adiabatic invariants, we derive generally applicable, readily computable formulas for the first several terms in the adiabatic invariant series. To demonstrate the utility of these formulas, we apply them to charged-particle dynamics in a strong magnetic field and magnetic field-line dynamics when the field lines are nearly closed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Fast correlation heating in moderately coupled electron–ion plasmas

If the electrons in a plasma are suddenly heated, the resulting change in Debye shielding causes the ion kinetic energy to quickly increase. For the first time, this correlation heating , which is much faster than collisional energy exchange, is rigorously derived for a moderately coupled, electron–ion plasma. The electron–ion mass ratio is taken to be the smallest parameter in the Bogoliubov–Born–Green–Kirkwood–Yvon hierarchy, smaller even than the reciprocal of the plasma parameter. This ordering differs from conventional kinetic theory by making the electron collision rates faster than the ion plasma frequency, which allows stronger coupling and makes the ion heating a function only of the total energy supplied to the electrons. The calculation uses known formulae for correlations in a two-temperature plasma, for which a new, elementary derivation is presented. Suprathermal ions may be created more rapidly by this mechanism than by ion–electron Coulomb collisions. This means that the use of a femtosecond laser pulse could potentially help to achieve ignition in certain fast ignition approaches to inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗