New SLED 3 system for Multi-mega Watt RF compressor. Chen Xu, Juwen Wang, Sami Tantawi

Similar documents
Conceptual Design Of An Ideal Variable Coupler For Superconducting Radiofrequency 1.3GHz Cavities. Chen Xu, Sami Tantawi

The Next Linear Collider Test Accelerator s RF Pulse Compression and Transmission Systems

Performance Measurements of SLAC's X-band. High-Power Pulse Compression System (SLED-II)

High Power, Magnet-free, Waveguide Based Circulator Using Angular-Momentum Biasing of a Resonant Ring

Christopher Nantista ISG-X SLAC June 17, 2003

Development of a 20-MeV Dielectric-Loaded Accelerator Test Facility

Projects in microwave theory 2017

Design, Development and Testing of RF Window for C band 250 kw CW Power Klystron

International Technology Recommendation Panel. X-Band Linear Collider Path to the Future. RF System Overview. Chris Adolphsen

RF AND MICROWAVE ENGINEERING

L-BAND COPLANAR SLOT LOOP ANTENNA FOR INET APPLICATIONS

Waveguides. Metal Waveguides. Dielectric Waveguides

ELEC4604. RF Electronics. Experiment 2

R.K.YADAV. 2. Explain with suitable sketch the operation of two-cavity Klystron amplifier. explain the concept of velocity and current modulations.

Projects in microwave theory 2009

DEVELOPMENT OF A BETA 0.12, 88 MHZ, QUARTER WAVE RESONATOR AND ITS CRYOMODULE FOR THE SPIRAL2 PROJECT

A Pin-Loaded Microstrip Patch Antenna with the Ability to Suppress Surface Wave Excitation

Evaluation of HOM Coupler Probe Heating by HFSS Simulation

Coupler Electromagnetic Design

RF simulations with COMSOL

On the RF system of the ILC

Design of a Novel Compact Cup Feed for Parabolic Reflector Antennas

Electromagnetics, Microwave Circuit and Antenna Design for Communications Engineering

Normal-conducting high-gradient rf systems

CAGE CAVITY: A LOW COST, HIGH PERFORMANCE SRF ACCELERATING STRUCTURE*

REVIEW OF HIGH POWER CW COUPLERS FOR SC CAVITIES. S. Belomestnykh

MICROWAVE MICROWAVE TRAINING BENCH COMPONENT SPECIFICATIONS:

Cavity BPMs for the NLC

THE MULTIPACTING STUDY OF NIOBIUM SPUTTERED HIGH-BETA QUARTER-WAVE RESONATORS FOR HIE-ISOLDE

30 GHz rf components for CTF3 (and CLIC) next part. I. Syratchev

arxiv:physics/ v1 [physics.optics] 28 Sep 2005

A DUAL-PORTED PROBE FOR PLANAR NEAR-FIELD MEASUREMENTS

RF Design of Normal Conducting Deflecting Cavity

RF STATUS OF SUPERCONDUCTING MODULE DEVELOPMENT SUITABLE FOR CW OPERATION: ELBE CRYOSTATS

3. (a) Derive an expression for the Hull cut off condition for cylindrical magnetron oscillator. (b) Write short notes on 8 cavity magnetron [8+8]

HYBRID ARRAY ANTENNA FOR BROADBAND MILLIMETER-WAVE APPLICATIONS

Physics Requirements Document Document Title: SCRF 1.3 GHz Cryomodule Document Number: LCLSII-4.1-PR-0146-R0 Page 1 of 7

ABSTRACT 1 CEBAF UPGRADE CAVITY/CRYOMODULE

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad

Review of New Shapes for Higher Gradients

Design of Rotman Lens Antenna at Ku-Band Based on Substrate Integrated Technology

2 Theory of electromagnetic waves in waveguides and of waveguide components

Christopher Nantista ISG8 SLAC June 25, 2002

High Power Over-Mode 90 Bent Waveguides for Circular TM 01 and Coaxial TEM Mode Transmission

A Millimeter Wave Center-SIW-Fed Antenna For 60 GHz Wireless Communication

The European Spallation Source. Dave McGinnis Chief Engineer ESS\Accelerator Division IVEC 2013

HIGH POWER INPUT COUPLERS FOR THE STF BASELINE CAVITY SYSTEM AT KEK

CHAPTER 2 MICROSTRIP REFLECTARRAY ANTENNA AND PERFORMANCE EVALUATION

Dual Band Dielectric Resonator Filter (DBDRF) with Defected Ground Structure (DGS)

HIGH POWER COUPLER FOR THE TESLA TEST FACILITY

Design and RF Measurements of an X-band Accelerating Structure for the Sparc Project

ACE3P and Applications to HOM Power Calculation in Cornell ERL

Rectangular Patch Antenna to Operate in Flame Retardant 4 Using Coaxial Feeding Technique

CIRCULARLY POLARIZED SLOTTED APERTURE ANTENNA WITH COPLANAR WAVEGUIDE FED FOR BROADBAND APPLICATIONS

Performance Analysis of Different Ultra Wideband Planar Monopole Antennas as EMI sensors

QUADRI-FOLDED SUBSTRATE INTEGRATED WAVEG- UIDE CAVITY AND ITS MINIATURIZED BANDPASS FILTER APPLICATIONS

SUPERCONDUCTING PROTOTYPE CAVITIES FOR THE SPALLATION NEUTRON SOURCE (SNS) PROJECT *

DESIGN OF A FABRY-PEROT OPEN RESONATOR AT RADIO FREQUENCIES FOR AN MgB2 TESTING PLATFORM

Possible High Power Limitations From RF Pulsed Heating *

OVERVIEW OF INPUT POWER COUPLER DEVELOPMENTS, PULSED AND CW*

APPLIED ELECTROMAGNETICS: EARLY TRANSMISSION LINES APPROACH

Energy Recovering Linac Issues

Multimoded RF Systems for Future Linear Colliders. Sami G. Tantawi

Citation Electromagnetics, 2012, v. 32 n. 4, p

Design of S-band re-entrant cavity BPM

RAJIV GANDHI COLLEGE OF ENGINEERING AND TECHNOLOGY Kirumampakkam,Puducherry DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING

A DUAL-PORTED, DUAL-POLARIZED SPHERICAL NEAR-FIELD PROBE

Development of a 20 MeV Dielectric-Loaded Test Accelerator

UNIVERSITI MALAYSIA PERLIS

Broadband Balanced Microstrip Antenna Fed by a Waveguide Coupler

Department of Electrical Engineering University of North Texas

A Multi-Moded RF Delay Line Distribution System (MDLDS) for the Next Linear Collider *

A HIGH-POWER LOW-LOSS MULTIPORT RADIAL WAVEGUIDE POWER DIVIDER

MAGNETO-DIELECTRIC COMPOSITES WITH FREQUENCY SELECTIVE SURFACE LAYERS

04th - 16th August, th International Nathiagali Summer College 1 CAVITY BASICS. C. Serpico

A Novel Dual-Band SIW Filter with High Selectivity

DESIGN OF COMPACT PLANAR RAT-RACE AND BRANCH- LINE HYBRID COUPLERS USING POLAR CURVES

SUPPRESSING ELECTRON MULTIPACTING IN TTF III COLD WINDOW BY DC BIAS

Couplers for Project X. S. Kazakov, T. Khabiboulline

TESLA RF POWER COUPLERS DEVELOPMENT AT DESY.

M. Y. Ismail and M. Inam Radio Communications and Antenna Design Laboratory (RACAD) Universiti Tun Hussein Onn Malaysia (UTHM) Batu Pahat, Malaysia

A Novel Meander Line Microstrip Log-Periodic Dipole Antenna for Dual-Polarized Radar Systems

Bandwidth Enhancement in Microstrip Rectangular Patch Antenna using Defected Ground plane

Generation and Absorption of the Untrapped Wakefield Radiation in the 3.9 GHz LCLS-II Cryomodule

Compact Triple-Band Monopole Antenna with Inverted-L Slots and SRR for WLAN/WiMAX Applications

Microwave cavities. Physics 401, Spring 2017 Eugene V. Colla

REVIEW ON SUPERCONDUCTING RF GUNS

DESIGN OF A NOVEL MICROSTRIP-FED DUAL-BAND SLOT ANTENNA FOR WLAN APPLICATIONS

Performance of the Prototype NLC RF Phase and Timing Distribution System *

High acceleration gradient. Critical applications: Linear colliders e.g. ILC X-ray FELs e.g. DESY XFEL

A Modified Elliptical Slot Ultra Wide Band Antenna

Circular Patch Antenna with CPW fed and circular slots in ground plane.

Index. Cambridge University Press Silicon Photonics Design Lukas Chrostowski and Michael Hochberg. Index.

Principles of Planar Near-Field Antenna Measurements. Stuart Gregson, John McCormick and Clive Parini. The Institution of Engineering and Technology

High-power multimode X-band rf pulse compression system for future linear colliders

STATE OF THE ART IN EM FIELD COMPUTATION*

Design of a full-band polariser used in WR-22 standard waveguide for satellite communications

APPLICATION OF HIGH FREQUENCY SYSTEM FOR IMPROVEMENT OF OUTPUT PROPERTIES OF STANDING WAVE ELECTRON LINEAR ACCELERATORS

Design of ESS-Bilbao RFQ Linear Accelerator

- reduce cross-polarization levels produced by reflector feeds - produce nearly identical E- and H-plane patterns of feeds

Transcription:

New SLED 3 system for Multi-mega Watt RF compressor Chen Xu, Juwen Wang, Sami Tantawi SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94309, USA Electronic address: chenxu@slac.stanford.edu Abstract A compact X band SLED is introduced for X band RF compressing application. This SLED compressor consists two major parts: a rectangular to circular waveguide converter and an overmoded spherical cavity. The RF compressor is designed to convert 50 magawatt X band RF power with pulse length 1.5 microseconds and deliver 200 megawatts with pulse length 106 nanoseconds to the X band accelerating structure. 1. Introduction SLAC has successfully designed and utilized X band and S band RF power compressor for NLC project and LCLS project. SLED is named as SLAC's Energy Doubling accelerating gradient. The conventional idea for SLED is that to store the first part energy of input RF pulse by cavities and releases the stored energy in cavities along with the second part of the pulse. Since the pulse power phase has been turned 180 degree, the output pulse is a combination in phase of the leakage power and input power. Thus, the pulse power is compressed and increased. Compression ratio, namely, the duration ratio between before and after compression is around 6. There are two critical points for this SLED technique: 1. The cavities must be over-coupled in order to obtain a output power higher than the original input power. 2. The pulse delivery network must be directional in order to combine leakage power and the flipped input pulse. The schematic system is shown in fig.1.

Fig.1. schematic drawing of the SLED microwave network. A 3db coupler is a four ports directional coupler. It is designed to convert a rectangular waveguide TE10 mode to equally feed two cavities without leakage to the fourth ports in fig.2. Fig.2. High power short slot hybrid 3db coupler. Moreover, when the power is leaked from these two cavities, the power from port 2 and 3 are canceling output to port 1 and combining output to port 4. This directional coupling characteristic is essential for the SLED technique. This conventional SLED utilized two single-moded pillbox cavities with quality factor around 10 5. The cavity mode at 2.8GHz is TE015 mode. The mode has a pattern showing in fig.3. The mode has

Q~10 5 and coupling coefficient ( ~5), therefore, the loaded Q L is 90KHz. 4 2 10 and mode bandwidth is The whole structure is utilized in the current LCLS accelerator main linac, and it helps boosting the LCLS accelerated electron energy from 6GeV to 12GeV. However, the two bulky cavities and a 3 db coupler have size of 1.5m1.5m 1.5m for S band operation. Our goal is reduce its size and increasing the compression ratio to around 10. This goal yields higher compressed power and shorter output pulse for X band operation. 2. Design concept We recently design three versions of dual modes circular polarized couplers. These couplers can convert TE10 rectangular waveguide into two TE11 modes with different polarization in a circular waveguide. This over-moded structure can replace the conventional 3db couplers even though it only has three physical ports. The over-moded output can feed an overmoded cavity to avoid the second cavity in the conventional SLED design. Similarly, the modes inside this cavity will have the same, Q and degeneration frequency. An ideal candidate for this application is a spherical cavity. Spherical cavity has analytical solution for TE and TM modes. A TE mnp or TM mnp modes are labelled for azimuthal, polar and radial directions. The TE mnp field scalar potential is formulated in following. The field potential patterns are shown in Fig. 4. The electromagnetic field for TE mnp can be derived the equation below : E E H k F r 2 1 Fr 1 Fr H r sin r r E 2 2 r 0 r ( ) 2 r 2 1 Fr 1 Fr H r r sin r

Fig.4. The electromagnetic field pattern for TE 114 mode. Note these three modes have identical field pattern but different polarizations. In a spherical structure, the mode orientation is determined by the excitation signal from a give port. For example, dipole input feeds on the spherical cavity wall will coupling TE114 dipole modes. However, the TE014 cannot be excited by any dipole mode in a circular waveguide, thus this TE014 would have very small coupling. In order to avoid close neighbor modes and obtain higher Q 0, TEn14 modes are chosen. In this Ten14 family, there are three degeneration modes: one TE014 and two TE114 modes. When spherical cavity has a radius of 5.8749cm, the TEn14 family has a degeneration frequency of 11.424GHz and the Q0 is around 10 5. We need to design a coupler to this spherical cavity, and make large enough to obtain the maximum output. Assumed the traveling back wave accelerating structure to feed has a feeding time around 100ns, the gradient reaches maximum when is around 9. The outputs of the circular polarizer are from one circular waveguide, and each polarization mode in this waveguide is used to feed the different modes in the spherical cavity. Two single mode cavities is replaced by two polarization in one cavity, and spherical cavity and polarizer supple two modes and feeds with identical RF characteristics. 3. Characteristic of the dual circular polarizer and spherical ball design 3.1. Dual modes Circular converter The polarizer has three physical ports but the output circular waveguide supports two modes. These two TE11 modes have the same power magnitude, different polarization and 90 degree phase delay. Using this two TE11 circular modes to feeds two TE114 modes in the cavity, the input power are split into two cavity modes which canceling each other at the output rectangular waveguide. The scattering matrix for this structure is shown below:

The electric field on the surface of this structure is shown in Fig. 4. Fig.4: The complex magnitude of E field on coupler surface. This structure has a fairly wide RF operation bandwidth. The output of this coupler is a TE11 circular mode rotating on port 3. The frequency spectrum of scattering parameters is shown in Fig 5, and it shows that the port 1 has neither reflection nor a transmission, and all power equally split into two modes in port 3.

Fig.5. Scattering parameters frequency spectrum. 3.2. Spherical cavity design. The radius of a spherical cavity is 5.8749cm when the resonating frequency is 11.424GHz. The unloaded Q is 5 110 and an ideal goal is around 10. It is convenient to connect the spherical cavity directly with the circular waveguide. An iris is introduced to control the coupling. The larger iris radius will bring larger. While the operating frequency 11.424GHz is beyond the cutoff frequency of the circular waveguide, the waveguide with an iris can cut off the 11.424GHz field within the iris region. However, the thickness of this iris aperture is small, thus the evanescent TE11 mode can still couple into the cavity. By optimizing this thin iris, one can obtain the coupling coefficient around 10. Fig.6 shows the cavity E field on the surface of a given snapshot.

Fig.6 The snapshot of E field pattern at different frequency inside of the spherical cavity. Left is the TE014 mode and right is one polarized TE114 mode. The change of the iris radius can vary the coupling coefficient, because the leakage to the circular waveguide is determined by the level of cutoff at the iris region. A larger iris can increase, but the change is not a linear variation. Meanwhile, the TE114 resonating frequency is also change with change of the iris radius. The scattering parameters are swept at different frequency in Fig. 7 and of each resonating modes are calculated with different iris radius. 5.8mm

6.2mm

6.6mm 7mm 8mm

9.5mm Fig.7. the scattering matrix and Smith chart of each frequency sweep. Since the simulation from HFSS is in a steady state, one can extract of each iris radius setups from the scattering parameter by the equation. 1 S 1 S 11 11 Smith charts are used to determine is under coupled or over coupled. The TE114 resonating frequency, scattering matrix and the dependency with the iris radius are shown in Fig8.

11.442 11.44 11.438 11.436 11.434 11.432 11.43 11.428 11.426 11.424 nim TE114 frequency vs Iris 3 4 5 6 7 8 9 10 0-0.5 Scattering vs Iris 3 4 5 6 7 8 9 10-1 -1.5-2 -2.5-3 -3.5-4 -4.5-5

Beta vs Iris 180 160 140 120 100 80 60 40 20 0 3 4 5 6 7 8 9 10 Fig.8, The resonating frequency, scattering matrix and dependency with the change of the iris radius. 3.3 Combination. A copper sheet with thickness of 5mm is machined to this structure. The geometry of this SLED is shown in Fig. 9. A snapshot of the surface E field is shown in Fig 10, and the input power from WR 90 is 1W. The peak E surface occurs at the waveguide, not on the surface of spherical cavity. When the input power is 50MW, the peak surface E field is 8.63MV/m.

Fig.9. the actual mechanical design of SLED system. 4. Frequency response, Transient response and Multi-physics concern. 4.1 Frequency response:

The cavity is simulated with frequency scan from GHz to GHz, and the frequency sweep is shown in Fig.7. Smith chart is also plotted in the fig.10. The center resonating peak suggests the two TE114 modes and the smaller peak on the right indicates the TE014 with very weak coupling. The TE114 modes have two polarizations, and the fields patterns are identical except the orientation. D Fig.10: The scattering parameters scan of driving frequency. 4.2 Transient response: Input pulse of this SLED compressor is a square pulse with a modulated rising time and at certain time duration, the input pulse phase is flipped 180 degree. The pulse is both plotted in time domain and frequency domain. The spectrum contains the center frequency and side band frequencies lobes and the cavity resonates at the center frequencies. After convoluting the cavity spectrum in Fig 7 with input pulse spectrum in fig.10, we inverse Fourier transform this output back to time domain and obtain an output with a compressed pulse in fig 11.

6 4 2 5600 5700 5800 5900 6000 2.0 1.5 1.0 0.5 2000 4000 6000 8000 10000 Fig.11: The scattering parameters scan of driving frequency.

1.0 0.5 0.0 5600 5700 5800 5900 6000 0.5 1.0

2.0 1.5 1.0 0.5 2000 4000 6000 8000 10000 2.0 1.5 1.0 0.5 2000 4000 6000 8000 10000 4.3 Thermal concern: The SLED spherical cavity has fairly high Q0, thus the loss on the cavity wall is small compared with the energy stored in the cavity. The major RF loss is from the waveguide, and the power loss on the convertor and the system within a RF cycle is shown in Fig.12. The input power is normalized to 1W. The integrated RF surface loss are W and W respectively.

Fig. 12. RF loss distribution on the converter and the system. By using this RF loss as input load, ANSYS workbench is used to calculate temperature rise on the copper internal surface. Note, the SLED is operated in a pulse mode, and the repetition rate of an X band klystron is 120hz and duty factor is 0.18. This means the HFSS result which is a CW calculation needs to multiply this duty factor to obtain averaged RF loss. A steady thermal calculation is simulated in ANSYS. The temperature distribution on the internal surface is illustrated in Fig. 13. Fig. 13. Temperature distribution on the system and the mode convertor. 4.4 Mechanical concern: The cavity has two resonating modes and the total Q0 is the same as the individual Q0 because the two modes are identical. The loaded Q is about one tenth of the Q0 and QL is around 10 4. The Lorentz force on the cavity wall can be calculated locally with equation: surface P E H 0 2 2 The electric field acts as an attractive force and magnetic field acts a repulsively force. The combination pressure on surface causes the deformation of the spherical cavity and leads to frequency detuned. To avoid this, a mechanical enhancement is added to support the exterior. The force and deformation on the surface are shown in the Fig. 14.

Fig.14. Surface force and deformation with input power 1W. The detune frequency of spherical cavity is formulated in equation below: 0 Um U U 0 e Assumed the is 10 and the input power is 50MW, and that make the power loss P d 4 1 2 P f on the wall is 16MW at steady state, thus the stored energy U PQ d L is around at 2.2 Joule at peak charge. The typical coefficient to quantify this frequency tuning is Hz/(MV/m) 2. As stated above, when the ratio between max E surface field over input power is about 3500V/mW. Thus, the tuning coefficient is Hz/W. We need a tuning range of MHz to operate 50MW on this SLED. 4.5 Tuning capability and Multipacting: The cavity requires a tuning system to accurately adjust resonating frequency to 11.424GHz. As stated above, the frequency will be tuned with stored energy, and we need a frequency tuner with a frequency range around 20Mhz. A waveguide with a plunger can be an ideal tuner. One can use either a circular waveguide with propagating or evanescent TE11 mode to tune this spherical cavity; however both tuners have different frequency variation range and tuning sensitivity. The setups of plunger tuners are demonstrated in Fig.15. A center rod with a stem is introduced to cut off the TE11 coaxial mode in the pipe for both tuners. It is necessary to implement a choke joint structure to prevent leakage in this coaxial for the propagating waveguide plunger. The space between two plates in Fig15b, is this choke, and the distance in between is a quarter of waveguide wavelength in the choke joint.

Fig.15. Setups of two waveguide plungers tuners. The frequency shift and changes are functions of the plunger location, and these dependencies are plotted in fig. 16 for both tuners. Fig.16. The resonating frequency of TE114 and coupling coefficient vs the waveguide plunger location. The x axis is the distance between the plunger head to the spherical cavity. For the evanescent waveguide plunger, a larger waveguide area will reduce more resonating TE114 frequency. When the plunger distance is larger than the evanescent delay length, the tuner approaches the lower tuning frequency limit. For the propagating waveguide case, tuner frequency shift resonates with the distance of plunger, and the oscillating period is the guided wavelength in the circular waveguide. Multipacting simulation is under further investigation for both plunger tuners.

5. Conclusion A novel compact SLED is introduced and designed. The RF and multi-physics simulation suggests that it can operate with input power of 50MW with pulse length of 1.5microseconds and output power of 150MW with pulse length 100ns. A prototype of this SLED is about to make and tested soon. Acknowledgement We would like to thank Dr David Farkas from SLAC national accelerating laboratory for some useful discussions. This work was supported by Department of Energy Contract No. DE-AC02-76SF00515. Reference 1. H. Padamsee, J. Knobloch, and T. Hays, RF Superconductivity for Accelerators. 2nd Edition (Wiley and Sons, New York,NY, 2008). 2. W.-D. Moeller for the TESLA Collaboration, High Power Coupler For The TESLA Test Facility, Proceedings of the 9th Workshop on the RF Superconductivity, 1999, Santa Fe, V.2, pp.577-581. 3. J. R. Delayen, L.R. Doolittle, T. Hiatt, J. Hogan, J. Mammosser. An R.F. Input Coupler System For The CEBAF Energy Upgrade Cryomodule. Proceedings of the 1999 Particle Accelerator Conference, New York, 1999. pp1462-1464. 4. S. Belomestnykh, et al., High Average Power Fundamental Input Couplers for the Cornell University ERL: Requirements, Design Challenges and First Ideas, Cornell LEPP Report ERL 02-8 (September 9, 2002). 5. RF windows 6. E. Snitzer. Cylindrical dielectric waveguide modes Journal of the Optical Society of America, Vol. 51, Issue 5, pp. 491-498 (1961) 7. Zaki, K.A. Atia, A.E. Modes in Dielectric-Loaded Waveguides and Resonators. Microwave Theory and Techniques, IEEE Transactions on V.31, I.12,1982. pp1039 1045.

8. Clarricoats, P.J.B. Properties of dielectric-rod junctions in circular waveguide, Electrical Engineers, Proceedings of the Institution of (Volume:111, Issue: 1 ). 1964. Pp.43 50 9. Clarricoats, P.J.B. ; Taylor, B.C. Evanescent and propagating modes of dielectric-loaded circular waveguide. Electrical Engineers, Proceedings of the Institution of (Volume:111, Issue: 12 ) 1964 pp.1951-1956. 10. Walter M. Elsasser. Attenuation in a Dielectric Circular Rod. J. Appl. Phys. 20, 1193 (1949) 11. Rothwell, E.J. and Frasch, L.L. Propagation characteristics of dielectric-rodloaded waveguides.microwave Theory and Techniques, IEEE Transactions on (Volume:36, Issue: 3 ) 1988. pp 594 600. 12. Optics book 13. S.F.Mahmoud. Electromagnetics waveguides theory and applications. The Institution of Engineering and Technology (December 1991) 14. A.A. Mostafaa, C.M. Krowneb, K.A. Zakic & S. Tantawid Hybrid-Mode Fields in Isotropic and Anisotropic Planar Microstrip Structures. Journal of Electromagnetic Waves and Applications Volume 5, Issue 6, 1991 Appendix:

i01 i01 -Cos[ 01] - e -Cos[ 01] + e Sin[ 01] 0 0 0 2 2 2 i02 i 02 -Cos[ 02] - e -Cos[ 02] + e Sin[ 02] 0 0 0 2 2 2 i01 i01 -Cos[ 01] + e -Cos[ 01] - e Sin[ 01] 0 0 0 2 2 S 2 i02 i02 -Cos[ 02] + e -Cos[ 02] - e Sin[ 02] 0 0 0 2 2 2 Sin[ 01] Sin[ 01] 0 0 Cos[ 01] 0 2 2 Sin[ 02] Sin[ 02] 0 0 0 Cos[ 02] 2 2 Table 5

2.0 1.5 1.0 0.5 5650 5700 5750 5800 5850 5900 0.5 1.0 1.5 1.0 0.5 5000 6000 7000 8000