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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Electron beam irradiated foil expansion simulated by a coupled HYDRA/MCNP6 application

Simulating the hydrodynamic evolution of a temporally pulsed electron beam with multiple metal foils requires coupling the nearly monoenergetic, relativistic electron beam to the foil target. This deposition may be effectively calculated by Monte Carlo transport. Once this energy is deposited on a computational mesh representing the foil target, highly developed radiation-hydrodynamic codes, such as HYDRA, can describe the subsequent evolution. Currently, HYDRA does not have an internal Monte Carlo electron transport capability. To partially rectify this deficiency, a PYTHON-based scheme to link MCNP6 to HYDRA has been developed for quadrilateral meshes. Results from a seven-foil tantalum target configuration are presented to demonstrate this capability.

42 ENGINEERING↗
DOE OSTI2024

Simulation study of ion beam used to produce Mo-99

An 820 mA CW positive ion source is being developed to produce Mo-99 using the fusion of deuterium and tritium ion beams on a rotating target to produce neutrons for use in the production of radiopharmaceuticals. The ion source consists of an RF plasma source, a multi-aperture extractor, and 300 kV accelerating column. This paper will describe a simulation study of the beam through the extractor grid and the accelerator to the target. The uniformity of beam distribution on the target is an important aspect of the simulation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Neural network model of neutral beam injection in the EAST tokamak to enable fast transport simulations

The neutral beam injection (NBI) system in EAST produces energetic neutral particles, which collide with electrons and ions in tokamak plasmas and heat the plasmas through Coulomb collisions. Moreover, it drives a non-inductive source of current, due to the charge-exchange collision between neutral particles and ions, and injects toroidal torque, which generates a toroidal rotation of the plasma. The effect caused by the NBI system, such as plasma heating, current drive, total neutron rate, momentum transfer, and shine-through, are modeled by a comprehensive module called NUBEAM. However, NUBEAM is computationally intensive since it relies on Monte Carlo methods. In this work, a neural network model has been developed as a surrogate model for NUBEAM in EAST. The database for neural-network model training, validation and testing is generated by running TRANSP for experimental discharges from recent EAST campaigns (after the latest NBI upgrade) while using the NUBEAM module. Simulation results illustrate that the trained neural network has the capability of replicating the predictions made by NUBEAM while demanding a significantly shorter execution time. Finally, these results indicate that surrogate models like the one proposed in this work could enable fast transport simulations for EAST after integrating them into a control-oriented predictive code such as COTSIM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fully Kinetic Simulations of Proton-beam-driven Instabilities from Parker Solar Probe Observations

The expanding solar wind plasma ubiquitously exhibits anisotropic nonthermal particle velocity distributions. Typically, proton velocity distribution functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by the Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam population that drift with respect to each other. This configuration triggers a proton beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing perpendicular heating at the expense of the parallel temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Accelerator Design Progress for EIC Strong Hadron Cooling

The Electron-Ion Collider will achieve a luminosity of 1034 cm-2 s-1 by incorporating strong hadron cooling to counteract hadron Intra-Beam Scattering, using a coherent electron cooling scheme. An accelerator will deliver the beams with key parameters, such as 1 nC bunch charge, and 1e-4 energy spread. The paper presents the design and beam dynamics simulation results. Methods to minimize beam noise, the challenges of the accelerator design, and the R&D topics being pursued are discussed.

Wang, E.↗
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 OSTI2022

CoSyR: A novel beam dynamics code for the modeling of synchrotron radiation effects

The self-consistent nonlinear dynamics of a relativistic charged particle beam interacting with its complete self-fields is a fundamental problem underpinning many of the accelerator design issues in high brightness beam applications, as well as the development of advanced accelerators. Particularly, synchrotron radiation induced effects in a magnetic dispersive beamline element can lead to collective beam instabilities and emittance growth. A novel beam dynamics code is developed based on a Lagrangian method for the calculation of the particles’ radiation near-fields using wavefront/wavelet meshes via the Green’s function of the Maxwell equations. These fields are then interpolated onto a moving mesh for dynamic update of the beam. This method allows radiation co-propagation and self-consistent interaction with the beam in 2D/3D simulations at greatly reduced numerical errors. Multiple levels of parallelisms are inherent in this method and implemented in our code CoSyR to enable at-scale simulations of nonlinear beam dynamics on modern computing platforms using MPI, multi-threading, and GPUs. Here, the current 2D implementation of CoSyR has been used to evaluate the transverse and longitudinal coherent radiation effects on the beam and to investigate beam optics designs proposed for mitigation of beam brightness degradation in a magnetic bunch compressor. In this paper, the design of CoSyR, as well as the benchmark with other coherent synchrotron radiation models, are described and discussed. Extension of the core algorithms to 3D is possible and planned.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Preliminary Look at the LBEG & MPEG Beam Transmissions between HPSim and Operation Data

This report summarizes recent work on estimating the transmission of LANSCE H- beams from the end of the present DTL to both the PSR stripper foil (LBEG) and WNR target 4 (MPEG). The LBEG beam might be considered a more typical LINAC beam and lends itself to continuous monitoring of transmission through the various stages of the accelerator. For the MPEG, however, the widely-spaced micropulses and lack of the requisite sensitivity current monitors throughout the accelerator make these measurements extremely difficult and more uncertain. Therefore, to assist in estimating the MPEG beam transmission, beam-dynamics simulations using HPSim were employed. These beam transmission estimates from the end of the Drift Tube Linac (DTL) to Target 4 (MPEG) and the PSR stripper foil (LBEG) are vital to determine the charge requirements for LAMP’s front end. In this technote, we demonstrate three major efforts in determining the transmission: (1) Convert the WNR beamline lattice from TRANSPORT to HPSim for use in the simulation; (2) Simulate the optimization process in the Central Control Room (CCR) that brings down the losses by up to a factor of 4000 between the end of the initial physics-based tuneup phase and production beam operation; (3) Analyze operational data to deduce measured transmissions. Table 1 shows the estimated losses with HPSim and operational analysis. Finally, a better measurement and other improvements to refine the results are also proposed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Thin-Foil Heating with Intense Relativistic Electron Beams

Thin-metal foils are used in intense relativistic beam experiments for various purposes. In planning these experiments and analyzing the results it is often desirable to have an estimate of the temperature reached during the beam pulse. This can be deduced from an estimate of the energy deposited by the beam, and the specific heat capacity of the material, which depends on the temperature. The methodology for doing this is explained herein, and a metric for foil survivability under heating by an intense relativistic electron beam is suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Reconstruction of Six-Dimensional Phase Space

A phase space is a mathematical representation of all possible physical states of a system. Particle beams at Fermilab exist within a six-dimensional (6D) phase space defined by three positional components, (x, y, z) and three momentum components, (px, py, pz). To reconstruct this space implies taking measurement data from detectors and mapping out particle behavior using computational methods. The beam detectors, however, are only able to detect spatial distribution among the events of the beam, therefore being limited to positional data. Also, due to the vast number of events in a particle beam, it is extremely difficult to analyze and differentiate every single one’s behavior. However, with Machine Learning (ML), which can distinguish between patterns and map out particle behavior more efficiently. We first used the particle beam software, G4beamline, to simulate a 10,000-event muon beam, adjusting parameters such as initial momentum magnitude (p¬0) and virtual detector position. Using ten virtual detectors, we analyzed p0 values such that minimum 9,990 events were analyzed by every detector. We then input the data from these beam simulations to a C++ program, that randomly selects 100 events, and creates a 2D histogram based on spatial distribution, detector position, and event intensity. This process is repeated 100 times to create 100 histograms per p0 value. These images were then input to a modified ResNet18 Convolutional Neural Network (CNN) for training, and to predict p0 from some unseen set of histograms. The model was accurate when trained on momentum increments of 5 MeV/c and provided with denser training samples around highly variable test values. These results displayed machine learning being able to accurately predict p0 from being trained on different particle behaviors.

Shirlee, Jermain [Fermilab]↗
DOE OSTI2021

Quasioptical modeling of wave beams with and without mode conversion. IV. Numerical simulations of waves in dissipative media

We report the first quasi-optical simulations of wave beams in a hot plasma using the quasi-optical code PARADE (PAraxial RAy DEscription) [K. Yanagihara, I. Y. Dodin, and S. Kubo, Phys. Plasmas 26, 072112 (2019)]. This code is unique in that it accounts for inhomogeneity of the dissipation-rate across the beam and mode conversion simultaneously. We show that the dissipation-rate inhomogeneity shifts beams relative to their trajectories in cold plasma and that the two electromagnetic modes are coupled via this process, an effect that was ignored in the past. Additionally, we propose a simplified approach to account for the dissipation-rate inhomogeneity. This approach is computationally inexpensive and simplifies the analysis of actual experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Beam and Spin Optics Simulation Tutorials, Using Zgoubi Exercises and Solutions 2 - RHIC Run 22 Spin Transparency Experiment

High electron and ion polarizations are paramount design requirements for the future Electron Ion Collider (EIC) at BNL. A “spin transparency” mode of operation, inherited from the “figure-8” ring design [1], allows preservation and control of electron and ion spin polarization in a collider or storage ring [2, 3]. It makes the ring lattice “invisible” to the spin and allows for polarization control by small quasi-static magnetic fields with practically no effect on beam orbital parameters. Beam polarization can then be fully controlled by small adjustments of snake axis orientation (a few degrees) and snake strengths (a few percent); polarization orientation at IPs can be controlled without spin rotators, spins can be adiabatically flipped without the need of ac magnets. The transparent spin mode has been demonstrated in simulations [2] for the purpose of being tested experimentally during RHIC Run 22. This experiment was aimed at providing data to validate the principles. Yellow RHIC ring (Fig. 1.1) was configured in the transparent spin mode by aligning the axes of its two helical snakes (located respectively in the 3 o’clock and 9 o’clock sextants) parallel to the longitudinal axis. Experiment goals and expected results, ad hoc procedures and hardware requirements, are addressed in Refs. [2, 3].

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Spot Size Optimization of the Scorpius Accelerator

The Scorpius Linear Induction Accelerator (LIA) is being developed by the Advanced Sources and Diagnostics (ASD) Project and will complement other U1a capabilities by providing a multi-pulse, DARHT-class, flash-radiography capability. Commissioning activities are expected to end in Spring of 2026. All of the diagnostics discussed in this chapter are relevant to LIAs that generate high-current relativistic electron beams, and some diagnostics will also be relevant to other accelerators. The Scorpius Accelerator will be the world’s most advance LIA. Although significant technological advancements have been incorporated throughout the machine, many performance aspects of Scorpius will be similar to those of FXR and DARHT Axis-I and -II. These LIAs offer the same challenges, whether it is suppressing beam instabilities, improving diagnostic accuracy, validating computer simulations and beam tunes, or reducing beam-target interaction effects. However, the ultimate figure of merit is consistent between these LIAs, which is spot size and dose.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Innovative approaches to neutron beam optimization: A case study of FP14 at LANSCE

This paper describes an experimental method to acquire high resolution energy- and spatially- resolved neutron beam spots using the time-gated neutron imaging system with Teledyne Pi-MAX4 camera. These experimental data offer a unique opportunity for benchmarking beam spot simulations. High-quality simulations depend significantly on a high-fidelity geometry model, which can be challenging for legacy facilities. We informed our MCNPX geometry model by latest metrology survey employing Leica laser tracker ATS600. It gave us a high fidelity description of our facility geometry. Such a robust integration of novel tools and methods yields a previously unattainable level of accuracy in both predicting and capturing neutron beam spots, marking a notable advancement over traditional methods reliant on static image plates. Here, to demonstrate the practical application of these tools, we are showing a non-uniform beam spot challenge at our Flight Path 14 (FP14) at the Los Alamos Neutron Science Center (LANSCE). Our precise MCNPX prediction of beam spot shifting as function of neutron energy was confirmed by experimental beam spot measured with extremely high level of detail. Results of this research demonstrate a significant leap in neutron beam optimization at LANSCE and set a new benchmark in beam spot characterization. The advanced methods presented here have potential for adoption at similar research facilities worldwide, aiming at substantial improvement in neutron beam delivery for experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Optimization of an electron injector for operating high-intensity polarized electron beam

In this work, the simulation and optimization of an electron injector that operates a high-intensity polarized electron beam are presented. The electron injector would provide a single-bunch electron beam at a repetition rate of 1 Hz with a bunch charge of 10 nC, kinetic energy of 400 MeV and an electron spin polarization of 85%. A direct current (DC) high-voltage photogun was employed to produce the 10 nC polarized electron beam with a kinetic energy of 350 keV, beam diameter of 1.44 cm (Gaussian distribution with σ of 0.36 cm) and full bunch length of 1.3 ns (Gaussian distribution with σ of 0.325 ns). The beam was compressed to 6 ps of RMS bunch length using one 114.24 MHz standing-wave sub-harmonic buncher, two 571.2 MHz standing-wave sub-harmonic bunchers and one 2.856 GHz traveling-wave buncher. The beam was accelerated finally to 400 MeV via eight 2.856 GHz traveling-wave Linacs. We performed the beam simulation, and the simulated results showed that the optimized RMS bunch length was 4.65 ps, RMS relative energy spread was 0.48% and normalized RMS transverse emittance was 30.13 mm · mrad at a beam energy of 400 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Beam Dynamics of the Muon $g\textrm{-}2$ Experiment

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab aims to measure the muon anomalous magnetic moment $a_{\mu}$ with unprecedented precision, potentially uncovering physics beyond the Standard Model of particle physics. The result based on Runs 1-3, released in 2023, achieved a precision of 0.20 ppm. The experiment circulates muons in a storage ring, measuring $a_{\mu}$ from decay positron time and energy measurements collected with calorimeters. To achieve the required accuracy, it is crucial to measure and control the magnetic field in the ring with high precision. Beam dynamics corrections are necessary for muons not orbiting exactly in the midplane, for their oscillations, and for electric field effects. Highly accurate beam dynamics simulations are instrumental for quantifying and validating the beam dynamics corrections, ultimately improving the precision of the $a_{\mu}$ measurement and facilitating the achievement of the ambitious $70\:\mathrm{ppb}$ systematic uncertainty goal. The measured field data was incorporated into models for simulations using three codes: \texttt{gm2ringsim} (an internal Geant4-based code), \textit{COSY INFINITY}, and \textit{BMAD}. The advantages of \texttt{gm2ringsim} include using CAD-based geometry and modelling the detector effects. \textit{COSY INFINITY} is a highly accurate and efficient code that uses high-order differential-algebraic transfer maps, precise fringe field calculations, and advanced symplectification methods. Symplectification is important for maintaining the physical correctness of the muon beam behaviour with high precision over the storage time, ensuring conservation of phase space volume and preventing artificial damping or excitation of particle motion. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous BNL experiment. Analyses of data from Runs 4-6 are ongoing, with results planned for release in 2025, potentially resolving the current tension between experiment and theory.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Experimental longitudinal emittance manipulation using laser-based photoionization in the Fermilab Linac

A series of simulations and beam studies were conducted at Fermilab’s linear accelerator to evaluate the effectiveness of longitudinal emittance control via laser-induced photoionization. While similar laser techniques have been employed at Fermilab to enhance injection and extraction efficiency into the Booster, the work presented here focuses on extending these methods to bunch-by-bunch manipulation. This approach utilizes fine-scale correction of the H- bunches’ longitudinal spatial distribution. In theory, loosely confined particles in longitudinal phase space contribute to emittance growth during acceleration. By selectively removing these outlying particles through laser scraping (H- + γ → H + e-), this growth can be reduced. This report presents experimental results from both symmetric and asymmetric longitudinal scraping of H⁻ bunches in the Fermilab linac, which were subsequently injected into Booster, and evaluates the broader applicability of this method for future high-intensity accelerator operations.

Landon, Parker [Boston U.; Fermilab]↗
DOE OSTI2021

Computational approaches to Coherent Synchrotron Radiation in two and three dimensions

Coherent Synchrotron Radiation (CSR) is an important and often detrimental effect in particle accelerators. While one-dimensional models have been successfully used to design and explain the behavior of modern machines, questions remain about their domain of validity. In recent years, two- and three-dimensional models have been developed that are amenable to efficient numerical computation. This article gives an overview of CSR computation from its discovery through the present state of the art.

43 PARTICLE ACCELERATORS↗