Accurate modelling of charged particle beams in linear accelerators

A higher order, energy conserving discretization technique for beam dynamics simulations is presented. The method is based on the discontinuous Galerkin (DG) formulation. It utilizes locally refined, non-conforming grids which are designed for high spatial resolution along the path of charged particle beams. Apart from this formulation, the paper introduces a class of general symplectic integrators which conserve discrete energy in a modified sense. Specialized split-operator methods with optimum dispersion properties in the direction of particle motion are, additionally, derived. The application examples given in the paper are performed in a high performance computing environment. They include the self-consistent simulation of the RF electron gun developed by the Photo Injector Test Facility at DESY Zeuthen (PITZ) project and the computation of short range wake fields for ultra-relativistic electron bunches.

[1]  Wojciech Rozmus,et al.  A symplectic integration algorithm for separable Hamiltonian functions , 1990 .

[2]  M. Suzuki,et al.  General theory of higher-order decomposition of exponential operators and symplectic integrators , 1992 .

[3]  L. Fezoui,et al.  Discontinuous Galerkin time‐domain solution of Maxwell's equations on locally‐refined nonconforming Cartesian grids , 2005 .

[4]  Thomas Weiland,et al.  A parallel 3D particle-in-cell code with dynamic load balancing , 2006 .

[5]  Thomas Weiland,et al.  TE/TM scheme for computation of electromagnetic fields in accelerators , 2005 .

[6]  Rolf Schuhmann,et al.  Algebraic Properties and Conservation Laws in the Discrete Electromagnetism , 1999 .

[7]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[8]  Rolf Schuhmann,et al.  Consistent material operators for tetrahedral grids based on geometrical principles , 2004 .

[9]  Rolf Schuhmann,et al.  Long-time numerical computation of electromagnetic fields in the vicinity of a relativistic source , 2003 .

[10]  Rolf Schuhmann,et al.  Consistent material operators for tetrahedral grids based on geometrical principles: Research Articles , 2004 .

[11]  G. Strang.  Accurate partial difference methods I: Linear cauchy problems , 1963 .

[12]  R. Ruth.  A Can0nical Integrati0n Technique , 1983, IEEE Transactions on Nuclear Science.

[13]  Rolf Schuhmann,et al.  The nonorthogonal finite integration technique applied to 2D- and 3D-eigenvalue problems , 2000 .

[14]  C. Birdsall,et al.  Plasma Physics via Computer Simulation , 2018 .

[15]  G. Rodrigue,et al.  High-order symplectic integration methods for finite element solutions to time dependent Maxwell equations , 2004, IEEE Transactions on Antennas and Propagation.

[16]  T. Weiland,et al.  TBCI and URMEL - New Computer Codes for Wake Field and Cavity Mode Calculations , 1983, IEEE Transactions on Nuclear Science.

[17]  Thomas Weiland,et al.  Detailed numerical studies of space charge effects in an FEL RF gun , 2002 .

[18]  Jan S. Hesthaven,et al.  High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids , 2006, J. Comput. Phys..

[19]  Loula Fezoui,et al.  A Nondiffusive Finite Volume Scheme for the Three-Dimensional Maxwell's Equations on Unstructured Meshes , 2002, SIAM J. Numer. Anal..

[20]  Daniel A. White,et al.  A high order mixed vector finite element method for solving the time dependent Maxwell equations on unstructured grids , 2005 .

[21]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[22]  Thomas Weiland,et al.  Time Integration Methods for Particle Beam Simulations with the Finite Integration Theory , 2005 .