C-2 Electromagnetic Field Computations
Paper Title Page
SAPAG04 HOM-Mitigation for Future SPS 33-Cell 200 MHz Accelerating Structures 35
  • P. Kramer, C. Vollinger
    CERN, Geneva, Switzerland
  The CERN SPS 200 MHz travelling wave (TW) accelerating structures pose an intensity limitation for the planned High Luminosity (HL-) LHC upgrade. Higher-order modes (HOMs) around 630 MHz have been identified as one of the main sources of longitudinal multi-bunch instabilities. Improved mitigation of these HOMs with respect to today’s HOM-damping scheme is therefore an essential part of the LHC injectors upgrade (LIU) project. The basic principles of HOM-couplers in cavities and today’s damping scheme are reviewed, before illustrating the numerous requirements an improved damping scheme for the future 33-cell structures must fulfil. These are, amongst others, the mitigation of HOMs situated in the lower part of the structure where there are no access ports for extraction, a sufficient overall damping performance and an acceptable influence on the fundamental accelerating passband (FPB). Different approaches tackling these challenges are investigated and their performance, advantages and pitfalls are evaluated by ACE3P and CST electromagnetic (EM) field solver suites.  
slides icon Slides SAPAG04 [2.184 MB]  
DOI • reference for this paper ※ https://doi.org/10.18429/JACoW-ICAP2018-SAPAG04  
About • paper received ※ 19 October 2018       paper accepted ※ 24 October 2018       issue date ※ 26 January 2019  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
Fast Multipole Methods for Multiparticle Simulations  
  • H. Zhang
    JLab, Newport News, Virginia, USA
  Funding: Authored by Jefferson Science Associates, LLC under U.S. DOE Contract No. DE-AC05-06OR23177. The U.S. Government retains a license to publish or reproduce this manuscript for U.S. Government purposes.
The fast multipole method (FMM) reduces the computation cost of the pairwise non-oscillating interaction between N particles from O(N2) to O(N). In the FMM, the contribution from a source particle is represented as a multipole expansion, while the contributions from multiple faraway sources can be combined into a local expansion around an objective particle. Without the dependence on a grid covering the whole domain under study, the FMM treats any charge distribution and geometry in a natural way. It avoids artificial smoothing due to the grid size and redundant computation on the free space grids. We will introduce the concept of the FMM using the Coulomb interaction as an example and then explain how the FMM can be extended to arbitrary non-oscillating interactions. Examples and discussions on how the FMM can be used in scientific simulations, especially in accelerator physics will also be provided.
slides icon Slides SUPAG02 [1.705 MB]  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
Recent Developments in Wake Field and Beam Dynamics Computation  
  • E. Gjonaj
    TEMF, TU Darmstadt, Darmstadt, Germany
  Funding: Work partially funded by DESY, Hamburg.
Wake potentials and beam coupling impedances can be calculated analytically only for simple structures and for special limiting cases. For the calculation of wake fields in ’real-world’ 3D accelerator structures, one has to rely on numerical electromagnetic field computation. Among the most successful numerical techniques for wake field calculations in the time domain are dispersion-free methods in the moving window. These techniques are particularly useful for short-range wake field calculations. Recently, this class of methods has been extended to include Surface Impedance Boundary Conditions (SIBC) based on the Auxiliary Differential Equation (ADE) technique. These boundary conditions allow the computation of resistive wall wake fields for 3D structures with arbitrary frequency dependent conductivity. An important application of this method is the calculation resistive wall wake fields in novel accelerator chambers with NEG and amorphous carbon coatings. Other developments to be discussed include the calculation of CSR-wakes in bunch compressors and undulator structures for x-ray sources. This task is computationally very difficult because of the curved bunch trajectory that leads to extremely high frequency and long-range wake fields. Time domain as well as frequency domain methods based on high order DG and FE discretization techniques for the electromagnetic fields computation in such structures will be presented.
slides icon Slides MOPLG02 [2.762 MB]  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
Exploring the Validity of the Paraxial Approximation for Coherent Synchrotron Radiation Wake Fields  
  • D. A. Bizzozero, H. De Gersem, E. Gjonaj
    TEMF, TU Darmstadt, Darmstadt, Germany
  Coherent synchrotron radiation (CSR) is an essential consideration in modern accelerators, yet is often computationally difficult to accurately model. A common approach used in simulating CSR effects uses the paraxial, or slowly-varying envelope approximation with a simple constant cross-section approximation of the geometry. While these approximations are often valid for the simulation of many accelerator components, we aim to more closely analyze the errors introduced by such approximations by comparing them with wake field solutions obtained by full-wave electromagnetic field simulations. The simulations are performed with CSRDG (Coherent Synchrotron Radiation with Discontinuous Galerkin), our GPU-enabled MATLAB code. Extended from earlier work [Coherent Synchrotron Radiation and Wake Fields With Discontinuous Galerkin Time Domain Methods, Proceedings of IPAC 2017, Copenhagen, Denmark], CSRDG evolves Maxwell’s equations the time domain after a curvilinear coordinate transformation and a Fourier series decomposition in a transverse direction.  
slides icon Slides MOPAG03 [2.149 MB]  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
Advances in Accelerator Modeling with Parallel Multi-Physics Code Suite ACE3P  
  • L. Xiao, L. Ge, Z. Li, C.-K. Ng
    SLAC, Menlo Park, California, USA
  ACE3P is a comprehensive set of parallel finite-element codes for multi-physics modeling of accelerator structures including integrated electromagnetic, thermal and mechanical effects. Recent advances of ACE3P have been focused on the development of multi-physics modeling capabilities, implementation of advanced numerical algorithms, and improvement of code performance on state-of-the-art high-performance computing (HPC) platforms for large-scale accelerator applications. A nonlinear eigensolver using the CORK algorithm [1] has been implemented in the eigensolver module Omega3P to enable accurate determination of damping factors of resonant modes above the beampipe cutoff frequency. It has enabled the first-ever direct calculation of trapped modes in the TESLA TTF cryomodules, providing reliable damping factors that were validated against measurements. A newly developed mechanical eigensolver in the multi-physics module TEM3P has allowed the determination of mechanical modes in Fermilab PIP-II high beta 650 MHz cryomodule, demonstrating mode coupling between the 6 cavities in the cryomodule. To exploit multi-core computer architectures on supercomputers, a hybrid MPI+OpenMP parallel programing has been developed in the particle tracking module Track3P to speed up dark current simulation in multiple cavities for the LCLS-II linac. Highlights of these developments and their impacts on accelerator modeling using HPC will be presented.
[1] R. Van Beeuman, Invited talk, this conference.
slides icon Slides TUPAF05 [1.355 MB]  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
TUPAG01 Computation of Eigenmodes in the BESSY VSR Cavity Chain by Means of Concatenation Strategies 253
  • T. Flisgen, A.V. Vélez
    HZB, Berlin, Germany
  • J. Heller, G. Zadeh, U. van Rienen
    Rostock University, Faculty of Computer Science and Electrical Engineering, Rostock, Germany
  Funding: The research leading to these results was supported by the German Bundesministerium für Bildungund Forschung, Land Berlin and grants of Helmholtz Association
Invited Talk: The computation of eigenmodes in chains of superconducting cavities with asymmetric couplers is a demanding problem. This problem typically requires the use of high-performance computers in combination with dedicated software packages. Alternatively, the eigenmodes of chains of superconducting cavities can be determined by the so-called State-Space Concatenation (SSC) approach that has been developed at the University of Rostock. SSC is based on the decomposition of the full chain into individual segments. Subsequently, the RF properties of every segment are described by reduced-order models. These reduced-order models are concatenated to a reduced-order model of the entire chain by means of algebraic side constraints arising from continuity conditions of the fields across the decomposition planes. The constructed reduced-order model describes the RF properties of the complete structure so that the field distributions, the coupling impedances and the external quality factors of the eigenmodes of the full cavity chain are available. In contrast to direct methods, SSC allows for the computation of the eigenmodes of cavity chains using desktop computers. The current contribution revises the scheme using the BESSY VSR cavity chain as an example. In addition, a comparison between a direct computation of a specific localized mode is described.
slides icon Slides TUPAG01 [3.483 MB]  
DOI • reference for this paper ※ https://doi.org/10.18429/JACoW-ICAP2018-TUPAG01  
About • paper received ※ 21 October 2018       paper accepted ※ 28 January 2019       issue date ※ 26 January 2019  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
High-Precision Lossy Eigenfield Analysis Based on the Finite Element Method  
  • W. Ackermann, H. De Gersem, V. Pham-Xuan
    TEMF, TU Darmstadt, Darmstadt, Germany
  A proper eigenanalysis of resonating particle accelerator components is particularly advantageous to characterize structures with high quality factors. While in former times eigenmode calculations have been concentrating on the lossless cases only, meanwhile also lossy structures with finite-conductive materials or with absorbing boundary conditions like PML or ports even with low quality factors are routinely available. In the lossless case where no damping is present, all eigenvalues are located along the real axis. If damping has to be modeled instead, the corresponding eigenvalues are distributed within the first quadrant of the complex plane that renders their determination much more expensive. One of the critical issues is that no resonance should be missed so that all desired eigenvalues in a given region of the complex plane can be precisely determined. We implemented two different eigenvalue solvers based on a distributed-memory architecture. While the first one is a classical Jacobi-Davidson eigenvalue solver which has been adopted to be used also within a complex-arithmetic environment, the second one is based on the contour-integral method which enables to determine all eigenvalues within a given closed contour in the complex plane. Both solvers are attached to a FEM processor with second-order edge elements on curved tetrahedra and can be used together in order to improve the computational efficiency. In the presentation a selection of successful real-world applications of the implemented parallel eigenvalue solvers will be given.  
slides icon Slides TUPAG03 [15.980 MB]  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
TUPAG04 Statistical Analysis of the Eigenmode Spectrum in the SRF Cavities with Mechanical Imperfections 265
  • A. Lunin, T.N. Khabiboulline, N. Solyak, A.I. Sukhanov, V.P. Yakovlev
    Fermilab, Batavia, Illinois, USA
  Funding: Work is supported by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the United States Department of Energy
The superconducting radio frequency (SRF) technology is progressing rapidly over last decades toward high accelerating gradients and low surface resistance making feasible the particle accelerators operation with high beam currents and long duty factors. However, the coherent RF losses due to high order modes (HOMs) excitation becomes a limiting factor for these regimes. In spite of the operating mode, which is tuned separately, the parameters of HOMs vary from one cavity to another due to finite mechanical tolerances during cavities fabrication. It is vital to know in advance the spread of HOM parameters in order to predict unexpected cryogenic losses, overheating of beam line components and to keep stable beam dynamics. In this paper we present the method of generating the unique cavity geometry with imperfections while preserving operating mode frequency and field flatness. Based on the eigenmode spectrum calculation of series of randomly generated cavities we can accumulate the data for the evaluation the HOM statistics. Finally, we describe the procedure for the estimation of the probability of the resonant HOM losses in the SRF resonators. The study of these effects leads to specifications of SC cavity and cryomodule and can significantly impact on the efficiency and reliability of the machine operation
slides icon Slides TUPAG04 [1.810 MB]  
DOI • reference for this paper ※ https://doi.org/10.18429/JACoW-ICAP2018-TUPAG04  
About • paper received ※ 15 October 2018       paper accepted ※ 28 January 2019       issue date ※ 26 January 2019  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
TUPAG07 Efficient Computation of Lossy Higher Order Modes in Complex SRF Cavities Using Reduced Order Models and Nonlinear Eigenvalue Problem Algorithms 270
  • H.W. Pommerenke, J. Heller, U. van Rienen
    Rostock University, Faculty of Computer Science and Electrical Engineering, Rostock, Germany
  Superconducting radio frequency (SRF) cavities meet the demanding performance requirements of modern accelerators and high-brilliance light sources. For the operation and design of such resonators, a very precise knowledge of their electromagnetic resonances is required. The non-trivial cavity shape demands a numerical solution of Maxwell’s equations to compute the resonant eigenfrequencies, eigenmodes, and their losses. For large and complex structures this is hardly possible on conventional hardware due to the high number of degrees of freedom required to obtain an accurate solution. In previous work it has been shown that the considered problems can be solved on workstation computers without extensive simplification of the structure itself by a combination of State-Space Concatenation (SSC) and Newton iteration to solve the arising nonlinear eigenvalue problem (NLEVP). First, SSC is applied to the complex, closed and thus lossless RF structure. SSC employs a combination of model order reduction and domain decomposition, greatly reducing the computational effort by effectively limiting the considered frequency domain. Next, a perturbation approach based on SSC is used to describe the resonances of the same geometry subject to external losses. This results in a NLEVP which can be solved efficiently by Newton’s method. In this paper, we expand the NLEVP solution algorithm by a contour integral technique, which increases the completeness of the solution set.  
slides icon Slides TUPAG07 [11.204 MB]  
DOI • reference for this paper ※ https://doi.org/10.18429/JACoW-ICAP2018-TUPAG07  
About • paper received ※ 18 October 2018       paper accepted ※ 24 October 2018       issue date ※ 26 January 2019  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
TUPAG22 Main and Fringe Field Computations for the Electrostatic Quadrupoles of the Muon g-2 Experiment Storage Ring 313
  • E.V. Valetov, M. Berz
    MSU, East Lansing, Michigan, USA
  Funding: This work was supported by the U.S. Department of Energy under Contract DE-FG02-08ER41546 and by Fermi Research Alliance for U.S. Department of Energy under Contract DE-AC02-07CH11359.
We consider semi-infinite electrostatic deflectors with plates of different thickness, including plates with rounded edges, and we calculate their electrostatic potential and field using conformal mappings. To validate the calculations, we compare the fringe fields of these electrostatic deflectors with fringe fields of finite electrostatic capacitors, and we extend the study to fringe fields of adjacent electrostatic deflectors with consideration of electrostatic induction, where field falloffs of semi-infinite electrostatic deflectors are slower than exponential and thus behave differently from most magnetic fringe fields. Building on the success with electrostatic deflectors, we develop a highly accurate and fully Maxwellian conformal mappings method for calculation of main fields of electrostatic particle optical elements. A remarkable advantage of this method is the possibility of rapid recalculations with geometric asymmetries and mispowered plates. We use this conformal mappings method to calculate the multipole terms of the high voltage quadrupole used in the storage ring of the Muon g-2 Experiment (FNAL-E-0989). Completing the methodological framework, we present a method for extracting multipole strength falloffs of a particle optical element from a set of Fourier mode falloffs. We calculate the quadrupole strength falloff and its effective field boundary (EFB) for the Muon g-2 quadrupole, which has explained the experimentally measured tunes, while simple estimates based on a linear model exhibited discrepancies up to 2%.
slides icon Slides TUPAG22 [3.780 MB]  
DOI • reference for this paper ※ https://doi.org/10.18429/JACoW-ICAP2018-TUPAG22  
About • paper received ※ 15 October 2018       paper accepted ※ 28 January 2019       issue date ※ 26 January 2019  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)  
WEPLG05 Review of Spectral Maxwell Solvers for Electromagnetic Particle-in-Cell: Algorithms and Advantages 345
  • R. Lehé, J.-L. Vay
    LBNL, Berkeley, California, USA
  Electromagnetic Particle-In-Cell codes have been used to simulate both radio-frequency accelerators and plasma-based accelerators. In this context, the Particle-In-Cell algorithm often uses the finite-difference method in order to solve the Maxwell equations. However, while this method is simple to implement and scales well to multiple processors, it is liable to a number of numerical artifacts that can be particularly serious for simulations of accelerators. An alternative to the finite-difference method is the use of spectral solvers, which are typically less prone to numerical artifacts. In this talk, I will review recent progress in the use of spectral solvers for simulations of plasma-based accelerators. This includes techniques to scale those solvers to large number of processors, extensions to cylindrical geometry, and adaptations to specific problems such as boosted-frame simulations.  
slides icon Slides WEPLG05 [2.861 MB]  
DOI • reference for this paper ※ https://doi.org/10.18429/JACoW-ICAP2018-WEPLG05  
About • paper received ※ 06 November 2018       paper accepted ※ 28 January 2019       issue date ※ 26 January 2019  
Export • reference for this paper using ※ BibTeX, ※ LaTeX, ※ Text/Word, ※ RIS, ※ EndNote (xml)