Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High efficiency RF sources developments

Calabazas Creek Research, Inc. (CCR) and its collaborators are developing high efficiency RF sources operating from a few hundred MHz to C-Band and power levels from tens to hundreds of kilowatts with the goal of providing MW-relevant sources. The efficiencies approach or exceed 80% with projected costs as low as $0.50/ Watt. Sources under development include magnetrons with phase and amplitude control, single and multi-beam klystrons, multi-beam power grid tubes, and multiple beam IOTs. A magnetron system achieved more than 80% efficiency with fast amplitude control using modulation of the phase locking signal. This would be a low cost, high efficiency RF source for superconducting accelerators. An L-Band, single beam klystron was built with simulated efficiency of 80%. The klystron has yet to be tested to confirm the simulation results. CCR is currently developing a multi-beam klystron to produce more than 200 kW CW at 80% efficiency. Also in development is a multiple beam triode to produce 200 kW CW from 300 MHz to approximately 1 GHz. Not only does the simulated efficiency exceed 75%, but it would be the lowest cost RF source in this frequency range. Finally, CCR recently concluded research for a multiple beam IOT at 700 MHz using third harmonic drive to boost efficiency toward 85%. Successful development and transition to production of these sources will significantly alter the cost/performance landscape for RF power generation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Implementation of compound refractive lenses for large field-of-view x-ray phase-contrast imaging during hypervelocity impact experiments

Synchrotron x-ray phase-contrast imaging (XPCI) offers time-resolved visualization of dynamic compression phenomena, but its intrinsically small field-of-view (FOV) limits the time that key features remain in frame. A novel approach to enlarge the FOV is achieved by positioning a two-dimensional parabolic compound refractive lens (CRL) upstream of the sample to deliberately defocus the white beam. Ray-tracing simulations and XPCI measurements show that this CRL configuration can expand the beam by ∼50% vertically and ∼15% horizontally based on the full width at half-maximum of the beam. Implementing the CRL, however, attenuates the photon flux and lowers signal-to-noise ratio (SNR). Task-based analysis using a calibration grid (30 μm dots) showed that both setups fail to consistently meet the Rose criterion (SNR ≥ 5) for features of this size in single-bunch imaging. Extrapolating the measured SNR Rose values suggests that the minimum consistently detectable feature lies closer to 30–40 μm for the standard XPCI setup and above 40 μm for CRL-XPCI. Despite this limitation, the CRL configuration nearly doubles the illuminated area, enabling simultaneous tracking of front and rear observations of boron carbide targets subjected to rod and sphere impacts at 1.0–2.6 km/s. Image tracking algorithms and photonic Doppler velocimetry were used to measure penetration and rear-surface velocity histories. Together, these measurements capture crack fronts, penetration, and material breakout, offering new benchmark data for validating high-strain-rate constitutive models of ceramic materials.

Ceramic materials↗
DOE OSTI2023

A simple, flexible technique for RF cavity wake-field calculations

It is typical in the accelerator field to model machine components, especially RF cavities, as parallel RLC resonators. In the interest of simulating and diagnosing beam instabilities, knowledge of the time-domain voltage waveform over an equivalent resonator by a bunch current often proves useful. This waveform may be found by convolving the bunch current with the RLC resonator impulse response. While analytical and quasi-analytical expressions are available in this regime for common distributions such as the Gaussian, analogous results for less standard distributions are difficult to obtain using direct methods, which opens the door for the development of a more generalized technique. In this paper, a formulation is created that allows for the simple computation of the time-domain voltage waveform of and RLC resonator. The formulation uses the Cauchy Residue Theorem to extract the convolution result from the Fourier Domain, and it only requires that the current distribution Fourier Transform be holomorphic and known at one specific evaluation point. This greatly simplifies the computation of the time domain voltage for a large amount of bunch distributions both common and uncommon.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Plasma-based particle sources

High-brightness particle beams generated by advanced accelerator concepts have the potential to become an essential part of future accelerator technology. In particular, high-gradient accelerators can generate and rapidly accelerate particle beams to relativistic energies. The rapid acceleration and strong confining fields can minimize irreversible detrimental effects to the beam brightness that occur at low beam energies, such as emittance growth or pulse elongation caused by space charge forces. Due to the high accelerating gradients, these novel accelerators are also significantly more compact than conventional technology. Advanced accelerators can be extremely variable and are capable of generating particle beams with vastly different properties using the same driver and setup with only modest changes to the interaction parameters. So far, efforts have mainly been focused on the generation of electron beams, but there are concepts to extend the sources to generate spin-polarized electron beams or positron beams. The beam parameters of these particle sources are largely determined by the injection and subsequent acceleration processes. Although, over the last decade there has been significant progress, the sources are still lacking a sufficiently high 6-dimensional (D) phase-space density that includes small transverse emittance, small energy spread and high charge, and operation at high repetition rate. This is required for future particle colliders with a sufficiently high luminosity or for more near-term applications, such as enabling the operation of free-electron lasers (FELs) in the X-ray regime. Major research and development efforts are required to address these limitations in order to realize these approaches for a front-end injector for a future collider or next-generation light sources. In particular, this includes methods to control and manipulate the phase-space and spin degrees-of-freedom of ultrashort plasma-based electron bunches with high accuracy, and methods that increase efficiency and repetition rate. These efforts also include the development of high-resolution diagnostics, such as full 6D phase-space measurements, beam polarimetry and high-fidelity simulation tools. A further increase in beam luminosity can be achieve through emittance damping. Emittance cooling via the emission of synchrotron radiation using current technology requires kilometer-scale damping rings. For future colliders, the damping rings might be replaced by a substantially more compact plasma-based approach. Here, plasma wigglers with significantly stronger magnetic fields are used instead of permanent-magnet based wigglers to achieve similar damping performance but over a two orders of magnitude reduced length.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Advancing Cross-Sectional Scanning Electron Microscopy of Perovskite Solar Cells

Organic–inorganic perovskites are an emerging class of photovoltaic materials. Despite achieving power conversion efficiencies surpassing 26%, the challenge of perovskite stability including degradation during exposure to operational conditions such as light, heat, humidity, water, oxygen, and electric fields is well known. Related, perovskite instability has limited high-resolution electron imaging and characterization techniques that can be used for understanding degradation mechanisms. Furthermore, we demonstrate perovskite device cross-section preparation using mechanical polishing in a water-free environment with cryogenic Ar ion milling. Scanning electron microscopy was then used in both backscattered electron and secondary electron imaging modes to obtain information about layer structure, grain aggregate structure, and compositional heterogeneity. Monte Carlo CASINO simulations inform optimum beam conditions and image acquisition parameters and the effects of accelerating voltage, dwell times, and frame averaging for practical image acquisition are reported.

14 SOLAR ENERGY↗
DOE OSTI2021

Stopping-power enhancement from discrete particle-wake correlations in high-energy-density plasmas

Three-dimensional (3D) simulations of electron beams propagating in high-energy-density plasmas using the quasistatic Particle-in-Cell (PIC) code QuickPIC demonstrate a significant increase in stopping power when beam electrons mutually interact via their wakes. Each beam electron excites a plasma wave wake of wavelength ~ 2 π c / ω p e , where c is the speed of light and ω p e is the background plasma frequency. Here we show that a discrete collection of electrons undergoes a beam-plasma-like instability caused by mutual particle-wake interactions that causes electrons to bunch in the beam, even for beam densities n b for which fluid theory breaks down. This bunching enhances the beam's stopping power, which we call “correlated stopping,” and the effect increases with the “correlation number” N b ≡ n b ( c / ω p e ) 3 . For example, a beam of monoenergetic 9.7 MeV electrons with N b = 1 / 8 , in a cold background plasma with n e = 10 26 cm – 3 (450 g cm – 3 DT), has a stopping power of 2.28 ± 0.04 times the single-electron value, which increases to 1220 ± 5 for N b = 64 . The beam also experiences transverse filamentation, which eventually limits the stopping enhancement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Distributed Stochastic Optimization of a Neural Representation Network for Time-Space Tomography Reconstruction

4D time-space reconstruction of dynamic events or deforming objects using X-ray computed tomography (CT) is an important inverse problem in non-destructive evaluation. Conventional back-projection based reconstruction methods assume that the object remains static for the duration of several tens or hundreds of X-ray projection measurement images (reconstruction of consecutive limited-angle CT scans). However, this is an unrealistic assumption for many in-situ experiments that causes spurious artifacts and inaccurate morphological reconstructions of the object. To solve this problem, we propose to perform a 4D time-space reconstruction using a distributed implicit neural representation (DINR) network that is trained using a novel distributed stochastic training algorithm. Our DINR network learns to reconstruct the object at its output by iterative optimization of its network parameters such that the measured projection images best match the output of the CT forward measurement model. Here, we use a forward measurement model that is a function of the DINR outputs at a sparsely sampled set of continuous valued 4D object coordinates. Unlike previous neural representation architectures that forward and back propagate through dense voxel grids that sample the object's entire time-space coordinates, we only propagate through the DINR at a small subset of object coordinates in each iteration resulting in an order-of-magnitude reduction in memory and compute for training. DINR leverages distributed computation across several compute nodes and GPUs to produce high-fidelity 4D time-space reconstructions. We use both simulated parallel-beam and experimental cone-beam X-ray CT datasets to demonstrate the superior performance of our approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

X-ray pulse generation with ultra-fast flipping of its orbital angular momentum

A method to temporally tailor the properties of X-ray radiation carrying Orbital Angular Momentum (OAM) is presented. In simulations, an electron beam is prepared with a temporally modulated micro-bunching structure which, when radiating at the second harmonic in a helical undulator, generates OAM light with a corresponding temporally modulated intensity. This method is shown to generate attosecond pulse trains of OAM light without the need for any additional external optics, making the wavelength range tunable. In addition to the OAM pulse train, the method can be adapted to generate radiation where the handedness of the OAM mode may also be temporally modulated (flipped).

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Initial Investigation of Coupled Electron-Photon Transport with SCEPTRE and RAPTURE (Rev.1.1)

The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Current-Regulated Arc Modulator Optimization

n this study, the performance of a current-regulated arc modulator was investigated with a focus on its role in initiating and sustaining plasma discharge within the Magnetron Body of the LINAC system. The analysis centered on how switching components, circuit topology, and feedback loop architecture influence critical factors such as energy efficiency, discharge stability, and long-term plasma containment. Particular attention was given to variations in pulse termination behavior, as observed through oscilloscope traces, which revealed inconsistencies affecting the duty factor and cathode temperature. These fluctuations have downstream effects on the cesium-coated cathode surface, thereby impacting H⁻ ion production and beam reliability. Simulation-based testing in LTspice was used to evaluate noise suppression techniques and arc current regulation schemes, revealing how optimized snubber networks, improved pulse shaping, and feedback stability can mitigate modulator-induced noise. The results ide

Campos, Nathan↗
DOE OSTI2020

The Analog Front-end for the LGAD Based Precision Timing Application in CMS ETL

The analog front-end for the Low Gain Avalanche Detector (LGAD) based precision timing application in the CMS Endcap Timing Layer (ETL) has been prototyped in a 65 nm CMOS mini-ASIC named ETROC0. Serving as the very first prototype of ETL readout chip (ETROC), ETROC0 aims to study and demonstrate the performance of the analog frontend, with the goal to achieve 40 to 50 ps time resolution per hit with LGAD (therefore reach about 30ps per track with two detector-layer hits per track). ETROC0 consists of preamplifier and discriminator stages, which amplifies the LGAD signal and generates digital pulses containing time of arrival and time over threshold information. This paper will focus on the design considerations that lead to the ETROC front-end architecture choice, the key design features of the building blocks, the methodology of using the LGAD simulation data to evaluate and optimize the front-end design. The ETROC0 prototype chips have been extensively tested using charge injection and the measured performance agrees well with simulation. The initial beam test results are also presented, with time resolution of around 33 ps observed from the preamplifier waveform analysis and around 41 ps from the discriminator pulses analysis. A subset of ETROC0 chips have also been tested to a total ionizing dose of 100 MRad with X-ray and no performance degradation been observed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Development and Testing of a new Current-Regulated Arc Modulator for the LINAC

In this study, the performance of a current-regulated arc modulator was investigated with a focus on its role in initiating and sustaining plasma discharge within the Magnetron Body of the LINAC system. The analysis centered on how switching components, circuit topology, and feedback loop architecture influence critical factors such as energy efficiency, discharge stability, and long-term plasma containment. Particular attention was given to variations in pulse termination behavior, as observed through oscilloscope traces, which revealed inconsistencies affecting the duty factor and cathode temperature. These fluctuations have downstream effects on the cesium-coated cathode surface, thereby impacting H⁻ ion production and beam reliability. Simulation-based testing in LTspice was used to evaluate noise suppression techniques and arc current regulation schemes, revealing how optimized snubber networks, improved pulse shaping, and feedback stability can mitigate modulator-induced noise. The results identified hardware level parameters that significantly enhance discharge repeatability and improve overall plasma performance under operational conditions.

Campos, Nathan [Unlisted, US; Fermilab]↗
DOE OSTI2023

Hybrid Simulation of Proton Cyclotron Waves Upstream of Mars Generated by Pickup Ion Beam Distribution

Linear instability analysis as well as a corresponding two‐dimensional hybrid simulation are performed to examine the excitation of the proton cyclotron waves observed upstream of Mars. The waves are believed to be excited by the pickup ions produced from the ionization of the Martian hydrogen exosphere. And previous statistical analysis of wave observations suggested that the waves are mostly related to pickup ion beam velocity distributions. While earlier linear instability analysis of pickup ion beam distributions has mainly been focused on the parallel unstable modes, our analysis reveals that the maximum growth rate occurs at very oblique propagation. The corresponding hybrid simulation confirms the linear analysis results and further demonstrates that the pickup ions are scattered toward an isotropic shell velocity distribution by the waves excited. Interestingly, the waves at oblique propagation gradually damp out and the system is eventually dominated by waves of quasi‐parallel propagation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Monte Carlo simulations demonstrating physics of equivalency of gamma, electron-beam, and X-ray for radiation sterilization

We report the sterilization of medical devices using the gamma rays from the decay of cobalt-60 has accumulated decades of experience of the performance of the materials and devices that are irradiated. The use of radiation using electron beams and x-rays has much less experience and this leads to questions of equivalency between these three technologies. Computer simulations were conducted to model the relevant physical processes of the interactions of each of the three forms of radiation in order to compare the spectra of electron energies at energies below 500 keV. It is predominantly the electrons below this threshold that produce the sterilization dose. No difference in energy spectra was seen between the three types of initial radiation. It is concluded that there is no energy dependent difference between gamma, e-beam, and x-ray for radiation sterilization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

3D Simulations Capture the Persistent Low-Mode Asymmetries Evident in Laser-Direct-Drive Implosions on OMEGA

Spherical implosions in Inertial Confinement Fusion (ICF) are inherently sensitive to perturbations that may arise from experimental constraints and errors. Control and mitigation of low-mode (long wavelengths) perturbations is a key milestone to improving implosion performances. Here, we present the first 3-D radiation-hydrodynamic simulations of directly driven ICF implosions with an inline package for polarized Crossed-Beam Energy Transfer (CBET). Simulations match bang times, yields (separately accounting for laser-induced high modes and fuel age), hot spot flow velocities and direction, for which polarized CBET contributes to the systematic flow orientation evident in the OMEGA implosion database. Current levels of beam mispointing, imbalance, target offset and asymmetry from polarized CBET degrade yields by more than 40%. The effectiveness of two mitigation strategies for low-modes is explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electron Energy Regression in High-Granularity Calorimeter Prototype

The dataset consists of simulations of calibrated reconstructed hits produced by a positron passing through the HGCAL test beam prototype. For the simulations, Monte Carlo method is used to produce the positrons with energy ranging from 10 to 350 GeV. The dataset contains the coordinates of the calibrated reconstructed hits in the prototype along with the calibrated energy in units of MIP. The HDF5 files can be extracted from the gzip files.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗