Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma models”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling nanoparticle charge distribution in the afterglow of non-thermal plasmas and comparison with measurements

Particle charging in the afterglows of non-thermal plasmas typically take place in a non-neutral space charge environment. We model the same by incorporating particle-ion collision rate constant models, developed in prior work by analyzing particle-ion trajectories calculated using Langevin Dynamics simulations, into species transport equations for ions, electrons and charged particles in the afterglow. A scaling analysis of particle charging and additional Langevin Dynamics calculations of the particle-ion collision rate constant are presented to extend the range of applicability to ion electrostatic to thermal energy ratios of 300 and diffusive Knudsen number (that scales inversely with gas pressure) up to 2000. The developed collision rate constant models are first validated by comparing predictions of particle charge against measured values in a stationary, non-thermal DC plasma from past PK-4 campaigns published in Phys. Rev. Lett. 93(8): 085001 and Phys. Rev. E 72(1): 016406). The comparisons reveal excellent agreement within ±35% for particles of radius 0.6,1.0,1.3 μm in the gas pressure range of ~20-150 Pa. The experiments to probe particle charge distributions by Sharma et al. (J. Physics D: Appl. Phys. 53(24): 245204) are modeled using the validated particle-ion collision rate constant models and the calculated charge fractions are compared with measurements. The comparisons reveal that the ion/electron concentration and gas temperature in the afterglow critically influence the particle charge and the predictions are generally in qualitative agreement with the measurements. Along with critical assessment of the modeling assumptions, several recommendations are presented for future experimental design to probe charging in afterglows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Resolution requirements for energy conservation in kinetic plasma simulations

The kinetic codes used to model the coupled dynamics of electromagnetic fields and charged particle transport have requirements for spatial, temporal, and charge resolution. These requirements may vary by the solution technique and scope of the problem. In this report, we investigate the resolution limits in the energy-conserving implicit particle-in-cell code CHICAGO. This report has the narrow aim of determining the maximum acceptable grid spacing for the dense plasmas generated in models of z-pinch target gases and power-flow electrode plasmas. In the 2D sample problem, the plasma drifts without external forces with velocity of 10 cm/µs. Simulations are scaled by plasma density to maintain uniform strides across the plasma and from the plasma to the boundaries. Additionally, the cloud-in-cell technique is used with 400 particles per cell and Δ t = 0.85× the Courant limit. For the linear cloud distribution, the criterion for conserving energy is Δ E/E tot < 0.01 for 50,000 time steps. The grid resolution is determined to crudely be Δx ≲ 3 l s , where l s is the electron collisionless skin depth. For the second-order cloud distribution the criterion is Δ E/E tot < 0.005 yielding Δx ≤ 15 l s . These scalings are functions of the chosen v d , Δ t , particles-per-cell, and number of steps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Generation of Large-Volume High-Pressure Plasma by Spatio-Temporal Control of Space Charge

Due to the fundamental limitations of scaling, very little progress has been made towards achieving a large-volume dense non-equilibrium plasma near atmospheric pressures. Commercially available state-of-the art glow-like atmospheric plasma devices for industrial applications are either narrow tubes or wide slits with narrow openings. Often multiple sources are put together in various configurations to process larger surfaces. The traditional methods of exciting electrodes create spatially fixed electric fields. As a result, the space charge at atmospheric pressure tends to be spatially confined resulting in non-uniformity which eventually leads to instability as the discharge is scaled. Theoretical work done under this project showed that it is possible to generate a spatially rotating electric field by exciting a set of electrodes with phase staggered sinusoidal waveforms. The modeling and simulations were done using plasma fluid models. It was shown that such a field can produce a uniform plasma. At the conclusion of the project, experimental proof of concept with an eight-electrode system in various gases (Air, Helium and Argon) was demonstrated. Power measurements and spectral investigation show that the concept can be used to generate a uniform stable plasma. This plasma source has the potential of opening new applications of nonthermal plasma including combustion of carbon-free fuel. The current limitation of the proposed method in scaling to higher volume and pressure is the need for multiple high voltage amplifiers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Simulating plasma wave propagation on a superconducting quantum chip

Quantum computers may one day enable the efficient simulation of strongly coupled plasmas that lie beyond the reach of classical computation in regimes where quantum effects are important and the scale separation is large. Here, in this article, we take a first step toward efficient simulation of quantum plasmas by demonstrating linear plasma wave propagation on a superconducting quantum chip. Using high-fidelity and highly expressive device-native gates, combined with an error-mitigation technique, we simulate the scattering of laser pulses from inhomogeneous plasmas. Our approach is made feasible by the identification of a suitable local spin model whose excitations mimic plasma waves, and whose circuit implementation requires a lower gate count than other proposed approaches that would require a future fault-tolerant quantum computer. This work opens avenues to study more complicated phenomena that cannot be simulated efficiently on classical computers, such as nonlinear quantum dynamics when strongly coupled plasmas are driven out of equilibrium.

general physics↗
DOE OSTI2023

Numerical and Experimental Study of an Aircraft Igniter Plasma Jet Discharge

The spark discharge of an aircraft plasma jet igniter is studied using high-fidelity numerical simulations and X-ray radiography measurements. The target problem here features the thermal expansion of hot gas introduced by the electric spark within a confined igniter cavity, which eventually evolves into a pulsed jet of a high-temperature kernel. A comprehensive set of models adapted from existing strategies for internal combustion engine spark plug discharge is extended to the target problem, including the modeling of energy deposition, plasma reactions, thermodynamic properties, and heat losses. A series of validation and parameter studies are performed and presented. The kernel size is found to be sensitive to heat losses arising from radiation and hot gas remained within the discharge cavity, rather than heat conduction to the wall in the discharge cavity. Depending on the enforced shape of the post-breakdown electric arc, the spark kernel can be off-centered, tilted, and considerably asymmetric. These features have been previously not considered when studying such igniter configurations and may have a first-order impact on the ignition process. Provided a proper setup of the heat loss models and electric arc shape, the numerical results are quantitatively comparable to the experimental results in terms of the kernel size, shape, and velocity throughout different stages after the spark discharge.

Tang, Yihao↗
DOE OSTI2021

Observation of ionization trends in a laboratory photoionized plasma experiment at Z

We report experimental and modeling results for the charge state distribution of laboratory photoionized neon plasmas in the first systematic study over nearly an order of magnitude range of ionization parameter ξ ∝ F/N e . The range of ξ is achieved by flexibility in the experimental platform to adjust either the x-ray drive flux F at the sample or the electron number density N e or both. Experimental measurements of photoionized plasma conditions over such a range of parameters enable a stringent test of atomic kinetics models used within codes that are applied to photoionized plasmas in the laboratory and astrophysics. From experimental transmission data, ion areal densities are extracted by spectroscopic analysis that is independent of atomic kinetics modeling. The measurements reveal the net result of the competition between photon-driven ionization and electron-driven recombination atomic processes as a function of ξ as it affects the charge state distribution. Results from radiation-hydrodynamics modeling calculations with detailed inline atomic kinetics modeling are compared with the experimental results. There is good agreement in the mean charge and overall qualitative similarities in the trends observed with ξ but significant quantitative differences in the fractional populations of individual ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of carbon pellets disruption mitigation in an NSTX-U plasma

Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted for an NSTX-U-like plasma to support the electromagnetic pellet injection concept (EPI). A carbon ablation model has been implemented in M3D-C1 and tested with available data. 2D simulations were conducted in order to estimate the amount of carbon needed to quench the plasma, finding that the content in a 1 mm radius vitreous carbon pellet (~3.2x10E20 atoms) would be enough if it is entirely ablated. 3D simulations were performed, scanning over pellet velocity and parallel thermal conductivity, as well as different injection directions and pellet concepts (solid pellets and shell pellets). The sensitivity of the thermal quench and other related quantities to these parameters has been evaluated. A 1 mm radius solid pellet only partially ablates at velocities of 300 m/s or higher, thus being unable to fully quench the plasma. To further enhance the ablation, approximations to an array of pellets and the shell pellet concept were also explored. 3D field line stochastization plays an important role in both quenching the center of the plasma and in heat flux losses, thus lowering the amount of carbon needed to mitigate the plasma when compared to the 2D case. This study constitutes an important step forward in `predict-first' simulations for disruption mitigation in NSTX-U and other devices, such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A hierarchal model for bacterial cell inactivation in solution by direct and indirect treatment using cold atmospheric plasmas

Cold atmospheric plasma devices have shown promise for a variety of plasma medical applications, including wound healing and bacterial inactivation often performed in liquids. In the latter application, plasma-produced reactive oxygen and nitrogen species (RONS) interact with and damage bacterial cells, though the exact mechanism by which cell damage occurs is unclear. Computational models can help elucidate relationships between plasma-produced RONS and cell killing by enabling direct comparison between dissimilar plasma devices and by examining the effects of changing operating parameters in these devices. In biological applications, computational models of plasma-liquid interactions would be most effective in design and optimization of plasma devices if there is a corresponding prediction of the biological outcome. In this work, we propose a hierarchal model for planktonic bacterial cell inactivation by plasma produced RONS in liquid. A previously developed reaction mechanism for plasma induced modification of cysteine was extended to provide a basis for cell killing by plasma-produced RONS. Results from the model are compared to literature values to provide proof of concept. Differences in time to bacterial inactivation as a function of plasma operating parameters including gas composition and plasma source configuration are discussed. Results indicate that optimizing gas-phase reactive nitrogen species production may be key in the design of plasma devices for disinfection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Achieving Higher Fidelity Pulsed-Power Simulations Using Advanced Gap Closure Modeling

This project has progressed several physics models in the EMPIRE plasma simulation code to achieve higher fidelity simulations of high-current diode operation in pulsed-power accelerators. In this report, we present details for the following major work products: (1) a set of verification problems covering all key processes involved in gap closure physics has been designed; this suite has facilitated feature vetting and overall model maturation, (2) a new EMPIRE exemplar has been developed: the Radiographic Integrated Test Stand 6 (RITS-6) diode, and (3) An exemplar for the Saturn accelerator exemplar was enabled by the models matured under this work to self-consistently simulate further into the diode pulse than previously possible (bipolar flow regime). These developments have lead to the highest confidence EMPIRE power flow predictions of the Saturn accelerator to date. Additionally, three modeling approaches for simulating electrode plasmas have been investigated. We report on these results and provide recommendations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Production and diffusion of H 2 O 2 during the interaction of a direct current pulsed atmospheric pressure plasma jet on a hydrogel

The interaction of cold atmospheric pressure plasma jets with hydrogels has been used as a model system to study the interaction of plasmas with tissues. In this study, we analyze the diffusion of reactive oxygen species (in particular H 2 O 2 ) and quantify the amount of plasma-produced H 2 O 2 species that penetrates into a gelatin hydrogel. We show that the diffusion constant of H 2 O 2 in 10% gelatin hydrogel is similar to its diffusion constant in water and that the production of H 2 O 2 in the hydrogel is significantly less than the production of H 2 O 2 in distilled water for the same plasma operation conditions suggesting that the scavenging of OH radicals at the plasma-gel interface significantly reduces the H 2 O 2 production.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

X-ray absorption spectroscopy measurements of radiatively ionized argon gas

X-ray absorption spectroscopy is a diagnostic tool that can characterize the temperature and ionization state of a plasma. This technique requires experiments to characterize the platform, careful data calibration, and comparison with atomic models to understand the plasma parameters. We performed ionizing radiation wave experiments at the OMEGA Laser Facility that used an ∼ 80 eV X-ray source to heat an Argon (Ar) gas cell at fill pressures of 3 atm. To diagnose the Ar plasma, we used a capsule backlighter offset 10 mm from the gas cell, to produce X-rays that were absorbed by the ionized Ar gas. The absorption analysis was calibrated using the significant line structure in the backlighter spectrum, which served as independent energy fiducials enabling definitive measurements of a 50-eV shift in the energy of the Ar K-edge due to ionization. We compare the measured absorption spectra to two independent atomic models, PrismSPECT and SCRAM, and show that the observed K-edge shift in the heated gas is consistent with ionization up to Ar 4+ and temperatures of 10 eV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

National Opacity Program: Analysis of Opacity-relevant X-ray Emission Spectra (Milestone Report ID 7121)

We have produced high energy density iron plasmas at temperatures above 1 keV and electron densities exceeding 1023 cm -3 (~ 1 g/cm 3 ) using the Orion laser at the Atomic Weapons Establishment. These plasmas were created by irradiating 50 µm diameter layered targets with frequency doubled (λ = 527 nm), 1 ps laser pulses focused to a 100 µm diameter producing an irradiance of ~ 2 x 10 18 W/cm 2 . The buried layer targets consist of 160 nm iron sulfide (FeS), 60 nm potassium chloride (KCl), and 15 nm carbon. The combined layers are tamped on both sides with 3 µm of parylene-N. The x-ray emission from the plasma was measured using two time-resolved, and four time-integrated Bragg crystal spectrometers, as well as one time-integrated imaging system. One time-resolved x-ray spectrometer measured emission from L-shell transitions in highly charged iron, the other from K-shell transitions in helium-like S 14+ and hydrogen-like S 15+ . The time-integrated spectrometers are intensity-calibrated and measured emission from K-shell transitions in sulfur, potassium, chlorine, and both K-shell and L-shell transitions in iron. The density and temperature of the plasma were determined by modeling the x-ray spectra using different spectral and hydrodynamic modeling packages. A brief overview of the uncertainties associated with the measurements and models are presented. We also give an overview of our 1-D HYDRA-DCA radiation-hydrodynamics model and improvements for future work. Our results aid in assessing experimental uncertainties associated with plasma uniformity and with x-ray emission employed as diagnostics in opacity experiments at temperatures and densities not achievable elsewhere and represent a significant step in creating and diagnosing plasmas near LTE. These results are summarized as part of the completion requirements for milestone 7121

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of high intensity orbital angular momentum beams for laser–plasma interactions

In this work, we explore the field of high intensity orbital angular momentum (OAM) beams, their generation with spiral phase mirrors, and the theory behind modeling both ideal and realistic beams. We explore OAM beam asymmetries introduced by aberrations in the beam, manufacturing defects, and bandwidth. A full three-dimensional description of the paraxial Laguerre–Gaussian (LG) modes suitable for modeling beams down to f/2 focal geometries is derived. A perturbative approach to modeling asymmetric OAM beams is introduced showing that only three LG modes are sufficient to model a wide variety of OAM asymmetries. The models are compared with experimental results followed by a discussion on the future of high intensity OAM beams in plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Symmetric Set of Transport Coefficients for Collisional Magnetized Plasma

Braginskii extended magnetohydrodynamics is used to model transport in collisional astrophysical and high energy density plasmas. We show that commonly used approximations to the α ⊥ and β ⊥ transport coefficients [e.g., Epperlein and Haines, Phys. Fluids 29 , 1029 (1986) ] have a subtle inaccuracy that causes significant artificial magnetic dissipation and discontinuities. This is because magnetic transport actually relies on β ∥ - β ⊥ and α ⊥ - α ∥ , rather than α ⊥ and β ⊥ themselves. We provide fit functions that rectify this problem and thus resolve the discrepancies with kinetic simulations in the literature. Furthermore, then implemented in the g orgon code, they reduce the predicted density asymmetry amplitude at laser ablation fronts. Recognizing the importance of α ⊥ - α ∥ and β ∥ - β ⊥ , we recast the set of coefficients. This makes explicit the symmetry of the magnetic and thermal transport, as well as the symmetry of the coefficients themselves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Controlled plasma–droplet interactions: a quantitative study of OH transfer in plasma–liquid interaction

Plasmas in contact with liquids are a rich source of OH radicals and have been extensively studied in the last decade to leverage the ability to generate chemically reactive species in gas phase plasmas to decompose organics. Multiphase transfer of OH radicals is highly transport limited and to overcome transport limits, the plasma activation of aerosols, small liquid droplets, interspersed in the plasma has been proposed. In this work, we report a combined experimental and modeling study of a controlled plasma–droplet interaction experiment using a diffuse RF glow discharge in He + 0.2% H 2 O with detailed plasma diagnostics, ex situ analysis of the plasma-induced chemistry in the droplet containing formate, droplet trajectory and size measurements. This enables a quantitative study of the reactivity transfer of OH from the gas phase plasma to the liquid phase and how its diffusion limitations impact formate decomposition in the water droplet. For a droplet with a diameter of 36 μm, we observed 50% reduction in formate concentration in the droplets after plasma treatment for droplet residence times in the plasma of ~10 ms. These short droplet residence times in the plasma allow in some cases for droplet size reductions of ~5% in spite gas temperatures of 360 K. A one-dimensional reaction–diffusion model was used to calculate the OH transport and formate oxidation inside the droplet and was able to predict the conversion of formate by plasma in a droplet without any fitting parameters. The model further shows that formate conversion is dominated by near-interfacial reactions with OH radicals and is limited by diffusion of formate in the droplet. Here, the results show that a controlled plasma–micro-droplet reactor as reported in this study might be an excellent tool for detailed quantitative plasma–liquid interaction studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Theory and simulations of cross-beam energy transfer between speckled laser beams

In this paper, we introduce a realistic modeling of laser beams smoothed by random phase plates to investigate cross-beam energy transfer (CBET). We compute the total field of the smoothed laser beams, which allows us to perform accurate simulations and to build a model accounting for the realistic structure of the speckles. This lets us confirm the results already discussed in a previous work given by Oudin et al. [Phys. Plasmas 29, 112112 (2022)]. When CBET is induced by a plasma flow, plane wave models can predict the exchange between smoothed laser beams. By contrast, when the exchange is caused by a wavelength shift between the laser fields of each beam, plane wave models overestimate the energy transfer, even for strongly Landau-damped ion acoustic waves. We also show that the laser angular spread, which results from its focusing, affects the resonance between the ponderomotive force induced by the laser fields and the ion acoustic wave. It increases the resonance width, which weakly depends on the Landau damping rate when there is a wavelength shift and no flow. In addition, the power transfer obtained with or without a phase plate is essentially the same in the perturbative regime and for small crossing angles.

Laser plasma interactions↗