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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Nonlinear δf particle in cell gyrokinetic simulations of the microtearing mode

Nonlinear simulations are carried out for the microtearing mode using particle-based δf gyrokinetic simulations for parameters relevant to spherical tokamaks. The present study finds that the microtearing mode can generate significant electron heat flux, which is predominantly carried out by the electromagnetic component of the heat flux with a negligible contribution from the electrostatic component. The mode sustains without the electrostatic component. We observe that the electron heat flux increases with the electron temperature gradient. The heat flux exhibits a very weak dependence on the collisions. It increases with electron β initially; however, at very high β, the electron heat flux is reduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Generation of ultrafast electron bunch trains via trapping into multiple periods of plasma wakefields

We demonstrate a novel approach to the generation of femtosecond electron bunch trains via laser-driven wakefield acceleration. We use two independent high-intensity laser pulses, a drive, and an injector, each creating their own plasma wakes. The interaction of the laser pulses and their wakes results in a periodic injection of free electrons in the drive plasma wake via several mechanisms, including ponderomotive drift, wake-wake interference, and pre-acceleration of electrons directly by strong laser fields. Electron trains were generated with up to four quasi-monoenergetic bunches, each separated in time by a plasma period. In conclusion, the time profile of the generated trains is deduced from an analysis of beam loading and confirmed using 2D particle-in-cell simulations

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code x-point gyrokinetic code (XGC), in order to study density pump-out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics code M3D-C1, taking into account the linear two-fluid plasma response. With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate and suppression of the electron heat flux in the steepest part of the edge pedestal. In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Laser-driven collisionless shock acceleration of ions from near-critical plasmas

Here, this paper overviews experimental and numerical results on acceleration of narrow energy spread ion beams by an electrostatic collisionless shockwave driven by 1 μm (Omega EP) and 10 μm (UCLA Neptune Laboratory) lasers in near critical density CH and He plasmas, respectively. Shock waves in CH targets produced high-energy ~50 MeV protons (ΔE/E of ≤30%) and 314 MeV C 6+ ions (ΔE/E of ≤10%). Observation of acceleration of both protons and carbon ions to similar velocities is consistent with reflection of particles off the moving potential of a shock front. For shocks driven by CO 2 laser in a gas jet, ~30 MeV peak in He ion spectrum was detected. Particle-in-cell simulations indicate that regardless of the target further control over its density profile is needed for optimization of accelerated ion beams in part of energy spread, yield and maximum kinetic energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Birefringence in thermally anisotropic relativistic plasmas and its impact on laser–plasma interactions

One of the paradigm-shifting phenomena triggered in laser–plasma interactions at relativistic intensities is the so-called relativistic transparency. As the electrons become heated by the laser to relativistic energies, the plasma becomes transparent to the laser light even though the plasma density is sufficiently high to reflect the laser pulse in the non-relativistic case. This paper highlights the impact that relativistic transparency can have on laser-matter interactions by focusing on a collective phenomenon that is associated with the onset of relativistic transparency: plasma birefringence in thermally anisotropic relativistic plasmas. Here, the optical properties of such a system become dependent on the polarization of light, and this can serve as the basis for plasma-based optical devices or novel diagnostic capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Three-dimensional cross-field flows at the plasma-material interface in an oblique magnetic field

This article describes experimental evidence that the magnetic presheath is a fully three-dimensional structure modified by ion–neutral collisions. Velocity distributions of both ions and neutrals, obtained via laser-induced fluorescence, show that cross field ion drifts do not result from entrainment of ions in a flowing neutral background. Ion flows parallel to E x B arise and accelerate to as much as 0.2c s within several ion gyroradii of the boundary surface, where c s is the sound speed. Within measurement resolution, the onset of the E x B aligned flow occurs at the same distance to the surface that ions begin to deflect from travel along magnetic field lines. Collisional fluid and particle-in-cell simulations of the boundary region are compared to the experimental measurements. We find that, in contrast to the classical collisionless Chodura model, collisional effects between the ions and the non-flowing neutral population are essential to quantitatively predict the observed ion drift velocities. No momentum coupling between ions and neutrals, separable from noise and other effects, is observed in either signal. Here, we discuss several explanations and implications of this observation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetic reconnection and kinetic waves generated in the Earth's quasi-parallel bow shock

Magnetic reconnection in quasi-parallel shocks, relevant to the Earth’s bow shock, is studied by means of two-dimensional full particle-in-cell simulations. As the Alfvénic Mach number increases, the propagation direction of the waves excited in the transition region changes, and the shock becomes more turbulent with more reconnection sites. In the higher Mach number shock, abundant electron-only reconnection sites are generated with scales on the order of the ion skin depth or less. Non-reconnecting current sheets can also generate electron jets and energy dissipation can occur there as well. However, non-reconnecting current sheets with the magnetic field reversal typically show a smaller energy dissipation rate than reconnecting current sheets. In the shock transition region, two types of waves are responsible for driving reconnection: one has a wavelength around three ion skin depths (d i ), and the other has a wavelength less than 1 d i . Here, electron and ion distribution functions show that in regions where the former type of waves is excited, there are two ion beams and a single-peaked electron distribution. In contrast, in regions where the latter type of waves is excited, there are multiple electron and ion beams. The waves propagating obliquely to the magnetic field bend the magnetic field lines, and magnetic reconnection occurs where oppositely directed field lines come into contact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Target normal sheath acceleration with a large laser focal diameter

In this work, the dependence of the laser-driven ion acceleration from thin titanium foils in the Target Normal Sheath Acceleration (TNSA) regime on target and laser parameters is explored using two dimensional particle-in-cell simulations. The oblique incidence (θ L = 45°) and large focal spot size (w 0 = 40μm) are chosen to take an advantage of quasi one-dimensional geometry of sheath fields and effective electron heating. This interaction setup also reveals low and achromatic angular divergence of a proton beam. It is shown that the hot electron temperature deviates from the ponderomotive scaling for short laser pulses and small pre-plasmas. This deviation is mainly due to the laser sweeping, as the short duration laser pulse each moment in time effectively heats only a fraction of a focal spot on the foil. This instantaneous partial heating results in an electron temperature deviation from the ponderomotive scaling and, thus, lower maximum proton energies than it could have been expected from the TNSA theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nonthermal electron and ion acceleration by magnetic reconnection in large laser-driven plasmas

Magnetic reconnection is a fundamental plasma process that is thought to play a key role in the production of nonthermal particles associated with explosive phenomena in space physics and astrophysics. Experiments at high-energy-density facilities are starting to probe the microphysics of reconnection at high Lundquist numbers and large system sizes. We have performed particle-in-cell (PIC) simulations to explore particle acceleration for parameters relevant to laser-driven reconnection experiments. We study particle acceleration in large system sizes that may be produced soon with the most energetic laser drivers available, such as at the National Ignition Facility. In these conditions, we show the possibility of reaching the multi-plasmoid regime, where plasmoid acceleration becomes dominant. Our results show the transition from X point to plasmoid-dominated acceleration associated with the merging and contraction of plasmoids that further extend the maximum energy of the power-law tail of the particle distribution for electrons. We also find for the first time a system-size-dependent emergence of nonthermal ion acceleration in driven reconnection, where the magnetization of ions at sufficiently large sizes allows them to be contained by the magnetic field and energized by direct X point acceleration. For feasible experimental conditions, electrons and ions can attain energies of ϵ max , e / k B T e > 100 and ϵ max , i / k B T i > 1000 . Using PIC simulations with binary Monte Carlo Coulomb collisions, we study the impact of collisionality on plasmoid formation and particle acceleration. The implications of these results for understanding the role reconnection plays in accelerating particles in space physics and astrophysics are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Similarity of capacitive radio-frequency discharges in nonlocal regimes

Similarity transformations are essential for correlating discharges at different scales, which are mostly utilized with local field or local energy approximations. In this work, we report the fully kinetic results from particle-in-cell/Monte Carlo collision simulations that unambiguously demonstrate the similarity of radio frequency (rf) discharges in nonlocal regimes where the electron energy relaxation length is much larger than the gap dimension. It is found that at a constant rf voltage amplitude, discharges will be similar if the gas pressure, inverse of gap distance, and rf driving frequency are all changed by the same scaling factor. The scaling relations of fundamental parameters are illustrated for rf discharges in the alpha-mode with secondary electron emission ignored, and the temporal electron kinetics are shown to have invariance in similar discharges. The results explicitly validate the scaling laws in nonlocal kinetic regimes, indicating promising application potentials of the similarity transformations across a wide range of kinetic regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Spatial core-edge coupling of the particle-in-cell gyrokinetic codes GEM and XGC

Two existing particle-in-cell gyrokinetic codes, GEM for the core region and XGC for the edge region, have been successfully coupled with a spatial coupling scheme at the interface in a toroidal geometry. Additionally, a mapping technique is developed for transferring data between GEM's structured and XGC's unstructured meshes. Two examples of coupled simulations are presented to demonstrate the coupling scheme. The optimization of GEM for graphics processing unit is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles

In this paper, we present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evolution of the self-injection process in long wavelength infrared laser driven LWFA

In this work, long wavelength infrared laser-driven plasma wakefield accelerators are investigated in the self-modulated laser wakefield acceleration (SM-LWFA) and blowout regimes using 3D particle-in-cell simulations. The simulation results show that in the SM-LWFA regime, self-injection arises with wave breaking, whereas in the blowout regime, self-injection is not observed under the simulation conditions. The wave breaking process in the SM-LWFA regime occurs at a field strength that is significantly below the 1D wave-breaking threshold. This process intensifies at higher laser power and plasma density and is suppressed at low plasma densities (≤1×10 17 cm –3 here). The produced electrons show spatial modulations with a period matching that of the laser wavelength, which is a clear signature of direct laser acceleration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of a chain of three plasma accelerator stages with the WarpX electromagnetic PIC code on GPUs

The fully electromagnetic particle-in-cell code WarpX is being developed by a team of the U.S. DOE Exascale Computing Project (with additional non-U.S. collaborators on part of the code) to enable the modeling of chains of tens to hundreds of plasma accelerator stages on exascale supercomputers, for future collider designs. The code is combining the latest algorithmic advances (e.g., Lorentz boosted frame and pseudo-spectral Maxwell solvers) with mesh refinement and runs on the latest computer processing unit and graphical processing unit (GPU) architectures. In this paper, we summarize the strategy that was adopted to port WarpX to GPUs, report on the weak parallel scaling of the pseudo-spectral electromagnetic solver, and then present solutions for decreasing the time spent in data exchanges from guard regions between subdomains. In Sec. IV, we demonstrate the simulations of a chain of three consecutive multi-GeV laser-driven plasma accelerator stages.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental and calculational investigation of laser-heated additive manufactured foams

Foam materials are starting to find application in laser-heated Hohlraums used to drive inertial confinement fusion implosions. Foams made using additive manufacturing (AM) techniques are now available and may have advantages over traditional chemical (aerogel) foams. Here, we present new experimental data on laser-heated AM foams. Samples of four different types of printed AM foams were heated using a single 527 nm laser beam at the Jupiter Laser Facility. The laser pulse was ∼180 J square pulse with an FWHM of 1.6 ns and a peak intensity of 3–4 × 1014 W/cm2. The foam densities ranged from 12 to 93 mg/cc (all supercritical for 527 nm light). We measured the backscattered light (power and spectrum), the transmitted light, side-on x-ray images, and the Ti K-shell emission that was used to infer the time-integrated temperature. The fraction of backscattered light was 6%–15% of the input laser energy. The pure carbon foam sample had less backscatter than a C8H9O3 foam of similar density, which was consistent with multi-fluid calculations that predicted less ion heating for the C8H9O3 foam. The level of backscatter and the thermal front speeds for the AM foams were similar to values measured for stochastic (aerogel) foams under similar conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transition characteristics and electron kinetics in microhollow cathode discharges

We demonstrate the transition characteristics and electron kinetics of microdischarges in a microgap with a cathode having microstructures using two-dimensional particle-in-cell/Monte Carlo collision (2d3v) simulations. It is found that the microdischarge is inside the hollow cavity at lower pressures, forming microhollow cathode discharges (MHCDs), while the dominant discharge moves outside the cavity at higher pressures, locating above the cathode rectangular protrusion tip. The spatial distributions of the microdischarge parameters (e.g., electric potential and electron density) are presented by capturing the transition characteristics. The electron kinetics of the microdischarges are examined based on the moment analysis of the electron Boltzmann equation. The collisional and non-collisional components of the electron power absorption are decomposed, which were seldomly investigated for microdischarges previously. The results provide fundamental insights into MHCD formation with a structured electrode at varying pressure conditions, which could be beneficial for the design and fabrication of microplasma devices in practical applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Laser-plasma acceleration beyond wave breaking

Laser wakefield accelerators rely on the extremely high electric fields of nonlinear plasma waves to trap and accelerate electrons to relativistic energies over short distances. When driven strongly enough, plasma waves break, trapping a large population of the background electrons that support their motion. Aside from limiting the maximum electric field, this trapping can lead to accelerated electron bunches with large energy spreads. Here we introduce a novel regime of plasma wave excitation and wakefield acceleration that allows for arbitrarily high electric fields while avoiding the deleterious effects of unwanted trapping. Here, the regime, enabled by spatiotemporal shaping of laser pulses, exploits the property that nonlinear plasma waves with superluminal phase velocities cannot trap charged particles and are therefore immune to wave breaking.

36 MATERIALS SCIENCE↗