High order time integrators for the simulation of charged particle motion in magnetic quadrupoles

02/22/2018
by   Abele Simona, et al.
0

Magnetic quadrupoles are essential components of particle accelerators like the Large Hadron Collider. In order to study numerically the stability of the particle beam crossing a quadrupole, a large number of particle revolutions in the accelerator must be simulated, thus leading to the necessity to preserve numerically invariants of motion over a long time interval and to a substantial computational cost, mostly related to the repeated evaluation of the magnetic vector potential. In this paper, in order to reduce this cost, we first consider a specific gauge transformation that allows to reduce significantly the number of vector potential evaluations. We then analyze the sensitivity of the numerical solution to the interpolation procedure required to compute magnetic vector potential data from gridded precomputed values at the locations required by high order time integration methods. Finally, we compare several high order integration techniques, in order to assess their accuracy and efficiency for these long term simulations. Explicit high order Lie methods are considered, along with implicit high order symplectic integrators and conventional explicit Runge Kutta methods. Among symplectic methods, high order Lie integrators yield optimal results in terms of cost/accuracy ratios, but non symplectic Runge Kutta methods perform remarkably well even in very long term simulations. Furthermore, the accuracy of the field reconstruction and interpolation techniques are shown to be limiting factors for the accuracy of the particle tracking procedures.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
07/10/2019

Uniformly accurate methods for three dimensional Vlasov equations under strong magnetic field with varying direction

In this paper, we consider the three dimensional Vlasov equation with an...
research
09/16/2019

A Particle Method without Remeshing

We propose a simple tweak to a recently developed regularisation scheme ...
research
12/19/2018

An arbitrary order time-stepping algorithm for tracking particles in inhomogeneous magnetic fields

The Lorentz equations describe the motion of electrically charged partic...
research
04/29/2022

Long term analysis of splitting methods for charged-particle dynamics

In this paper, we rigorously analyze the energy, momentum and magnetic m...
research
10/15/2021

New applications for the Boris Spectral Deferred Correction algorithm for plasma simulations

The paper investigates two new use cases for the Boris Spectral Deferred...
research
03/31/2022

An implicit symplectic solver for high-precision long term integrations of the Solar System

Compared to other symplectic integrators (the Wisdom and Holman map and ...
research
05/15/2020

Performance of the BGSDC integrator for computing fast ion trajectories in nuclear fusion reactors

Modelling neutral beam injection (NBI) in fusion reactors requires compu...

Please sign up or login with your details

Forgot password? Click here to reset