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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Plasma gradient effect on direct laser acceleration

The transfer of a high-intensity laser pulse energy to a high-energy electron beam via the direct laser acceleration mechanism is shown to be significantly enhanced through control of the plasma density gradient. Experiments performed using the OMEGA EP facility's high-intensity beams altered the plasma density and gradients by changing the Mach number and the angle of the gas-jet nozzle to the laser axis. When a long density gradient at the rear of the target is used, the total high-energy electron number measured was enhanced by 4.5 times compared to a shorter rear gradient. Complementary two-dimensional simulations, which follow the laser field evolution and the corresponding electron dynamics, strongly support the key trends observed in the experiment. The effect is twofold, the long density gradient provides the longest acceleration distance while it minimizes the formation of the sheath field as the electron beam exits into the vacuum. This study shows the importance of tailoring the plasma density.

Laser plasma interactions↗
DOE OSTI2025

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗
DOE OSTI2025

Neutron generation dynamics inside a MA-class dense plasma focus Z-pinch

Dense plasma focii (DPFs) are appealing as energy efficient sources of short pulses of ions, neutrons, and x rays. The output of these sources is expected to scale with input current (I4), but has been shown to drop at the MA level [S. K. H. Auluck, “On the failure of neutron yield scaling in the dense plasma focus,” Phys. Plasmas 30, 080701 (2023)]. New results on the MegaJOuLe Neutron Imaging Radiography DPF showed neutron yield production in agreement with the input current scaling beyond the previously observed drop. This work provides insight into the pinch formation on a DPF and reports on the two different mechanisms leading to neutron generation inside a DPF using a combination of kinetic simulations and experimental data. A combination of particle-in-cell (PIC) and 1D shock theory results are used to describe the pinch formation and disassembly and the corresponding thermonuclear and beam-target mechanisms. The temporal evolution of the pinch column predicted by the PIC simulations shows qualitative agreement with the experimental data from plasma photon emission as well as temporal neutron pulse shapes. In MJ-class DPFs, both thermonuclear and beam-target mechanisms can occur over the course of the implosion and contribute to the total neutron production. Hence the neutron source size of a DPF will change throughout the implosion. Experimental neutron radiographs show the increase in source size as the pinch breaks apart, in agreement with simulation's prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Nonlocal, diamagnetic electromagnetic effects in magnetically insulated transmission lines

We identify the time-dependent physics responsible for the critical reduction of current losses in magnetically insulated transmission lines (MITLs) due to uninsulated space charge-limited currents of electrons emitted by field stress. A drive current of sufficiently short pulse length introduces a strong enough time dependence that steady-state results alone become inadequate for the complete understanding of current losses. The time-dependent physics can be described as a nonlocal, diamagnetic electromagnetic response of space charge limited currents. As the pulse length is increased or equivalently, the MITL length reduced, these time-dependent effects diminish and current losses converge to those predicted by the well-known Child–Langmuir law in the external (vacuum) fields. We present a simple one-dimensional (1D) model that encapsulates the essence of this physics. We find excellent agreement with 2D particle-in-cell simulations for two MITL geometries, Cartesian parallel plate and azimuthally symmetric straight coaxial. Based on the 1D model, we explore various scaling dependencies of MITL losses with relevant parameters, e.g., peak current, pulse length, geometrical dimensions, etc. We propose an improved physics model of magnetic insulation in the form of a Hull curve, which could also help improve predictions of current losses by common circuit element codes, such as BERTHA. Finally, we describe how to calculate the temperature rise due to electron impact within the 1D model.

Computer simulation↗
DOE OSTI2025

Increased electron, positron, and x-ray production from high intensity laser interactions using micro-wire targets

We report increases in energetic electrons, positrons, and x-rays emitted from high-intensity laser interactions (10 18−20 W/cm 2 ) with structured silicon micro-wires on the surface of a 1 mm gold converter target using a 10 ps laser pulse. A total of four different wire configurations are tested, where the gaps (7–28 μm) between the wires and the thicknesses (3–6 μm) of the wires are varied, while the height remains constant (⁠ ~25 μm). We observe the largest enhancement in electrons, positrons, and x-rays with the sparsest wire configurations. The electron temperature (T e ≈6 MeV) remains consistent across all shots, regardless of whether wires or planar targets are used. This suggests that the observed enhancement is due to increased laser light absorption by the accelerated electrons over a long scale length. Two-dimensional particle-in-cell simulations confirm that absorption is significantly enhanced with the wire target. Additionally, specific simulations examining laser pointing on different parts of the wire structure reveal that, while the final electron spectrum remains largely insensitive, the angular distribution is highly sensitive to these variations.

Bremsstrahlung↗
DOE OSTI2025

On the role of self-generated magnetic fields in focused ion beams driven by intense picosecond laser pulses

We investigate the physics of ion beam focusing in the interaction of intense picosecond laser pulses with concave foils, using two-dimensional particle-in-cell simulations. Curved targets have been observed to deliver brighter multi-MeV ion beams than flat targets, and they could be used to ignite pre-compressed ICF fusion capsules in a scheme called ion fast ignition. Focusing is achieved by shaping the target and relying on an ion acceleration mechanism in which the ions are accelerated normal to the target surface; subsequently, the plasma expansion grows magnetic fields that advect with the plasma and deflect ions away from the target symmetry axis, causing the focal spot to shift downstream from its geometric location. The drop in both electric and magnetic field strengths in the expanding plasma makes the ion beam envelope asymmetric with respect to the focusing plane.

Focused ion beam↗
DOE OSTI2025

Multimodal azimuthal oscillations in electron beam generated E×B plasma

Electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E×B) fields are promising for low-damage material processing. However, these plasmas can be subject to the formation of azimuthally propagating structures that enhance the radial transport of energetic charged species, which can harm the gentle processing capability of the plasma. In this work, we investigate the azimuthal structure formation in an e-beam generated E×B plasma using experimental diagnostics and 2D3V particle-in-cell simulations. Our findings demonstrate the formation of multiple simultaneously occurring azimuthally propagating modes that exhibit nontrivial radial dependence. It is suggested that the multimodal azimuthal spectrum is caused in part by the radial outflow of ions in the plasma.

electron beams↗
DOE OSTI2025

Modeled sensitivity of multi-MA accelerator performance to electrode contaminant inventory

Significant particle-in-cell code development has enabled simulations of power flow in multi-MA accelerators to include the desorption of surface contaminants, their ionization into surface plasmas, and the impact of these plasmas on efficiency. The simulations base desorption on an Arrhenius equation, whose most significant unknown is the surface contaminant inventory. The sensitivity of power-flow simulations to this inventory is studied here using Sandia National Laboratories' Z accelerator with a 7-nH MagLIF load [Phys. Plasmas 17, 056303 (2010)]. Simulations are conducted in 3D cylindrical coordinates for the current-adder, or “convolute,” region of Z and in 2D for the final feed only. Simulated contaminant inventories are varied from 1 to 32 monolayers (MLs) in 2D, and 2 to 4 ML in 3D. The results reveal sensitivities to the local ratio of E/B⁠. The high B-field, low E-field region near the short-circuit load is insensitive to the contaminant inventory, where assumed values of 4–32 ML change the load current by ≤ 2%, and agree with experiment to within 2% at peak current. A 1-ML value is the outlier, increasing the load current by 5%, but still within measurement uncertainty. In contrast, the relatively higher E-field, lower B-field convolute region has slower contaminant desorption and higher-magnitude E-field penetration of the surface plasmas. The current loss in the convolute region does increase with contaminant inventory. The loss assuming 4 ML is 12% larger than for 2 ML, with 4 ML being the better match to experiment.

Arrhenius equation↗
DOE OSTI2025

Laser wavelength dependence of particle acceleration mechanisms in high intensity laser–solid density plasma interactions

We investigate the generation of relativistic electrons and the subsequent ion acceleration due to target-normal sheath acceleration when ultra-intense (⁠ I > 10 18 W/cm 2 ⁠) short pulse (⁠ τ L < 10ps⁠) lasers are incident onto solid density targets as laser wavelength is varied. Scaling laws for the hot electron temperature, T hot ⁠, and the maximum ion energy, E max ⁠, are recast as a function of laser wavelength. These predictions are compared to results from particle-in-cell computer simulations in a variety of geometries, including cases where realistic plasma density profiles as determined by a radiation hydrodynamics code are used. It is found that the wavelength dependence observed in simulation is less pronounced than what is predicted from the well-established scaling laws. An assessment of how switching to longer laser wavelengths, specifically 2 μm Tm:YLF technology, would impact current high energy density science applications and diagnostics is made.

Electromagnetism↗
DOE OSTI2026

Scalings of laser-driven MeV electron sources with relativistic multi-picosecond laser pulses

We investigate the generation of MeV electrons in relativistic multi-picosecond laser–solid interactions, using two-dimensional particle-in-cell simulations. Scanning a parameter space of intensities from 10 18 to 10 20 W/cm 2 and pulse durations between 100 fs and 10 ps, we discuss scalings of the effective hot electron temperature and MeV electron dose relevant to applications, as well as the direct laser acceleration mechanism that drives electrons to energies much higher than the ponderomotive potential of the laser pulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Saturation of the kinetic ballooning instability due to the electron parallel nonlinearity

The electron parallel nonlinearity (EPN) is implemented in the gyrokinetic particle-in-cell turbulence code GEM [Y. Chen and S. E. Parker, J. Comp. Phys. 220, 839 (2007)]. Application to the Cyclone Base Case reveals a strong effect of EPN on the saturated heat transport above the kinetic ballooning mode (KBM) threshold. Evidence is provided to show that the strong effect is associated with the electron radial motion due to magnetic fluttering, which turns fine structures of the KBM eigenmode in radius into fine structures in velocity and increases the magnitude of the EPN term in the kinetic equation.

Gyrokinetic simulations↗
DOE OSTI2026

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks↗
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 OSTI2024

Study of synchronous RF pulsing in dual frequency capacitively coupled plasma

Low-pressure multi-frequency capacitively coupled plasmas (CCPs) are used for numerous etch and deposition applications in the semiconductor industry. Pulsing of the radio-frequency (RF) sources enables control of neutral and charged species in the plasma on a millisecond timescale. The synchronous (i.e. simultaneous, in-phase) pulsing of both power sources in a dual frequency CCP is examined in this article. Due to the low gas pressure, modeling has been done using the electrostatic particle-in-cell/Monte Carlo collision method. Here, the objective of this work is to investigate the sensitivity of the plasma properties to small changes in timing during synchronous pulsing of the two RF sources. It is demonstrated that small deviations in the on and off times of the two RF sources can lead to major changes in the plasma characteristics. This high sensitivity is of concern for process repeatability but can be utilized to enable better control of the dynamics of plasma-surface interaction. In the simulations, the pulsing parameters (on and off times and ramp rates) are varied and the temporal evolution of plasma characteristics such as electron density (n e ), species current at the electrode, and electron temperature are examined. It is demonstrated that if the low-frequency (LF) source is turned off a few μs before (or after) the high-frequency source, n e during the off-state is significantly higher (or lower) due to the frequency coupling effect. Similarly, turning on the LF source with a small delay results in a sharp increase in the plasma density when the RF sources are turned on.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Imaging nuclei by smashing them at high energies: how are their shapes revealed after destruction?

High-energy nuclear collisions have recently emerged as a promising ‘imaging-by-smashing’ approach to reveal the intrinsic shapes of atomic nuclei. Here, I outline a conceptual framework for this technique, explaining how nuclear shapes are encoded during quark–gluon plasma (QGP) formation and evolution, and how they can be decoded from final-state particle distributions. Furthermore, I highlight the method’s potential to advance our understanding of both nuclear structure and QGP physics.

imaging by smashing↗
DOE OSTI2023

Search for photons above 10 19 eV with the surface detector of the Pierre Auger Observatory

We use the surface detector of the Pierre Auger Observatory to search for air showers initiated by photons with an energy above 10 19 eV. Photons in the zenith angle range from 30° to 60° can be identified in the overwhelming background of showers initiated by charged cosmic rays through the broader time structure of the signals induced in the water-Cherenkov detectors of the array and the steeper lateral distribution of shower particles reaching ground. Applying the search method to data collected between January 2004 and June 2020, upper limits at 95% CL are set to an E -2 diffuse flux of ultra-high energy photons above 10 19 eV, 2 × 10 19 eV and 4 × 10 19 eV amounting to 2.11 × 10 -3 , 3.12 × 10 -4 and 1.72 × 10 -4 km -2 sr -1 yr -1 , respectively. While the sensitivity of the present search around 2 × 10 19 eV approaches expectations of cosmogenic photon fluxes in the case of a pure-proton composition, it is one order of magnitude above those from more realistic mixed-composition models. The inferred limits have also implications for the search of super-heavy dark matter that are discussed and illustrated.

ultra high energy photons and neutrinos↗
DOE OSTI2024

Comparison of chamber beam geometry robustness to mispointing, imbalance and target offset for direct-drive laser fusion facilities

This study focuses on the optimization of beam chamber geometry designs for future direct-drive laser facilities. It provides a review of leading target chamber geometries, with a particular emphasis on random errors. Through comprehensive solid-sphere illuminations and analysis, we identify an optimized beam geometry design, highlighting its robustness and performance under realistic experimental conditions. Three major sources of random errors are evaluated, closely linked to experimental evaluations at OMEGA. The findings underscore the importance of optimizing the irradiation system alongside beam pattern considerations to enhance the efficiency and reliability of inertial confinement fusion experiments. We conclude that for a desired illumination uniformity of 1% in the presence of system errors, the split icosahedron design is the most robust. However, for a 0.3% uniformity goal, the charged-particle, icosahedron, and t-sphere methods exhibit similar performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Flying focus with arbitrary directionality for spatiotemporal control of laser intensity

Flying focus techniques produce laser pulses whose focal points travel at arbitrary, controllable velocities. While this flexibility can enhance a broad range of laser-based applications, existing techniques constrain the motion of the focal point to the propagation direction of the pulse. Here, we introduce a flying focus configuration that decouples the motion of the focus from the propagation direction. Here, a chirped laser pulse focused and diffracted by a diffractive lens and grating creates a focal point that can move both along and transverse to the propagation direction. The focal length of the lens, grating period, and chirp can be tuned to control the direction and velocity of the focus. Simulations demonstrate this control for a holographic configuration suited to high-power pulses, in which two off-axis pump beams with different focal lengths encode the equivalent phase of a chromatic lens and grating in a gas or plasma. For low-power pulses, conventional solid-state or adaptive optics can be used instead. Multi-dimensional control over the focal trajectory enables new configurations for applications, including laser wakefield acceleration of ions, nonlinear Thomson scattering, and surface-plasmon emission of THz radiation.

Classical optics↗