Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electron beam simulations”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Light-Matter Interaction near the Schwinger Limit Using Tightly Focused Doppler-Boosted Lasers

Strong-field quantum electrodynamics (SF QED) is a burgeoning research topic dealing with electromagnetic fields comparable to the Schwinger field (≈1.32×10 18 V/m). While most past and proposed experiments rely on reaching this field in the rest frame of relativistic particles, the Schwinger limit could also be approached in the laboratory frame by focusing to its diffraction limit the light reflected by a plasma mirror irradiated by a multipetawatt laser. Here, we explore the interaction between such intense light and matter with particle-in-cell simulations. We find that the collision with a relativistic electron beam would enable the study of the nonperturbative regime of SF QED, while the interaction with a solid target leads to a profusion of SF QED effects that retroact on the interaction. In both cases, relativistic attosecond pair jets with high densities are formed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1979

Charging rates of metal-dielectric structures

Metal plates partially covered by 0.01-centimeter-thick fluorinated ethylene-propylene (FEP) Teflon were charged in the Lewis Research Center's geomagnetic substorm simulation facility using 5-, 8-, 10-, and 12-kilovolt electron beams. Surface voltage as a function of time was measured for various initial conditions (Teflon discharged or precharged) with the metal plate grounded or floating. Results indicate that both the charging rates and the levels to which the samples become charged are influenced by the geometry and initial charge state of the insulating surfaces. The experiments are described and the results are presented and discussed. NASA charging analyzer program (NASCAP) models of the experiments have been generated, and the predictions obtained are described. Implications of the study results for spacecraft are discussed.

Purvis, C. K.↗
NASA NTRS1987

Effect of wave localization on plasma instabilities

The Anderson model of wave localization in random media is involved to study the effect of solar wind density turbulence on plasma processes associated with the solar type III radio burst. ISEE-3 satellite data indicate that a possible model for the type III process is the parametric decay of Langmuir waves excited by solar flare electron streams into daughter electromagnetic and ion acoustic waves. The threshold for this instability, however, is much higher than observed Langmuir wave levels because of rapid wave convection of the transverse electromagnetic daughter wave in the case where the solar wind is assumed homogeneous. Langmuir and transverse waves near critical density satisfy the Ioffe-Reigel criteria for wave localization in the solar wind with observed density fluctuations -1 percent. Numerical simulations of wave propagation in random media confirm the localization length predictions of Escande and Souillard for stationary density fluctations. For mobile density fluctuations localized wave packets spread at the propagation velocity of the density fluctuations rather than the group velocity of the waves. Computer simulations using a linearized hybrid code show that an electron beam will excite localized Langmuir waves in a plasma with density turbulence. An action principle approach is used to develop a theory of non-linear wave processes when waves are localized. A theory of resonant particles diffusion by localized waves is developed to explain the saturation of the beam-plasma instability. It is argued that localization of electromagnetic waves will allow the instability threshold to be exceeded for the parametric decay discussed above.

Levedahl, William Kirk↗
DOE OSTI2024

CEBAF Injector for K Long Beam Conditions

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab concurrently operates four experimental Halls with distinct bunch charge specifications and repetition rates. Numerous critical beam parameters within CEBAF are configured in the injector, some remaining unchanged throughout the accelerator. Consequently, the injector plays a crucial role in determining final beam characteristics, including bunch structure, beam sizes, bunch lengths, energy spread, and beam transmission. The Jefferson Lab KL experiment is scheduled to take place at CEBAF in Hall D, featuring a much lower bunch repetition rate of 7.80 MHz or 15.59 MHz, below the nominal values of 249.5 MHz or 499 MHz. Although the proposed average current of 5 ?A or 10 ?A is low compared to the maximum CEBAF cur- rent of approximately 180 ?A, the corresponding bunch charge is unusually high for CEBAF injector operation. This study focuses on the behavior of low-repetition-rate, high-bunch- charge (0.32 to 0.64 pC) beams in the CEBAF injector. We investigated the evolution and transmission of low-charge beams to space-charge dominated high-charge beams in the front end of the CEBAF injector for two configurations: the pre-existing CEBAF Phase 1 injector upgrade, operated at 130 kV, and the existing CEBAF Phase 2 injector upgrade, operated at 140 kV, 180 kV, and to be operated at 200 kV. The electron beam through the CEBAF injector is characterized using beam dynamics simulations and comparisons with the available measurements performed at 130 kV. Multi-objective genetic optimizations of the CEBAF injector were performed to determine the operating magnetic elements and RF settings for the evolution and transmission of low, moderate, and high charge beams in the CEBAF injector at 180 kV and 200 kV DC gun voltages. Subsequently, simulations at the same voltages were conducted to obtain the beam characteristics at the front end of the CEBAF injector. The laser spot size and laser pulse length at the cathode were varied to observe their effects on beam transmission in the injector at different voltages (130 kV, 180 kV, and 200 kV). Experimental studies at 130 kV, 140 kV, and 180 kV validate the simulations. Beam study measurements are carried out using EPICS tools, while optimizations and simulations are facilitated by General Particle Tracer. Based on the findings, optimal parameters for the upcoming Jefferson Lab KL experiment are proposed, utilizing a lower repetition rate and higher bunch charge

Pokharel, Sunil↗
DOE OSTI2026

Collider-quality electron bunches from an all-optical plasma photoinjector

We present an approach for generating collider-quality electron bunches using a plasma photoinjector. The approach leverages recently developed techniques for the spatiotemporal control of laser pulses to produce a moving ionization front in a nonlinear plasma wave. The moving ionization front generates an electron bunch with a current profile that balances the longitudinal electric field of an electron beam-driven plasma wave, creating a uniform accelerating field across the bunch. Particle-in-cell (PIC) simulations of the ionization stage show the formation of an electron bunch with 220 pC charge and low emittance (ɛ 𝑥 = 171 nm rad, ɛ 𝑦 = 76 nm rad). Quasistatic PIC simulations of the acceleration stage show that the bunch is efficiently accelerated to 24 GeV over 2 m with a final energy spread of less than 1% and emittances of ɛ 𝑥 = 189 nm rad and ɛ 𝑦 = 80 nm rad. This high-quality electron bunch meets the requirements outlined by the Snowmass process for intermediate-energy colliders and compares favorably to the beam quality of proposed and existing accelerator facilities. The results establish the feasibility of plasma photoinjectors for future collider applications making a significant step toward the realization of high-luminosity, compact accelerators for particle physics research.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

Particle accelerators are among the most versatile and important tools of scientific discovery. The Nation's accelerators are responsible for a wealth of advances in materials science, chemistry, the biosciences, particle physics, and nuclear physics. They also have important applications to national security, the environment, energy, medicine, and on the quality of people's lives. LANL has a long history of making pioneering contributions to Accelerator Science including key contributions to the field of Computational Accelerator Physics. These include the development of early beam dynamics codes with space charge (such as PARMILA and PARMELA), the development of rf cavity codes and magnet codes (including Poisson and Superfish), and the development and distribution of codes to the accelerator community through the Los Alamos Accelerator Code Group. LANL researchers also helped pioneer the development of massively parallel space-charge codes. In project t22_accelsim we have moved beyond electrostatic models of collective effects (i.e., solving the Poisson equation in the bunch frame) to fully electromagnetic models based on the Lienard-Wiechert formalism. This approach enables the large-scale simulation of radiation production and collective effects in high brightness electron beams. This is highly relevant to LANL given its future goal of developing an X-ray Free Electron Laser (XFEL). It also directly impacts a LANL LDRD project to develop an undulator-based non-invasive beam profile monitor for beams created in laser-plasma accelerator systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Directed gas phase preparation of ethynylallene (H 2 CCCHCCH; X 1 A') via the crossed molecular beam reaction of the methylidyne radical (CH; X 2 Π) with vinylacetylene (H 2 CCHCCH; X 1 A')

The gas-phase bimolecular reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A') was conducted at a collision energy of 20.3 kJ mol –1 under single collision conditions exploiting the crossed molecular beam experimental results merged with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations. The laboratory data reveal that the bimolecular reaction proceeds barrierlessly via indirect scattering dynamics through long-lived C 5 H 5 reaction intermediate(s) ultimately dissociating to C 5 H 4 isomers along with atomic hydrogen with the latter predominantly originating from the vinylacetylene reactant as confirmed by the isotopic substitution experiments in the D1-methylidyne–vinylacetylene reaction. Combined with ab initio calculations of the potential energy surface (PES) and statistical Rice–Ramsperger–Kassel–Marcus (RRKM) calculations, the experimental determined reaction energy of –146 ± 26 kJ mol –1 along with the distribution minimum of T(θ) at 90° and isotopic substitution experiments suggest ethynylallene (p1; ΔrG = –230 ± 4 kJ mol –1 ) as the dominant product. The ethynylallene (p1) may be formed with extensive rovibrational excitation, which would result in a lower maximum translational energy. Further, AIMD simulations reveal that the reaction dynamics leads to p1 (ethynylallene, 75%) plus atomic hydrogen with the dominant initial complex being i1 formed by methylidyne radical addition to the double C[double bond, length as m-dash]C bond in vinylacetylene. Altogether, combining the crossed molecular beam experimental results with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations, ethynylallene (p1) is expected to represent the dominant product in the reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A').

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Excitation of surface waves in 3D ion beam neutralization

Neutralization of beams with 2D and 3D geometries by the electrons emitted from an external source is studied using particle-in-cell simulations. Our work reveals that the high-energy electrons excite Trivelpiece–Gould (TG) surface waves in the beams with 3D axisymmetric geometries. These high-energy electrons are generated because of a large amplitude electrostatic solitary wave (ESW) that forms near the electron source and has an electric potential amplitude more than three times the electron thermal energy. Here, we also find that surface wave excitation only happens when the beam radius is large enough at the ion source to attract enough electrons that could form the large amplitude ESW. A comparison of the 3D TG surface wave dispersion relation with an expression for 2D surface waves reveals that they become excited in 3D axisymmetric but not in 2D planar beam because of a higher phase speed requirement in the latter case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Theory and simulation of electromagnetic beam modes and whistlers

Using particle-in-cell simulations and analytical methods, a study of the nonlinear evolution of electromagnetic instabilities driven by an anisotropic electron beam in an external magnetic field is performed. The unstable waves are either whistlerlike or beam-mode-like depending on the external field strength and beam velocity. The evolution of the particle distribution differs significantly in the two regimes. Even in the presence of a faster electrostatic instability, the electromagnetic waves grow to a significant amplitude. In certain cases, an energetic tail is formed, resulting in enhanced Cerenkov emission of electrostatic waves. The initial evolution of the particle distribution is explained in terms of the interaction of a given linearly unstable wave with the self-consistent perturbed distribution.

Newman, David L.↗
NASA NTRS2018

Numerical Studies of Hall Thruster Acceleration Region Electron Transport

Azimuthal direction 1-D and azimuthal and axial direction 2-D Particle in Cell (PIC) simulations of the Hall thruster acceleration and near-plume regions are presented. The axial extent of the 2-D computational domain is based on the ExB drift parameter, 𝑫𝑬𝒙𝑩, the number of times around the channel an electron ExB drifts before being scattered. With only high energy main beam ions, the simulations produced a broad potential profile in the acceleration region. The addition of a small amount of ionization in the acceleration domain at a rate consistent with Hall thruster parameters causes the potential profile to develop two distinct regions. In the upstream portion of the computational domain, the potential profile is very steep, typical of Hall thruster acceleration regions. In the downstream portion, the potential profile is almost flat, like potentials in Hall thruster near plume regions. This bifurcation happens naturally; there is no other change in code parameters between the two results. Without the addition of a small amount of ionization in the acceleration region, the axial electron transport is very low. With ionization, the calculated ratio of electron current to ion current is the same order of those measured in Hall thrusters. This suggests that in regions with lower energy electrons and ions, plasma waves are effective in allowing crossfield transport, and axial electric fields are small. In the acceleration region, the plasma waves are less effective at scattering the energetic electrons, and high axial electric fields are needed to produce the observed cross-field electron transport.

Ortega, Alejandro Lopez↗
NASA NTRS1989

Coherent Cerenkov radiation from the Spacelab 2 electron beam

The plasma environment of the Spacelab 2 mission was investigated through the deployment of the Plasma Diagnostics Package (PDP) by the Space Shuttle Orbiter and the Orbiter's ejection of a continuous 1-keV/50-mA electron beam along a field line. As the PDP flew by the beam, its plasma-wave instrument detected intense whistler-mode radiation originating from the beam. A detailed model has been developed of the coherent Cerenkov emission process, using a one-dimensional computer simulation of the beam to model the expected phase space structure of the electrons. The power calculated for the modeled 200-m beam segment can easily account for the measured whistler mode wave power.

Farrell, W. M.↗
NASA NTRS2024

Additively Manufactured Plastic Plasma Spectrometer

We report results in development and testing of a low resource top hat electrostatic analyzer (ESA) for space plasma measurements. This device has been additively manufactured (3D-printed) using fused deposition modeling. The classic top hat design is composed of four plastic pieces, without any surface coatings. The three conducting electrodes are printed from carbon nanotube infused polyether ether ketone (CNT-PEEK). The fourth piece, an insulating electrode support, uses pure PEEK. This ESA is designed to detect electrons in space from 10eV up to 30 keV. We demonstrate that the printed CNT-PEEK is sufficiently electrically conductive to support the fast high voltage slewing often required for high time resolution measurements. The plastic ESA has been successfully vibrated beyond standard pre-flight levels, tested under keV electron beam illumination over a wide range of temperatures, and tested under UV illumination, simulating the solar Ly-α flux. In comparison with an identical machined aluminum ESA, our CNT-PEEK ESA provides nominal energy/angle band passes, closely matching simulation. These band passes imply minimal impact from surface charging at beam energies of 2-3 keV, although more investigation is needed. We also find that the CNT-PEEK ESA provides far superior out-of-band electron rejection and UV photon rejection compared to the machined aluminum ESA. We do not detect any problems with trapped gases or outgassing. This development offers potential for significant mass savings, implementation of otherwise unattainable geometric configurations, and dramatic simplification in manufacturing and assembly processes required for development of space plasma instruments.

Quetzal Larrick↗
DOE OSTI2021

Particle-In-Cell Computer Simulations of Beam Spill in the Scorpius Injector

Computer simulations of beam transport in the Scorpius injector were performed using the particle-in-cell computer code LSP. The simulations were used to predict spot size growth due to electron stimulated emission of ions from beam spill on the anode stalk and the entrance transport region.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Physics↗
DOE OSTI2022

Ultrabright Electron Bunch Injection in a Plasma Wakefield Driven by a Superluminal Flying Focus Electron Beam

Here, we propose a new method for self-injection of high-quality electron bunches in the plasma wakefield structure in the blowout regime utilizing a “flying focus” produced by a drive beam with an energy chirp. In a flying focus the speed of the density centroid of the drive bunch can be superluminal or subluminal by utilizing the chromatic dependence of the focusing optics. We first derive the focal velocity and the characteristic length of the focal spot in terms of the focal length and an energy chirp. We then demonstrate using multidimensional particle-in-cell simulations that a wake driven by a superluminally propagating flying focus of an electron beam can generate GeV-level electron bunches with ultralow normalized slice emittance (~30 nm rad), high current (~17 kA), low slice energy spread (~0.1%), and therefore high normalized brightness (> 10 19 A/ m 2 /rad 2 ) in a plasma of density ~10 19 cm –3 . The injection process is highly controllable and tunable by changing the focal velocity and shaping the drive beam current. Nearterm experiments at FACET II where the capabilities to generate tens of kA, < 10 fs drivers are planned, could potentially produce beams with brightness near 10 20 A / m 2 / rad 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Physics of heavy ions (1989-1990)

The results from studies on polar wind ion heating due to kinetic ion beam instabilities and the effects of such ion heating on the outflow of O(+) in the polar wind are presented and discussed. First, the linear instabilities associated with an O(+) and H(+) polar wind plasma in the presence of O(+) and H(+) beams for a range of O(+)/H(+) beam densities, T(sub e)/T(sub i), and ion beam speeds were examined. Then, nonlinear heating of the polar wind ions was studied, using numerical simulations. The O(+) and H(+) polar wind ions were modeled by isotropic Maxwellian distributions, and the electrons, O(+) beams, and H(+) beams were modeled by drifting Maxwellian distributions. The effects of the kinetic ion heating on the outflow of the polar wind ions were examined from the ionosphere, using a time-dependent hydrodynamic model. A numerical code to solve the O(+) and H(+) continuity and momentum equations in a flux tube from ionospheric to magnetospheric altitudes were developed. The effects of ion heating were included by allowing for the altitudinal variation of the ion temperatures in the momentum equation. The ion temperature profiles were specified based on the ion heating characteristics found from previous kinetic simulations. It was assumed that heating occurred above 1500 km and increased to a saturated value of temperature that was obtained directly from the kinetic simulation study. The characteristics of the dynamical polar wind without ion heating were studied, and a flux tube on closed field lines that suddenly became open at t = 0 was simulated. Then, the effects of ion heating were included. To gain some physical insight, two limiting cases were considered: preferential H(+) heating and preferential O(+) heating. How O(+) heating can lead to enhanced polar wind O(+) fluxes in the polar magnetosphere is shown.

Ashour-Abdalla, M.↗
DOE OSTI2023

Refined interpretation of electron temperature response to neutral beam injection at DIII-D

Accurate particle and power deposition profiles of neutral beam injection (NBI) are essential to transport studies, and that information is usually acquired through Monte Carlo simulations with a given collisional model. The deposition process of the energetic beam particles leads to the informative electron temperature (T e ) evolution trajectory, which can be captured by electron cyclotron emission (ECE) system due to its good spatial and temporal resolution. Previously, some work has been done to interpret the T e responses to the pulsed NBI as a linear heating source with Fourier-based techniques, although that approach fell short when the fast ion slowing-down time becomes significant (~100 ms). It has been observed in DIII-D that the modulated NBI pulses (10–50 Hz) reduce local core T e values ~0.1 keV through cold electron dilution in high-T e (>2 keV) plasmas alongside accumulative heating. Here, a novel approach to interpret the T e response to NBI was developed by linearizing and modeling the detailed T e evolution trajectory using coherently averaged ECE data based on the different time scales of the terms in the local power and particle balance equations. The technique does not require absolute calibrations of ECE and is independent of collisional models. The resulting beam deposition profiles show good consistency and reasonable agreement with Monte Carlo calculations based on the atomic data from the Atomic Data and Analysis Structure (ADAS). Local electron density response measured by Thomson scattering (TS) also suggests the same features when the beam pulse is large enough for that diagnostic to resolve. The remaining discrepancies are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

PUFFIn – A user friendly fast interface for calculating and visualizing the dose distribution in materials

This article describes a new graphical user interface that uses the PENELOPE Monte Carlo code to calculate dose distributions in materials. It is named PUFFIn, for the Penelope User-Friendly Fast Interface and was developed as an educational and scoping analysis tool for non-experts in radiation modeling. PUFFIn enables the user to visualize and compare the dose distributions in objects (e.g., sterilized healthcare products) irradiated with cobalt-60 gamma-rays, electron beam (E-beam) or X-rays. From such comparisons, the user can determine the most efficient product and/or packaging designs for any given radiation field – whether for conceptual or existing products and packaging. PUFFIn is distributed in a complete, self-contained package of software, including the PENELOPE radiation transport code, a standard graphics package, and a set of simple exercises. The package is available to any user at no cost and requires a minimal amount of training compared to other similar software. PUFFIn's capabilities are described, as well as validation measurements performed at Texas A&M University E-beam facility and the Aerial E-beam Facility in France.

97 MATHEMATICS AND COMPUTING↗