Engineering Papers⌕ Search

Engineering topics

Mitchell, Chad

Publications and source records attributed to Mitchell, Chad.

12 records
All sourcesNASA NTRSDOE OSTI
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

Chromatic Transport Models for ImpactX

This note describes the basic levels of Hamiltonian approximation (for straight-axis elements) that are used in the symplectic beam dynamics modeling code ImpactX. These include: purely linear models, models based on a paraxial (chromatic) expansion of the Hamiltonian, and models based on the exact nonlinear Hamiltonian. The focus here is on paraxial (chromatic) models.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Poisson Equation for a (General) Homogeneous d-Dimensional Ellipsoid with Applications to Beam Envelope Tracking

This note describes the solution of the free-space Poisson equation in the interior of a $d$-dimensional homogeneous ellipsoid, and the associated space charge fields. An explicit formula (\ref{Sformula}) is provided that relates the $d\times d$ matrix describing the space charge (quadratic) potential to the $d\times d$ covariance matrix of the ellipsoid. For the cases $d=2$ and $d=3$, this result is used to determine the linear map corresponding to a space charge kick, that may be used to push the beam $6\times 6$ covariance matrix during envelope tracking. The treatment of upright ellipsoids for $d=2$ and $d=3$ is well-represented in the literature. However, the approach taken here emphasizes a general ellipsoid with arbitrary correlations in any dimension. The Appendix provides a general solution of the free-space Poisson equation in dimension $d$ for a source distribution with ellipsoidal symmetry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Elliptic multipoles and the modeling of narrow-gap bend magnets in accelerators

We highlight the virtues of 2D elliptic-multipole field expansions in modeling the magnetic fields of narrow-aperture, straight-axis bending magnets with parallel faces, addressing the limitations of the conventional circular multipole series when the beam-orbit sagitta exceeds the magnet's vertical half-gap. The elliptic multipoles provide a convenient way to represent the field in all aspects of the magnet development (design, particle-tracking simulations, measurements). We propose a numerically robust method of data analysis to determine the elliptic (or circular) multipoles from stretched-wire measurements with the wire moving on an arbitrary path.

Venturini, Marco↗
DOE OSTI2024

Integrated Simulation of PIP-II at Fermilab

We describe progress towards a community software ecosystem for efficient modeling of the Fermilab PIP-II complex for design validation and virtual test stands.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

ImpactX v0.1

ImpactX is the next generation of the IMPACT-Z code. It is a s-based simulation code for modeling intense beams in particle accelerators using symplectic tracking methods and includes collective effects. It is multi-node parallel and supports modern compute hardware such as GPUs, modern algorithms such as mesh-refinement and realistic geometries (embedded boundaries).

Huebl, Axel↗
DOE OSTI2021

Design of double- and multi-bend achromat lattices with large dynamic aperture and approximate invariants

A numerical method to design nonlinear double- and multi-bend achromat (DBA and MBA) lattices with approximate invariants of motion is investigated. The search for such nonlinear lattices is motivated by Fermilab’s Integrable Optics Test Accelerator (IOTA), whose design is based on an integrable Hamiltonian system with two invariants of motion. While it may not be possible to design an achromatic lattice for a dedicated synchrotron light source storage ring with one or more exact invariants of motion, it is possible to tune the sextupoles and octupoles in existing DBA and MBA lattices to produce approximate invariants. In our procedure, the lattice is tuned while minimizing the turn-by-turn fluctuations of the Courant-Snyder actions J x and J y at several distinct amplitudes, while simultaneously minimizing diffusion of the on-energy betatron tunes. The resulting lattices share some important features with integrable ones, such as a large dynamic aperture, trajectories confined to invariant tori, robustness to resonances and errors, and a large amplitude-dependent tune-spread. Compared to the nominal NSLS-II lattice, the single- and multi-bunch instability thresholds are increased and the bunch-by-bunch feedback gain can be reduced.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Multipass Simulations of Space Charge Compensation using Electron Columns at IOTA

Defocusing repulsive forces due to self space charge fields lead to degradation of high-intensity particle beams. Being of particular concern for low- and medium-energy proton beams, they result in emittance growth, beam halo formation, and beam loss. They set stringent limits on the intensity of frontier accelerators; therefore, the mitigation of space charge effects is a crucial challenge to improve proton beam intensity. The space charge effects in a positively charged proton beam can be effectively compensated using negatively charged electron columns. In this paper, we present the results of simulations using Synergia of the Electron Column lattice for IOTA. Beam loss due to space charge effects and aperture restrictions have been studied, as well as bunch formation and matching using an adiabatic ramp of the RF cavity. The results show the need for space charge compensation, and provide the basis for integration of the Synergia and Warp codes in order to form a complete simulation of space charge compensation using an Electron Column in IOTA.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Weak convergence to equilibrium of statistical ensembles in integrable Hamiltonian systems

This article explores the long-time behavior of the bounded orbits associated with an ensemble of initial conditions in a nondegenerate integrable Hamiltonian system. Such systems are inherently nonlinear and subject to highly regular phase space filamentation that can drive the ensemble of orbits toward a stationary state. Describing the statistical ensemble by a probability density on a neighborhood of a family of invariant tori, it is proved that the probability density describing the ensemble at time t converges weakly to an invariant density as t → ∞. More generally, we provide sufficient conditions for convergence to equilibrium of a multiphase system in action-angle form. These ideas are applied to an illustrative exactly soluble example. Finally, this work is relevant for understanding the statistical mechanics of integrable and near-integrable Hamiltonian systems

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Simulation study on the emittance compensation of off-axis emitted beam in RF photoinjector

To make full use of the photocathode material and improve its quantum efficiency lifetime, it can be necessary to operate the laser away from the cathode center in photoinjectors. In RF guns, the off-axis emitted beam will see a time-dependent RF effect, which would generate a significant growth in transverse emittance. It has been demonstrated that such an emittance growth can be almost completely compensated by orienting the beam on a proper orbit in the downstream RF cavities along the injector [1]. In this paper we analyze in detail the simulation techniques used in reference [1] and the issues associated with them. The optimization of photoinjector systems involving off-axis beams is a challenging problem. To solve this problem, one needs advanced simulation tools including both genetic algorithms and an efficient algorithm for 3D space charge. In this paper, we report on simulation studies where the two codes ASTRA and IMPACT-T are used jointly to overcome these challenges, in order to optimize a system designed to compensate for the emittance growth in a beam emitted off axis.

36 MATERIALS SCIENCE↗
DOE OSTI2014

An Adaptive Unified Differential Evolution Algorithm for Global Optimization

In this paper, we propose a new adaptive unified differential evolution algorithm for single-objective global optimization. Instead of the multiple mutation strate- gies proposed in conventional differential evolution algorithms, this algorithm employs a single equation unifying multiple strategies into one expression. It has the virtue of mathematical simplicity and also provides users the flexibility for broader exploration of the space of mutation operators. By making all control parameters in the proposed algorithm self-adaptively evolve during the process of optimization, it frees the application users from the burden of choosing appro- priate control parameters and also improves the performance of the algorithm. In numerical tests using thirteen basic unimodal and multimodal functions, the proposed adaptive unified algorithm shows promising performance in compari- son to several conventional differential evolution algorithms.

43 PARTICLE ACCELERATORS↗