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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2021

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A study of space charge induced non-linearity in the Single Line Of Sight camera

We report a new generation of gated x-ray detectors at the National Ignition Facility has brought faster, enhanced imaging capabilities. Their performance is currently limited by the amount of signal they can be operated with before space charge effects in their electron tube start to compromise their temporal and spatial response. We present a technique to characterize this phenomenon and apply it to a prototype of such systems, the Single Line Of Sight camera. Results of this characterization are used to benchmark Particle In Cell simulations of the electrons drifting inside the detector, which are found to well reproduce experimental data. These simulations are then employed to predict the optimum photon flux to the camera, with the goal to increase the quality of the images obtained on an experimental campaign while preventing the appearance of deleterious effects. They also offer some insight into some of the improvements that can be brought to the new pulse-dilation systems being built at Lawrence Livermore National Laboratory

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Theory, observations, and simulations of kinetic entropy in a magnetotail electron diffusion region

In this work, we examine velocity-space kinetic entropy, a spatially local measure of entropy for systems out of thermal equilibrium, during an encounter of an electron diffusion region at a magnetic reconnection site in Earth’s magnetotail by the Magnetospheric Multiscale (MMS) mission. We start by generalizing the theory of kinetic entropy to the case of non-uniform velocity space grids and transforming the equations into spheri- cal energy coordinates useful to experimental plasma detectors. The theory is then applied to MMS data and compared to particle-in-cell simulations of reconnection. We demonstrate that the entropy-based non-Maxwellianity measure from the MMS data is of sufficiently high precision to reliably identify non-Maxwellian distributions and therefore the measurements when kinetic effects are most significant. By comparing two different non-Maxwellian measures, we show that total entropy density suffers from “information loss” because it lacks a dependence on the velocity space grid, and so has lost information about how well a distribution function is resolved. Local velocity-space kinetic entropy density recovers this information. We quantify information loss and argue that the considerations needed to minimize it are crucial for instruments designed to measure distribution functions in situ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hybrid particle-in-cell simulations of electromagnetic coupling and waves from streaming burst debris

Various systems can be modeled as a point-like explosion of ionized debris into a magnetized, collisionless background plasma-including astrophysical examples, active experiments in space, and laser-driven laboratory experiments. Debris streaming from the explosion parallel to the magnetic field may drive multiple resonant and non-resonant ion-ion beam instabilities, some of which can efficiently couple the debris energy to the background and may even support the formation of shocks. We present a large-scale hybrid (kinetic ions + fluid electrons) particle-in-cell simulation, extending hundreds of ion inertial lengths from a 3D explosion, that resolves these instabilities. We show that the character of these instabilities differs notably from the 1D equivalent by the presence of unique transverse structure. Additional 2D simulations explore how the debris beam length, width, density, and speed affect debris-background coupling, with implications for the generation of quasi-parallel shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Strong reconnection electric fields in shock-driven turbulence

Turbulent magnetic reconnection in a quasi-parallel shock under parameters relevant to the Earth's bow shock is investigated by means of a two-dimensional particle-in-cell simulation. The addressed aspects include the reconnection electric field, the reconnection rate, and the electron and the ion outflow speeds. In the shock transition region, many current sheets are generated in shock-driven turbulence, and electron-only reconnection and reconnection where both ions and electrons are involved can occur in those current sheets. The electron outflow speed in electron-only reconnection shows a positive correlation with the theoretical speed, which is close to the local electron Alfvén speed, and a strong convection electric field is generated by the large electron outflow. As a result, the reconnection electric field becomes much larger than those in the standard magnetopause or magnetotail reconnection. In shock-driven reconnection that involves ion dynamics, both electron outflows and ion outflows can reach of the order of 10 times the Alfvén speed in the X-line rest frame, leading to a reconnection electric field the same order as that in electron-only reconnection. An electron-only reconnection event observed by the magnetospheric multiscale mission downstream of a quasi-parallel shock is qualitatively similar to those in the simulation and shows that the outflow speed reaches approximately half the local electron Alfvén speed, supporting the simulation prediction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Toward the core-edge coupling of delta-f and total-f gyrokinetic models

The coupling of core reduced delta-f and edge total-f gyrokinetic models could enable a significant speed-up of the whole device modeling simulations. It is thus interesting to compare these models in the core where they will be coupled. The main difference between these models is the presence of the zeroth order term on the right-hand side of the total-f gyrokinetic equation. This term is responsible for driving a radial electric field that is also associated with Geodesic acoustic mode-like oscillations. To investigate the coupling between these models, the subsequent large-scale gyrokinetic 3D turbulence simulations have been performed starting from an already saturated gyrokinetic axisymmetric equilibria (collisions are not included). This work has also been extended to couple different axisymmetric and turbulent models available in XGC. It is shown that the radial electric field and its drive have to be taken care of consistently while coupling different models together.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. I. Nonlinear and kinetic effects

Results are presented from a series of simulations examining the susceptibility of the cross-beam energy transfer (CBET) instability to nonlinear processes in the context of direct-drive inertial confinement fusion experiments on the OMEGA laser facility. These form the basis for the second paper of this series [A. G. Seaton, L. Yin, R. Follett, B. J. Albright, and A. Le, “Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth,” Phys. Plasmas 29, 042707 (2022)], where we examine the efficacy of increases in laser bandwidth at suppressing CBET. We choose laser and plasma conditions for the simulations that are favorable to CBET and promote nonlinearity. Through a comparison of outputs from the particle-in-cell code vector particle in cell (VPIC) and the linearized fluid code laser-plasma simulation environment (LPSE), a series of nonlinear effects have been identified in the kinetic simulations that include particle trapping, the two-ion wave decay, and ion-acoustic wave self-focusing. These effects produce time-dependent energy transfer, in contrast to the linearized fluid simulations in which a steady state is reached after an initial transient. Ion trapping is shown to allow for increased energy transfer relative to fluid simulations, with the remaining nonlinear processes acting to reduce the energy transfer. Nonlinear dynamics is contrasted for low- and high-intensity beams as well as between speckled and planar beams. For the parameters under consideration, beam profile has a significant effect on nonlinear dynamics, though the greatest sensitivity is to beam intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth

The response of the cross-beam energy transfer instability (CBET) to laser bandwidth is investigated through a combination of theory and simulation. Existing linear theory is generalized to treat broadband lasers, demonstrating that CBET is most effectively suppressed when the bandwidth exceeds the ion-acoustic wave (IAW) frequency. It is shown that for such bandwidths, reverse (seed to pump) transfer becomes possible, which reduces the net energy transfer rapidly as bandwidth is increased. The CBET gain exponent in this regime scales with bandwidth (Δω) as Δω−3 for Gaussian or Lorentzian laser spectra with different scalings possible for other spectra. Comparison of our theory with linearized fluid and particle-in-cell simulations, performed with the laser-plasma simulation environment (LPSE) and vector particle in cell (VPIC) codes, respectively, finds that the model is accurate in the absence of nonlinear processes. However, linear analysis also finds that the IAW energy density scales as Δω−1, implying that nonlinear effects may be more difficult to control than the CBET scaling would suggest. Indeed, nonlinear effects are found to be present in VPIC simulations with high-intensity lasers, despite minimal apparent CBET. Nonlinear processes in the VPIC cases include particle trapping, the two-ion wave decay, and ion wave self-focusing. In some high intensity VPIC cases, these effects lead to net energy transfer from seed to pump and increases to backscatter stimulated Brillouin scattering reflectivities. Finally, for a given bandwidth, we show that improved control of nonlinear processes can be achieved via smoothing by spectral dispersion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A scintillator attenuation spectrometer for intense gamma-rays

A new type of compact high-resolution high-sensitivity gamma-ray spectrometer for short-pulse intense gamma-rays (250 keV to 50 MeV) has been developed by combining the principles of scintillators and attenuation spectrometers. The first prototype of this scintillator attenuation spectrometer (SAS) was tested successfully in Trident laser experiments at LANL. Later versions have been used extensively in the Texas Petawatt laser experiments in Austin, TX, and more recently in OMEGA-EP laser experiments at LLE, Rochester, NY. The SAS is particularly useful for high-repetition-rate laser applications. Furthermore, we give a concise description of the design principles, capabilities, and sample preliminary results of the SAS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Electron-only reconnection and associated electron heating and acceleration in PHASMA

In this work, using incoherent Thomson scattering, electron heating and acceleration at the electron velocity distribution function (EVDF) level are investigated during electron-only reconnection in the PHAse Space MApping (PHASMA) facility. Reconnection arises during the merger of two kink-free flux ropes. Both push and pull type reconnection occur in a single discharge. Electron heating is localized around the separatrix, and the electron temperature increases continuously along the separatrix with distance from the X-line. The local measured gain in enthalpy flux is up to 70% of the incoming Poynting flux. Notably, non-Maxwellian EVDFs comprised of a warm bulk population and a cold beam are directly measured during the electron-only reconnection. The electron beam velocity is comparable to, and scales with, electron Alfvén speed, revealing the signature of electron acceleration caused by electron-only reconnection. The observation of oppositely directed electron beams on either side of the X-point provides “smoking-gun” evidence of the occurrence of electron-only reconnection in PHASMA. 2D particle-in-cell simulations agree well with the laboratory measurements. The measured conversion of Poynting flux into electron enthalpy is consistent with recent observations of electron-only reconnection in the magnetosheath [Phan et al., Nature 557, 202 (2018)] at similar dimensionless parameters as in the experiments. The laboratory measurements go beyond the magnetosheath observations by directly resolving the electron temperature gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Understanding the surface wave characteristics using 2D particle-in-cell simulation and deep neural network

Here, the characteristics of the surface waves along the interface between a plasma and a dielectric material have been investigated using kinetic particle-in-cell simulations. A microwave source of GHz frequency has been used to trigger the surface wave in the system. The outcome indicates that the surface wave gets excited along the interface of plasma and the dielectric tube and appears as light and dark patterns in the electric field profiles. The dependency of radiation pressure on the dielectric permittivity and supplied input frequency has been investigated. Further, we assessed the capabilities of neural networks to predict the radiation pressure for a given system. The proposed deep neural network model is aimed at developing accurate and efficient data-driven plasma surface wave devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Laser-driven, ion-scale magnetospheres in laboratory plasmas. II. Particle-in-cell simulations

Ion-scale magnetospheres have been observed around comets, weakly magnetized asteroids, and localized regions on the Moon and provide a unique environment to study kinetic-scale plasma physics, in particular in the collision-less regime. In this work, we present the results of particle-in-cell simulations that replicate recent experiments on the large plasma device at the University of California, Los Angeles. Using high-repetition rate lasers, ion-scale magnetospheres were created to drive a plasma flow into a dipolar magnetic field embedded in a uniform background magnetic field. The simulations are employed to evolve idealized 2D configurations of the experiments, study highly resolved, volumetric datasets, and determine the magnetospheric structure, magnetopause location, and kinetic-scale structures of the plasma current distribution. We show the formation of a magnetic cavity and a magnetic compression in the magnetospheric region, and two main current structures in the dayside of the magnetic obstacle: the diamagnetic current, supported by the driver plasma flow, and the current associated with the magnetopause, supported by both the background and driver plasmas with some time-dependence. From multiple parameter scans, we show a reflection of the magnetic compression, bounded by the length of the driver plasma, and a higher separation of the main current structures for lower dipolar magnetic moments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inviscid damping of an elliptical vortex subject to an external strain flow

Inviscid spatial Landau damping is studied experimentally for the case of oscillatory motion of a two-dimensional vortex about its elliptical equilibrium in the presence of an applied strain flow. Here, the experiments are performed using electron plasmas in a Penning–Malmberg trap. They exploit the isomorphism between the two-dimensional Euler equations for an ideal fluid and the drift-Poisson equations for the plasma, where plasma density is the analog of vorticity. Perturbed elliptical vortex states are created using E x B strain flows, which are generated by applying voltages to electrodes surrounding the plasma. Measurements of spatial Landau damping (also called critical-layer damping) are in agreement with previous studies in the absence of an applied strain, where the damping is due to a resonance between the local fluid motion and the vortex oscillations. Interestingly, the damping rate does not change significantly over a wide range of applied strain rates. This can be accurately predicted from the initial vorticity profile, even though the resonant frequency is reduced substantially due to the applied strain. For higher amplitude perturbations, nonlinear trapping oscillations also exhibit behavior similar to the strain-free case. In principle, higher-order effects of the applied strain, such as separatrix crossing of peripheral vorticity and interactions with harmonics of the fundamental resonance, are expected to change the damping rate. However, this occurs only for conditions that are not realized in the experiments described here. Vortex-in-cell simulations are used to investigate the possible roles of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stable electron beam propagation in a plasma column

The stability of plasma-based accelerators against transverse misalignments and asymmetries of the drive beam is crucial for their applicability. Without stabilizing mechanisms, even small initial offsets of the drive beam centroid can couple coherently to the plasma wake, grow, and ultimately lead to emittance degradation or beam loss for a trailing witness beam. In this work, we demonstrate the intrinsic stability of a beam propagating in a plasma column. This result is relevant in the context of plasma-based positron acceleration, where a wakefield suitable for the transport and acceleration of a positron witness beam is generated in a plasma column by means of an electron drive beam. The stable propagation of the drive beam is a necessary condition for the experimental implementation of this scheme. The differences and similarities of stabilizing mechanisms in a plasma column compared to a homogeneous plasma are identified via theory and particle-in-cell simulations. Experimental tolerances are given, demonstrating the experimental feasibility of the scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗