De Prisco, Renato; Comunian, Michele; Grespan, Franscesco; Piscent, Andrea; Karlsson, Anders

Similar documents
Preliminary design studies of a 100 MeV H /H + LINAC as injector for SNS synchrotron/ads LINAC

Developmental Studies of High Current Proton Linac for ADS Program

Generating the optimal magnetic field for magnetic refrigeration

Analytical Studies of the Influence of Regional Groundwater Flow by on the Performance of Borehole Heat Exchangers

Design of a Double Spoke Cavity for C-ADS

Resistivity survey at Stora Uppåkra, Sweden

Sound absorption properties of a perforated plate and membrane ceiling element Nilsson, Anders C.; Rasmussen, Birgit

Linac JUAS lecture summary

COMPARISON OF THERMAL HAND MODELS AND MEASUREMENT SYSTEMS

From the Wideröe gap to the linac cell

On the physical limit of radar absorbers

ELECTROMAGNETIC SIMULATION CODES FOR DESIGNING CAVITIES -1

Compact Linac for Deuterons (Based on IH Structures with PMQ Focusing)

Aalborg Universitet. Lecture 14 - Introduction to experimental work Kramer, Morten Mejlhede. Publication date: 2015

Beam Dynamics for CSNS Linac. Jun Peng September. 18, 2012

RF LINACS. Alessandra Lombardi BE/ ABP CERN

Resistivity imaging for mapping of quick clays

Comparative Assessment of HIPPI Normal Conducting Structures

The dipole moment of a wall-charged void in a bulk dielectric

A numerical study of vortex-induced vibrations (viv) in an elastic cantilever

The Wien Bridge Oscillator Family

Evaluation of In-Vacuum Wiggler Wakefield Impedances for SOLEIL and MAX IV

Transmuted distributions and extrema of random number of variables

ESS Linac Beam Physics Design Update

Geoelectrical and IP Imaging Used for Pre-investigation at a Tunnel Project. Danielsen, Berit Ensted; Arver, Henrik; Karlsson, T; Dahlin, Torleif

Lattice Boltzmann Modeling From the Macro- to the Microscale - An Approximation to the Porous Media in Fuel Cells -

Design possibilities for impact noise insulation in lightweight floors A parameter study

A beam theory fracture mechanics approach for strength analysis of beams with a hole

Neutronics calculations for the ITER Collective Thomson Scattering Diagnostics

A comparison study of the two-bladed partial pitch turbine during normal operation and an extreme gust conditions

Higher-order spectral modelling of the diffraction force around a vertical circular cylinder

NEW DEVELOPMENT IN HIGH POWER RFQ ACCELERATORS*

Validity of PEC Approximation for On-Body Propagation

CERN Accelerator School. RF Cavities. Erk Jensen CERN BE-RF

Lecture 5: Photoinjector Technology. J. Rosenzweig UCLA Dept. of Physics & Astronomy USPAS, 7/1/04

Reflection and absorption coefficients for use in room acoustic simulations

A Note on Powers in Finite Fields

Benchmark of Femlab, Fluent and Ansys

BEAM DYNAMICS ISSUES IN THE SNS LINAC

Nicolas PICHOFF. France CEA-Direction d Île de France Service de Physique et Applications des Accélérateurs (SP2A)

Shielding calculations with MCNPX at the European spallation source

ANALYSIS OF HIGH ORDER MODES IN 1.3 GHZ CW SRF ELECTRON LINAC FOR A LIGHT SOURCE

Embankment dam seepage assessment by resistivity monitoring

Reflection of sound from finite-size plane and curved surfaces

Lidar calibration What s the problem?

Sorption isotherms for textile fabrics and foam used in the indoor environment

A Procedure for Classification of Cup-Anemometers

Graphs with given diameter maximizing the spectral radius van Dam, Edwin

Application of the MCNPX-McStas interface for shielding calculations and guide design at ESS

Design of an RF Photo-Gun (PHIN)

Compressor Lattice Design for SPL Beam

Realtime 3D stress measurement in curing epoxy packaging

2D PIM Simulation Based on COMSOL

Physics design of a CW high-power proton Linac for accelerator-driven system

Novel, Hybrid RF Injector as a High-average. Dinh Nguyen. Lloyd Young

Spoke and other TEM-class superconducting cavities. J.L. Muñoz, ESS-Bilbao Academy-Industry Matching Event CIEMAT, Madrid, 27.May.

On numerical evaluation of two-dimensional phase integrals

Geophysical techniques for investigating shallow and deep landslides: results in the framework of PREVIEW project (FP6 UE)

On machine measurements of electrode wear in micro EDM milling

Unknown input observer based scheme for detecting faults in a wind turbine converter Odgaard, Peter Fogh; Stoustrup, Jakob

Citation for published version (APA): Burcharth, H. F. (1988). Stress Analysis. Aalborg: Department of Civil Engineering, Aalborg University.

Optimizing Reliability using BECAS - an Open-Source Cross Section Analysis Tool

Modelling Salmonella transfer during grinding of meat

Lectures on accelerator physics

Short Introduction to CLIC and CTF3, Technologies for Future Linear Colliders

Geometry explains the large difference in the elastic properties of fcc and hcp crystals of hard spheres Sushko, N.; van der Schoot, P.P.A.M.

Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)

Design of a 10 MeV normal conducting CW proton linac based on equidistant multi-gap CH cavities *

Multi (building)physics modeling

RF cavities (Lecture 25)

Determination of packaging induced 3D stress utilizing a piezocoefficient mapping device

Simulating wind energy resources with mesoscale models: Intercomparison of stateof-the-art

An Equivalent Source Method for Modelling the Lithospheric Magnetic Field Using Satellite and Airborne Magnetic Data

Aalborg Universitet. Comparison between Different Air Distribution Systems Nielsen, Peter Vilhelm. Publication date: 2006

Convergence analysis for a class of LDPC convolutional codes on the erasure channel

Radiation field inside the tunnel of the Linear Collider TESLA

Higher order hierarchical H(curl) Legendre basis functions applied in the finite element method: emphasis on microwave circuits

Benchmarking of optimization methods for topology optimization problems

Physics 610. Adv Particle Physics. April 7, 2014

TULIP all-linac: a high gradient linac for proton therapy

Matrix representation of a Neural Network

CEA Saclay Codes Review for High Intensities Linacs Computations

Numerical modeling of the conduction and radiation heating in precision glass moulding

Aalborg Universitet. The number of independent sets in unicyclic graphs Pedersen, Anders Sune; Vestergaard, Preben Dahl. Publication date: 2005

Kommerciel arbejde med WAsP

On the Unique D1 Equilibrium in the Stackelberg Model with Asymmetric Information Janssen, M.C.W.; Maasland, E.

Straightforward and accurate technique for post-coupler stabilization in drift tube linac structures

The introduction to learning Poisson/Superfish & Astra codes. Sara Zarei Winter 96

Steady State Comparisons HAWC2 v12.2 vs HAWCStab2 v2.12

Feasibility of non-linear simulation for Field II using an angular spectrum approach

Longitudinal Dynamics

IFMIF High energy beam line design and beam expansion using non-linear multipole lenses and "step-like" magnet

Proposal to convert TLS Booster for hadron accelerator

Spatial decay rate of speech in open plan offices: the use of D2,S and Lp,A,S,4m as building requirements Wenmaekers, R.H.C.; Hak, C.C.J.M.

Numerical modeling of magnetic induction and heating in injection molding tools

A comparison between two different method of moments codes for periodic problems

Varying accelerating fields

Gao, Xiang; Zhu, Meifang; Rusek, Fredrik; Tufvesson, Fredrik; Edfors, Ove

Transcription:

ESS DTL RF MODELIZATION: FIELD TUNING AND STABILIZATION De Prisco, Renato; Comunian, Michele; Grespan, Franscesco; Piscent, Andrea; Karlsson, Anders 2013 Link to publication Citation for published version (APA): de Prisco, R., Comunian, M., Grespan, F., Piscent, A., & Karlsson, A. (2013). ESS DTL RF MODELIZATION: FIELD TUNING AND STABILIZATION. Paper presented at 4th International Particle Accelerator Conference (IPAC 2013), Shanghai, China. General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal Take down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. L UNDUNI VERS I TY PO Box117 22100L und +46462220000

ESS DTL RF MODELIZATION: FIELD TUNING AND STABILIZATION R. De Prisco, ESS, Lund, Sweden and Lund University, Lund, Sweden M. Comunian, F. Grespan, A. Pisent, INFN/LNL, Legnaro, Italy A. Karlsson, Lund University, Lund, Sweden Abstract The Radio Frequency (RF) design of Drift Tube Linac (DTL) of the European Spallation Source, ESS, has be defined in order to satisfy the accelerating field requirements of beam dynamic studies and to reduce peak field levels in the critical areas. The electro-magnetic field is stabilized with post-couplers. The cells geometries of the DTL are optimized to accommodate permanent magnet quadrupoles (PMQ), to get maximum shunt impedance, to meet the Moretti criterion at the low energy part and to facilitate the mechanical construction. INTRODUCTION The linac of the ESS will accelerate 63 ma of H -ions at 352.21 MHz up to 2.5 GeV. The pulses are 2.86 ms long with duty cycles of 4%. The average beam power on the target will be 5 MW and the peak power 125 MW. The linac will have a normal-conducting front-end up to 77.7 MeV followed by three families of superconducting cavities and a high-energy beam transport to the spallation target [1]. In particular the DTL of the ESS accelerator will accelerate the beam from 3.0 to 77.7 MeV. Permanent Magnet Quadrupoles (PMQs) are used as focusing elements in a FODO lattice. The main properties of the ESS DTL are summarized in Table 1 [2]. Table 1: Tanks properties of the ESS DTL Parameter/Tank 1 2 3 4 Cells 66 36 29 25 E 0 [MV/m] 2.8 to 3.2 3.16 3.16 3.16 E MAX /E Kilp 1.4 1.43 1.39 1.37 φ S -35 to -24-24 -24-24 L T [cm] 7.95 7.62 7.76 7.72 R b [mm] 10 10 11 12 PCs [#] 23 25 28 24 Tuners [#] 24 24 24 24 Modules [#] 4 4 4 4 P Cu,MAX [MW] 0.91 0.91 0.92 0.95 E OUT [MeV] 21.4 41.0 60.0 77.7 P T OT [MW] 2.06 2.12 2.10 2.07 Slug tuners in a DTL should compensate the frequency effects due to construction errors. The evaluation of the static frequency error is done by applying realistic tolerances on the important dimension of the cavity (tank renato.deprisco@esss.se diameter, drift-tube lengths, drift tube diameter, face angles). The tuners have diameter of 90 mm, are distributed uniformly every 30 cm along the tank and are located at 45 with respect to the stem axis in order not to influence the frequency of the PC 0-mode by tuner penetration. The tuner sensitivity is 6.98 (khz/mm) m, linear around 50 mm of penetration. TUNING OF ACCELERATING FIELD Let the z-axis be the symmetry axis of the tank, its origin be on the low energy end wall and the pulsation, or angular resonance frequency, ω (z) be constant along z and equal to ω 0. The effects of a variation of the cell geometries on the pulsation in a DTL can be described by waveguide theory. Then the DTL is replaced by an empty equivalent tank in which the geometric variation of the cells changes the resonance frequency, ω 0, along z such that ω (z) = ω 0 + ω (z). From continuos to discrete tuning Even if it is possible to define a continuous function ω (z), z [0, L T ], where L T is the tank length, it is convenient to define ω(z) as piecewise constant function. Let suppose that the tank is composed of N cells, each of length L n, n [1, N], n N. Let define ω n = L 1 n L n ω (z) dz ω 0. It is possible to show that if ω n ω 0 in each cell, then the electric field along z-axis, E z (z), must satisfy the following equation: d 2 E z (z) dz 2 2ω 0 c 2 ω ne z (z), z [0 n, L n ], (1) where 0 n is the relative origin, along z-axis, of the n-cell. Boundary and continuity conditions L 1 n Let define the piecewise constant electric fields E 0n = Ln 0 n E z (z) dz. Assume that the beam dynamic requirements imply the constraints E 0n, n [1, N]. It is, now, assumed that the total stored energy in RF fields is constant. The knowledge of E 0n gives the possibility to define E z (0), E z (L T ) and de z (z) /dz, z [0, L T ]. By knowing E z (z), utilizing the boundary conditions for TM modes (de z /dz = 0 in 0 and L T ) and by imposing the continuity conditions between adjacent cells it is possible to solve the equation (1), cell by cell, and to find ω n, n [1, N]. In the following the z-dependence of E z and its derivatives are omitted.

The implicit slope As a first tuning step the RF cavity computer code, SU- PERFISH (DTLfish), has been used to change the face angle, α f, in order to reach the target pulsation ω 0 in all cells of the tank. The details of the drift tube section is shown in the Fig. 1. The second step is to assembly all cells of Figure 1: Representative half unit cell geometry longitudinal section: Tank Diameter, D; cell Length, L; gap length, g; drift tube diameter, d; bore Radius, R b ; Flat length, F; inner nose Radius, R i ; outer nose Radius, R o ; face angle, α f. It is useful to divide the tank into three parts, macrocells, each one composed of one or more cells. The first macrocell satisfies the boundary conditions on the end wall. It ramps the accelerating field and compensates for the implicit slope caused by the F cells with total length L F = F n=1 L n. The second macrocell adjusts the slope to match to the boundary conditions at the other end wall. It is composed by L cells with total length L L = N n=n L+1 L n. The third one is located between the first two and composed by (N L) (F + 1) cells with total length L L = N L n=f +1 L n. Instead to solve the equation (1), cell by cell, it is convenient to solve this equation for the macrocells. To do this the the constant C = 2ω 0 /c is introduced. For the first macrocell the ω F must satisfy the equation (2): C ω F sinh (C ω F L F ) E z (0) = S T. (2) For the second macrocell the ω L must satisfy the equation (3): C ω L sin (C ) ω L L L E z (L T L L ) + + cos (C ) ω L L L S T = 0. Tuning of the ESS DTL For the ESS DTL the desired electric field E 0n goes from 2.8 MV/m to 3.2 MV/m in the first tank. It is constant and equal to 3.16 MV/m in the last three tanks. In Fig. 2 the desired E 0n are shown. The accelerating field is calculated by MDTfish, after the tuning of the two macrocells by using equations (2) and (3). Specifically the first three and the last three cells have been used to reach the total slope in the first tank. The first two and the last two cells are used to compensate the implicit slope in the second tank. The first and the last cells are used to compensate the implicit slope in the third and fourth tanks. (3) the tank and simulate the entire structure in SUPERFISH (MDTfish). Even if all the cells, different in length, have the same pulsation, ω 0, the accelerating field is not constant as we can expect from equation (1). This is due to the fact that there is not a perfect mode matching between the adjacent cells built, individually, by DTLfish. The mismatch produces a natural tilt of the accelerating field that must be compensated; from now this tilt is denoted implicit slope. In a ramped gradient tank the implicit slope must be properly added to the desired one to compensate it. In a constant gradient tank the implicit slope is the only one to compensate since the desired one is zero. We denote the sum of the implicit slope and the desired one the total slope, S T. Macrocells Figure 2: Desired E 0, E 0 and percentage Error. Compensation of stems effect Stems connect drift tubes with the inner lateral wall of each tank. They break the z-axial cylindrical symmetry. SUPERFISH can only handle axially symmetric structures. However, the code treats the stem as a volume perturbation and by using the Slater theorem it is possible to calculate, cell by cell, the pulsation shift, ω Stem,n, n [1, N] induced by them. To compensate this perturbations one strategy is to introduce volume perturbations (in the ESS DTL case the face angles have been changed) that induce a pulsation shift ω Comp,n = ω Stem,n, n [1, N]. This strategy has been validated by COMSOL 3D-simulation

in a representative tank with 20 cells, from the first to the twentieth cell of the second tank, and long 1.921 m. In the Fig. 3 the stems effect and the face angles compensation are shown. the inner neighborhood of the cavity lateral wall respectively. The field stabilization is achieved by utilizing Post Couplers (PC). COMSOL has been used to design the PC. PCs can be introduced as series of inductance, L P C, and capacitance, C P C, in parallel with the capacitance C P positioned in the longitudinal center of the drift tubes. The values of L P C and C P C must be choosen in order to stabilize the accelerating field. It is possible to outdistance the PCs until the E 0 is within the 1% of the desired value. The results are confirmed by COMSOL 3D-simulations in a representative tank of 21 cells. In the Fig. 5 the transmission line model and COMSOL result are shown with their relative error in the case of stabilization achieved by using one PC every two cells. Figure 3: Effects of stems and compensative face angles. DTL TRANSMISSION LINE MODEL In the neighborhood of TM 010 mode the electromagnetic field, in each tank of the DTL, can be approximate by the electromagnetic field of a circular cavity filled by anisotropic media [3] (µ is constant and ɛ = ẑẑɛ z + ˆrˆrɛ t ). The equivalent cell of the transmission line is shown in the Fig. 4 with: Figure 4: Equivalent DTL cell of the transmission line DTL model. L S = πω 011 E t (D/2) MAX z (ω011 2 ω2 010 ) H t (D/2) MAX L ; (4) T C P = 1 C S = ω010 2 L ; (5) S π H t (D/2) MAX z ω 011 L T E t (D/2) MAX, (6) where ω 010 and ω 011 are the pulsations of the TM 010 and the TM 011 mode respectively, E t (D/2) and H t (D/2) are the electric and the magnetic transverse fields evaluated in Figure 5: Stabilized field with one PC every two cells obtained by transmission line model and COMSOL. Fixed the number of PCs and their radius in the case of constant gradient tanks the post length are then adjusted such that the frequency of the PC 0-mode is close to the operation frequency of 352.21 MHz (confluence). After reaching of the confluence, the posts in the first tank have been bent, without changing the PC length along its cylindrical symmetry axis, to equalize the product of capacitive coupling (to drift tube) times the field level on both sides of each PC [4]. The compensation of the PC is done with the same strategy that was used for the compensation of the stems. CONCLUSION To design the ESS DTL it was implemented a software, RENdtl, to automatically tuning the accelerating field, compensate the stems and the PCs effect, determine the PCs numbers and their optimum length, define the minimum numbers of tuners, that allows to obtain a really exhaustive RF DTL design. REFERENCES [1] H. Danared, Design of ESS accelerator, IPAC 12, New orleans, May 2012, THPPP071, p. 3904.

[2] M. Comunianet et al., DTL Design for ESS, LINAC 12, Tel-Aviv, September 2012. [3] M.K. Ke et al., Continuous transmission line modeling of a DTL including beam loading and other perturbations, Particle Accelerators, p. 109 (1994). [4] J.H. Billen et al., A 3λ/4 post coupler for drift-tube linacs, IEEE Transaction on Nuclear Science, p. 3184 (1985).