Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Kinetic modeling of solitary wave dynamics in a neutralizing ion beam

In this work, we characterize the formation and evolution of electrostatic solitary waves (ESWs) in the space-charge neutralization of ion beams using particle-in-cell simulations. These waves become excited when the electrons emitted from an external filament source initiate a two-stream instability in the beam. We show that such electrostatic waves become excited in both two-dimensional (2D) and three-dimensional (3D) beams with different shapes and sizes. Through a 1D Bernstein–Greene–Kruskal (BGK) analysis of the 2D beam, we find that the non-Maxwellian nature of the beam electrons gives rise to large-sized ESWs that are not predicted by BGK theory since it assumes a Maxwellian electron velocity distribution in the beam. Finally, we show that a 1D BGK theory is inadequate to describe ESWs in 3D beams because of complex electron trajectories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The effects of pre-plasma scale length on the relativistic electron beam directionality

The effects of pre-plasmas on the electron beam directionality was experimentally and numerically investigated. Single material and layered targets made of Ti and/or CH were used to simultaneously measure high-energy (≥3 MeV) electrons along two directions, pre-pulse energy and pre-plasma density. The electron directionality is quantified by using a new parameter, the electron energy ratio of the total kinetic energies along the two directions. Measurements and radiation–hydrodynamic (RH) simulations show that a large (≥3.5 μm) plasma scale length at the critical surface enhances electrons along the laser axis, and such pre-plasma conditions could only be achieved with the CH targets. Particle-in-cell simulations were performed on the RH generated pre-plasmas from Ti and CH targets, and the results show that the CH target provided conditions for higher forward momentum gains by electrons. First, the CH target allowed longer distances for electrons to interact with laser. Second, the intense laser pulse modified the critical surface, but the resulting surface differed. The CH target resulted in a smooth surface where a retro-reflection was observed while the Ti target resulted in a rippled surface that scattered the reflected light. As results, the CH electrons gained higher forward momentum via a direct-laser-acceleration in the counter propagating laser fields. So the results presented in this article show a way of controlling the high-energy electron directionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Proton acceleration in an overdense hydrogen plasma by intense CO 2 laser pulses with nonlinear propagation effects in the underdense pre-plasma

We report on proton acceleration from intense CO 2 laser-irradiated hydrogen plasmas at near-critical densities, with the density gradient steepened by Nd:YAG laser ablation-driven hydrodynamic shocks. While the experimental results, such as the quasi-monoenergetic proton spectra and their scaling with respect to the laser energy, are generally in agreement with the simulations, certain laser shots produced significantly higher proton energies than anticipated during the experiment. The increased proton energy may be linked to nonlinear propagation effects in the steepened plasma density ramp before the critical surface, including relativistic self-focusing and, for the case of temporally-structured laser pulses observed in the experiment, focusing of the trailing pulse through the plasma channel formed by the leading pulse 25 ps ahead. The occurrence of channel focusing in the underdense hydrogen plasma is supported by a subsequent pump-probe experiment with a dark-field imaging technique, where the formation of ion channels was observed after the passage of an intense CO 2 laser pulse

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Investigation of the significance of excited species in wave-heated dense non-ideal microplasmas using particle-in-cell Monte Carlo-collision modeling

A computational model for the interaction of a high-pressure microplasma with an electromagnetic wave is presented. A one-dimensional particle-in-cell Monte Carlo collision model is used to investigate the plasma non-ideality effects in a second-stage laser-heated xenon plasma with a comprehensive chemistry mechanism, including excited species. A xenon microplasma at a temperature of 300 K and pressure of 10 bar becomes non-ideal after the laser heating stage with lower ionization and electronic excitation barriers, which makes the presence of excited species important to the plasma generation process. For these conditions, two-step ionization via excited species becomes a major contributor to plasma ionization. We find that full ionization of the plasma is obtained in ~2 ps when excited species are included in the chemistry model, as opposed to ~2.5 ps for a chemistry mechanism that does not consider their effect. With excited species, the ionization mode transitions from direct ionization via electron-neutral collisions to direct and stepwise ionization as the plasma generation progresses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electron-field instability: Excitation of electron plasma waves by an electric field

Electric fields are commonplace in plasmas and affect transport by driving currents and, in some cases, instabilities. The necessary condition for instability in collisionless plasmas is commonly understood to be described by the Penrose criterion, which quantifies a sufficient relative drift between different populations of particles that must be present for wave amplification via inverse Landau damping. For example, electric fields generate drifts between electrons and ions that can excite the ion-acoustic instability. Here, we use particle-in-cell simulations and linear stability analysis to show that the electric field can drive a fundamentally different type of kinetic instability, named the electron-field instability. This instability excites electron plasma waves with wavelengths ≳30λ De , has a growth rate that is proportional to the electric field strength, and does not require a relative drift between electrons and ions. The Penrose criterion does not apply when accounting for the electric field. Furthermore, the large value of the observed frequency, near the electron plasma frequency, further distinguishes it from the standard ion-acoustic instability, which oscillates near the ion plasma frequency. The ubiquity of macroscopic electric fields in quasineutral plasmas suggests that this instability is possible in a host of systems, including low-temperature and space plasmas. In fact, damping from neutral collisions in such systems is often not enough to completely damp the instability, adding to the robustness of the instability across plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Cross-beam energy transfer in conditions relevant to direct-drive implosions on OMEGA

In cross-beam energy transfer (CBET), the interference of two laser beams ponderomotively drives an ion-acoustic wave that coherently scatters light from one beam into the other. This redirection of laser beam energy can severely inhibit the performance of direct-drive inertial confinement fusion (ICF) implosions. To assess the role of nonlinear and kinetic processes in direct-drive-relevant CBET, the energy transfer between two laser beams in the plasma conditions of an ICF implosion at the OMEGA laser facility was modeled using particle-in-cell simulations. For typical laser beam intensities, the simulations are in excellent agreement with linear kinetic theory, indicating that nonlinear processes do not play a role in direct-drive implosions. At higher intensities, CBET can be modified by pump depletion, backward stimulated Raman scattering, or ion trapping, depending on the plasma density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Investigation of boosted proton energies through proton radiography of target normal sheath acceleration fields in the multi-ps regime

Multi-kilojoule, multi-picosecond short-pulse lasers, such as the National Ignition Facility-Advanced Radiographic Capability laser and the OMEGA-Extended Performance laser, which have been constructed over the last two decades, enable exciting opportunities to produce high-brightness, high-energy laser-driven proton sources for applications in high-energy-density science like proton fast ignition for inertial fusion energy, particle radiography, and materials science studies. Results on these platforms have demonstrated enhanced accelerated proton energies and electron temperatures when compared to established scaling laws. Recent work has developed a new scaling for proton TNSA in the multi-ps regime. However, this new physics in the multi-ps regime motivates the need to understand the origin of the enhancement in proton energies. Toward this goal, here this work presents the first measurements of the TNSA accelerating sheath field in the multi-ps regime for pulse durations of 0.6, 5, and 10 ps. This measurement was achieved by using a separate TNSA proton source to radiograph the spatiotemporal profile of the accelerating sheath that is responsible for proton acceleration. The use of stacked radiochromic film detectors allows for a discrete time profile of the radiographs, thus enabling the measurement of the temporal and spatial evolution of the accelerating field. In performing this measurement, we extract quantities such as the sheath strength as a function of time and pulse duration, which shows that longer pulse durations sustain a stronger electric field for a longer duration when compared to sub-ps laser pulses, which may enable the observed boosted proton energies and proton conversion efficiencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Kinetic simulation of a 50 mTorr capacitively coupled argon discharge over a range of frequencies and comparison to experiments

The effect of driving frequency in the range of 13.56–73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct “temperatures” for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2–10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25% of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model, and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A gas dynamic perspective on particle lifting in electrostatic discharge-like devices using multiphase particle-in-cell approach

Although particle–laden electrostatic discharges are widely used in laboratories as well as in industrial applications, the mechanism of particle lifting for particles initially at rest in such highly unsteady systems is not well understood. A multiphase gas–particle solver is developed using the multiphase particle-in-cell (MP-PIC) approach to emulate the interaction of a compressible shock-dominated gas phase with the dense particle phase. First, the two-way coupled solver is initially used to study the interaction of a planar traveling shock with a vertical curtain of particulates. The gas and particle phase evolution was found to be in good agreement with a similar experimental study in Ling et al. [Phys. Fluids 24, 113301 (2012)]. Second, the MP-PIC code is used to study the interaction of an expanding blast wave with a thick bed of particles. The simulation considered forces such as quasi-steady drag, pressure-gradient, added-mass, Saffman, and Magnus forces. We observe that the vertical liftoff particles close to the shock impingement point in this configuration are associated with the quasi-steady drag, pressure gradient, and added-mass forces. Also, the Saffman lift and Magnus forces contribute to lifting particles located radially farther away from the shock impingement point. In addition, the study finds a decrease in particle lifting efficiency with decreasing plasma kernel length and shock strength.

42 ENGINEERING↗
DOE OSTI2023

Faraday effect in collisional magnetized plasmas

Faraday rotation is a valuable diagnostic tool for investigating laboratory and astrophysical plasmas, but it has almost exclusively been treated in the collisionless limit (despite the non-negligible role of Coulomb collisions in many laboratory plasmas). Here, we show that Faraday rotation can occur in collisional plasmas and that the usual effect is curiously altered by collisions. Namely, an initially linearly polarized light wave—propagating parallel to a uniform magnetic field—becomes elliptically polarized in the collisional plasma, and collisional absorption of the wave also occurs. Moreover, the ellipticity and rotation angle are quantifiably sensitive to the fidelity of the collisional transport coefficients. As we will demonstrate with particle-in-cell EPOCH simulations, these effects offer a stringent diagnostic tool for benchmarking multi-scale plasma simulation codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Evolution and hot electron generation of laser–plasma instabilities in direct-drive inertial confinement fusion

A series of 2D in-plane plane wave particle-in-cell simulations find distinctive paths of laser-plasma instability evolution in OMEGA-scale implosions, depending on the initial electron temperature. At low temperatures, two-plasmon decay (TPD) dominates in both initial growth and the steady state. At high temperatures, the initial dominant modes switch to stimulated Raman scattering, but TPD still dominates a steady state characterized by cavitation and Langmuir turbulence. A hot electron scaling is also obtained from the simulations that, when combined with laser/plasma conditions from hydro simulations, can predict hot electron generation in implosions that do not employ smoothing-by-spectral-dispersion (SSD). It also shows that under the same laser/plasma conditions, SSD can reduce hot electron generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Similarity-based scaling networks for capacitive radio frequency discharge plasmas

We demonstrate similarity-based scaling networks for capacitive radio frequency (RF) plasmas, which extensively correlate discharge characteristics under varied conditions, incorporating the transition from original to similarity states. Based on fully kinetic particle-in-cell simulations, similar RF discharges in argon are demonstrated with three external control parameters (gas pressure, gap distance, and driving frequency) simultaneously tuned. A complete set of scaling pathways regarding fundamental discharge parameters is obtained, from which each plasma state finds its neighboring node with only one control parameter tuned. The results from this study provide a promising strategy for plasma multi-parameter mapping, enabling effective cross-comparisons, prediction, and manipulation of RF discharge plasmas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the force exerted on a non-spherical asymmetric dust grain from homogeneous, stationary, isotropic, non-magnetized plasma

Whereas the conventional wisdom suggests that the force between non-magnetized homogeneous, stationary, isotropic plasma, and the dust grain is only possible for the case of relative plasma–grain velocity, it is shown that stationary non-spherical asymmetric dust grain immersed in stationary, non-magnetized, isotropic plasma can experience a force caused by the grain–plasma interactions. The component of the force due to scattering of plasma particles in the limit of infinite Debye length is considered analytically. Both the particle scattering and absorption force components are modeled numerically in the limits of infinite and finite Debye length using a newly developed 2D3V Aspherical Particle-in-Cell code. The code simulates interactions of dust grain of selected non-spherical asymmetric shape with plasmas using dust shape conforming coordinates. Here, the simulations confirm the existence of the force on non-spherical asymmetric grain in stationary non-magnetized plasma and show that the plasma screening effects can lead to reversal of the force direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On parallel laser beam merger in plasmas

Self-focusing instability is a well-known phenomenon of nonlinear optics, which is of great importance in the field of laser–plasma interactions. Self-focusing instability leads to beam focusing and, consequently, breakup into multiple laser filaments. The majority of applications tend to avoid a laser filamentation regime due to its detrimental role on laser spot profile and peak intensity. In our work, using nonlinear Schrödinger equation solver and particle-in-cell simulations, we address the problem of interaction of multiple parallel beams in plasmas. We consider both non-relativistic and moderately relativistic regimes and demonstrate how the physics of parallel beam interaction transitions from the familiar self- and mutual-focusing instabilities in the non-relativistic regime to a moderately relativistic regime, where an analytical description of filament interaction is not available.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Laser-to-proton conversion efficiency studies for proton fast ignition

We assess the conversion efficiency from intense picosecond laser pulses to multi-MeV ion beams for a wide range of laser and target parameters, using 2D kinetic particle-in-cell simulations. Scalings are addressed in a quasi-one-dimensional geometry, leaving out beam divergence. Then, we study the conversion efficiency into a narrow spatial band along the laser axis for flat and hemispherical targets in large-scale 2D simulations. Combining these findings allows us to calculate the energy requirements for ignition of a compressed ICF target with an intense proton beam in a fast-ignition scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Review and meta-analysis of electron temperatures from high-intensity laser–solid interactions

The accelerated electron spectrum from high-intensity laser–solid interaction is often conveniently described using a Boltzmann distribution, whose temperature is known within the field as the hot-electron temperature. The importance of the electron temperature is highlighted by the sheer number of experimental and simulation studies on the subject over the past three decades. Recently, multi-kJ, multi-ps pulses have yielded electron spectra with temperatures far beyond the expected ponderomotive result. Expressions that predict the electron temperature considering laser parameters beyond intensity and wavelength have been developed, albeit using small datasets. In this review, we present what is, to the best of our knowledge, the largest dataset of electron temperatures gathered from experimental measurements and particle-in-cell simulations. This dataset allows us to compare existing analytical and empirical hot-electron temperature scaling models over a wide parameter range. We also develop new scaling models that incorporate the laser pulse duration of the laser and the plasma scale length. Three models that include pulse-duration and scale length dependence are especially successful at predicting both simulated and experimental data. The dataset will soon be made publicly available to encourage further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Intermittency, bursty turbulence, and ion and electron phase-space holes formation in collisionless current-carrying plasmas

In the previous studies of nonlinear saturation of the Buneman instability caused by high electron drift velocity relative to ions, the phase-space holes and the plateau on the electron velocity distribution function were identified as features of the saturation stage of instability [notably in the paper by Omura et al., J. Geophys. Res. 108, 1197 (2003)]. We have performed a much longer simulation of the Buneman instability and observed a secondary instability. This secondary instability generates fast electron-acoustic waves. By analyzing the phase-space plot of ions and electrons, we show that the fast electron heating and the formation of the plateau of electron velocity distribution function are not due to the quasi-linear diffusion but due to the nonlinear interaction of ion- and electron-acoustic solitary waves (phase-space holes) by exchange of trapped electrons in each wave. We also report the details on the intermittent and bursty nature of turbulence driven by this instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗