Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle method”

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 937 records · Page 52

Measuring characteristic differences between high- and low-performing discharges on the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) DPF

A dense plasma focus (DPF) is a compact coaxial plasma gun, which completes its discharge as a Z-pinch, producing short (<100 ns) pulses of ions, x rays, and/or neutrons. Lawrence Livermore National Laboratory recently constructed and began operating a new device, the MJOLNIR (MegaJOuLe Neutron Imaging Radiography) DPF, which is designed for single-pulse flash neutron radiography. This device has achieved neutron yields of up to 4.1×10 11 neutrons/pulse at 3.3 MA peak current, and higher-current commissioning is under way. Like most DPFs, MJOLNIR exhibits variable yields in some configurations. Here, we present evidence of the role of parasitic current paths within the gun in stochastically influencing the yield. First through “conditioning shots,” where new hardware has been introduced, we show that increased run-down and run-in speeds correlate with higher yields. These observations are consistent with current being delivered to the electrodes but not to the main plasma sheath, degrading the implosion-driving force. Once nominal conditions are established, we correlate low-performing discharges with smaller current dip and associated voltage spike for a fixed machine configuration. A snow-plow model is able to recreate small-magnitude current dips through the introduction of a parasitic current path, and particle-in-cell simulations establish how parasitic current paths lower the ion beam energy available to produce neutrons. Finally, we observe an increased likelihood of shots with low yield and smaller current dip with increasing fill pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗

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↗

Investigating the effects of electron bounce-cyclotron resonance on plasma dynamics in capacitive discharges operated in the presence of a weak transverse magnetic field

Recently, Patil et al. [Phys. Rev. Res. 4, 013059 (2022)] have reported the existence of an enhanced operating regime when a low-pressure (5 mTorr) capacitively coupled discharge (CCP) is driven by a very high radio frequency (60 MHz) source in the presence of a weak external magnetic field applied parallel to its electrodes. Their particle-in-cell simulations show that a significantly higher bulk plasma density and ion flux can be achieved at the electrode when the electron cyclotron frequency equals half of the applied radio frequency for a given fixed voltage. In the present work, we take a detailed look at this phenomenon and further delineate the effect of this “electron bounce-cyclotron resonance (EBCR)” on the electron and ion dynamics of the system. We find that the ionization collision rate and stochastic heating are maximum under resonance condition. The electron energy distribution function also indicates that the population of tail-end electrons is highest for the case where EBCR is maximum. Formation of electric field transients in the bulk plasma region is also seen at lower values of applied magnetic field. Finally, we demonstrate that the EBCR-induced effect is a low-pressure phenomenon and weakens as the neutral gas pressure increases. The potential utility of this effect to advance the operational performance of CCP devices for industrial purposes is discussed in this report.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High repetition rate mapping of the interaction between a laser plasma and magnetized background plasma via laser induced fluorescence

In this work, the laminar coupling of energy between a laser-produced plasma and a background magnetized plasma was investigated via planar laser induced fluorescence diagnostic and magnetic flux probes. Experiments performed on the Large Plasma Device at the University of California, Los Angeles, mapped out the two-dimensional spatiotemporal evolution of the laser-plasma (debris) ion velocity distribution function (VDF) to assess debris-background coupling in a sub-Alfvénic regime. The acquisition of these data necessitates high repetition rate (1 Hz) as each dataset is the accumulation of thousands of laser shots, which would not be feasible in single-shot experiments. Fully kinetic, three-dimensional particle-in-cell simulations are compared to the measured VDFs to provide a framework in which we can understand the coupling of a sub-Alfvénic plasma flow through a preformed, magnetized plasma. The simulations display the same departure from the expected gyromotion of the debris plasma as observed in the experimental data, and in conjunction with the measured magnetic field traces, have led to the direct observation of the collisionless coupling via laminar fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of collimated vortex gamma-rays from intense Poincaré beam–plasma interaction

We report on numerical calculations in which a multi-petawatt γ-ray beam is generated using a novel configuration based on fully structured light irradiating an overdense plasma waveguide. We analyze how the relativistic laser pulse efficiently confines and accelerates plasma electrons to GeV-scale energies and drives a quasi-static field that induces magneto-bremsstrahlung radiation. Multiphoton Compton scattering of electrons in the intense part of the laser also occurs although the radiated energy-density is comparatively lower. The emitted γ-rays carry orbital angular momentum, are highly collimated, and account for upwards of 15% of the incident field energy in one particular case. A comparison of the laser-to-particle angular momentum and energy transfer efficiencies is made between the cases of irradiation by a circularly polarized Laguerre–Gauss mode and one type of full Poincaré beam, and it is found that the latter yields an order-of-magnitude enhancement. In conclusion, the essential characteristics of the interaction are validated with three-dimensional particle-in-cell simulations that include quantum electrodynamical effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma electron contribution to beam emittance growth from Coulomb collisions in plasma-based accelerators

Coulomb collisions with background plasma can cause emittance growth in plasma accelerators. This paper extends the theory to consider collisions with not only motionless plasma ions but also plasma electrons with relativistic motion, based on the Frankel cross section. The theory is verified by particle-in-cell simulations with a Monte Carlo collision module. Further, it is shown that the electron contribution has the same amount as that of ions in linear acceleration regime and may not be negligible in nonlinear regime depending on the plasma electron density and its relativistic bulk velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excitation of long electrostatic solitary waves in ion beam neutralization process

Unusually long electrostatic solitary waves (ESWs) are discovered in a particle-in-cell simulation study of the process of ion beam neutralization by electron emission from a filament. These ESWs are long because the density perturbation responses to the potential wells created by the ESWs are very small. Notably, the density perturbation is small because the trapped (positive) and untrapped (negative) electron density perturbations nearly compensate each other because of a non-Maxwellian electron velocity distribution in the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pressure–strain interaction. III. Particle-in-cell simulations of magnetic reconnection

How energy is converted into thermal energy in weakly collisional and collisionless plasma processes, such as magnetic reconnection and plasma turbulence, has recently been the subject of intense scrutiny. The pressure–strain interaction has emerged as an important piece, as it describes the rate of conversion between bulk flow and thermal energy density. In two companion studies, we presented an alternate decomposition of the pressure–strain interaction to isolate the effects of converging/diverging flow and flow shear instead of compressible and incompressible flow, and we derived the pressure–strain interaction in magnetic field-aligned coordinates. Here, we use these results to study pressure–strain interaction during two-dimensional anti-parallel magnetic reconnection. We perform particle-in-cell simulations and plot the decompositions in both Cartesian and magnetic field-aligned coordinates. We identify the mechanisms contributing to positive and negative pressure–strain interaction during reconnection. Furthermore, this study provides a roadmap for interpreting numerical and observational data of the pressure–strain interaction, which should be important for studies of reconnection, turbulence, and collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced electron acceleration by high-intensity lasers in extended (confined) preplasma in cone targets

Here we report on experimental results from a high-intensity laser interaction with cone targets that increase the number (×3) and temperature (×3) of the measured hot electrons over a traditional planar target. This increase is caused by a substantial increase in the plasma density within the cone target geometry, which was induced by 17 ± 9 mJ prepulse that arrived 1.5 ns prior to the main high intensity (>10 19 W/cm 2 ). Three-dimensional hydrodynamic simulations are conducted using hydra which show that the cone targets create substantially longer and denser plasma than planar targets due to the geometric confinement of the expanding plasma. The density within the cone is a several hundred-micron plasma “shelf” with a density of approximately 10 20 n e /cc. The HYDRA simulated plasma densities are used as the initial conditions for two-dimensional particle-in-cell simulations using EPOCH. These simulations show that the main acceleration mechanism is direct-laser-acceleration, with close agreement between experimentally measured and simulated electron temperatures. Further analysis is conducted to investigate the acceleration of the electrons within the long plasma generated within a compound parabolic concentrator by the prepulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗