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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Positron Injector for Ce+BAF

A baseline concept for a continuous wave (CW) polarized positron injector was developed for the CEBAF at Jefferson Lab. This concept is based on the generation of CW longitudinally polarized positrons by a high-current, polarized electron beam 1 mA, 130-370 MeV, and 90% longitudinal polarization that passes through a rotating, water-cooled, tungsten target. The simulations for the Ce+BAF injector at the Low Energy Recirculator Facility (LERF) are performed, including positron beam generation, capture, energy selection, and acceleration to 123 MeV. The positron yield (or positron current) and longitudinal polarization are calculated considering the longitudinal and transverse CEBAF acceptances (<1% energy spread, 1.2 mm bunch length and normalized emittance of 100 mm·mrad). The impact of target thickness, electron beam energy, and transverse size on positron yield within the required emittance limit is evaluated. Target thermal and structural FEM analyses (ANSYS Fluent) are performed to determine the maximum electron beam current and minimum transverse beam size on the target.

Benesch, Jay [Thomas Jefferson National Accelerato↗
DOE OSTI2026

The Memory Scaling of Reverse-Mode Differentiation in Particle Accelerator Simulations with Space Charge

The recent development of differentiable simulation codes for particle accelerators has enabled gradient-based workflows that promise finer control and more realistic modeling of accelerator facilities. However, when using reverse-mode automatic differentiation, the memory usage continuously increases during the simulation, and can potentially exceed the available hardware memory - especially when costly space charge computation is included. To study the memory requirements for differentiable simulations, we have implemented space charge in Cheetah, a PyTorch-based beam tracking code that supports reverse-mode differentiation. We find that the memory usage for reverse-mode differentiation grows linearly with the number of macroparticles and cells, and that it is proportional to the number of space charge kicks involved in the simulation. This general scaling can be used to evaluate whether a given differentiable simulation is feasible given hardware memory constraints.

Dhamrait, Arjun↗
DOE OSTI2025

Monte Carlo simulation analysis for radiation damage on Glidcop Al-15 caused by 17-20 MeV/u heavy ions

The Facility for Rare Isotope Beams (FRIB) SRF linear heavy-ion accelerator is designed to accelerate all ions up to uranium to a maximum beam power of 400 kW. Several beam intercepting devices (BIDs) are essential to the successful operation of the accelerator, including a low power charge selector (LPCS) made of copper containing 0.15% precipitated aluminum oxide by weight called Glidcop Al-15. As FRIB ramps up the primary beam power beyond the current 20 kW level, the charge selector must withstand higher radiation damage rates, typically measured in displacements per atom (dpa). Significant beam induced radiation damage including significant deformation, like swelling, blistering, and cracking has been observed on a recently removed LPCS. We also observed physical features up to about 5.5 mm wide and appear to be deeper than the projected range of any ions. This paper presents the results of Monte Carlo simulations carried out using the Particle and Heavy Ion Transport code System (PHITS), quantifying the total damage dose and ion concentration, etc. Simulating an accurate irradiation history is essential to determining the scope of post irradiation examination (PIE) work.

Fletcher, Ethan [Michigan State U.]↗
DOE OSTI2021

Advanced assessment of beam-induced background at a muon collider

Renewed international interest in muon colliders motivates the continued investigation of the impacts of beam-induced background on detector performance. This continues the effort initiated by the Muon Accelerator Program and carried out until 2017. The beam-induced background from muon decays directly impacts detector performance and must be mitigated by optimizing the overall machine design, with particular attention paid to the machine detector interface region. In order to produce beam-induced background events and to study their characteristics in coordination with the collider optimization, a flexible simulation approach is needed. To achieve this goal we have chosen to utilize the combination of LineBuilder and Monte Carlo FLUKA codes. We report the results of beam-induced background studies with these tools obtained for a 1.5 TeV center of mass energy collider configuration. Good agreement with previous simulations using the MARS15 code demonstrates that our choice of tools meet the accuracy and performance requirements to perform future optimization studies on muon collider designs.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Optimization of a High Bunch Charge ERL Injection Merger for PERLE

Delivery of high charge electron bunches into the main loop of an ERL (energy recovery linac) while preserving the emittance is challenging. This is because at the typical injection momentum, space charge forces still have a significant effect on the beam dynamics. In this work we consider the design of the merger for PERLE, an ERL test facility to be based at IJCLab in France. Previous simulations have shown that the baseline DC gun based injector can achieve the required emittance at the booster linac exit. The quality of the 500 pC bunches must then be preserved with space charge through the merger at total beam energy of 7 MeV keeping the emittance below 6 mm mrad. The beam dynamics in the merger were simulated using the code OPAL and optimised using a genetic algorithm. Three possible merger schemes were investigated. The goal of the optimisation was to minimise the emittance growth while also achieving the required Twiss parameters to match onto the spreader at the main linac exit. A three dipole solution is then examined in more detail.

Hounsell, B.↗
DOE OSTI2023

An empirical method for geometric calibration of a photon counting detector-based cone beam CT system

BACKGROUND: Geometric calibration is essential in developing a reliable computed tomography (CT) system. It involves estimating the geometry under which the angular projections are acquired. Geometric calibration of cone beam CTs employing small area detectors, such as currently available photon counting detectors (PCDs), is challenging when using traditional-based methods due to detectors’ limited areas. OBJECTIVE: This study presented an empirical method for the geometric calibration of small area PCD-based cone beam CT systems. METHODS: Unlike the traditional methods, we developed an iterative optimization procedure to determine geometric parameters using the reconstructed images of small metal ball bearings (BBs) embedded in a custom-built phantom. An objective function incorporating the sphericities and symmetries of the embedded BBs was defined to assess performance of the reconstruction algorithm with the given initial estimated set of geometric parameters. The optimal parameter values were those which minimized the objective function. The TIGRE toolbox was employed for fast tomographic reconstruction. To evaluate the proposed method, computer simulations were carried out using various numbers of spheres placed in various locations. Furthermore, efficacy of the method was experimentally assessed using a custom-made benchtop PCD-based cone beam CT. RESULTS: Computer simulations validated the accuracy and reproducibility of the proposed method. The precise estimation of the geometric parameters of the benchtop revealed high-quality imaging in CT reconstruction of a breast phantom. Within the phantom, the cylindrical holes, fibers, and speck groups were imaged in high fidelity. The CNR analysis further revealed the quantitative improvements of the reconstruction performed with the estimated parameters using the proposed method. CONCLUSION: Apart from the computational cost, we concluded that the method was easy to implement and robust.

Instruments & Instrumentation↗
DOE OSTI2021

Beam physics research with the IOTA electron lens

The electron lens in the Fermilab Integrable Optics Test Accelerator (IOTA) will enable new research in nonlinear integrable optics, space-charge compensation, electron cooling, and the stability of intense beams. This research addresses scientific questions on high-brightness beams and operational challenges of high-power accelerators for nuclear and particle physics. We review the roles that electron lenses play in this field and the physical principles behind their applications. The design criteria and specifications for the IOTA storage ring and electron lens are then discussed. We conclude with a description of the components of the apparatus.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modulation of intensity emerging from zero effort (MIEZE) with extended Fourier time at large scattering angle

Modulation of Intensity Emerging from Zero Effort (MIEZE) is a neutron resonant spin echo technique that allows one to measure time correlation scattering functions in materials by implementing radio-frequency (RF) intensity modulation at the sample and the detector. The technique avoids neutron spin manipulation between the sample and the detector and, thus, could find applications in cases where the sample depolarizes the neutron beam. However, the finite sample size creates a variance in the path length between the locations where scattering and detection happen, which limits the contrast in intensity modulation that one can detect, in particular, toward long correlation times or large scattering angles. We propose a modification to the MIEZE setup that will enable one to extend those detection limits to longer times and larger angles. Furthermore, we use Monte Carlo simulations of a neutron scattering beamline to show that by tilting the RF flippers in the primary spectrometer with respect to the beam direction, one can shape the wave front of the intensity modulation at the sample to compensate for the path variance from the sample and the detector. The simulation results indicate that this change enables one to operate a MIEZE instrument at much increased RF frequencies, thus improving the effective energy resolution of the technique. For the MIEZE instrument simulated, it shows that for an incident beam with the maximum divergence of 0.33°, the maximum Fourier time can be increased by a factor of 3.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Effects of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer

The influences of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer (CBET) are examined using two- and three-dimensional particle-in-cell simulations in parameter regimes relevant to recent CBET experiments at the OMEGA laser facility. In mid-Z plasma irradiated by an intense pump beam and weaker probe beam, ion trapping, collisional de-trapping, and plasma flow induced by thermal effects are shown to affect the CBET gain. Ion trapping can enhance or detune the CBET resonance [Nguyen et al., Phys. Plasmas 28, 082705 (2021)]. Collisional de-trapping can affect the CBET gain at low seed beam intensity near the onset threshold for ion trapping. Thermal-effects-induced flow can also detune the CBET resonance at a level comparable to that from trapping at low seed beam intensity. As a consequence, the CBET gain is sensitive to collisions and dimensionality at low seed beam intensity where ion trapping is weak but is insensitive to collisions and dimensionality at high seed beam intensity where ion trapping is strong.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

NASA’s first ground-based Galactic Cosmic Ray Simulator: Enabling a new era in space radiobiology research

With exciting new NASA plans for a sustainable return to the moon, astronauts will once again leave Earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic radiation (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPEs with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include carcinogenesis, central nervous system (CNS) effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including circulatory and heart disease, as well as potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy-ion exposures and the validation of countermeasures in a relevant space environment. Historically, most research on understanding space radiation-induced health risks has been performed using acute exposures of monoenergetic single-ion beams. However, the space radiation environment consists of a wide variety of ion species over a broad energy range. Using the fast beam switching and controls systems technology recently developed at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory, a new era in radiobiological research is possible. NASA has developed the “GCR Simulator” to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (approximately 65%–75%) and helium ions (approximately 10%–20%) with heavier ions (Z ≥ 3) contributing the remainder. The GCR simulator exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the 5 heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader system is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low-energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires approximately 75 minutes. To more closely simulate the low-dose rates found in space, sequential field exposures can be divided into daily fractions over 2 to 6 weeks, with individual beam fractions as low as 0.1 to 0.2 mGy. In the large beam configuration (60 × 60 cm 2 ), 54 special housing cages can accommodate 2 to 3 mice each for an approximately 75 min duration or 15 individually housed rats. On June 15, 2018, the NSRL made a significant achievement by completing the first operational run using the new GCR simulator. This paper discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NSRL to accelerate our understanding and mitigation of health risks faced by astronauts. Ultimately, the GCR simulator will require validation across multiple radiogenic risks, endpoints, doses, and dose rates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Beam Hardening Artifact Reduction in X-Ray CT Reconstruction of 3D Printed Metal Parts Leveraging Deep Learning and CAD Models

Nondestructive evaluation (NDE) of additively manufactured (AM) parts is important for understanding the impacts of various process parameters and qualifying the built part. X-ray computed tomography (XCT) has played a critical role in rapid NDE and characterization of AM parts. However, XCT of metal AM parts can be challenging because of artifacts produced by standard reconstruction algorithms as a result of a confounding effect called “beam hardening.” Beam hardening artifacts complicate the analysis of XCT images and adversely impact the process of detecting defects, such as pores and cracks, which is key to ensuring the quality of the parts being printed. In this work, we propose a novel framework based on using available computer-aided design (CAD) models for parts to be manufactured, accurate XCT simulations, and a deep-neural network to produce high-quality XCT reconstructions from data that are affected by noise and beam hardening. Using extensive experiments with simulated data sets, we demonstrate that our method can significantly improve the reconstruction quality, thereby enabling better detection of defects compared with the state of the art. We also present promising preliminary results of applying the deep networks trained using CAD models to experimental data obtained from XCT of an AM jet-engine turbine blade.

Ziabari, Amir↗
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↗
DOE OSTI2021

Identification of localized radiation damage in power MOSFETs using EBIC imaging

The rapidly increasing use of electronics in high-radiation environments and the continued evolution in transistor architectures and materials demand improved methods to characterize the potential damaging effects of radiation on device performance. Here, electron-beam-induced current is used to map hot-carrier transport in model metal-oxide semiconductor field-effect transistors irradiated with a 300 KeV focused He + beam as a localized line spanning across the gate and bulk Si. By correlating the damage to the electronic properties and combining these results with simulations, the contribution of spatially localized radiation damage on the device characteristics is obtained. Furthermore, this identified damage, caused by the He + beam, is attributed to localized interfacial P b centers and delocalized positive fixed-charges, as surmised from simulations. Comprehension of the long-term interaction and mobility of radiation-induced damage are key for future design of rad-hard devices.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Generation of focusing ion beams by magnetized electron sheath acceleration

Abstract We present the first 3D fully kinetic simulations of laser driven sheath-based ion acceleration with a kilotesla-level applied magnetic field. The application of a strong magnetic field significantly and beneficially alters sheath based ion acceleration and creates two distinct stages in the acceleration process associated with the time-evolving magnetization of the hot electron sheath. The first stage delivers dramatically enhanced acceleration, and the second reverses the typical outward-directed topology of the sheath electric field into a focusing configuration. The net result is a focusing, magnetic field-directed ion source of multiple species with strongly enhanced energy and number. The predicted improvements in ion source characteristics are desirable for applications and suggest a route to experimentally confirm magnetization-related effects in the high energy density regime. We additionally perform a comparison between 2D and 3D simulation geometry, on which basis we predict the feasibility of observing magnetic field effects under experimentally relevant conditions.

42 ENGINEERING↗
DOE OSTI2022

Whistler waves generated by nongyrotropic and gyrotropic electron beams during asymmetric guide field reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

EIC hadron ring longitudinal stability simulations at injection

At the EIC injection energy (23.8 GeV for protons), the hadron beam experiences a relatively strong space charge field. The corresponding wakefield, in conjunction with the wakes from the RF cavities (i.e. the geometrical wakes), may become the cause of single-bunch instabilities, both longitudinal and transverse. In this report, we focus on the longitudinal stability aspect of this picture, using a simplified model for the ring.

43 PARTICLE ACCELERATORS↗