Search NASA⌕ Search

SEARCH · Search NASA

Results for “PIC simulation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Multiscale Studies of the Three-Dimensional Dayside X-Line

We review recent experience from the Cluster, Double Star, and THEMIS missions for lessons that apply to the upcoming Magnetospheric Multiscale Mission (MMS) being developed for launch in 2014. On global scales, simulation and statistical studies lead to mean congurations of dayside reconnection, implying specific relative alignments of the inflow magnetic fields and X-line, with implications for MMS operations designed to maximize the number of close encounters with the diffusion region. At intermediate MHD-to-ion scales, reconstruction of features created by one or two X-lines have developed to the point where data from a cluster of spacecraft can determine their temporal trends and the approximate three-dimensional X-line structure. Recent petascale particle-in-cell (PIC) simulations of reconnection encompass three spatial dimensions with excellent resolution, and make striking predictions of electron scale physics that creates complex interacting flux ropes under component reconnection. High time resolution measurements from MMS will determine the detailed electron scale kinetics embedded within the global and MHDion scale contexts. These developments will lead to the refinement of our three-dimensional multiscale picture of reconnection, yielding improved understanding of the global, MHD, and local physics controlling the onset or quenching, variability, and mean rate of reconnection. This in turn will enable improved predictability of the structural features created by transient reconnection, and their space weather consequences.

Reconnection↗
DOE OSTI2021

Measuring Saturn's Electron Beam Energy Spectrum using Webb's Wedges

It is very difficult to measure the voltage of the load on the Saturn accelerator. Time-resolved measurements such as vacuum voltmeters and V-dot monitors are impractical at best and completely change the pulsed power behavior at the load at worst. We would like to know the load voltage of the machine so that we could correctly model the radiation transport and tune our x-ray unfold methodology and circuit simulations of the accelerator. Step wedges have been used for decades as a tool to measure the end - point energies of high energy particle beams. Typically, the technique is used for multi-megavolt accelerators, but we have adapted it to Saturn's modest <2 MV end-point energy and modified the standard bremsstrahlung x-ray source to extract the electron beam without changing the physics of the load region. We found clear evidence of high energy electrons >2 MV. We also attempted to unfold an electron energy spectrum using a machine learning algorithm and while these results come with large uncertainties, they qualitatively agree with PIC simulation results.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Numerical simulation of terahertz radiation by laser-driven plasma dipole oscillation

Laser-plasma interactions have many theoretical and technological applications. One is the use of coherently accelerated electrons to provide novel sources of THz radiation. Our research focuses on simulating the cross/self-interactions between two high intensity, ultra-short, counter propagating and detuned laser pulses and an initial neutral target for controlled ionization. Unlike our previous studies of laser-matter interaction over preformed plasma, we explore the injection and collision of laser pulses to induce background plasma driven by the self-guided laser wakefield mechanism, which is then used to perturb the plasma resulting in induced dipole oscillations leading to transverse radiation. Inducing a cylindrical spatial plasma column within the laser beam radius regime provides a stable, spatially localized plasma channel. The emitted radiation from the plasma dipole oscillation (PDO) will not be affected by surrounding plasma absorption, resulting in effective radiation distribution. Finally, results include 3D EM-PIC simulations and a comparison of the self- ionizing plasma against the preformed plasma to assess the efficiency of the mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Laser Plasma Interaction in the MG Magnetic Field (Final Report)

Strong magnetic fields change properties of plasma. Employment of the magnetic fields can improve conditions for inertial confinement fusion. Compression and heating of the magnetized plasma is a base of the MagLIF pulsed power approach to the controlled fusion. Study of expansion and heating of plasma in the external magnetic fields, development of plasma instabilities, parametric and other effects are important for basic physics and applications. The main goal of the research for the grant was experimental investigation of laser produced plasma in the well characterized and controlled MG magnetic fields. We developed a platform for investigation of high intensity laser-plasma interaction (LPI) in the MG fields using a university-scale pulsed power machine. A Zebra machine at UNR generated transverse magnetic fields of 3 MG on the surface of the rod load and longitudinal fields of 1.4 MG in the coil loads. The magnetic fields were characterized with two-color Faraday rotation diagnostics. The pulsed power machine provided almost static magnetic fields for LPI compared to the much faster laser driven magnetic fields. We found that eddy currents in the rising magnetic field generate plasma on the metal targets. Eddy currents did not impact Si and CH dielectric targets. A Leopard laser with intensity in the focus up to 10 19 W/cm 2 produced plasma in the MG magnetic fields. The modified reflected Faraday diagnostics was used to study laser driven B-field at the MTW laser at Laboratory for Laser Energetics, UR. We demonstrated a slow dependence of the magnetic field in the coil on the applied laser intensity. The magnetic field increased by a factor of 2 if laser intensity increased by a factor of 30. This scaling is important for comparison of experiments at different laser facilities. Expansion of the plasma in the 2-3 MG external B-field was studied. Laser-produced plasma in the azimuthal magnetic field took the unique form of a thin disc expanding radially with a velocity of 250 km/s and confined in the vertical direction. A HYDRA MHD program at the University of Rochester was applied for simulations of LPI in the magnetic field B >3 MG. Simulations for the plasma disc parameters and the expansion dynamics were in a good agreement with experiments. Generation of narrow dense plasma jets in the longitudinal magnetic field of 0.6-0.8 MG was studied. Narrow plasma jets reached a length of 3-4 mm with the electron plasma density of (0.2-1.2)x10 20 cm -3 . A jet tip propagated with the velocity of 160-200 km/s. MHD simulations showed a good agreement of the dynamics of the formation of plasma jets with experiments. These jets are relevant to astrophysical jets. Simulations of K-shell spectra of Si plasma was performed for LPI in the B-field. The laser prepulse generated plasma near the laser target. A PrismSPECT model with a MeV electron beam produced by the laser was used. The spectral modelling showed the increased plasma density by a factor by 2-3 in the magnetic field. This was in agreement with the observed dynamics of the plasma plume confined by the axial magnetic field. The two-plasmon decay was studied in the 2.5-3 MG transversal magnetic fields. Strong 2-3.5 nm spectral widening and a 2-4 nm shifts of “red” and “blue” 3/2ω 0 spectral components were observed. The large shift and widening exceed the expected temperature and magnetic shift. PIC simulations are performing to clarify physics of this effect. The experimental research program for the grant is completed. New plasma effects in the MG fields were studied at the university-scale pulsed power generator. Two graduate and two undergraduate UNR and UR students carried experiments at the Zebra generator, supported and developed plasma diagnostics, provided data processing, and performed MHD simulations of LPI in the strong magnetic field. Research results are published in 9 referred papers and presented in 12 conferences. Program manager: Kramer U. Akli; Collaborators: R. C. Mancini and H. Sawada, University of Nevada, Reno (UNR); R. Betti and A. V. Maximov, University of Rochester (UR).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm

A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different -space operators in Faraday’s and Ampere’s law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. Here, the new algorithm was implemented into Osiris by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, eA / mc 2 , exceeding for typical cell sizes and time steps.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Improved C 1 shape functions for simplex meshes

Using continuously-differentiable (C 1 ) interpolation methods previously developed for computer visualization and other applications, a new class of shape functions is constructed and applied to electrostatic (ES) Particle-In-Cell (PIC) simulations on simplex (triangle or tetrahedral) meshes. The resulting shape functions are non-negative and partition unity,and lead to a reduction of shot noise and aliasing associated with particle interpolation to and from the mesh. For an energy-conserving implementation as chosen, C 1 interpolation of the potential leads to a continuous (C 0 ) electric field, with associated improved fidelity of particle motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Electron energization dynamics in interaction of self-generated magnetic vortices in upstream of collisionless electron/ion shocks

Abstract Relativistic collisionless shocks are considered responsible for particle energization mechanisms leading to particle acceleration. While electron energization in shock front region of electron/ion collisionless shocks are the most studied, the mechanism of electron energization in interaction with self-generated magnetic vortices (MVs) in the upstream region is still unclear. We investigate electron energization mechanism in the upstream region of electron/ion relativistic collisionless shocks, using two dimensional particle-in-cell (PIC) simulations. We discuss mechanism of electron energization which takes place in the upstream region of the shock, where the counter stream particles interact with incoming flow. The energy gain of electrons happens during their interaction with evolving fields of self-generated magnetic vortices in this region. Three Fermi-like electron energization scenarios are discussed. Stochastic acceleration of electrons in interaction with fields of MV leads to anisotropic heating of fast electrons due to diffusion in the momentum space of electrons and, finally, synergetic effect of evolving fields of MVs leads to the formation of a power-law tail of supra-thermal particles.

Naseri, N.↗
DOE OSTI2025

Laser based 100 GeV electron acceleration scheme for muon production

High energy muons, due to their unique ability to penetrate deeply into matter, can enable radiography of structures that cannot be probed by other forms of radiation. Current terrestrial sources of muons require conventional GeV-TeV particle accelerators which are hundreds to thousands of meters in size. Laser wakefield acceleration (LWFA) can achieve acceleration gradients of two-to-three orders of magnitude greater than conventional accelerators, thus shrinking the accelerator to a number of meters. We propose a concept for a compact muon source based on the first self-consistent PIC simulations of an all optical LWFA that uses a guiding channel to achieve electron energies of 100 GeV in a distance of 6 m with a driving laser energy of 300 J in a single stage. From the resulting electron energy spectrum we estimate muon production for this source. We show that this accelerator, coupled with high average power laser driver technology, provides the basis for a high energy and high flux muon source.

laser-produced plasmas↗
DOE OSTI2020

Self-acceleration and energy channeling in the saturation of the ion-sound instability in a bounded plasma

A novel regime of the saturation of the Pierce-type ion-sound instability in bounded ion-beam-plasma system is revealed in 1D PIC simulations. It is found that the saturation of the instability is mediated by the oscillating virtual anode potential structure. The periodically oscillating potential barrier separates the incoming beam ions into two groups. One component forms a supersonic beam which is accelerated to an energy exceeding the energy of the initial cold ion beam. The other component is organized as a self-consistent phase space structure of trapped ions with a wide energy spread – the ion hole. The effective temperature (energy spread) of the ions trapped in the hole is lower than the initial beam energy. In the final stage the ion hole expands over the whole system length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Suppression of power losses during laser pulse propagation in underdense plasma slab

For current state-of-the-art terawatt lasers, the primary laser scattering mechanisms in plasma include Forward Raman Scattering (FRS), excitation of plasma waves, and the self-modulational instability (SMI). Using 2D PIC simulations, we demonstrate the dominance of the FRS in the regime with medium-to-low density plasma and non-relativistic laser fields. However, the use of multi-colored lasers with frequency detuning exceeding the plasma frequency, Δω > ω pe , suppresses the FRS. Finally, the laser power can then be transmitted efficiently.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Proton beam enhancements from ultrafast laser interactions with compound parabolic concentrators

Target normal sheath acceleration (TNSA) is among the most widely studied laser–plasma ion acceleration mechanisms. In this Letter, we report on studies of proton acceleration from flat Cu targets with cone-like focusing structures called compound parabolic concentrators (CPCs) coupled to their front surface. The CPC acts as a non-imaging focusing optic that enhances the laser intensity at the target's front surface by ∼2×, from 9 × 10 18 to 2 × 10 19 W/cm 2 . This effect drives a reduction in the particle source size from the typical laser spot size of 120 μm to the CPC tip size of 65 μm and significant enhancements in the peak energy and temperature (2.6 × increase) of the resulting TNSA proton beam. We also observe increased opening angle of the beam. 2D PIC simulations have been conducted and replicate the experimental behavior. These increases have implications for long focal length facilities that require higher energy proton beams, which now can be achieved without significant infrastructure changes.

3D printing↗
DOE OSTI2020

Transport of kJ-laser-driven relativistic electron beams in cold and shock-heated vitreous carbon and diamond

We report experimental results on relativistic electron beam (REB) transport in a set of cold and shock-heated carbon samples using the high-intensity kilojoule-class OMEGA EP laser. The REB energy distribution and transport were diagnosed using an electron spectrometer and x-ray fluorescence measurements from aCu tracer buried at the rear side of the samples. The measured rear REB density shows brighter and narrower signals when the targets were shock-heated. Hybrid PIC simulations using advanced resistivity models in the target warm-dense-matter (WDM) conditions confirm this observation. We show that the resistivity response of the media, which governs the self generated resistive fields, is of paramount importance to understand and correctly predict the REB transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dynamic focusing of laser driven positron jets by self-generated fields

Abstract Focusing effect of laser-driven positron jets by self-generated target sheath fields has been observed for the first time experimentally and the results are supported by the computational studies. In the experiment, OMEGA EP short-pulse (0.7 ps, 500 J) irradiates mm-size gold targets with a concave back surface and reference flat-surface targets. Both targets exhibited positrons with quasi-monoenergetic energy peaks while targets with concave curvature also showed increased number of positrons at the detector. The data is consistent with hybrid-PIC simulations confirming that the time-varying electric fields driven by electrons escaping from the target significantly change the trajectories of positrons. These simulations show a small radius of curvature on the rear side increases the relative focusing effect and the positrons to electrons ratio in the escaping plasma. For the smallest radius of curvature, positron jets that are up to 10 times denser can be achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Plasma instabilities dominate radioactive transients magnetic fields: the self-confinement of leptons in Type Ia and core-collapse supernovae, and kilonovae

The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $\beta ^+$ and $\beta ^-$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material.

instabilities↗
DOE OSTI2021

Electron acceleration at oblique angles via stimulated Raman scattering at laser irradiance >10 16 W cm –2 μm 2

The generation of hot, directional electrons via laser-driven stimulated Raman scattering (SRS) is a topic of great importance in inertial confinement fusion (ICF) schemes. Little recent research has been dedicated to this process at high laser intensity, in which back, side, and forward scatter simultaneously occur in high energy density plasmas, of relevance to, for example, shock ignition ICF. In this paper, we present an experimental and particle-in-cell (PIC) investigation of hot electron production from SRS in the forward and near-forward directions from a single speckle laser of wavelength λ 0 = 1.053μm, peak laser intensities in the range I 0 = 0.2–1.0×10 17 Wcm –2 and target electron densities between ne = 0.3–1.6% n c , where n c is the plasma critical density. As the intensity and density are increased, the hot electron spectrum changes from a sharp cutoff to an extended spectrum with a slope temperature T = 34 ± 1 keV and maximum measured energy of 350 keV experimentally. Multidimensional PIC simulations indicate that the high energy electrons are primarily generated from SRS-driven electron plasma wave phase fronts with k vectors angled ~50° with respect to the laser axis. These results are consistent with analytical arguments that the spatial gain is maximized at an angle which balances the tendency for the growth rate to be larger for larger scattered light wave angles until the kinetic damping of the plasma wave becomes important. The efficiency of generated high energy electrons drops significantly with a reduction in either laser intensity or target electron density, which is a result of the rapid drop in growth rate of Raman scattering at angles in the forward direction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Transport of an intense proton beam from a cone-structured target through plastic foam with unique proton source modeling

Laser-accelerated proton beams are applicable to several research areas within high-energy density science, including warm dense matter generation, proton radiography, and inertial confinement fusion, which all involve transport of the beam through matter. Here, we report on experimental measurements of intense proton beam transport through plastic foam blocks. The intense proton beam was accelerated by the 10 ps, 700 $\textit{J}$ OMEGA EP laser irradiating a curved foil target, and focused by an attached hollow cone. The protons then entered the foam block of density 0.38 g/cm 3 and thickness 0.55 or 1.00 mm. At the rear of the foam block, a Cu layer revealed the cross section of the intense beam via proton- and hot electron-induced Cu-K α emission. Images of x-ray emission show a bright spot on the rear Cu film indicative of a forward-directed beam without major breakup. 2D fluid-PIC simulations of the transport were conducted using a unique multi-injection source model incorporating energy-dependent beam divergence. Along with postprocessed calculations of the Cu – $K_α$ emission profile, simulations showed that protons retain their ballistic transport through the foam and are able to heat the foam up to several keV in temperature. The total experimental emission profile for the 1.0 mm foam agrees qualitatively with the simulated profile, suggesting that the protons indeed retain their beamlike qualities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗