particle accelerator conference, pac09, vancouver, british columbia, canada, 4 th -8 th may 2009
DESCRIPTION
- PowerPoint PPT PresentationTRANSCRIPT
THE INFLUENCE OF CELL MISSALIGNMENTS AND CAVITY PERTURBATIONS ON LARGE ACCELERATING
STRUCTURES INVESTIGATED USING MODE MATCHING AND THE GLOBABLISED SCATTERING
MATRIX TECHNIQUE I.R.R Shinton, R.M. Jones and V. Khan; The Cockcroft institute of accelerator science, Daresbury, WA4 4AD, UK; The school of physics and astronomy, The University of Manchester, Manchester, M13 9PL, UK
1 1.5 2 2.5 3 3.5 4 4.50
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
frequency: GHz
Sca
tter
ing
coef
ficie
nt
Full Structure Simulation
Cascaded Simulation
Particle Accelerator Conference, PAC09, Vancouver, British Columbia, Canada, 4 th-8th May 2009
1 1.5 2 2.5 3 3.5 4 4.50
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
Frequency: GHz
S21
2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 4.4
0.1
0.2
0.3
0.4
0.5
0.6
0.7
frequency: Ghz
Scatt
ering c
oeff
icie
nt
Tesla 9cell Sliced
Tesla 9cell
GSM A Powerful Electromagnetic Technique
Fast, efficient, accurate simulations
Efficient incorporation of defects into the calculation
Large scale simulations Accurate simulations for large scale structures
The generalised scattering matrix technique is a mature RF method [1] requiring little in the way of computational resources or time allowing large structures, beyond the means of the traditionally employed numerical techniques, to be modelled. The generalised scattering matrix technique (and similar methods such as the coupled scattering calculation CSC [2], [3], [4]) has been shown to be capable of rapidly and accurately simulating structures [5]. The technique is exact as demonstrated in [5] and is capable of incorporating misalignments and defects into the calculation in an efficient manner, refer to [6], allowing rapid RMS calculations to be preformed. From the GSM calculation (or a similar scattering matrix technique) the e.m. field may be re-derived in by using GSM in conjunction with a number of techniques. Presented here is a GSM formulation which allows rapid field calculation across an accelerating structure in which we have employeed a GSM technique in conjunction with a mode matching method.
Comparison between the cascaded and fully simulated structure for the dominant TE11 mode as a function of the TE11 mode for the S21 matrix for a structure constructed from 3 TESLA middle cells.
The cascaded results in very good agreement with the fully simulated results with an average error of 0.014% but the amount of computational time saved through cascading is considerable; the total time required to calculate the fully simulated structure was 29hrs 25mins and 20sec, while the cascading calculation (using the referenced unit cell S matrix values) took a mere 1.5sec to run on a single processor machine.
GSM result for the dominant TE11 mode as a function of the TE11 mode for the S21 matrix for the FLASH accelerator – 6 TESLA cryomodules – 432 cells.
Large scale simulations beyond the means of traditionally employed numerical techniques such as the finite element method can be achieved using GSM.
The generalised cascaded S21 matrix for a 9-cell TESLA structure comparison between an unperturbed structure and the RMS of a randomly perturbed structure for the TE11 mode scattered into the TE11 mode. Here the middle cells in the structure were randomly perturbed using slice sizes of 0mm, 1mm, 2mm and 3mm with the RMS taken over 20 different simulations (representing 20 possible fabricated cavities).
The inclusion of realistic defects for an RMS calculation in traditionally used numerical techniques is prohibitively time consuming as it will require re-meshing of the problem domain. A GSM calculation circumvents this issue by only altering the specific regions (where the defects are located) in its calculation to include these effects directly.
1 1.5 2 2.5 3 3.5 40
0.5
1
1.5
2
2.5
3
3.5
4
4.5
5
frequency: Ghz
S21
1
x10
-3 (
%)
GSM compilation error across a hypothetical symmetric 100m structure cascading from any cell in such a structure should result in the same answer – showing that the error obtained from cascading from any cell is less that 5x10-3%.GSM is capable of consistent accurate answers for very large structures, as seen in the above hypothetical example
TH5PFP044
a b
A Mode Matching GSM Technique for field determination YaZaYaZUS tt ˆˆ2
1
21
USaS 2111
UaYaZUS t 1
22ˆ2
2212 SUaS
22
12
022
mnn
nnnm
XbaXXJb
baXJXaa
212
011
021 n
N
nzz eSTST
2212
011
021|| nn
N
nzz eYSTST
2111
1
2211011 STSTSTSTUTS ggggz
21
1
1122021 STSTSUTS ggz
Region 2 electric field
3
13
2321
31,3 n
N
n
zzikn YeenS n
N
n
zzikn
neenS13
2321||
31,3
Region 3 electric field
N
n
zzikn
zzikn
nn eenSee11
1111
11||
11 1,1
1
11
1111
11
11 1,1 n
N
n
zzikn
zzikn YeenSee nn
Region 1 electric field
The S matrix at the transition is obtained by mode matching the fields in terms of the modal amplitudes C, in which we consider the waves to be propagating from one port to the other, in this manner we obtain:
In the above equations is the identity matrix and is the impedance and is the admittance. Using the above equations the S matrix of any structure may be constructed using a series of WN steps and GSM. The e.m... field in any structure can be subdivided into one of three different regions, as depicted in diagram 2. Region 1 lies between z=0 and z=1, here the S matrix consists of the GSM formed z=z1 and z=z3. Region 2 lies between z=z1 and z=z2, the derivation of this S matrix was previously determined by the authors of this paper in [3].
UZ Y
2
3
4
5
6
7
8
9
10
11
12
13
Modematching is a well established technique [4], [5] and has been applied to general accelerator problems in the past for sharp transitions consisting of adjoining WN sections [6], [7], [8], [9], [10], [11]. The mode matching procedure presented here differs from previous studies. Herein we will consider the case for a propagating monopole mode being launched form the starting port. We analytically derive the S matrix for a WN or a NW junction. Let us consider the circular WN junction in diagram 1 with n modes in region 1 and m modes in region 2 where m<=n. The analytical normalised e.m. fields in a circular beam pipe have been derived by N. Marcuvitz [12]. By mode matching the electric field at the interface between the regions i.e. we obtain what shall be referred to as the scalar product in equation 1 (where the subscripts represent the modes) is the electric field, the ^ superscript represents region II quantities.
aee ˆ
a
e
Any subsection of a structure will be classified as a region 2 section, only the first and last sections are classified as region 1 and 3 respectively. The S matrix of any region 2 subsection consists of a propagating matrix representing the infinite propagating wave series…..and a reflected matrix… representing the infinite reflected wave series both which are decaying within the subsection [3]. Here the superscripts I and II represent the GSM’s left and right of the subsection respectively, ……..………...(……………. ) and g is the length of the subsection. Region 3 lies between z=z2 and z=z3, here the S matrix consists of the GSM between z=0 and z=z2. Since the S matrices are directly proportional to the modal amplitudes then after applying the mode matching procedure we obtain the following formulae for the three regions:
021S
011S
ziknmz
neT 2 0,1 mnmn
Diagram 1: Diagram of a WN transition of radius b in the Wide section (region I) and of radius a in the Narrow section (region II)
Diagram 2: Diagram representing the three classes of regions that are encountered in a structure. Here a WNW structure is used as an example.
1
An exact analytical formulation for cavities
cos12
k
15
The relationship between phase and suseptance can be derived using either transmission line theory or from considering the ABCD matrix [5]:
lkBlk zz sin
2coscos
klB
21cos
Bethe perturbation theory allows an analytical formulation, which gives the constant as:
Combining GSM and mode matching in conjunction with analytical analysis of Bethe [13] allows exact analytical expressions to be obtained for the suseptance …..of a structure (which is equivalent to …., the coupling constant. Let us consider the case for a monopole accelerating mode; form circuit model theory [14] we can express the relationship between angular frequency and phase for the thin iris approximation as:
kB
XJlb
ak
12
3
3
4
23
21
4
2
3
Xa
kXJbB z
In equation 17 is the iris radius. Rewriting and comparing equations 14 and 16 in terms of …allows the following equation to be derived, which describes the underlying physics of the problem (for the thin iris approximation).
cosa
14
In equation 2 is the cavity length and is the propagation constant: (where b is designated as the equator radius and is the zero order Bessel function root of the dominant mode). If …………then equation 15 is simplified to:
zkl 22 bXckz
X
1kl
16
17
18
19
The dispersion curve obtained from the combination of equations 15 and 18 is limited due to the thin iris approximation. An alternative method which circumvents this limitation is to consider the S matrix of a cell in an infinite periodic structure (obtained from GSM):
1,12
1,11,11cos
21
211
221
S
SS
2
1
2221
1211
2
1
a
a
SS
SS
b
b
11,1
1,12
11
11
S
iSB
20
21
After some rearranging we obtain the following equation:
An exact formulation for the dominant mode can also be obtained by considering a load suseptance in parallel with a transmission line.
iB
In equation 20 a represent the incident waves and b the reflected waves on each port, the subscripts refer to the ports. As this cell is in an infinite periodic structure: 21 aeb i 12 aeb i
References[1] I.R.R. Shinton and R.M. Jones, PAC, 2007, p2218; [2] I.R.R. Shinton and R.M. Jones, 13th International workshop on RF Superconductivity, 2007, WEP80; [3] I.R.R. Shinton and R.M. Jones, EPAC, 2008, p619; [4] R.E. Collin, Foundations for microwave engineering,1966; [5] R.E. Collin, Field theory of guided waves, 1991 ; [6] U.van. Rienen, AIP Conf. Proc, 1997, Vol 391, p89; [7] S.A. Heifets and S.A. Kheifets, SLAC-PUB-5907; [8] S.A. Heifets and S.A. Kheifets, SLAC-PUB-5562; [9] V.A.Dolgashev, ICAP, 1998, p149; [10] V.Dolgashev, T.Higo,Proceedings of the 8th international workshop on linear colliders, LC97,1997.; [11] L. Carin, Computational analysis of cascaded coaxial and circular waveguide discontinuities, Masters thesis, University of Maryland, 1986, 12] N.Marcuvitz, The wave guide handbook, 1986; [13] H. Henke, TET-Note 96/04, August 1996; [14] K.L.F. Bane SLAC-PUB-5783; [15] L.F. Wang, Y.Z. Lin, T.Higo and K. Takata, PAC, 2001, pg 3045; [15] R.M.Jones and V.F.Khan, LINAC08, 2008, MOP105
Number of cells HFSS R/Q GSM R/Q Percentage difference Time HFSS:s Time GSM:s1 51.76 51.18 1.12 55 1.133 97.57 98.93 1.37470939 559 2.936 203 197.8 2.628918099 1658 5.15
Here we see a slice of the results for a WNW transition for the variation of iris thickness and iris radius for a particular equator radius b=9.88mm at 120deg phase – the technique allows a full parameter sweep. The colour axes are the frequency in GHz and the R/Q in Ohms respectively
0 20 40 60 80 100 120 140 160 18011.84
11.86
11.88
11.9
11.92
11.94
11.96
11.98
12
12.02
12.04
Phase:deg
f:GH
z
GSM Analytical
HFSSv11 eigen infinity
0 20 40 60 80 100 120 140 160 18011.84
11.86
11.88
11.9
11.92
11.94
11.96
11.98
12
Phase:deg
f:GH
z
GSM analytical
KN7cHFSSv11 eigen infinity
Rapid cavity optimisation and large scale simulations
In this example we are using a simple sharp WNW transition to demonstrate the usefulness of the techniqueHere we are comparing analytical dispersion curve obtained from using KN7c (a longitudinal mode matching code), The analytical S matrix calculation and the eigen mode solution from HFSSv11 using linear elements.Note eigen mode values have been extrapolated to infinite mesh valuesMaximum error between all three codes is less than 2MHz
The technique can easily be applied to cavities with curved surfaces in which the curvature can be approximated by a series of conjoined sharp transitions.Here we see the technique applied to one of the alternate CLIC geometries [8]. Using 0.001GHz steps 30 modes per port and 50 segments to approximate the iris in the GSM mode matching method (equations 1-5 and 21)Note eigen mode values have been extrapolated to infinite mesh values for HFSSv11Maximum error between the codes is less than 3MHz
The method utilising equations 1-13 and equation 21 allows an analytical dispersion curve to be obtained as in a single calculation, calculations of the electromagnetic fields within the cell are preformed concurrently. This allows investigation of the parameter space by varying aspects of the geometry of the cell in question iteratively and recording the effect on the desired design parameters. Any axisymmetric cavity geometry can be represented by the technique previously outlined, in which curved geometry may be approximated as a series of conjoined sharp transitions. The versatility and rapid nature of the method make it an ideal for rapid cavity optimization and large scale calculations.
The time taken to generate the entire Dispersion curve using equation 21 (the GSM method) was less than 1min, while the meshed based numerical FEM method of HFSSv11 (eigen solver) required 45mins or mode per point!
The time taken to generate the entire Dispersion curve using equation 21 (the GSM method) was less than 1min, while the meshed based numerical FEM method of HFSSv11 (eigen solver) required 30mins or mode per point!
Each R/Q calculation required less than 1s (using a PC equipped with a 3GHz CPU and 4GHz of RAM ) – thus allowing large parameter sweeps/ large scale simulations to be undertaken
Compared RQ’s calculated using GSM mode matching to those obtained using HFSSv11 for an idealised sharp WNW structure based on the design from [15] i.e. a 2/3 phased structure.The driven modal solver in HFSSv11 was used, all simulations were conducted using the discrete solver at 11.93GHz (since 1, 3 and 6 cells all share this resonant frequency), simulations were run until 0.2% convergence was achieved. Meshing for the cells were 7401, 65056 and 154215 quadratic mesh cells.