Long-Term Simulations of Beam-Beam Dynamics on GPUs

Future machines such as the electron-ion colliders (JLEIC), linac-ring machines (eRHIC) or LHeC are particularly sensitive to beam-beam effects. This is the limiting factor for long-term stability and high luminosity reach. The complexity of the non-linear dynamics makes it challenging to perform such simulations which require millions of turns. Until recently, most of the methods used linear approximations and/or tracking for a limited number of turns. We have developed a framework which exploits a massively parallel Graphical Processing Units (GPU) architecture to allow for tracking millions of turns in a sympletic way up to an arbitrary order and colliding them at each turn. The code is called GHOST for GPU-accelerated High-Order Symplectic Tracking. As of now, there is no other code in existence that can accurately model the single-particle non-linear dynamics and the beam-beam effect at the same time for a large enough number of turns required to verify the long-term stability of a collider. Our approach relies on a matrix-based arbitrary-order symplectic particle tracking for beam transport and the Bassetti-Erskine approximation for the beam-beam interaction

Accurate and efficient spin integration for particle accelerators

Accurate spin tracking is a valuable tool for understanding spin dynamics in particle accelerators and can help improve the performance of an accelerator. In this paper, we present a detailed discussion of the integrators in the spin tracking code gpuSpinTrack. We have implemented orbital integrators based on drift-kick, bend-kick, and matrix-kick splits. On top of the orbital integrators, we have implemented various integrators for the spin motion. These integrators use quaternions and Romberg quadratures to accelerate both the computation and the convergence of spin rotations. We evaluate their performance and accuracy in quantitative detail for individual elements as well as for the entire RHIC lattice. We exploit the inherently data-parallel nature of spin tracking to accelerate our algorithms on graphics processing units.Comment: 43 pages, 17 figure

Coupling of an SPH-based solver with a multiphysics library

Financiado para publicación en acceso aberto: Universidade de Vigo/CISUGA two-way coupling between the Smoothed Particle Hydrodynamics-based (SPH) code with a multiphysics library to solve complex fluid-solid interaction problems is proposed. This work provides full access to the package for the use of this coupling by releasing the source code, completed with guidelines for its compilation and utilization, and self-contained template setups for practical uses of the novel implemented features, is provided here. The presented coupling expands the applicability of two different solvers allowing to simulate fluids, multibody systems, collisions with frictional contacts using either non-smooth contact (NSC) or smooth contact (SMC) methods, all integrated under the same framework. The fluid solver is the open-source code DualSPHysics, highly optimised for simulating free-surface phenomena and structure interactions, uniquely positioned as a general-purpose Computational Fluid Dynamics (CFD) software with a GPU-accelerated solver. Mechanical systems that comprise collision detection and/or multibody dynamics are solved by the multiphysics library Project Chrono, which uses a Discrete Element Method (DEM). Therefore, this SPH-DEM coupling approach can manage interactions between fluid and complex multibody systems with relative constraints, springs, or mechanical joints.Ministerio de Ciencia e Innovación | Ref. PID2020-113245RB-I00Xunta de Galicia | Ref. ED431C 2021/44Xunta de Galicia | Ref. ED481A-2021/337Ministerio de Ciencia e Innovación, Xunta de Galicia con fondos de la Unión Europea NextGenerationEU y el Fondo Europeo Marítimo y de Pesca | Ref. PRTR-C17.I

DualSPHysics modelling to analyse the response of Tetrapods against solitary wave

Financiado para publicación en acceso aberto: Universidade de Vigo/CISUGThe stability of Tetrapod armour units against solitary waves using Smoothed Particle Hydrodynamics (SPH) method is analysed in this work. To this purpose, the SPH-based DualSPHysics code was coupled with the multiphysics library Project Chrono. Tetrapod units are placed above a submerged mound. DualSPHysics solves the fluid-solid interaction, while Project Chrono solves the Tetrapod-mound interactions based on the contact and material properties of the block surface. The motion of the units during the simulation was compared with the physical model experiments where Tetrapods are made of mortar, and the mound is in PVC. The numerical results expressed as displacements of Tetrapods and damage ratio under different solitary waves are in reasonable agreement with the experiments, proving the capability of the DualSPHysics code to simulate challenging environments under the same numerical framework. The validated tool is then applied to study the stability for different coefficients of friction between mound and Tetrapods aiming at simulating the effects of different materials and surface roughness.Ministerio de Ciencia e Innovación | Ref. PID2020-113245RB-I00Ministerio de Ciencia e Innovación | Ref. TED2021-129479A-I00Xunta de Galicia | Ref. ED431C 2021/44Xunta de Galicia | Ref. ED481A-2021/ 33

Computational modeling of low-density ultracold plasmas

2017 Summer.Includes bibliographical references.In this dissertation I describe a number of different computational investigations which I have undertaken during my time at Colorado State University. Perhaps the most significant of my accomplishments was the development of a general molecular dynamic model that simulates a wide variety of physical phenomena in ultracold plasmas (UCPs). This model formed the basis of most of the numerical investigations discussed in this thesis. The model utilized the massively parallel architecture of GPUs to achieve significant computing speed increases (up to 2 orders of magnitude) above traditional single core computing. This increased computing power allowed for each particle in an actual UCP experimental system to be explicitly modeled in simulations. By using this model, I was able to undertake a number of theoretical investigations into ultracold plasma systems. Chief among these was our lab's investigation of electron center-of-mass damping, in which the molecular dynamics model was an essential tool in interpreting the results of the experiment. Originally, it was assumed that this damping would solely be a function of electron-ion collisions. However, the model was able to identify an additional collisionless damping mechanism that was determined to be significant in the first iteration of our experiment. To mitigate this collisionless damping, the model was used to find a new parameter range where this mechanism was negligible. In this new parameter range, the model was an integral part in verifying the achievement of a record low measured UCP electron temperature of 1.57 ± 0.28K and a record high electron strong coupling parameter, Γ, of 0.35 ± 0.08. Additionally, the model, along with experimental measurements, was used to verify the breakdown of the standard weak coupling approximation for Coulomb collisions. The general molecular dynamics model was also used in other contexts. These included the modeling of both the formation process of ultracold plasmas and the thermalization of the electron component of an ultracold plasma. Our modeling of UCP formation is still in its infancy, and there is still much outstanding work. However, we have already discovered a previously unreported electron heating mechanism that arises from an external electric field being applied during UCP formation. Thermalization modeling showed that the ion density distribution plays a role in the thermalization of electrons in ultracold plasma, a consideration not typically included in plasma modeling. A Gaussian ion density distribution was shown to lead to a slightly faster electron thermalization rate than an equivalent uniform ion density distribution as a result of collisionless effects. Three distinct phases of UCP electron thermalization during formation were identified. Finally, the dissertation will describe additional computational investigations that preceded the general molecular dynamics model. These include simulations of ultracold plasma ion expansion driven by non-neutrality, as well as an investigation into electron evaporation. To test the effects of non-neutrality on ion expansion, a numerical model was developed that used the King model of the electron to describe the electron distribution for an arbitrary charge imbalance. The model found that increased non-neutrality of the plasma led to the rapid expansion of ions on the plasma exterior, which in turn led to a sharp ion cliff-like spatial structure. Additionally, this rapid expansion led to additional cooling of the electron component of the plasma. The evaporation modeling was used to test the underlying assumptions of previously developed analytical expression for charged particle evaporation. The model used Monte Carlo techniques to simulate the collisions and the evaporation process. The model found that neither of the underlying assumption of the charged particle evaporation expressions held true for typical ultracold plasma parameters and provides a route for computations in spite of the breakdown of these two typical assumptions