Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accelerator Simulation”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

RF Gap Transformation in PARMILA

The PARMILA (Phase And Radial Motion in Ion Linear Accelerators) code is used to design ion linacs and simulate the beam dynamics therein. The original code was developed in the 1960’s and later the Los Alamos Accelerator Code Group developed a PC version that has been used worldwide. Recently, in an effort to modernize the code and use it as the basis for a high-performance beam dynamics simulator with the LANSCE linac, the details of the code have come under additional scrutiny. PARMILA is a “z code” that uses various transformations to approximately represent the effect of the beam line element on the beam. In the case of the RF gap, it treats the particle dynamics inside the gap as a combination of drift + kick at the mid-plane + drift. Recently, we uncovered an undocumented keyword “off-exact” in PARMILA that enables the code to use the full modified-Bessel functions instead of low-order approximations in the RF gap transformation. This change improved the accuracy of the transformation for off-axis particles. However, a discrepancy between the analytical formula and the code was found for the transformation of the transverse coordinates. The purpose of this tech note is to document the derivation of the gap transformation and the related DTL cell length calculation used in the PARMILA code.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Reconstruction of implanted marker trajectories from cone-beam CT projection images using interdimensional correlation modeling

Purpose: Cone-beam CT (CBCT) is a widely used imaging modality for image-guided radiotherapy. Most vendors provide CBCT systems that are mounted on a linac gantry. Thus, CBCT can be used to estimate the actual 3-dimensional (3D) position of moving respiratory targets in the thoracic/abdominal region using 2D projection images. The authors have developed a method for estimating the 3D trajectory of respiratory-induced target motion from CBCT projection images using interdimensional correlation modeling. Methods: Because the superior–inferior (SI) motion of a target can be easily analyzed on projection images of a gantry-mounted CBCT system, the authors investigated the interdimensional correlation of the SI motion with left–right and anterior–posterior (AP) movements while the gantry is rotating. A simple linear model and a state-augmented model were implemented and applied to the interdimensional correlation analysis, and their performance was compared. The parameters of the interdimensional correlation models were determined by least-square estimation of the 2D error between the actual and estimated projected target position. The method was validated using 160 3D tumor trajectories from 46 thoracic/abdominal cancer patients obtained during CyberKnife treatment. The authors’ simulations assumed two application scenarios: (1) retrospective estimation for the purpose of moving tumor setup used just after volumetric matching with CBCT; and (2) on-the-fly estimation for the purpose of real-time target position estimation during gating or tracking delivery, either for full-rotation volumetric-modulated arc therapy (VMAT) in 60 s or a stationary six-field intensity-modulated radiation therapy (IMRT) with a beam delivery time of 20 s. Results: For the retrospective CBCT simulations, the mean 3D root-mean-square error (RMSE) for all 4893 trajectory segments was 0.41 mm (simple linear model) and 0.35 mm (state-augmented model). In the on-the-fly simulations, prior projections over more than 60° appear to be necessary for reliable estimations. The mean 3D RMSE during beam delivery after the simple linear model had established with a prior 90° projection data was 0.42 mm for VMAT and 0.45 mm for IMRT. Conclusions: The proposed method does not require any internal/external correlation or statistical modeling to estimate the target trajectory and can be used for both retrospective image-guided radiotherapy with CBCT projection images and real-time target position monitoring for respiratory gating or tracking.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Nonthermally Dominated Electron Acceleration during Magnetic Reconnection in a Low-beta Plasma

This work was motivated by electron acceleration during solar flares. After some introductory remarks on proposed particle acceleration mechanisms and questions needing answers, dynamic simulations and simulation results are presented including energy spectra and the formation of the power law distribution. In summary, magnetic reconnection is highly efficient at converting the free magnetic energy stored in a magnetic shear and accelerating electrons to nonthermal energies in low-β regime. The nonthermal electrons have a dominant fraction and form power-law energy spectra with spectral index p ~ 1 in low-β regime. Electrons are preferentially accelerated along the curvature drift direction along the electric field induced by the reconnection outflow. The results can be applied to explain the observations of electron acceleration during solar flares.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. I. A Systematic Study

Blazars are relativistic magnetized plasma outflows from supermassive black holes that point very close to our line of sight. Their emission is nonthermal-dominated and highly variable across the entire electromagnetic spectrum. Relativistic magnetic reconnection has been proposed as the driver of particle acceleration during blazar flares. While recent particle-in-cell (PIC) simulations have self-consistently studied the evolution of magnetic reconnection and particle acceleration therein, the resulting radiation signatures have not been systematically explored. In particular, the polarization signatures, which directly reflect the characteristic strongly dynamical evolution of magnetic field during reconnection, have not been carefully investigated. In this paper, we present a systematic study of radiation and polarization signatures arising from magnetic reconnection in blazars, based on combined PIC and polarized radiation transfer simulations with various physical parameters. We identify a harder-when-brighter trend in the spectral evolution. Moreover, higher-frequency bands (ultraviolet to X-ray) tend to flare earlier than lower-frequency bands (infrared to optical) in the synchrotron spectral component. Most importantly, polarization signatures appear more variable with higher frequencies. We find that the variation in temporal polarization depends strongly on the guide field strength. Specifically, reconnection with a significant guide field component leads to a very high polarization degree that contradicts typical blazar observations, while large polarization angle rotations are unique signatures of magnetic reconnection between nearly antiparallel magnetic field lines. These rotations are at least 90° and can extend to >180°, and they may be in either direction. Furthermore, these results imply that blazars that have shown large polarization angle rotations intrinsically have more nearly antiparallel magnetic field morphology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Simulation of radiation damping in rings, using stepwise ray-tracing methods

The ray-tracing code Zgoubi computes particle trajectories in arbitrary magnetic and/or electric field maps or analytical field models. It includes a built-in fitting procedure, spin tracking many Monte Carlo processes. The accuracy of the integration method makes it an efficient tool for multi-turn tracking in periodic machines. Energy loss by synchrotron radiation, based on Monte Carlo techniques, had been introduced in Zgoubi in the early 2000s for studies regarding the linear collider beam delivery system. However, only recently has this Monte Carlo tool been used for systematic beam dynamics and spin diffusion studies in rings, including eRHIC electron-ion collider project at the Brookhaven National Laboratory. Some beam dynamics aspects of this recent use of Zgoubi capabilities, including considerations of accuracy as well as further benchmarking in the presence of synchrotron radiation in rings, are reported here.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

TU-H-CAMPUS-TeP1-02: Seated Treatment: Setup Uncertainty Comparable to Supine

Purpose: For some head and neck patients, positioning in the supine position is not well tolerated. For these patients, treatment in a seated position would be preferred. We have evaluated inter- and intra- fraction uncertainty of patient set-up in a novel treatment chair which is compatible with modern linac designs. Methods: Five head-and-neck cancer patients were positioned in the chair, fitted with immobilization devices, and imaged with orthogonal X-rays. The couch (with chair attached) was rotated to simulate delivery (without actual treatment), another set of images were acquired, providing a measure of intra-fraction displacement. The patient then got off of and back onto the chair and the process was repeated, thus providing a measure of inter-fraction set-up uncertainty. Six sub-regions in the head-and-neck were rigidly registered to evaluate local intra- and interfraction displacement. Image guidance was simulated by first registering one sub-region; the residual displacement of other sub-regions was then measured. Additionally, a patient questionnaire was administered to evaluate tolerance of the seated position. Results: The chair design is such that all advantages of couch motions may be utilized. Average inter- and intrafraction displacements of all sub-regions in the seated position were less than 2 and 3 mm, respectively. When image guidance was simulated, interfraction displacements were reduced by an average of 4 mm, providing comparable setup to the supine position. The enrolled patients, who had no indication for a seated treatment position, reported no preference for the seated or the supine position. Conclusion: The novel chair design provides acceptable inter- and intra-fraction displacement, with reproducibility similar to that observed for patients in the supine position. Such a chair will be utilized for patients who cannot tolerate the supine position and use with CBCT images for planning, in a fixed-beam linac system, and for other treatment sites is under investigation. Funding: Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

The Hamiltonian Mechanics of Stochastic Acceleration

We show how to end the physical Langevin equation describing the trajectories of particles undergoing collisionless stochastic acceleration. These stochastic dierential equations retain not only one-, but two-particle statistics, and inherit the Hamiltonian nature of the underlying microscopic equations. This opens the door to using stochastic variational integrators to perform simulations of stochastic interactions such as Fermi acceleration. We illustrate the theory by applying it to two example problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Reconstruction of beam parameters and betatron radiation spectra measured with a Compton spectrometer

The photon flux resulting from high-energy electron beam interactions with high-field systems, such as those found in the upcoming FACET-II experiments at the SLAC National Accelerator Laboratory, yields deep insight into the electron beam’s underlying dynamics during the interaction. However, extracting this information is an intricate process. To demonstrate how to approach this challenge using modern methods, this paper utilizes simulated data that models plasma wakefield acceleration-derived betatron radiation in experiments to determine reliable methods of reconstructing key beam and beam-plasma interaction properties. For betatron radiation measurements, translating the observed 200⁢ keV to 30⁢ MeV photon double-differential energy-angle spectra obtained from an advanced Compton spectrometer requires testing multiple methods to optimize the pipeline from its response to incident electron beam information. The paper compares maximum likelihood estimation and machine learning to refine the translation of photon spectra into precise electron beam metrics, such as spot size, energy, and emittance, enhancing the understanding of beam behavior within these dense, high-field environments. We also introduce machine learning and the expected maximization algorithm to reconstruct the primary photon spectrum, employing a multilayer neural network for regression analysis of the energy and angle spectra. With appropriate modifications, the advanced methods reproduce relevant incident beam parameters with high accuracy, even for beam sizes in the <10 μ⁢m range. This capacity is critical to understanding intense beam propagation and its optimization in plasma.

Beam code development & simulation techniques↗
DOE OSTI2013

Increasing energy coupling into plasma waves by tailoring the laser radial focal spot distribution in a laser wakefield accelerator

By controlling the focal spot quality with a deformable mirror, we are able to show that increasing the fraction of pulse energy contained within the central part of the focal spot, while keeping the total energy and central spot size constant, significantly increases the amount of energy transferred to the wakefield: Our measurements show that the laser loses significantly more laser energy and undergoes greater redshifting and that more charge is produced in the accelerated beam. Three dimensional particle in cell simulations performed with accurate representations of the measured focal spot intensity distribution confirm that energy in the wings of the focal spot is effectively wasted. Even though self-focusing occurs, energy in the wings of the focal spot distribution is not coupled into the wakefield, emphasising the vital importance of high quality focal spot profiles in experiments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Parallel heat flux and flow acceleration in open field line plasmas with magnetic trapping

The magnetic field strength modulation in a tokamak scrape-off layer (SOL) provides both flux expansion next to the divertor plates and magnetic trapping in a large portion of the SOL. Previously, we have focused on a flux expander with long mean-free-path, motivated by the high temperature and low density edge anticipated for an absorbing boundary enabled by liquid lithium surfaces. Here, the effects of magnetic trapping and a marginal collisionality on parallel heat flux and parallel flow acceleration are examined. The various transport mechanisms are captured by kinetic simulations in a simple but representative mirror-expander geometry. The observed parallel flow acceleration is interpreted and elucidated with a modified Chew-Goldberger-Low model that retains temperature anisotropy and finite collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

MO-F-CAMPUS-T-03: Data Driven Approaches for Determination of Treatment Table Tolerance Values for Record and Verification Systems

Purpose: To determine appropriate couch tolerance values for modern radiotherapy linac R&V systems with indexed patient setup. Methods: Treatment table tolerance values have been the most difficult to lower, due to many factors including variations in patient positioning and differences in table tops between machines. We recently installed nine linacs with similar tables and started indexing every patient in our clinic. In this study we queried our R&V database and analyzed the deviation of couch position values from the acquired values at verification simulation for all patients treated with indexed positioning. Mean and standard deviations of daily setup deviations were computed in the longitudinal, lateral and vertical direction for 343 patient plans. The mean, median and standard error of the standard deviations across the whole patient population and for some disease sites were computed to determine tolerance values. Results: The plot of our couch deviation values showed a gaussian distribution, with some small deviations, corresponding to setup uncertainties on non-imaging days, and SRS/SRT/SBRT patients, as well as some large deviations which were spot checked and found to be corresponding to indexing errors that were overriden. Setting our tolerance values based on the median + 1 standard error resulted in tolerance values of 1cm lateral and longitudinal, and 0.5 cm vertical for all non- SRS/SRT/SBRT cases. Re-analizing the data, we found that about 92% of the treated fractions would be within these tolerance values (ignoring the mis-indexed patients). We also analyzed data for disease site based subpopulations and found no difference in the tolerance values that needed to be used. Conclusion: With the use of automation, auto-setup and other workflow efficiency tools being introduced into radiotherapy workflow, it is very essential to set table tolerances that allow safe treatments, but flag setup errors that need to be reassessed before treatments.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Novel Materials R&D for Next-Generation Accelerator Target Facilities

High-Entropy Alloys and Electrospun Nanofiber materials are two novel classes of materials that can offer improved resistance to beam-induced radiation damage and thermal shock. Research to develop these new materials specifically for multi-megawatt accelerator target applications, such as beam windows and particle-production targets, has recently begun. The research program will combine in-beam experiments with complementary simulations to tailor the microstructures of these novel materials for use in next-generation accelerator target facilities. Iterative simulations to optimize the material composition, physics performance and beam-induced thermomechanical response will guide the material design and fabrication processes based on established figures of merit. Ensuing material irradiation experiments using low-energy ions and prototypic high-energy protons, followed by extensive post-irradiation material characterization, will then assess and qualify the selected novel materials. This talk will provide an overview of the novel materials development research program initiated at Fermilab through my DOE Early Career Research Program award.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Comparing Measured Driver Behavior Distributions to Results from Car-Following Models using SUMO and Real-World Vehicle Trajectories from Radar: SUMO Default vs. Radar-Measured CF model Parameters

In this study, the physical principles governing car-following (CF) behavior and their impact on traffic flow at signalized intersections are investigated. High temporal-resolution radar data is used to provide valuable insights into actual CF behavior, including acceleration, deceleration, and time headway distribution. Demand-calibrated SUMO simulations are run using empirical CF parameter distributions, and three CF models are evaluated: IDM, EIDM, and Krauss. By emulating radar data in SUMO and processing simulated vehicle traces, discrepancies between empirical and simulated parameter distributions are identified. Further analysis includes comparisons with default SUMO CF model parameters. The findings reveal that measured accelerations differ from CF model parameter accelerations and using the empirical value ($\mu = 0.89m/s^2$) leads to unrealistic simulations that fail volume-based calibration. Default parameters for all three models reasonably approximate the mean and median of measured parameters, but fail to capture the true distribution shape, partly due to homogeneity when using default parameters. The results show that it is more effective to simulate with the default parameters provided by SUMO rather than using measurements of real-world distributions without additional calibration. Future work will investigate closing the loop between the measured real-world and SUMO distributions using traditional calibration tactics, as well as assess the impact of calibrated vs. default CF parameters on simulation outputs like fuel consumption.

Schrader, Max (ORCID:0000000336737672)↗
DOE OSTI2022

Direct Measurement of Ice-Ablator Interface Motion for Instability Mitigation in Indirect Drive ICF Implosions

In indirect drive inertial confinement fusion (ICF) implosions hydrodynamic instability growth at the imploding capsule ablator-DT fuel interface can reduce fuel compressibility and inject ablator into the hot spot hence reducing hot spot pressure and temperature. As a mitigation strategy, a gentle acceleration of this interface is predicted by simulations and theory to significantly reduce this instability growth in the early stage of the implosion. We have performed high-contrast, time-resolved x-ray refraction enhanced radiography (RER) to accurately measure the level of acceleration as a function of the initial laser drive time history for indirect-drive implosions on the National Ignition Facility. Here, we demonstrate a transition from no acceleration to 20 ± 1.8μm ns -2 acceleration by tweaking the drive that should reduce the initial instabilities by an order of magnitude at high modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

The ray-tracing code Zgoubi - latest developments

The reader is referred to the ICFA-BD-NL No.43, and to a more recent status article in NIMA, regarding the technical details of the Zgoubi tracking engine and its accuracy, which make it "near-symplectic" - sufficiently so to model millions of turns in large hadron colliders, - and its potential for pushing particles in complex accelerator structures and accelerator optical elements or their field maps. Extended documentation can be found in the regular FFA workshops as well as in the SourceForge code development site where in particular many examples of beam lines, FFAs and other ring simulations are maintained. There is also a comprehensive user's guide. Further documentation, including an attempt to demonstrate 6D simulation in all possible types of particle accelerators, can be found in the "Computer Laboratory" Stony Brook University "PHY689" semester series.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Comparison between simulated and observed LHC beam backgrounds in the ATLAS experiment at E beam =4 TeV

Results of dedicated Monte Carlo simulations of beam-induced background (BIB) in the ATLAS experiment at the Large Hadron Collider (LHC) are presented and compared with data recorded in 2012. During normal physics operation this background arises mainly from scattering of the 4 TeV protons on residual gas in the beam pipe. Methods of reconstructing the BIB signals in the ATLAS detector, developed and implemented in the simulation chain based on the \textsc{Fluka} Monte Carlo simulation package, are described. The interaction rates are determined from the residual gas pressure distribution in the LHC ring in order to set an absolute scale on the predicted rates of BIB so that they can be compared quantitatively with data. Through these comparisons the origins of the BIB leading to different observables in the ATLAS detectors are analysed. The level of agreement between simulation results and BIB measurements by ATLAS in 2012 demonstrates that a good understanding of the origin of BIB has been reached.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Space-Charge Simulation of Integrable Rapid Cycling Synchrotron

Integrable optics is an innovation in particle accelerator design that enables strong nonlinear focusing without generating parametric resonances. We use a Synergia space-charge simulation to investigate the application of integrable optics to a high-intensity hadron ring that could replace the Fermilab Booster. We find that incorporating integrability into the design suppresses the beam halo generated by a mismatched KV beam. Our integrable rapid cycling synchrotron (iRCS) design includes other features of modern ring design such as low momentum compaction factor and harmonically canceling sextupoles. Experimental tests of high-intensity beams in integrable lattices will take place over the next several years at the Fermilab Integrable Optics Test Accelerator (IOTA) and the University of Maryland Electron Ring (UMER).

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Beam-induced backgrounds measured in the ATLAS detector during local gas injection into the LHC beam vacuum

Inelastic beam-gas collisions at the Large Hadron Collider (LHC), within a few hundred metres of the ATLAS experiment, are known to give the dominant contribution to beam backgrounds. These are monitored by ATLAS with a dedicated Beam Conditions Monitor (BCM) and with the rate of fake jets in the calorimeters. These two methods are complementary since the BCM probes backgrounds just around the beam pipe while fake jets are observed at radii of up to several metres. In order to quantify the correlation between the residual gas density in the LHC beam vacuum and the experimental backgrounds recorded by ATLAS, several dedicated tests were performed during LHC Run 2. Local pressure bumps, with a gas density several orders of magnitude higher than during normal operation, were introduced at different locations. The changes of beam-related backgrounds, seen in ATLAS, are correlated with the local pressure variation. In addition the rates of beam-gas events are estimated from the pressure measurements and pressure bump profiles obtained from calculations. Using these rates, the efficiency of the ATLAS beam background monitors to detect beam-gas events is derived as a function of distance from the interaction point. These efficiencies and characteristic distributions of fake jets from the beam backgrounds are found to be in good agreement with results of beam-gas simulations performed with theFluka Monte Carlo programme.

43 PARTICLE ACCELERATORS↗