Beam Dynamics Simulation of the Solenoid Sextupole Error in the LCLS-II Injector
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Several beam simulation codes are used to help gain a better understanding of beam dynamics in the DARHT LIAs. The most notable of these fall into the following categories: for beam production – Tricomp Trak orbit tracking code, LSP Particle in cell (PIC) code, for beam transport and acceleration – XTR static envelope and centroid code, LAMDA time-resolved envelope and centroid code, LSP-Slice PIC code, for coasting-beam transport to target – LAMDA time-resolved envelope code, LSP-Slice PIC code. These codes are also being used to inform the design of Scorpius.
Purpose: To evaluate the Monte Carlo simulated beam data with the measured commissioning data for the Mevion S250 proton therapy system. Method: The Mevion S250 proton therapy system utilizes a passive double scattering technique with a unique gantry mounted superconducting accelerator and offers effective proton therapy in a compact design concept. The field shaping system (FSS) includes first scattering foil, range modulator wheel (RMW), second scattering foil and post absorber and offers two field sizes and a total of 24 treatment options from proton range of 5 cm to 32 cm. The treatment nozzle was modeled in detail using TOPAS (TOolkit for PArticle Simulation) Monte Carlo code. The timing feathers of the moving modulator wheels were also implemented to generate the Spread Out Bragg Peak (SOBP). The simulation results including pristine Bragg Peak, SOBP and dose profiles were compared with the data measured during beam commissioning. Results: The comparison between the measured data and the simulation data show excellent agreement. For pristine proton Bragg Peaks, the simulated proton range (depth of distal 90%) values agreed well with the measured range values within 1 mm accuracy. The differences of the distal falloffs (depth from distal 80% to 20%) were also found to be less than 1 mm between the simulations and measurements. For the SOBP, the widths of modulation (depth of proximal 95% to distal 90%) were also found to agree with the measurement within 1 mm. The flatness of the simulated and measured lateral profiles was found to be 0.6 % and 1.1 %, respectively. Conclusion: The agreement between simulations and measurements demonstrate that TOPAS could be used as a viable platform to proton therapy applications. The matched simulation results offer a great tool and open opportunity for variety of applications.
Independent Modal Space Control (IMSC) method avoids control spillover generated by conventional control schemes such as Coupled Modal Control by decoupling the large flexible space structure into independent subsystems of second order and controlling each mode independently. The IMSC implementation requires that the number of actuators be equal to that of modeled modes, which is in general very huge. Consequently the number of required actuators is unrealizable. Two methods are proposed for the implementation of IMSC with reduced number of actuators. In the first method, the first m modes are optimized, leaving the last (n-m) modes unchanged. In the second method, generalized inverse matrices are employed to design the feedback controller so that the control scheme is suboptimal with respect to IMSC. The performance of the proposed methods is tested by performing computer simulation on a simply support beam. Simulation results are presented and discussed.
This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.
Two computed tomography techniques are explored to reconstruct beam transverse phase space using both simulated beam and multi-wire profile data in the Fermilab Muon Test Area ("MTA") beamline. Both Filtered Back-Projection ("FBP") and Simultaneous Algebraic Reconstruction Technique ("SART") algorithms [2] are considered and compared. Errors and artifacts are compared as a function of each algorithm’s free parameters, and it is shown through simulation and MTA beamline profiles that SART is advantageous for reconstructions with limited profile data.
This article presents the characterization of the in-flight beams, the beam window functions, and the associated uncertainties for the Planck Low Frequency Instrument (LFI). The structure of the paper is similar to that presented in the 2013 Planck release; the main differences concern the beam normalization and the delivery of the window functions to be used for polarization analysis. The in-flight assessment of the LFI main beams relies on measurements performed during observations of Jupiter. By stacking data from seven Jupiter transits, the main beam profiles are measured down to -25 dB at 30 and 44 GHz, and down to -30 dB at 70 GHz. It has been confirmed that the agreement between the simulated beams and the measured beams is better than 1% at each LFI frequency band (within the 20 dB contour from the peak, the rms values are 0.1% at 30 and 70 GHz; 0.2% at 44 GHz). Simulated polarized beams are used for the computation of the effective beam window functions. The error budget for the window functions is estimated from both main beam and sidelobe contributions, and accounts for the radiometer band shapes. The total uncertainties in the effective beam window functions are 0.7% and 1% at 30 and 44 GHz, respectively (at ℓ ≈ 600); and 0.5% at 70 GHz (at ℓ ≈ 1000).
A fully three-dimensional (3D) model of an electron beam focused by a periodic permanent magnet (PPM) stack has been developed. First, the simulation code MAFIA was used to model a PPM stack using the magnetostatic solver. The exact geometry of the magnetic focusing structure was modeled; thus, no approximations were made regarding the off-axis fields. The fields from the static solver were loaded into the 3D particle-in-cell (PIC) solver of MAFIA where fully 3D behavior of the beam was simulated in the magnetic focusing field. The PIC solver computes the time-integration of electromagnetic fields simultaneously with the time integration of the equations of motion of charged particles that move under the influence of those fields. Fields caused by those moving charges are also taken into account; thus, effects like space charge and magnetic forces between particles are fully simulated. The electron beam is simulated by a number of macro-particles. These macro-particles represent a given charge Q amounting to that of several million electrons in order to conserve computational time and memory. Particle motion is unrestricted, so particle trajectories can cross paths and move in three dimensions under the influence of 3D electric and magnetic fields. Correspondingly, there is no limit on the initial current density distribution of the electron beam, nor its density distribution at any time during the simulation. Simulation results including beam current density, percent ripple and percent transmission will be presented, and the effects current, magnetic focusing strength and thermal velocities have on beam behavior will be demonstrated using 3D movies showing the evolution of beam characteristics in time and space. Unlike typical beam optics models, this 3D model allows simulation of asymmetric designs such as non- circularly symmetric electrostatic or magnetic focusing as well as the inclusion of input/output couplers.
Future high current proton accelerators will need to minimize beam loss due to space-charge in order to achieve safe operation while achieving the desired physics goals. One method of space-charge compensation to be tested at the Integrable Optics Test Accelerator (IOTA) at Fermilab is the Electron Column. The concept for this device is to allow a circulating beam to ionize a small region of relatively high pressure residual gas, while using electric and magnetic fields to confine and shape the resulting plasma electrons. If the profile of the electrons is matched to the beam profile transversely and longitudinally, the electrons should counteract the space-charge force of the proton beam. Simulations of the IOTA proton beam circulating through the Electron Column have been performed, with the evolution of the electron plasma and its effect on the beam studied.
Abstract Simulation of the impact of charge-exchange (CX) reactions on beam ions in the Mega Amp Spherical Tokamak (MAST) Upgrade was compared to measurements carried out with a fission chamber (neutron fluxes) and a fast ion deuterium-alpha (FIDA) diagnostic. A simple model was developed to reconstruct the outer-midplane neutral density based on measurements of deuterium-alpha emission from edge neutrals, and on Thomson scattering measurements of electron density and temperature. The main computational tools used were the ASCOT orbit-following code and the FIDASIM code for producing synthetic FIDA signals. The neutral density reconstruction agrees qualitatively with SOLPS-ITER modelling and yields a synthetic passive FIDA signal that is consistent with measurement. When CX losses of beam ions are accounted for, predicted neutron emission rates are quantitatively more consistent with measurement. It was necessary to account for CX losses of beam ions in simulations to reproduce the measured passive FIDA signal quantitatively and qualitatively. The results suggest that the neutral density reconstruction is a good approximation, that CX with edge neutrals causes significant beam-ion losses in MAST Upgrade, typically 20% of beam power, and that the ASCOT fast-ion CX model can be used to accurately predict the redistribution and loss of beam ions due to CX.
Beam-beam interactions pose substantial challenges in the design and operation of circular colliders, significantly affecting their performance. In particular, the weak-strong simulation approach is pivotal for investigating single-particle dynamics during the collider design phase. This paper evaluates the limitations of existing models in weak-strong simulations, noting that while they accurately account for energy changes due to slingshot effects, they fail to incorporate longitudinal coordinate changes ( z variation). To address this gap, we introduce two novel transformations that enhance Hirata’s original framework by including both z variation and slingshot effect-induced energy changes. Through rigorous mathematical analysis and extensive weak-strong simulation studies, we validate the efficacy of these enhancements in achieving a more precise simulation of beam-beam interactions. Our results reveal that although z variation constitutes a higher-order effect and does not substantially affect the emittance growth rate within the specific design parameters of the Electron-Ion Collider, the refined model offers improved accuracy, particularly in scenarios involving the interaction between beam-beam effects and other random diffusion processes, as well as in simulations incorporating realistic lattice models. Published by the American Physical Society 2024
The topic of validating structural optimization methods by use of experimental results is addressed. The need for validating the methods as a way of effecting a greater and an accelerated acceptance of formal optimization methods by practicing engineering designers is described. The range of validation strategies is defined which includes comparison of optimization results with more traditional design approaches, establishing the accuracy of analyses used, and finally experimental validation of the optimization results. Examples of the use of experimental results to validate optimization techniques are described. The examples include experimental validation of the following: optimum design of a trussed beam; combined control-structure design of a cable-supported beam simulating an actively controlled space structure; minimum weight design of a beam with frequency constraints; minimization of the vibration response of helicopter rotor blade; minimum weight design of a turbine blade disk; aeroelastic optimization of an aircraft vertical fin; airfoil shape optimization for drag minimization; optimization of the shape of a hole in a plate for stress minimization; optimization to minimize beam dynamic response; and structural optimization of a low vibration helicopter rotor.
In the energy recovery linac (ERL) based eRHIC design, the electron beam, accelerated in a superconducting ERL, collides with the proton beam circulating in the RHIC ring ones, undergoing a very strong beam-beam interaction with the protons, which warrants careful examination. Following the analysis of the transverse disruption and linear mismatch effects in the electron beam caused by beam-beam collision discussed in [2], we adopt the beam parameters listed in Table I for the eRHIC linac-ring design. After a review of the electron disruption process, we perform long term strong-strong beam-beam simulations with a modification of K. Ohmi’s strong-strong beam-beam code BBSS, with the modification consisting of the regeneration of the electron beam after every collision with the proton beam. With the parameters listed in Table I, a kink instability is observed in the proton beam. By reducing the proton bunch length, thus decreasing the electron disruption parameter, the kink instability is observed to disappear for a proton bunch length σ z,p < 10 cm, or electron disruption parameter d e < 1.5. The effect of random offsets in the electron bunch centroids on the long term stability of the proton beam is studied next, with a proton bunch length σ z,p = 6.3 cm (d e = 0.9), thus below the kink instability threshold. The random offsets of the electron bunch centroids are assumed to be Gaussian distributed, with standard deviation offs, and "delta-correlated" in time (white noise). Studies to take into account finite time-correlation effects and random offsets generated by non-Gaussian random processes will be considered in the future. The effect from fluctuations in the electron bunch population has been studied analytically and will be discussed numerically in a companion paper.
To simulate the distributed spectra of plasmas in space, the differential current density spectrum of the plasma was divided into a number of energy bands and the beam energy and current were calculated for each band to provide a piecewise reproduction of the distributed spectrum. Beam energies and current densities were chosen to match the velocity moments of the plasma distribution function. The velocity moments are averages related to physical quantities such as particle density, flux, pressure, and energy flux, and have been used extensively to characterize the measured properties of plasmas in space. Combinations of one, two, and three beams were found to match two to six velocity moments of Maxwellian distributions. A computational model was used to compare the charging of a spacecraft by plasmas with distributed spectra and by monoenergetic beams.
Purpose: Monte Carlo (MC) simulation is typically regarded as the most accurate dose calculation method for proton therapy. Yet for real clinical cases, the overall accuracy also depends on that of the MC beam model. Commissioning a beam model to faithfully represent a real beam requires finely tuning a set of model parameters, which could be tedious given the large number of pencil beams to commmission. This abstract reports an automatic beam-model commissioning method for pencil-beam scanning proton therapy via an optimization approach. Methods: We modeled a real pencil beam with energy and spatial spread following Gaussian distributions. Mean energy, and energy and spatial spread are model parameters. To commission against a real beam, we first performed MC simulations to calculate dose distributions of a set of ideal (monoenergetic, zero-size) pencil beams. Dose distribution for a real pencil beam is hence linear superposition of doses for those ideal pencil beams with weights in the Gaussian form. We formulated the commissioning task as an optimization problem, such that the calculated central axis depth dose and lateral profiles at several depths match corresponding measurements. An iterative algorithm combining conjugate gradient method and parameter fitting was employed to solve the optimization problem. We validated our method in simulation studies. Results: We calculated dose distributions for three real pencil beams with nominal energies 83, 147 and 199 MeV using realistic beam parameters. These data were regarded as measurements and used for commission. After commissioning, average difference in energy and beam spread between determined values and ground truth were 4.6% and 0.2%. With the commissioned model, we recomputed dose. Mean dose differences from measurements were 0.64%, 0.20% and 0.25%. Conclusion: The developed automatic MC beam-model commissioning method for pencil-beam scanning proton therapy can determine beam model parameters with satisfactory accuracy.
The use of artificial electron beams as probes of the distant magnetosphere is discussed. The accelerators, the method of injecting and analyzing the beams using television techniques and particle counters, and how the beams simulate the natural trapped electrons are summarized. It is found that when an electron beam is injected from a sounding rocket in the ionosphere, the rocket potential rises to a positive value and collects a current of ionospheric plasma electrons. It is proposed to use optical methods on the Space Shuttle to analyze electron beams injected in orbit, and in view of the results of the ECHO 5 experiments, such optical techniques may prove ineffective due to background luminosity from the ground and other sensitivity limitations.