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

Suppressing the enhancement of stimulated Raman scattering in inhomogeneous plasmas by tuning the modulation frequency of a broadband laser

The stimulated Raman scattering (SRS) instability can inhibit the performance of laser driven inertial confinement fusion (ICF) implosions by scattering light into unwanted directions or by generating hot electrons that preheat the target fuel. In principle, ICF target designs can avoid parameter regimes conducive to large, linear SRS gains. In practice, kinetic inflation—the nonlinear enhancement of SRS due to electron trapping in the excited plasma wave—makes this difficult. Here we show that laser bandwidth in the form of frequency modulation can either decrease or increase the inflationary SRS (iSRS) threshold in inhomogeneous plasmas depending on the maximum chirp of the laser pulse. The threshold, mapped out by a series of particle-in-cell simulations, exhibits a minimum when the frequency change within the pulse cancels the spatial detuning due to density inhomogeneities along the trajectory of the scattered light. In conclusion, by tuning the pump laser parameters away from this minimum, the iSRS threshold can be larger than at zero bandwidth, providing a path to mitigating kinetic inflation in ignition relevant plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Enhancing positron production using front surface target structures

We report a target design which produced a substantial gain in relativistic electron-positron pair production using high-intensity lasers and targets with large-scale micro-structures on their surface. Comparing to an unstructured target, a selected Si microwire array target yielded a near 100% increase in the laser-to-positron conversion efficiency and produced a 10MeV increase in the average emitted positron energy under nominally the same experimental conditions. We had established a multi-scale particle-in-cell simulation scheme to simulate both the laser absorption and the subsequent pair productions in a thick metal target. Here, the experimental results are supported by the simulations demonstrating the performance increase is due to a higher conversion efficiency of laser energy into electrons with kinetic energies greater than 10 MeV due to enhanced direct laser acceleration of electrons enabled by the microwire array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

How Hall electric fields intrinsically chaotize and heat ions during collisionless magnetic reconnection

The origin of extreme and anomalous ion heating during magnetic reconnection has been a longstanding mystery. It is shown in this paper that, at the core of the process, the stochastic heating mechanism is responsible for this energization. The in-plane Hall electric fields that intrinsically develop during collisionless reconnection not only accelerate ions, but also thermalize them by inducing chaotic motion. The mechanism explains various observed features of ion heating, e.g., its relation to the guide field strength and the ion mass-to-charge ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Plasma physics in strong-field regimes: Theories and simulations

In strong electromagnetic fields, unique plasma phenomena and applications emerge whose description requires recently developed theories and simulations [Y. Shi, Ph.D. thesis, Princeton University (2018)]. In the classical regime, to quantify effects of strong magnetic fields on three-wave interactions, a convenient formula is derived by solving the fluid model to the second order in general geometry. As an application, magnetic resonances are exploited to mediate laser pulse compression, whereby higher intensity pulses can be produced in wider frequency ranges, as confirmed by particle-in-cell simulations. In even stronger fields, relativistic-quantum effects become important, and a plasma model based on scalar quantum electrodynamics (QED) is developed which unveils observable corrections to Faraday rotation and cyclotron absorption in strongly magnetized plasmas. Beyond the perturbative regime, lattice QED is extended as a numerical tool for plasma physics, using which the transition from wakefield acceleration to electron-positron pair production is captured when laser intensity exceeds the Schwinger threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Liner implosion experiments driven by a dynamic screw pinch

This paper expands upon recent experimental results [Campbell et al., Phys. Rev. Lett. 125, 035001 (2020)], where thin-foil liner implosions were driven by a dynamic screw pinch (DSP) and found to have magneto-Rayleigh–Taylor instability (MRTI) amplitudes up to three times smaller than in implosions driven by a standard z-pinch (SZP). Here, the expanded discussion presented herein includes: (1) a detailed comparison of the MRTI growth measured in the experiment with that calculated from theory; (2) measurements of axial magnetic field injection into the liner interior prior to the implosion, as well as the subsequent compression of this field during the implosion; (3) an in-depth description of how the helical geometry of the DSP can result in earlier implosion and stagnation times relative to the SZP; and (4) particle-in-cell simulations showing different electron drift behavior in the anode–cathode gap of the DSP relative to the SZP, and how this difference may be related to the different current waveforms recorded during the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reflectivity and spectrum of relativistic flying plasma mirrors

Flying plasma mirrors induced by intense lasers via laser wake field acceleration scheme have been proposed as a promising way to generate few-cycle EUV or x-ray lasers. In addition, if such a relativistic plasma mirror can accelerate, then it would serve as an analog black hole to investigate the information loss paradox associated with the black hole Hawking evaporation. Among these applications, the reflectivity, which is usually frequency-dependent, would affect the outgoing photon spectrum and, therefore, impact on the analysis of the physics under investigation. In this paper, these two issues are investigated analytically and numerically with one-dimensional particle-in-cell simulations. Based on our simulation results, we propose a new model that provides a better estimate of the reflectivity than those studied previously. Moreover, we found that the peak frequency of the reflected spectrum of a Gaussian incident wave deviates from the expected value, 4γ 2 ω, due to the dependence of reflectivity on the frequency of the incident wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Energy-angle distribution of the ions in the RF sheath of ICRH antennas

Radio-frequency sheaths forming at the surface of ion cyclotron resonance heating (ICRH) actuators in fusion experiments are associated with enhanced impurity sputtering from ICRH plasma-facing components (PFCs). The minimization of impurity fluxes from the ICRH PFC is a vital task for the usage of ICRH systems. Capturing the ion kinetics is important for an accurate description of plasma–material interactions, because the ion dynamics plays a crucial role in RF sheaths. Here, we present a hybrid particle-in-cell (hPIC) model able to capture the kinetic behavior of the ions. We analyze the kinetic ion energy-angle distributions (IEADs) impacting the RF antenna and its dependence on different plasma and RF sheath parameters. In particular, the IEAD dependence on RF frequency and magnetic field alignment is investigated. Using hPIC, we simulated a case emulating the latest experimental campaign from JET. The simulation showed that under specific plasma and RF parameters, the kinetic motion of ions results in a cusp formation in the IEAD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2021

A particle-in-cell code comparison for ion acceleration: EPOCH, LSP, and WarpX

There are now more Particle-in-Cell (PIC) codes than ever before that researchers use to simulate intense laser-plasma interactions. To date, there have been relatively few direct comparisons of these codes in the literature, especially for relativistic intensity lasers interacting with thin overdense targets. To address this we perform a code comparison of three PIC codes: EPOCH, LSP, and WarpX for the problem of laser-driven ion acceleration in a 2D(3v) geometry for a 10 20 W cm -2 intensity laser. We examine the plasma density, ion energy spectra, and laser-plasma coupling of the three codes and find strong agreement. We also run the same simulation 20 times with different random seeds to explore statistical fluctuations of the outputs. We then compare the execution times and memory usage of the codes (without “tuning” to improve performance) using between 1 and 48 processors on one node. We provide input files to encourage larger and more frequent code comparisons in this field.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cross-beam energy transfer saturation by ion trapping-induced detuning

The performance of direct-drive inertial confinement fusion implosions relies critically on the coupling of laser energy to the target plasma. Cross beam energy transfer (CBET), the resonant exchange of energy between intersecting laser beams mediated by ponderomotively driven ion-acoustic waves (IAW), inhibits this coupling by scattering light into unwanted directions. Furthermore, the variety of beam intersection angles and varying plasma conditions in an implosion results in IAWs with a range of phase velocities. Here we show that CBET saturates through a resonance detuning that depends on the IAW phase velocity and that results from trapping-induced modifications to the ion distribution functions. For smaller phase velocities, the modifications to the distribution functions can rapidly thermalize in the presence of mid-Z ions, leading to a blueshift in the resonant frequency. For larger phase velocities, the modifications can persist, leading to a redshift in the resonant frequency. Ultimately, these results may reveal pathways towards CBET mitigation and inform reduced models for radiation hydrodynamics codes to improve their predictive capability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Lower-hybrid drift waves and their interaction with plasmas in a 3D symmetric reconnection simulation with zero guide field

We investigate lower-hybrid drift waves (LHDW) in symmetric magnetic reconnection with zero guide field using three-dimensional particle-in-cell simulations. The long-wavelength mode with kρiρe∼1 develops in the bifurcated electron current layer around the X-line within the width of the electron meandering motion from the mid-plane, where ρi(e) is the ion (electron) gyroradius. The short-wavelength mode with kρe∼1 develops in the separatrix region downstream of the electron outflow jet, producing electron vortices in the background flow frame. Electrons follow the E × B drift with corrections from the diamagnetic drift and are heated inside the vortices with diverging electric fields. In the vortices, ions have comparable E × B and inertia drifts, which together mostly cancel the diamagnetic drift. Toward the center of diverging field vortices, ions are decelerated, leading to a decrease in the perpendicular temperature, while the loss of low-energy ions results in an increase in the parallel temperature. Parallel electric fields exist as a combination of the LHDW wave field projected to the magnetic field direction and the penetration of whistler waves that are mainly outside of the LHDW layer. The magnetic flux tube is twisted in the vortices. The twist may potentially lead to slippage reconnection, as indicated by the non-uniform parallel potential variation across field lines, while the periodic variations of the twisting directions are a limiting factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evidence for Trapping-Induced Nonlinear Frequency Shifts in Langmuir Waves Driven via Stimulated Raman Scattering

Thomson scattering is used to detect the spectra of Langmuir waves driven through the backward stimulated Raman scattering process in a diffraction limited laser focal spot. Measured Langmuir wave spectral frequencies are found to vary in time and have broadened spectral power, consistent with a nonlinear frequency shift of the driven Langmuir wave due to electron-trapping. Broadening of the Langmuir wave spectral power is observed to decrease in time, consistent with measured variations in the frequency shift of the driven Langmuir waves. Furthermore, the observed spectral broadening is consistent with the temporally short (ps), bursty nature of backward stimulated Raman scattered light observed in simulations that cannot be resolved by the Thomson scattering diagnostic. Comparison of the broadened spectrum with time integrated spectra from two-dimensional particle-in-cell simulations shows favorable comparison in the broadened spectral widths, supporting the supposition of electron-trapping induced, nonlinear shifting of daughter Langmuir wave frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Whistler waves generated by nongyrotropic and gyrotropic electron beams during asymmetric guide field reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fermi acceleration in relativistic collisionless plasma shocks correlates with anisotropic energy gains

Collisionless shocks generated by two colliding relativistic electron-positron plasma shells are studied using particle-in-cell simulations. Shocks are mediated by the Weibel instability (WI), and the kinetic energy of the fastest accelerated particles is found to be anisotropically modified by WI-induced electric fields. Specifically, we show that all particles interacting with the shock bifurcate into two groups based on their final relativistic Lorentz factor γ: slow (γ < γ bf ) and fast (γ < γ bf ), where γ bf is the bifurcation Lorentz factor that was found to be approximately twice the initial (upstream) Lorentz factor γ0. We have found that the kinetic energies of the slow particles are equally affected by the longitudinal and transverse components of the shock electric field, whereas the fast particles are primarily accelerated by the transverse field component.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Absorption of relativistic multi-picosecond laser pulses in wire arrays

We study the interaction of intense multi-picosecond laser pulses with arrays of carbon wires attached to solid substrates. Here, we find that laser absorption in wire arrays resembles that in flat targets with very large uniform plasma density gradients. Performing two-dimensional particle-in-cell simulations, we optimize target parameters like wire thickness and -distance for energy absorption of a 2picosecond laser pulse with a large focal spot; this has implications for x-ray- and charged particle source development.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Simulations of ion heating due to ion-acoustic instabilities in presheaths

Particle-in-cell, direct simulation Monte Carlo simulations reveal that ion-acoustic instabilities excited in presheaths can cause significant ion heating. Ion-acoustic instabilities are excited by the ion flow toward a sheath when the neutral gas pressure is small enough and the electron temperature is large enough. A series of 1D simulations were conducted in which neutral plasma (electrons and ions) was uniformly sourced with an ion temperature of 0.026 eV and different electron temperatures (0.1 eV–50 eV). Ion heating was observed when the electron-to-ion temperature ratio exceeded the minimum value predicted by linear response theory to excite ion-acoustic instabilities at the sheath edge (T e /T i ≈ 28). When this threshold was exceeded, the temperature equilibration rate between ions and electrons rapidly increased near the sheath so that the local temperature ratio did not significantly exceed the threshold for instability. This resulted in significant ion heating near the sheath edge, which also extended back into the bulk plasma; presumably due to wave reflection from the sheath. Furthermore, this ion-acoustic wave heating mechanism was found to decrease for higher neutral pressures, where ion-neutral collisions damp the ion-acoustic waves and ion heating is instead dominated by inelastic collisions in the presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗