Differential algebraic fast multipole-accelerated boundary element method for nonlinear beam dynamics in arbitrary enclosures
Not Available
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.
Not Available
Not Available
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
Each of the two proton beams in the High-Luminosity Large Hadron Collider (HL-LHC) will carry a total energy of 720 MJ. One concern for machine protection is the energy stored in the transverse beam tails, estimated to potentially reach up to 5% of the total stored energy. Several failure scenarios could drive these tails into the collimators, potentially causing damage and therefore severely affecting operational efficiency. Hollow Electron Lenses (HEL) were integrated in the HL-LHC baseline to mitigate this risk by depleting the tails in a controlled way. A hollow-shaped electron beam runs co-axially to the hadron beam over about 3 m, such that halo particles at large amplitudes become unstable, while core particles ideally remain undisturbed. Residual fields from e-beam asymmetries can, however, induce emittance growth of the beam core. Various options for the pulsing of the HEL are considered and are compared using two figures of merit: halo depletion efficiency and core emittance growth. This contribution presents simulations for these two effects with different HEL pulsing modes using the final HL-LHC optics, that was optimized at the location of the lenses.
Not Available
We are commissioning a 2.5-MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MADX simulation codes. Our findings reveal agreement within a factor of 2 in core phase-space density predictions up to 100 synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.
We are commissioning a 2.5~MeV proton beam for the Integrable Optics Test Accelerator at Fermilab, allowing experiments in the strong space-charge regime with incoherent betatron tune shifts nearing 0.5. Accurate modelling of space-charge dynamics is vital for understanding planned experiments. We compare anticipated emittance growth and beam loss in the bare IOTA configuration using transverse space-charge models in Xsuite, PyORBIT, and MAD-X simulation codes. Our findings reveal agreement within 30\% in core phase-space density predictions up to 100~synchrotron periods at moderate beam currents, while tail distributions and beam loss show significant differences.
Nonlinear kickers (NLKs) enable off-axis injection in ultralow-emittance storage rings by providing a strong kick to the injected beam while remaining nearly transparent to the stored beam. In hybrid schemes, a conventional four-kicker bump defines the injected trajectory, and the NLK reduces the first-turn action under constrained beam offset and optics conditions. Effective operation additionally requires stable and reproducible first-turn injection trajectories. We develop a compact action–angle framework that expresses NLK dynamics in terms of Courant–Snyder invariants and yields an analytical bound on achievable action reduction. This formulation provides direct design rules for NLK placement, phase advance, injected-beam offset, and kicker field profile. Within this framework, we identify the 8-wire NLK as a practical baseline and extend its design by relaxing the square-geometry constraint, enabling inward shifting of the off-axis field peak while preserving on-axis field and gradient cancellation. Application to the NSLS-II lattice shows how aperture, pulsed-power, and mechanical constraints combine to determine a coupled design solution. Multi-turn tracking confirms that candidate NLK locations maintain sufficient stay-clear (aperture-clearance) margin, while the optimized wire geometry reduces the required current and Lorentz force load. The results establish a unified approach for NLK-assisted injection design and provide a practical pathway for upgrades in diffraction-limited storage rings.
The Integrable Optics Test Accelerator (IOTA) ring is a new Fermilab facility dedicated to beam physics experiments, currently operating with 150 MeV electrons. Space charge effects are expected to be significant when it operates with 2.5 MeV protons. In this contribution, we present results of a suite of validation tests of PyORBIT, a PICstyle space charge code. Single particle dynamics of quasiintegrable optics using an octupole string in IOTA is compared with MADX, and shown to be in good agreement. Requirements for the convergence of space charge computations are systematically established and when possible, tests involving space charge are compared with theoretical predictions.
A computational procedure suitable for the solution of equations of motions for flexible multibody systems has been developed. The flexible beams are modeled using a fully nonlinear theory which accounts for both finite rotations and large deformations. The present formulation incorporates physical measures of conjugate Cauchy stress and covariant strain increments. As a consequence, the beam model can easily be interfaced with real-time strain measurements and feedback control systems. A distinct feature of the present work is the computational preservation of total energy for undamped systems; this is obtained via an objective strain increment/stress update procedure combined with an energy-conserving time integration algorithm which contains an accurate update of angular orientations. The procedure is demonstrated via several example problems.
A computational procedure suitable for the solution of equations of motions for flexible multibody systems has been developed. A fully nonlinear continuum approach capable of accounting for both finite rotations and large deformations has been used to model a flexible beam component. The beam kinematics are referred directly to an inertial reference frame such that the degrees of freedom embody both the rigid and flexible deformation motions. As such, the beam inertia expression is identical to that of rigid body dynamics. The nonlinear coupling between gross body motion and elastic deformation is contained in the internal force expression. Numerical solution procedures for the integration of spatial kinematic systems can be directily applied to the generalized coordinates of both the rigid and flexible components. An accurate computation of the internal force term which is invariant to rigid motions is incorporated into the general solution procedure.
The Fermilab Accelerator Science and Technology (FAST) facility has a unique configuration of a photocathode rf gun beam injecting two TESLA-type single cavities (CC1 and CC2) in series prior to the cryomodule. Beam propagation off axis in these cavities can result in both long-range and short-range transverse wakefields which can lead to emittance dilution within the macropulses and micropulses, respectively. Two configurations of a Hamamatsu C5680 streak camera viewing a downstream OTR screen were utilized to track centroid shifts during the macropulse (framing mode) for the long-range case and during the micropulse for the short-range case (~10-micron spatial resolution and 2-ps temporal resolution). Steering off axis before CC1, resulted in a 100-kHz bunch centroid oscillation within the macropulse that was detected by the downstream rf BPMs and the streak camera*. At 500 pC/b, 50b, and 4-mrad off-axis vertical steering into CC2, we observed an ~ 100-micron head -tail centroid shift in the streak camera image y(t) profiles which we attributed to a short-range wakefield effect. Additional results for kick-angle compensations and model results will be presented.
Parameters of the post-disruption runaway electron (RE) beam in the collisionless background plasma achieved after deuterium injection are investigated in DIII-D. The spatially resolved RE energy distribution function is measured for the first time during the RE plateau stage by inverting hard X-ray bremsstrahlung spectra. It has maximum energy up to 20 MeV and a non-monotonous feature at 5{6 MeV observed only in the core of the beam supporting the possibility of kinetic instabilities. The RE current profile is reconstructed for the first time using the spatially resolved RE energy distribution. It is found to be more peaked than the pre-disruption plasma current, with higher internal inductance, suggesting preferential formation of REs in the core plasma or potentially a radially inward motion of the REs. The accessed relatively low-current (180 kA) RE beam is found stable, likely due to its elevated q profile. From this base stable equilibrium, instability is accessed by ramping the solenoid current. Under a large applied accelerating loop voltage an internal kink mode is observed. The kink mode leads to sawtooth-like relaxation of the RE current profile, but drives no RE loss. Under conditions of a large decelerating voltage, chirping Alfvénic instabilities driven by REs during the RE plateau stage are observed for the first time. The low-frequency (0.1{7 MHz) Alfvénic instabilities, supposedly compressional Alfvén eigenmodes, are found to be correlated with a RE loss signal. Finally, this provides novel plasma conditions to study these instabilities and investigate opportunities to utilize them for RE control. Bifurcated electron cyclotron emission (ECE) spectrum is measured during the collisionless RE plateau, with a break point at ≈100 GHz suggesting resonant absorption of the ECE at low frequencies.
The Electron Ion Collider aims to produce luminosities of 10³⁴ cm⁻²s⁻¹ . The machine will operate over a broad range of collision energies with highly polarized beams. The coexistence of highly radiative electrons and nonradiative ions produce a host of unique effects. Strong hadron cooling will be employed for the final factor of 3 luminosity boost.
Several linac driven X-ray Free Electron Lasers (XFELs) are being developed to provide high brightness photon beams with very short, tunable wavelengths. In this paper, three XFEL configurations are proposed that achieve LCLS-like performance using X-band linac drivers. These linacs are more versatile, efficient and compact than ones using S-band or C-band rf technology. For each of the designs, the overall accelerator layout and the shaping of the bunch longitudinal phase space are described briefly.
Abstract not provided.
The outline for this presentations is as follows: I. Introduction. II. Injector beamline and streak camera viewing optical transition radiation (OTR) screen at X121. -Strategy of beam steering off axis into TESLA Cavities to generate wakefields and beam effects. III. Previous long long-range wakefield test, Higher Higher-order Modes (HOMs). IV. Initial observations of short short-range wakefield effects. V. Summary.
The Strong Hadron Cooler Energy Recovery Linac (SHC-ERL) for the Electron-Ion Collider (EIC) is required to produce an electron beam with a bunch charge of 1 nC, an average current of 100 mA, and a beam energy of both 150 and 55 MeV, with strict requirements for the transverse emittance, slice energy spread, and longitudinal distribution. This paper covers the design in detail, including preliminary considerations of higher order effects and diagnostics.