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

DOI10.1016/J.CPC.2019.01.018MaRDI QIDQ6100274zbMATH OpenWikidataFDO

Authors Abele Simona, Luca Bonaventura, Thomas Pugnat, Barbara Dalena

Publication date 12 May 2023

Published in Computer Physics Communications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1802.08157

zbMATH Keywords

Hamilton equationshigh order ODE methodsmagnetic fieldsparticle acceleratorssymplectic methods

Mathematics Subject Classification ID

Numerical analysis (65-XX) Optics, electromagnetic theory (78-XX)

Abstract: 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.

Recommendations

• Explicit high-order symplectic integrators for charged particles in general electromagnetic fields
• A high-order integrator for the Schrödinger equation with time-dependent, homogeneous magnetic field
• Large-stepsize integrators for charged-particle dynamics over multiple time scales
• Arbitrarily high-order energy-preserving methods for simulating the gyrocenter dynamics of charged particles
• High order Gauss-Seidel schemes for charged particle dynamics
• Long-term analysis of a variational integrator for charged-particle dynamics in a strong magnetic field
• High-order integration scheme for relativistic charged particle motion in magnetized plasmas with volume preserving properties
• High-order energy-conserving line integral methods for charged particle dynamics
• Explicit symmetric exponential integrators for charged-particle dynamics in a strong and constant magnetic field
• On a large-stepsize integrator for charged-particle dynamics

Cites work

• Geometric Numerical Integration
• scientific article; zbMATH DE number 3197857 (Why is no real title available?)
• The mimetic finite difference method for the 3D magnetostatic field problems on polyhedral meshes

Cited in

(2)

• IsoGeometric approximations for electromagnetic problems in axisymmetric domains
• BEM-based magnetic field reconstruction by ensemble Kálmán filtering

This page was built for publication: High order time integrators for the simulation of charged particle motion in magnetic quadrupoles

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6100274)