Search NASA⌕ Search

SEARCH · Search NASA

Results for “tracking code”

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 37 records · Page 2

Machine learning surrogate for charged particle beam dynamics with space charge based on a recurrent neural network with aleatoric uncertainty

In this work, we develop a machine learning (ML) model with aleatoric uncertainty for the low energy beam transport (LEBT) region of the LANSCE linear accelerator in which we model the transport of a space-charge-dominated 750 keV proton beam through a lattice of 22 quadrupole magnets. Our ML model is developed based on data generated by a Kapchinsky–Vladimirsky (KV) envelope model of beam transport. We show that a recurrent neural network can be used as a dynamical surrogate model for fast prediction of the LEBT beam envelope. Furthermore, we endow the model with the prediction of aleatoric uncertainty and compare three different approaches. We demonstrate that the ML-based uncertainty quantification models are well calibrated and produce good estimates of the regions where the model is less certain about its predictions. This ML framework is a necessary step in the development of a real-time virtual diagnostic tool with uncertainty quantification that can be integrated into more complex downstream tasks (e.g., adaptive control or learning flexible control policies via reinforcement learning) for improved efficiency in beam operations. In future work, we plan to expand on this preliminary study by considering more realistic envelope models that include longitudinal momentum spread and dispersive effects in bending magnets, as well as particle tracking codes with 3D space charge (such as and ). Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

LiTrack

LiTrack is a fast 1D longitudinal phase-space tracking code written in MATLAB, which includes standard cavities, chicanes, wakefields, steady-state CSR, and longitudinal space charge.

Lou, William↗

PDC VOL Connector (PDCVOL) v0.1

PDCVOL is an extension software originated from the Proactive Data Containers (PDC) software (LBNL Track Code 2017-098), it is an HDF5 VOL connector that allows HDF5 applications to store and access data with the PDC system without the need to modify their code. The application can take advantage of the PDC optimization features and achieve better I/O performance.

Byna, Suren↗

Transverse BBU Suppression with Feedback for Energy-Recovery Linacs

The Electron–Ion Collider (EIC) achieves its design luminosity by cooling the ion beam with a high-current electron beam generated in an energy-recovery linac (ERL). The baseline ERL lattice employs a BNL five-cell cavity that is frequency-scaled to 197 MHz, 591 MHz, and 1.773 GHz, raising concerns about multibunch beam-breakup (BBU) instabilities. Threshold currents for each frequency option are established with two independent BBU tracking codes, providing cross-validated operating margins that guide cavity selection. To further increase the current limit, we incorporate a digital transverse feedback (FB) system that targets the dominant higher-order modes (HOMs). Simulations show the FB raises the BBU threshold by roughly an order of magnitude: power-spectral-density analysis of the beam centroid at the linac exit confirms strong suppression of the aliased HOM peak, although some spectral growth appears at secondary frequencies. This combined study quantifies baseline BBU limits, demonstrates effective active mitigation, and charts a practical path toward robust, high-current operation of the EIC ERL.

Accelerator Physics↗

3D Tracking of a Single Electron in IOTA

High-resolution observations of single-particle dynamics have potential as a powerful tool in the diagnostics, tuning and design of storage rings. We are presenting the results of experiments with single electrons that were conducted at Fermilab’s IOTA ring to explore the feasibility of this approach. A set of sensitive, high-resolution digital cameras was used to detect the synchrotron radiation emitted by an electron, and the resulting images were used to reconstruct the time evolution of oscillation amplitudes in all three degrees of freedom. From the evolution of the oscillation amplitudes, we deduce transverse emittances, momentum spread, damping times, beam energy and estimated residual-gas density and composition. To our knowledge, this is the first time that the dynamics of a single particle in a storage ring has been tracked in all three dimensions. We discuss farther development of a single particle diagnostics that may allow reconstruction of its turn-by-turn coordinates over macroscopic periods of time facilitating ultra-precise lattice diagnostics and direct benchmarking of tracking codes.

43 PARTICLE ACCELERATORS↗

Improving the Operational Lifetime of the CEBAF Photo-Gun by Anode Biasing

The operating lifetime of GaAs-based photocathodes in DC high voltage electron photo-guns is dominated by the ionization rate of residual beamline gas molecules. In this work, experiments were performed to quantify the improvement in photocathode charge lifetime by biasing the photo-gun anode with a positive voltage, which repels ions generated downstream of the anode. The photo-cathode charge lifetime improved by almost a factor of two when the anode was biased compared to the usual grounded configuration. Simulations were performed using the particle tracking code General Particle Tracer (GPT) with a new custom element. The simulation results showed that both the number and energy of ions play a role in the pattern of QE degradation. The experiment results and conclusions supported by GPT simulations will be presented.

Yoskowitz, J. T.↗

FFA@CEBAF beam transport error and tolerance simulation studies

The Continuous Electron Beam Accelerator Facility (CEBAF) is a 12 GeV recirculating electron accelerator at the Thomas Jefferson National Accelerator Facility (JLAB). Major upgrades to the accelerator are being investigated which include a new 650 MeV injection beamline and state-of-the-art fixed-field alternating (FFA) gradient recirculation arcs. The upgrade will extend the energy of the electron beam to over 20 GeV. In this paper, we provide an error and tolerance simulation study of the amended beam optics transport of the existing accelerator tuned for 22 GeV operation. The study is conducted with the particle tracking codes elegant and Bmad in two parts. In the first part, we treat each section of the accelerator (electromagnetic arcs and linacs) modularly with ideal conditions at the beginning. The second part is a pseudo start-to-end (S2E) simulation with accumulated errors propagating from one beamline to the next.

Accelerator Physics↗

SINGLE BUNCH TRACKING ON THE TEN-PASS ER@CEBAF ENERGY RECOVERY BEAMLINE

The proposed ten-pass energy recovery linac (ERL) demonstration (five accelerating, five decelerating) at the CEBAF accelerator, ER@CEBAF, involves a multi-GeV en- ergy range of a continuous electron beam. New CEBAF transverse optics were designed for this ERL demonstration. This redesign incorporates additional components in Arc A, including a path length chicane and new quadrupoles to ensure proper dispersion localization. The new five energy recovery passes with a shared arc transport scheme chal- lenge the overall beamline optics design, including large beta functions in the CEBAF spreaders and recombiners. Here we discuss results of bunch tracking performed using the elegant tracking code for the full ER@CEBAF beamline

Neththikumara, I.↗

ParFlow Sand Tank: A tool for groundwater exploration

The ParFlow Sand Tank model is an open source application designed to allow users to interactively simulate and visualize groundwater movement through the subsurface. The app is designed for both research and education; teaching hydrogeology concepts and making it easy explore and run sophisticated groundwater simulations. Our goal is to support increased accessibility and usability of research grade hydrology tools for research and teaching. The Sand Tank application simulates groundwater and surface water fluxes as well as contaminant transport in real time using the integrated physical hydrology model ParFlow (Kollet & Maxwell, 2006; Maxwell & Miller, 2005; Osei-Kuffuor et al., 2014) and the particle tracking code EcoSlim (Maxwell et al., 2019). ParFlow is a numerical hydrology model that simulates spatially distributed groundwater and surface water flow. It is a well established research tool with more than 90 publications documenting its development use to advance our understanding of groundwater dynamics and groundwater surface water interactions from the hillslope to the continental scale e.g. (Condon et al., 2020; Condon & Maxwell, 2019; Maxwell & Condon, 2016). It is designed for efficient parallel computation and has been run on many platforms spanning from laptops to supercomputers. However, one of the challenges of ParFlow is that it requires significant training and hydrologic expertise to develop simulations. The Sand Tank application makes this model accessible to anyone for education and exploration. Our application uses ParFlow for its simulation backend and ParaView for the data loading and processing. The communication infrastructure relies on the ParaViewWeb framework. We use model templates deployed in Docker images to setup the Sand Tank framework. Users can build the application locally or interact with it through our web deployment. When interacting with a template users can interactively change model parameters like subsurface processes or pump/inject water into the subsurface and watch the system respond to their changes in real time as the simulation runs. Additionally, our template setup will allow more advanced users to build custom templates of increasing complexity for both research and educational purposes.

54 ENVIRONMENTAL SCIENCES↗

Instabilities driven by the fundamental crabbing mode

The Electron Ion Collider (EIC) requires crab cavities to correct for crossing angle effects. These cavities produce a time varying horizontal kick that makes the collision head on in the center of mass frame of the bunches. This results in a luminosity that is very close to the the luminosity obtained without a crossing angle. Producing the kick requires a transverse RF field. This field requires a superconducting resonator and reasonable power levels require loaded quality factors of order a million. This results is a very large transverse impedance for the crabbing mode. The horizontal beta function is large and this, coupled with the large transverse impedance of the resonator can lead to strong instabilities. In the first section of the paper a simple model is used to describe the situation without feedback. There are large ranges of tune leading to strong instabilities in both the electron and hadron storage rings of the EIC. These estimates are confirmed with simulations. In the second section, RF feedback is used to reduce the apparent impedance of the crabbing mode. Simulating the actual feedback would require both a significant rewrite of the tracking code and a significant reduction in computational speed. Instead, the apparent impedance is fitted using one pole filters, which allows use of the existing simulation code. In the third section a transverse damping system is outlined.

43 PARTICLE ACCELERATORS↗

Acceptance Testing of the Scorpius Sweep Magnet

A beam-dump magnet for Scorpius has been defined by LANL and designed and built by Radiabeam Technologies. Initial verification of the Radiabeam design used a magnetic field model provided by Radiabeam and analyzed by LANL tools during September 2021. Simulations indicated that the design was sufficient to meet deflection requirements to protect sensitive equipment in the case of a failed target. In January 2022, the magnet was fabricated and Radiabeam conducted measurements of the magnetic field as a function of position. The measurements deviate somewhat from the simulation. In particular, the peak magnetic field is 10-15% smaller in magnitude than expected from simulation. This note compares the magnetic field measurements to simulation and to measurements of the existing DARHT magnet, which was designed for the same purpose. We show a comparison of the measured beam deflection performance between Scorpius and DARHT and show updated results of LANL’s electron beam tracking code using the Scorpius measurements. Using field measurements we find that the DARHT magnet will deflect a 19.6 MeV electron beam by 46.6° and the Scorpius magnet will deflect by 44.6°. These results show that the as-built Scorpius magnet has met design requirements for beam deflection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Beam Dynamics Simulations of Transient Beam Loading Effects in the 10GeV EIC Electron Storage Ring

We report on beam dynamics studies of transient beam loading effects in the 10 GeV EIC electron storage ring [1]. The studies are carried out with time-dependent Vlasov-Fokker-Planck simulations performed with the parallel, particle tracking code SPACE [2], which allows to follow self-consistently the dynamics of h bunches, where h in the number of RF buckets, in arbitrary multi-bunch configurations. The specific goal of the numerical simulations is to determine stable RF cavity settings under heavy beam loading. We also study the option to operate with a passive, third-harmonic cavity (3HC) system for bunch lengthening, addressing both stability and the performance limitation due to a gap in the uniform filling pattern for ion clearing.

43 PARTICLE ACCELERATORS↗

Higher Order Multipole Analysis of Beam Crabbing Mode at 197 MHz for Electron-Ion Collider

In the local crabbing scheme for Electron Ion Collider interaction region, beam rotation is initiated by one set of crab cavities at a high beta region with near 90 deg phase advance to the interaction point (IP), and ceased by an identical set of cavities on the symmetrically opposite side. The fundamental mode of all crab cavities is designed to apply a momentum kick to the particle bunches to complete this half cycle of rotation in the horizontal plane, meanwhile the electromagnetic (EM) field can be expanded into a series of harmonics. Despite the lowest dipole mode, higher harmonics cause instabilities during beam operation and loss of particles. We used SimTrack, which is the 6-d symplectic element-by-element tracking code developed at Brookhaven National Laboratory in 2015, to study the threshold of each higher-order harmonic multipolar mode with the criteria of maintaining a dynamic aperture of 6 σ or greater for the Hadron Storage Ring (HSR). In this paper, we showed the threshold for each higher harmonic multipole and their comparison with the multipole strength for the current crab cavity design.

43 PARTICLE ACCELERATORS↗

Uncertainty quantification in MELCOR Safety analysis of ARIES reactor designs

MELCOR-TMAP is a combined thermal-hydraulics and tritium tracking code developed to simulate severe accident scenarios in fission and fusion power plants. Here, we demonstrate the results of MELCOR-TMAP analyses on historical ARIES program reference designs. By coupling MELCOR-TMAP with the open source RAVEN probabilistic risk analysis framework’s Bayesian UQ capabilities, we also demonstrate key uncertainties in material properties with the highest impact on tritium inventory and plant risk.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling of standing wave RF cavities for tracking through multi-pass energy recovery linac

Short bunches, high current and multiple linac pass are all characteristics of Energy Recovery Linacs (ERLs), which may result in collective effects. They in turn, may affect the beam, degrading its quality, or even yield to instabilities causing a beam loss. To study and mitigate these effects one needs a numerical simulation code, that can take into account both the collective effects, as well as, particular ERL features, such as a multi-turn design that does not reach a steady state or the multiple passages of the beam through Radio-Frequency (RF) cavities at different energies. CODAL [1], a code developed by SOLEIL in collaboration with IJCLab, enables such studies. It is a 6 dimensional (6D) tracking code applying 'kicks' based on the integration of the local Hamiltonian for each element of the lattice. It is also capable of simulating space charge, wakefields and coherent synchrotron radiation. However, to correctly take into account the ERL dynamics, an upgrade had to be made to include the effect of a standing wave RF cavity in 6D. In this paper, we will concentrate on the implementation and benchmarking (with DESY’s tracking code ASTRA [2]) of both the longitudinal and the transverse models (by J.B. Rosenzweig and L. Serafini [3]), which we use to carry out tracking of fully analytical 6D RF cavity.

Bogacz, A.↗

Extracting symplectic maps for space-charge dominated beams

Symplecticity of transfer maps is important for reliable evaluation of space-charge dominated beams in accelerators. Unfortunately, most simulation codes that include collective effects, such as space charge, do not use canonical phase-space variables and therefore are not symplectic in the presence of electromagnetic fields. In this paper, we present a numerical method to extract local linear symplectic transfer maps using particle tracking simulation code for space-charge dominated beams. We demonstrate this method for the photoinjector (113 MHz SRF gun) section of the Coherent electron Cooling (CeC) Proof of Principle (POP) experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation of transfer map calculation for electrostatic deflectors in the code COSY INFINITY

The code COSY INFINITY uses a beamline coordinate system with a Frenet–Serret frame relative to the reference particle, and calculates differential algebra-valued transfer maps by integrating the ODEs of motion in the respective vector space over a differential algebra (DA). In our work, we described and performed computation of the DA transfer map of an electrostatic spherical deflector in a laboratory coordinate system using two conventional methods: (1) by integrating the ODEs of motion using a numerical integrator and (2) by computing analytically and in closed form the properties of the respective elliptical orbits from Kepler theory. We compared the resulting transfer maps with (3) the DA transfer map of COSY INFINITY’s built-in electrostatic spherical deflector element ESP and (4) the transfer map of the electrostatic spherical deflector computed using the program GIOS, which uses analytic formulas from a paper1 by Hermann Wollnik regarding second-order aberrations. In addition to the electrostatic spherical deflector, we studied an electrostatic cylindrical deflector, where the Kepler theory is not applicable. We computed the DA transfer map by the ODE integration method (1), and we compared it with the transfer maps by (3) COSY INFINITY’s built-in electrostatic cylindrical deflector element ECL and (4) GIOS. The transfer maps of electrostatic spherical and cylindrical deflectors obtained using the direct calculation methods (1) and (2) are in excellent agreement with those computed using (3) COSY INFINITY. On the other hand, we found a significant discrepancy with (4) the program GIOS.

43 PARTICLE ACCELERATORS↗

Nonuniform Absorption and Scattered Light in Direct-Drive Implosions Driven by Polarization Smoothing

Laser-direct-drive symmetric implosions on OMEGA illuminate a target with 60 laser beams and are designed to produce spherical implosions. Each beam is smoothed using orthogonal polarizations obtained by passing through distributed polarization rotators (DPR’s). Observations of the scattered light from OMEGA implosions do not show the expected symmetry and have much larger variation than standard predictions. For the first time, we have quantified the scattered-light nonuniformity from individual beams and identified the DPR’s as the source of the enhanced nonuniformity. An instrument was invented that isolated and measured the variation in the intensity and polarization of the light scattered from each OMEGA beam. The asymmetric intensity and polarization measurements are explained when the on-target offsets between the two orthogonal polarizations produced by the DPR’s are modeled using a 3-D cross-beam energy transfer (CBET) code that tracks the polarizations of each beam. The time-integrated nonuniformity in laser absorption and scattered light due to CBET and the DPR polarization offsets during high-performance OMEGA implosions is predicted to be significant and dominated by low spherical harmonic mode numbers. Further, the nonuniformity is predicted to be greatly reduced by replacing the DPR’s with new optics that create smaller offsets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗