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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Numerical simulations of colliding jets in an external wind: application to 3C 75

ABSTRACT The radio galaxy 3C 75 is remarkable because it contains a pair of radio-loud active galaxies, each of which produces a two-sided jet, with the jet beams appearing to collide and merge to the west of the galaxies. Motivated by 3C 75, we have conducted three-dimensional hydrodynamic simulations of jet collisions. We have extended previous studies by modelling the physical properties of the cluster atmosphere, including an external wind, and using realistic jet powers obtained from observational data. We are able to produce a morphology similar to that of 3C 75. The simulations imply that direct contact between the bulk jet flows on the west of the source is required to produce a morphology consistent with 3C 75. We quantify how the merging jets decelerate, how the wind deflects the jets and cocoons, the entrainment of intra-cluster material into the cocoons, the cocoon energetics, and how the jet interactions generate enstrophy. By comparing simulations of pairs of two-sided jets with those of single two-sided sources, we determine how the interaction between two bipolar jets changes their evolution. The unprecedented sensitivity and angular resolution of upcoming observatories will lead to the detection of many more complex sources at high redshift, where interacting jets are expected to be more numerous. The morphology of these complex sources can provide significant insight into the conditions in their environments.

Musoke, G.↗

DOE OSTI2023

Beam delivery and beamstrahlung considerations for ultra-high energy linear colliders

As part of the Snowmass'21 community planning excercise, the Advanced Accelerator Concepts (AAC) community proposed future linear colliders with center-of-mass energies up to 15 TeV and luminosities up to 50 × 1034 cm -2 s -1 in a compact footprint. In addition to being compact, these machines must also be energy efficient. We identify two challenges that must be addressed in the design of these machines. First, the Beam Delivery System (BDS) must not add significant length to the accelerator complex. Second, beam parameters must be chosen to mitigate beamstrahlung effects and maximize the luminosity-per-power of the machine. In this paper, we review advances in plasma lens technology that will help to reduce the length of the BDS system, and we detail new Particle-in-Cell simulation studies that will provide insight into beamstrahlung mitigation techniques. We apply our analysis to both $e$ + $e$ - and $γγ$ colliders. The challenges and solutions described in this paper are considered independently. In conclusion, a unified, self-consistent concept for a BDS system for a 15 TeV linear collider will be the subject of future work.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

First evaluation of meson and τ lepton spectra and search for heavy neutral leptons at ILC beam dump

A beam dump experiment can be seamlessly added to the proposed Inter- national Linear Collider (ILC) program because the high energy electron beam should be dumped after the collision point. The ILC beam dump experiment will provide an excellent opportunity to search for new long-lived particles. Since many of them can be produced by a rare decay of standard model particles, we evaluate spectra of the mesons and τ lepton at the decay based on the PHITS and PYTHIA8 simulations. As a motivated physics case, we study the projected sensitivity of heavy neutral leptons at the ILC beam dump experiment. The heavy neutral leptons can also be produced via deep inelastic scattering and Z boson decay at the ILC main detector, which we include in the projection. With the multi-track signal, the reach would be greatly extended in mass and coupling, even compared with the other proposed searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

The effect of laser entrance hole foil thickness on MagLIF-relevant laser preheat

The magnetized liner inertial fusion (MagLIF) scheme relies on coupling laser energy into an underdense fuel raising the fuel adiabat at the start of the implosion. To deposit energy into the fuel, the laser must first penetrate a laser entrance hole (LEH) foil which can be a significant energy sink and introduce mix. In this paper, we report on experiments investigating laser energy coupling into MagLIF-relevant gas cell targets with LEH foil thicknesses varying from 0.5 μm to 3 μm. Two-dimensional (2D) axisymmetric simulations match the experimental results well for 0.5 μm and 1 μm thick LEH foils but exhibit whole-beam self-focusing and excessive penetration of the laser into the gas for 2 μm and 3 μm thick LEH foils. Additionally, better agreement for the 2 μm-thick foil is found when using a different thermal conductivity model in 2D simulations, while only 3D Cartesian simulations come close to matching the 3 μm-thick foil experiments. The study suggests that simulations may over-predict the tendency for the laser to self-focus during MagLIF preheat when thicker LEH foils are used. This effect is pronounced with 2D simulations where the azimuthally symmetric density channel effectively self-focuses the rays that are forced to traverse the center of the plasma. The extra degree of freedom in 3D simulations significantly reduces this effect. The experiments and simulations also suggest that, in this study, the amount of energy coupled into the gas is highly correlated with the laser propagation length regardless of the LEH foil thickness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Atomic understanding of structural deformations upon ablation of graphene

We investigate the atomic rearrangement in graphene under femtosecond pulse illumination with reactive molecular dynamics simulations and compare with ultra-fast laser ablation experiments. To model the impact of the laser pulse irradiation, heat is locally applied to a selected area of the graphene layer and the resulting structural deformation is simulated as a function of time, providing a detailed understanding of the bond breaking process under laser illumination and subsequent re-equilibration afterAQ3 the pulse is turned off. Analysis of the atomic dynamics indicates that the types of defects formed depend on the pulse energy and exposure duration. By varying the exposed area, we determine that the shape of the ablated area is not only a function of the pulse energy, but also of the beam spot size and pulse repetition. Furthermore, we apply a machine learning approach to extrapolate our simulated data to experimental length scales and reproduce the trends in ablated area as a function of temperature. Furthermore, our study provides a first step towards understanding the design parameters for graphene nano-patterning.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

DOE OSTI2021

Evaluation of the accuracy of stopping and range of ions in matter simulations through secondary ion mass spectrometry and Rutherford backscattering spectrometry for low energy heavy ion implantation

The freely available “Stopping and Range of Ions in Matter” (SRIM) code is used for evaluating ion beam ranges and depth profiles. We present secondary ion mass spectrometry and Rutherford backscattering experimental results of Si samples implanted with low energy Sb ions to evaluate the accuracy of SRIM simulations. We show that the SRIM simulation systematically overestimates the range by 2–6 nm and this overestimation increases for larger ion implantation energy. For the lowest energy implantation investigated, here we find up to a 25% error between the SRIM simulation and the measured range. Finally, the ion straggle shows excellent agreement between simulation and experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Searching for the Most Harmful Field Errors in the HSR IR Superconducting Magnets

In this project, we improve beam stability for the Electron-Ion Collider. Magnetic field errors can reduce beam stability, making it essential to identify the field errors that have the greatest impact on accelerator performance. However, this is particularly challenging because beam stability depends on the complex interactions of many magnetic field errors, resulting in a high-dimensional and nonlinear optimization problem. We determine which field errors are the most influential for the large physical aperture superconducting magnet B2PF, a critical magnet in the Interaction Region (IR) in the Hadron Storage Ring (HSR). We complete and analyze nearly 30,000 simulations on the Brookhaven National Laboratory Linux Cluster by varying 18 nonlinear magnetic field errors. We evaluate beam stability using the dynamic aperture and the tune diffusion. We identify the field errors that most strongly influence beam stability and establish quantitative field error tolerances that improve accelerator performance.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Celeritas R&D Report: Accelerating Geant4

Celeritas is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a Graphics Processing Unit (GPU)-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.4, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling straightforward integration of Celeritas into the high energy physics (HEP) experimental frameworks CMSSW and ATLAS FullSimLight. On the Perlmutter supercomputer, Celeritas performs EM physics between 3× and 18× faster using the machine’s Nvidia GPUs compared to using only CPUs, corresponding to an electrical power efficiency up to a factor of 5. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.7–2.2× on GPU and 1.2× on CPU. In a CMS test application using tt¯ events and the prototype Run 4 configuration, compared to Geant4 CPU, Celeritas with a Nvidia A100 improves overall throughput up to a factor of 2.7× but cannot be efficiently shared with more than 8 cores.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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 OSTI2022

Prospects for Machine Learning and Pulse Shaping on the Scorpius Accelerator [Poster]

The Advanced Sources and Detectors (ASD) project aims to build Scorpius, a multi pulse linear induction accelerator capable of delivering a 1.4 kA electron beam at energies up to 24 MeV. One of the primary advancements of Scorpius is the use of solid state pulsed power (SSPP) to provide flexibility in pulse shaping by independently triggering 45 individual stages stacked in each of 984 line replaceable units (LRU), with 168 LRUs dedicated to the injector. By leveraging circuit modeling of each LRU stage, a machine learning model of the SSPP will be developed to allow for optimization of the pulse shape, including pulse flattening and reflection mitigation. Particle-in-cell simulations of Scorpius have, for example, demonstrated that reducing reflections during multi-pulse operation mitigates beam spill by preventing the production of off-energy electrons between pulses, thereby abating stimulated ion desorption from the wall and beam charge neutralization. This machine learning model will be validated and tuned with experimental data collected from the Scorpius injector and Integrated Test Stand

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

The positron arm of a plasma-based linear collider

Plasma-based acceleration of electrons has produced high-energy beams at high accelerating gradients with a narrow energy spread and high efficiency both in experiments and simulations. Furthermore, it is now being considered as a complementary approach to the use of radiofrequency cavities in next-generation lepton accelerators. However, compared with electrons, plasma-based positron acceleration is at the present time much less advanced. Although high-gradient positron acceleration in a plasma has been achieved, we are one to three orders of magnitude away from delivering the high-quality positron beams needed for a future high-energy linear collider. Here we review the status of plasma-based acceleration of electrons and positrons and discuss the prospects for substantial progress towards developing the positron arm of a plasma-based electron–positron linear collider in the next decade.

Plasma-based accelerators↗

DOE OSTI2022

A Non-Cooperative Game-Based Distributed Beam Scheduling Framework for 5G Millimeter-Wave Cellular Networks

Here, this paper studies the problem of distributed beam scheduling for 5G millimeter-Wave (mm-Wave) cellular networks where base stations (BSs) belonging to different operators share the same spectrum without centralized coordination among them. Our goal is to design efficient distributed scheduling algorithms to maximize the network utility, which is a function of the achieved throughput by the user equipment (UEs), subject to the average and instantaneous power consumption constraints of the BSs. We propose a Media Access Control (MAC) and a power allocation/adaptation mechanism utilizing the Lyapunov stochastic optimization framework and non-cooperative games. In particular, we first decompose the original utility maximization problem into two sub-optimization problems for each time frame, which are a convex optimization problem and a non-convex optimization problem, respectively. By formulating the distributed scheduling problem as a non-cooperative game where each BS is a player attempting to optimize its own utility, we provide a distributed solution to the non-convex sub-optimization problem via finding the Nash Equilibrium (NE) of the game whose weights are determined optimally by the Lyapunov optimization framework. Finally, we conduct simulation under various network settings to show the effectiveness of the proposed game-based beam scheduling algorithm in comparison to that of several reference schemes.

42 ENGINEERING↗

DOE OSTI2023

Improved calibration of rf cavities for relativistic electron beams: Effects of secondary corrections and experimental verification

In the aspect of longitudinal beam bunching, the bunching strength can be controlled by the rf cavity phase and voltage. However, these machine parameters are different from those that interact with the beam itself. In order to gain control of the beam-cavity interaction, cavity calibration must be performed. Furthermore, it relies on fitting the beam energy gain versus cavity phase to a calibration function. Under the conventional assumption of relativistic beam conditions, the calibration function is a first harmonic sinusoidal function (a sinusoidal function with a period of 2π). However, this expression is insufficient for a high-voltage bunching cavity. Due to beam acceleration inside the cavity, an energy bias and a second harmonic function should be included to modify the conventional calibration function, even for a relativistic electron beam. In this paper, we will derive this modification and provide a comparison to both the Coherent Electron Cooling Experiment and the IMPACT T simulation, respectively.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Full-Cycle Simulations of the Fermilab Booster

The Proton Improvement Plan phase II (PIP-II) project currently under construction at FNAL will replace the existing 400 MeV normal conducting linac with a new 800 MeV superconducting linac. The beam power in the downstream rapid-cycling Booster synchrotron will be doubled by raising the machine cycle frequency from 15 to 20 Hz and by increasing the injected beam intensity by a factor 1.5. This has to be accomplished without raising uncontrolled losses beyond the administrative limit of 500 W. In addition, slip-stacking efficiency in the Recycler, the next machine in the accelerator chain, sets an upper limit on the longitudinal emittance of the beam delivered by the Booster. As part of an effort to better understand potential losses and emittance blow-up in the Booster, we have been conducting full cycle 6D simulations using the code PyORBIT. The simulations include space charge, wall impedance effects and transition crossing. In this paper, we discuss our experience with the code and present representative results for possible operational scenarios.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Full Cycle Simulations of The Fermilab Booster

43 PARTICLE ACCELERATORS↗

DOE OSTI2026

Fermilab PIP II machine protection system digitized data noise elimination scheme and its FPGA implementation

In Fermilab's PIP-II machine protection system, beam loss signals from various detectors are digitized at 125 MS/s. Noise from both high-frequency sources and low-frequency 60 Hz AC power equipment can contaminate the data. To suppress noise across these ranges especially 60 Hz and its harmonics, which overlap with beam loss signal frequencies advanced digital processing beyond standard filtering is required. Several real-time functional blocks were simulated and tested on an FPGA: (1) a dual time-constant discharging integrator filter, (2) a de-ripple baseline extraction and storage block, and (3) a fast-recovery discharging integrator. The nonlinear IIR integrator filter removes high-frequency noise and feeds into the baseline extractor. Upon detecting abrupt beam loss, it switches to a longer time constant to prevent baseline distortion. The de-ripple block calculates a valid baseline by averaging over multiple 60 Hz periods, storing results in a 4096-word FPGA RAM. This baseline is subtracted from raw data before integration by the fast-recovery block, which resets quickly after use. All blocks achieved expected performance.

Wu, J. [Fermilab] (ORCID:0000000344329521)↗

DOE OSTI2025

Benchmarking fast-neutron response of 7 Li-enriched C 7 LYC scintillators: Simulations and experiment

New insights into the fast neutron response of 7 Li-enriched C 7 LYC scintillators have been obtained through detailed Geant4 simulations of 35 Cl(n,p) and 35 Cl(n, α) reactions that are primarily responsible for the fast neutron response of C 7 LYC. The light output of a 1” x1” C 7 LYC detector was simulated using the Birks equation for light quenching in scintillators. The simulations were benchmarked against neutron scattering data from a 12 C target obtained with a white neutron beam at the Los Alamos LANSCE facility and fourteen 1”x1” C 7 LYC detectors. Incident neutron energies extracted from time-of-flight measurements, and appropriate recoil corrections for the different detector angles were used to map the response of C 7 LYC to mono-energetic neutrons. Resonances in the 35 Cl(n,p) reaction cross-section below 3 MeV incident neutron energy were directly observed in C 7 LYC. The experimental data and simulations are compared with ENDF databases, as well as recent direct cross-section measurements of neutrons scattered from a 35 Cl target. In conclusion, the results suggest a scrutiny of the 35 Cl(n,p) and 35 Cl(n, α) cross sections at low energy in the ENDF library.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗