Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Beam dynamics”

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

A Novel Two Stage Collimation Unit for Fermilab Booster

A new two-stage collimation unit (2SC) for Fermilab Booster will be installed during 2024 summer shutdown. It is a supplementary collimator for existing single stage Booster collimators. Design details of this 2SC adapted to Booster conditions are described. Results of beam dynamics simulations on collimation efficiency and evaluation of collimator shielding with MARS code on this new system are presented. The analysis on beam dynamics shows that the system will help to collimate up to about 1% of the beam in a confined way and, prompt and residual activation meet Fermilab Radiological Control limits in current as well as during PIP-II era high intensity operation of the Booster.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Extending JuTrack’s capabilities to the FRIB accelerator to enhance online modeling

JuTrack is a Julia-based accelerator modeling and tracking package that utilizes compiler-level automatic differentiation (AD) to enable fast and accurate derivative calculations. While JuTrack provides a solid foundation for beam dynamics simulations, its capabilities must be extended to support the Facility for Rare Isotopes (FRIB) linac. This includes modeling heavy-ion linac accelerator components such as the liquid-lithium charge stripper, which facilitates efficient acceleration by remove electrons from heavy isotopes, and incorporating multi-charge state acceleration tracking, which allows for charge-dependent beam dynamics. These extensions address challenges such as the beam matching and optimization of multi charge state through various accelerating structures and beam-material interaction modeling while maintaining the auto differentiation capability. This work focuses on adapting JuTrack to incorporate these elements, enhancing its online modeling abilities. We present modifications to JuTrack’s framework and demonstrate their performance in FRIB simulations.

Accelerator Physics↗
DOE OSTI2025

On Final Results From the Muon $g\mathrm{-}2$ Experiment at Fermilab

The Muon g-2 Experiment at Fermilab has measured the muon anomalous magnetic moment, a_mu, with unprecedented precision, leveraging a dataset from Runs 1 6 that is 21 times larger than its Brookhaven predecessor. This talk will present the experiment s final result, which serves as a benchmark for testing the Standard Model with high precision. The measurement was performed in a storage ring with a highly uniform magnetic field, where precise beam dynamics understanding is vital to determine the muon anomalous precession frequency. Key to this effort were advanced simulation tools including COSY INFINITY, precise fringe field modeling, and calculations of beam dynamics characteristics like tunes and chromaticity. We will conclude by exploring how the experiment's findings reshape our understanding of particle physics.

Valetov, Eremey [Michigan State U.]↗
DOE OSTI2025

Leveraging prior mean models for faster Bayesian optimization of particle accelerators

Tuning particle accelerators is a challenging and time-consuming task that can be automated and carried out efficiently using suitable optimization algorithms, such as model-based Bayesian optimization techniques. One of the major advantages of Bayesian algorithms is the ability to incorporate prior information about beam physics and historical behavior into the model used to make control decisions. In this work, we examine incorporating prior accelerator physics information into Bayesian optimization algorithms by utilizing fast executing, neural network models trained on simulated or historical datasets as prior mean functions in Gaussian process models. We show that in ideal cases, this technique substantially increases convergence speed to optimal solutions in high-dimensional tuning parameter spaces. Additionally, we demonstrate that even in non-ideal cases, where prior models of beam dynamics do not exactly match experimental conditions, the use of this technique can still enhance convergence speed. Finally, we demonstrate how these methods can be used to improve optimization in practical applications, such as transferring information gained from beam dynamics simulations to online control of the LCLS injector, and transferring knowledge gained from experimental measurements across different operating modes, such as accelerating different ion species at the ATLAS heavy ion accelerator.

43 PARTICLE ACCELERATORS↗
NASA NTRS1992

A piecewise continuous Timoshenko beam model for the dynamic analysis of tapered beam-like structures

Distributed parameter modeling offers a viable alternative to the finite element approach for modeling large flexible space structures. The introduction of the transfer matrix method into the continuum modeling process provides a very useful tool to facilitate the distributed parameter model applied to some more complex configurations. A uniform Timoshenko beam model for the estimation of the dynamic properties of beam-like structures has given comparable results. But many aeronautical and aerospace structures are comprised of non-uniform sections or sectional properties, such as aircraft wings and satellite antennas. This paper proposes a piecewise continuous Timoshenko beam model which is used for the dynamic analysis of tapered beam-like structures. A tapered beam is divided into several segments of uniform beam elements. Instead of arbitrarily assumed shape functions used in finite element analysis, the closed-form solution of the Timoshenko beam equation is used. Application of the transfer matrix method relates all the elements as a whole. By corresponding boundary conditions and compatible conditions a characteristic equation for the global tapered beam has been developed, from which natural frequencies can be derived. A computer simulation is shown in this paper, and compared with the results obtained from the finite element analysis. While piecewise continuous Timoshenko beam model decreases the number of elements significantly; comparable results to the finite element method are obtained.

Shen, Ji Yao↗
DOE OSTI2025

Study of Quadrupole Fringe Fields in the Interaction Region of the Hadron Storage Ring of the Electron Ion Collider

Fringe fields in quadrupole magnets are usually neglected in studies of beam dynamics at accelerators. However, the extreme optical parameters present in the final focus of a collider such as the Electron–Ion Collider (EIC) may give rise to effects that should not be overlooked. The calculation of quadrupole fringe fields presented in this study follows the procedure outlined in Ref. [1], specialized to the case of a straight reference orbit (i.e., with no dipole field component). A right-handed Cartesian coordinate system is employed, with the $z$-axis aligned with the quadrupole axis and $x$ and $y$ denoting the horizontal and vertical transverse coordinates, respectively. The magnetic quadrupole field gradient, $G = \partial B_y / \partial x$, transitions from its peak value inside the quadrupole—where it is nearly independent of the longitudinal coordinate $z$—to zero at some distance beyond the magnet edge. Consequently, $G$ is treated as a function of $z$. The region over which this variation occurs is defined as the quadrupole fringe field region. The study begins with the development of a description of the magnetic field in the fringe region using a power-series expansion in the transverse coordinates $x$and $y$, consistent with the longitudinally varying gradient. A model for the $z$-dependence of the gradient is then proposed and adjusted to reproduce magnetic field data obtained from three-dimensional field calculations. To evaluate the impact of the fringe fields on beam dynamics, the corresponding vector potential is derived and incorporated into a Hamiltonian formulation of particle motion. The significance of the fringe fields is quantified by calculating the amplitude-dependent tune shift from the Hamiltonian. Using linear beam optics parameters of the Hadron Storage Ring (HSR) of the EIC, the tune shift due to the fringe fields of all quadrupole magnets in the IR-6 interaction region is evaluated. Finally, the resulting tune shifts are compared with those arising from other nonlinear field components present in the HSR.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Digital twin framework for PIP-II linac: AI-driven multi-scale modeling from ion source to 800 MeV

The PIP-II superconducting linac at Fermilab is designed to deliver multi-megawatt proton beams for neutrino physics and other high-intensity applications. To expedite commissioning and enhance operational reliability, we have developed an EPICS-based data flow framework that seamlessly integrates digital twins (DT) with physical twins (PT). These digital twins comprise high-fidelity beam dynamics models or data-driven surrogate models connected to their physical counterparts through real-time diagnostics and advanced machine-learning algorithms.Central to this framework is Linac_Gen, an accelerated simulation tool that incorporates convolutional neural networks, random forests, and genetic algorithms to provide up to a tenfold speedup in optimizing the accelerator geometry model. An EPICS translator layer ensures interoperability by efficiently mapping lattice parameters across diverse simulation platforms.Our EPICS-based framework supports multiple operational modes—monitoring, passive learning, closed-loop control, and online learning—covering the entire machine lifecycle. By leveraging HPC resources and multi-objective optimization techniques, the digital twin enables adaptive trajectory correction, real-time fault detection, and predictive modeling of beam stability. This comprehensive approach paves the way for robust, high-intensity operation and data-driven accelerator R&D at Fermilab.

Pathak, Abhishek [Fermilab]↗
DOE OSTI2026

Analytic Chromaticity Formulas for the Muon g -2 experiment at fermilab

The Muon g-2 Experiment (E989) at Fermilab measured the muon anomalous magnetic moment aμ with unprecedented precision of 127 ppb using a storage ring. The experiment pushed the limits in terms of accuracy and systematic errors. Achieving the 78-ppb total systematic uncertainty required precise beam dynamics modeling and corrections for effects on the measured muon precession frequency. Here, we derived the first analytic aberration formulas up to the second order for the muon g-2 storage ring’s combined-function electrostatic quadrupoles (superimposed magnetic dipole and electric quadrupole fields) using an order-by-order perturbation method. From these, we obtained the exact chromaticity formulas for three ring models of different granularity and validated them against numerical calculations using COSY INFINITY, achieving analytic-numerical agreement to ⁠$\mathscr{O}(10^{-10})$. This work resolved discrepancies between previous approximate derivations and provided essential beam dynamics results for Runs 4-6 analyses. We also calculated nonlinear chromaticities up to ninth order. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous muon g-2 experiment at Brookhaven National Laboratory, with the final result announced in June 2025.

aberrations↗
DOE OSTI2026

Status of 22 GeV CEBAF Upgrade Using FFA Arcs

Jefferson Lab is exploring an upgrade to extend CEBAF’s energy to ~22 GeV within the existing tunnel by incorporating Fixed-Field Alternating-Gradient (FFA) arcs. In this scheme, the two lowest-energy electromagnetic arcs are removed, the remaining arcs are reassigned to lower energies, and a new pair of high-energy FFA arcs enables six additional recirculations using the existing SRF linacs. The non-scaling FFA lattice employs Halbach-derived permanent magnets with provisions for dipole, quadrupole, and higher-order components, offering large momentum acceptance and reduced operating cost. Implementing this concept requires updated linac optics, modifications to the remaining arcs, and a redesigned electromagnetic switchyard. Ongoing beam-dynamics studies address synchrotron-radiation–induced energy loss and emittance growth, along with strategies for their mitigation.

Ogur, S. [Thomas Jefferson National Accelerator Fa↗
DOE OSTI2026

Dynamic Multipoles in the RCS

The Rapid Cycling Synchrotron (RCS) of the Electron--Ion Collider (EIC) will accelerate polarized electrons up to 18 GeV at a repetition rate of 1 Hz. Dynamic multipoles driven by eddy currents in the vacuum chamber can affect beam dynamics and polarization preservation during the energy ramp. Building on the analytical treatment of S.Y. Lee for a ramped normal-conducting dipole magnet, we use conformal transformations to derive compact expressions for dynamic multipoles in higher-order ramped magnets with circular beam pipes. We further derive analytical frequency-domain transfer functions for both self-response and cross-response multipoles, providing compact dynamical models for their evolution during arbitrary ramp waveforms. The resulting analytical predictions are benchmarked against finite-element simulations, including the effects of practical pole geometry and beam-pipe misalignment. Finally, the expected dynamic multipoles in the RCS magnets are estimated and their implications for beam dynamics and polarization preservation are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

Efficient accelerator modeling and particle tracking are key for the design and configuration of modern particle accelerators. In this work, we present JuTrack, a nested accelerator modeling package developed in the Julia programming language and enhanced with compiler-level automatic differentiation (AD). With the aid of AD, JuTrack enables rapid derivative calculations in accelerator modeling, facilitating sensitivity analyses and optimization tasks. Here we demonstrate the effectiveness of AD-derived derivatives through several practical applications, including sensitivity analysis of space-charge-induced emittance growth, nonlinear beam dynamics analysis for a synchrotron light source, and lattice parameter tuning of the future Electron-Ion Collider (EIC). Through the incorporation of automatic differentiation, this package opens up new possibilities for accelerator physicists in beam physics studies and accelerator design optimization.

43 PARTICLE ACCELERATORS↗
NASA NTRS1978

Experiments on current-driven three-dimensional ion sound turbulence. I - Return-current limited electron beam injection. II - Wave dynamics

Pulsed electron beam injection into a weakly collisional magnetized background plasma is investigated experimentally; properties of the electron beam and background plasma, as well as the low-frequency instabilities and wave dynamics, are discussed. The current of the injected beam closes via a field-aligned return current of background electrons. Through study of the frequency and wavenumber distribution, together with the electron distribution function, the low-frequency instabilities associated with the pulsed injection are identified as ion acoustic waves driven unstable by the return current. The frequency cut-off of the instabilities predicted from renormalized plasma turbulence theory, has been verified experimentally.

Stenzel, R. L.↗
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

Simulation and measurement of beam-induced heating of ceramic vacuum chambers

In this article, we summarize recent theoretical and experimental studies of the impedance and beam-induced heating of titanium-coated ceramic vacuum chambers used in the NSLS-II injection kickers. The impedance was calculated using the field matching theory assuming planar approximation and compared with the mpedanceake2 code. For the coating thickness of a few microns, we demonstrated that the beam-induced power is dissipated in the titanium coating and that the longitudinally averaged two-dimensional power density is approximated by an analytical expression, thus allowing the use of a simplified model of the power density as input for the code to simulate the temperature distribution with realistic nonuniform thickness of the Ti coating. For a few values of the NSLS-II beam current, we measured the beam-induced heating of two ceramic chambers using thermal sensors installed along the chamber and compared the measurement results with the simulations. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
NASA NTRS1991

Hierarchic Extensions in the Static and Dynamic Analysis of Elastic Beams

Approximate solutions of static and dynamic beam problems by the p-version of the finite element method are investigated. Within a hierarchy of engineering beam idealizations, rigorous formulations of the strain and kinetic energies for straight and circular beam elements are presented. These formulations include rotating coordinate system effects and geometric nonlinearities to allow for the evaluation of vertical axis wind turbines, the motivating problem for this research. Hierarchic finite element spaces, based on extensions of the polynomial orders used to approximate the displacement variables, are constructed. The developed models are implemented into a general purpose computer program for evaluation. Quality control procedures are examined for a diverse set of sample problems. These procedures include estimating discretization errors in energy norm and natural frequencies, performing static and dynamic equilibrium checks, observing convergence for qualities of interest, and comparison with more exacting theories and experimental data. It is demonstrated that p-extensions produce exponential rates of convergence in the approximation of strain energy and natural frequencies for the class of problems investigated.

Watson, Robert A.↗
DOE OSTI2024

3D small-gain formula allowing strong focusing and harmonic lasing for a ring-based x-ray free electron laser oscillator

We present a detailed derivation of a formula for the small-gain calculation for an x-ray free electron laser oscillator (XFELO) based on a medium-energy (3–4 GeV) storage ring. We found harmonic lasing and strong focusing are essential for this beam energy range. Taking the small-signal low-gain formula developed by Kim and his colleagues, we modified it in such a way that the gain can be calculated without the “no focusing approximation,” and a strong focusing can be applied, as well as harmonic lasing. In this formula, the gain is represented as a product of two factors with one of them depending only on the harmonic number, undulator period, and gap. Using this factor, we show that it is favorable to use harmonic lasing to achieve hard x-ray FEL working in the small-signal low-gain regime with the medium-energy electron beam. Our formula also allows FEL optimization by varying the vertical gradient of the undulator, the vertical dispersion, and the horizontal and vertical focusing, independently. As an example, we applied this formula to study the feasibility of an XFELO option for the National Synchrotron Light Source II (NSLS-II) upgrade. Since a quite high peak current is required for the FEL, collective effects of beam dynamics in medium-energy synchrotrons significantly affect the electron beam parameters. We carried out a multiparameter optimization taking collective effects into account. Note, even though our example is for a ring-based XFELO at 3 to 4 GeV, the formula and, in particular, the approach developed here may be applied to other types of FELs. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2024

Using a particle-in-cell model for accelerator control room applications

Many accelerator control rooms rely on envelope models to simulate beam dynamics because they are fast and accurate at tracking the beam core. Particle-in-Cell (PIC) models, however, can track particles inside and outside the core and, with the improvements of computers, are now fast enough to be used in control rooms. The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is currently developing a tool to use a Particle-in-Cell model for control room applications. This report covers the progress so far and the future goals of using PyORBIT, a Particle-in-Cell simulation model, in the SNS control room.

Cathey, Brandon↗