Spectral-domain approach (SDA) is applied to modeling propagation characteristics of multiconductor structure on superconducting strip lines with signal strips and ground planes of finite thickness in lossy layered media. Equivalent surface impedance is used to transform the superconducting strip of finite thickness into an infinitely thin strip. A novel procedure of a calculation of the equivalent surface impedance is proposed. Two sets of basis functions (Chebyshev and Legendre polynomials) are compared by their ability to accurate modeling of microwave losses in the structure. The model applicability is illustrated by numerical results compared with experimental and simulated data

Deleniv, A.

Gashinova, M.

Vendik, I. B.

English

Spectral-domain approach (SDA) is applied to modeling propagation characteristics of multiconductor structure on superconducting strip lines with signal strips and ground planes of finite thickness in lossy layered media. Equivalent surface impedance is used to transform the superconducting strip of finite thickness into an infinitely thin strip. A novel procedure of a calculation of the equivalent surface impedance is proposed. Two sets of basis functions (Chebyshev and Legendre polynomials) are compared by their ability to accurate modeling of microwave losses in the structure. The model applicability is illustrated by numerical results compared with experimental and simulated data.</p

Gashinova, Marina

University of Birmingham Research Portal

SDA full-wave analysis of boxed multistrip lines of finite thickness embedded in a layered lossy medium

AMS Acta

SDA FULL-WAVE ANALYSIS OF BOXED MULTISTRIP LINES OF FINITE THICKNESS EMBEDDED IN 
A LAYERED LOSSY MEDIUM 
 
A.Deleniv*, M. Gashinova**, and I. B. Vendik** 
 
* Chalmers University of Technology, Kemivägen 9, SE-412 96 Göteborg, Sweden, 
Phone: +4631-772 3605, Fax: +4631-164513, Email:anatoli@ep.chalmers.se 
**Electrotechical University, St.Petersburg,  197376, Russia, Phone/Fax: +7(812)3460867, Email: mwlab@eltech.ru 
 
ABSTRACT 
 
Spectral-domain approach (SDA) is applied to modeling propagation characteristics of multiconductor structure on 
superconducting strip lines with signal strips and ground planes of finite thickness in lossy layered media. Equivalent 
surface impedance is used to transform the superconducting strip of finite thickness into an infinitely thin strip. A novel 
procedure of a calculation of the equivalent surface impedance is proposed. Two sets of basis functions (Chebyshev and 
Legendre polynomials) are compared by their ability to accurate modeling of microwave losses in the structure. The 
model applicability is illustrated by numerical results compared with experimental and simulated data. 
 
THEORETICAL ANALYSIS 
 
A CAD tool for design of filters based on sections of multicoupled planar transmission lines Vendik  et al. (1) is 
developed and improved. An efficient full-wave method is derived to calculate complex propagation parameters of the 
boxed multi-conductor planar lines embedded in a layered isotropic media as shown in Fig. 1. The method can be 
applied to the analysis of superconducting and normal metal lines. 
 
The dielectric layers are characterized by the thickness hiUp ,hjLow and the complex dielectric permittivity 
 
 ( ) ( ) ωσεε jiji)j(i )Low(Up j−= .     (1) 
  
For a superconducting layer the dielectric constant is determined as  
 
ωσε )j(isc
)j(i
)Low(Up j−= ,     (2) 
 
where σsc  is the complex conductivity of a superconductor, Vendik (2). 
 
The goal of the analysis is to define a vector of the complex propagation constants ( ) ( ) ( )ωαωβω jk z −=  in the wave 
propagation direction z and the characteristic impedance for each mode propagating in the structure. 
 
With the commonly used spectral domain notation, the linear relationships between fields and currents on the 
conductor interface can be written as 
 
( )
( )
( )
( )










=





nx
nz
xxxz
zxzz
nx
eq
s
nz
eq
s
J
J
ZZ
ZZ
JZ
JZ
α
α
α
α
~
~
~~
~~
~
~
,    (3) 
where eqsZ  is the equivalent surface impedance, ijZ
~
 is the component of Green’s dyad, jiJ ,
~
 is the component of the 
surface current density, αn is the Fourier variable. We use the analytical procedure for evaluation of the spectral Green’s 
dyad based on so-called immitance approach Itoh (3). An extension to the lossy case is made by use of complex 
permittivity ( )jiε  for i-th (j-th) dielectric layer and the complex propagation constant zk . Contribution of the 
superconducting strip conductivity to the loss factor can be accurately estimated with a modified set of boundary 
conditions applied to an infinitely thin strip. The equivalent surface impedance is used to transform a superconducting 
strip of finite thickness into an infinitely thin strip (Fig. 2). The real superconducting strip is characterized by the 
tangential electric and magnetic field components on the upper and lower interfaces related by impedance matrix of a 
thin conducting layer Van den Berghe et al.(4): 
 


















−
×=





2
1
2221
1211
2
1 ~
~
10
01
~
~
H
H
ZZ
ZZ
n
E
E r     (4) 
)csc(
)cot(
2112
2211
tjZZZ
tjZZZ
scsc
scsc
β
β
==
−== ,    (5) 
 
where Zsc is the surface impedance and βsc is the complex propagation constant of the superconductor. 
We  define the electric field in the plane 0≤y≤t of the original strip as 
 
21
~
)sinh(
))(sinh(~
)sinh(
)sinh()(~ E
t
tyE
t
yyE
sc
sc
sc
sc
β
β
β
β −−=     (6) 
 
and use the following approximate ratio of the magnetic fields on the upper and lower interfaces for lossless dielectric 
media 
 
1
1
2
1
~
~
Low
Up
H
H
ε
ε
= .      (7) 
 
Keeping the equality between the power stored and dissipated within the volume of a real superconducting strip and the 
equivalent infinitely thin one, the following equation for the eqivalent surface impedance of infinitely thin strip was 
obtained 
 
( )( )
( )
( )
( )( )
( ) dykt
tyk
t
yZ
h
sc
scsc
sc
scsc
UpLowsc
eq
s
2
0
211
2
sinh
sinh1
sinh
sinh
1
1
∫
−−
+
=
β
σβ
β
σβ
εεσ
,  (8) 
 
where  
 
,1 1
1
1211
Up
LowZZk
ε
ε
−=  121
1
112 ZZk
Up
Low −=
ε
ε .    (9) 
 
In performing solution for planar transmission lines by SDA, we consider two different sets of basis functions ϕ1 and 
ϕ2 to expand the unknown electric currents of the system (3). The set ϕ1 is the first type of Chebyshev polynomials Tp 
with Maxwell weighting function. The set ϕ2 is based on orthogonal Legendre polynomials Pp and is suitable for 
modeling finite edge current. In our codes, we use these sets to model both the axial current density and the derivative 
of the transversal one 
 
).()('    ),()( xbxJxaxJ
k p
k
p
k
px
k p
k
p
k
pz ∑∑∑∑ == ϕϕ     (10) 
 
The system of linear equations (3) has to be arranged to a homogeneous one. We formulate two different systems to be 
solved using the sets ϕ1 and ϕ2: 
 




















−
−
−
−
−
=





)('~
)(~
~~
~
~
)(~
)(~
2
nx
nz
n
eq
sxx
n
xz
n
zxeq
szz
nx
nz
J
J
ZZ
j
Z
j
Z
ZZ
V
E
α
α
αα
α
α
α     (11) 
 
for ϕ1 and 




















−
−
−
−
−
=




Π
)('~
)(~
~~
~~
)(~
)(~
2
32
nx
nz
n
eq
sxx
n
xz
n
zx
n
eq
szz
nx
nz
J
J
ZZ
j
Z
j
ZZZ
V
E
α
α
αα
αα
α
α ,    (12) 
 
for ϕ2, with 
∫
−
=
x
L
zx dxxExV ')'()( ,     (13) 
t
L    x0
y
w1
D
w2 wn
m1 mn
s1
εUp
I
εLow
I hLow
I
hLow
I
hUp
I
hUp
I
∫ ∫− −=Π
x
L
x
L zz
dxdxxExE '')'()( .    (14) 
    
The both systems can be presented in general form as 
[ ] ( )[ ][ ].)(~ ,,~)(~ nznn JkFE αωαα =      (15) 
 
Solving (15) by Galerkin’s method leads to a homogeneous system of linear equations 
 
[ ][ ] 0),( =XkA zω      (16) 
 
with [X] being unknown expansion coefficients of the current density. The entries of [A(ω,kz)] are numerical series , 
whose general term in (16) is  
 
[ ] [ ]*~),,(
~)(~
*)(~)(~
1,,
,,
j
qznn
n
k
p
n
k
p
n
n
k
p
jk
qp kFA ϕωααϕαϕαϕ µν
µ
ν ∑∑
=   (17) 
 
where µ,ν stand for x- or z- components of spectral dyad, k, j=1..n are the strip numbers, p,q are the numbers of basis 
function. 
In order to reduce numerical efforts for solving the dispersion equation, the effective speed-up algorithm in the case of 
using the basis function set ϕ1 was implemented Gano et al. (5), while for the set ϕ2 the large argument approximation 
for Legendre polynomial basis functions was used. 
 
NUMERICAL RESULTS. 
 
The convergence behavior of Galerkin's method for both sets of basis function was studied. The propagation constants 
versus number of applied basis function are shown in Fig. 3 for superconducting (YBCO) microstrip line. Evidently, the 
Chebyshev polynomial-based set demonstrates quite poor convergence (typically, more than 15 basis functions have to 
be applied to provide 1% accuracy). It makes this basis set computationally inefficient. The basis set based on the 
Legendre polynomials shows better behavior (6-7 basis functions are enough to obtain the result close to converged 
value). All data presented below are computed using 7 basis functions of this set. 
 
A comparison between results of calculation of the HTS multistrip structure characteristics by the method developed 
with the experimental characteristics was done for the multistrip structure of two coupled S-shape YBCO resonators on 
LaAlO3 substrates (Fig. 4). Simulated and measured transmission coefficients at two different temperatures (T = 20 K 
and 77 K) Vendik et al. (6) are presented in Fig. 5. A good agreement gives a reliable confirmation of correctness of the 
method developed. 
 
CONCLUSIONS 
 
A computationally efficient full-wave method of analysis of HTS multistrip multilayered structure is proposed and 
discussed. The method allows taking into account real characteristics of the dielectric layers and conducting strips of a 
finite thickness. Numerous examples of simulation of different structures revealed a high efficiency and reliability of 
the method. As result, the complex propagation constants and related modal impedances of the multiconductor structure 
containing superconducting (normal metal) strip and lossy multilayer dielectric can be calculated. The method can be 
effectively used for a simulation and design of planar filter structures including tunable devices using ferroelectric layer. 
 
 
~ 
~ 
~ ~ 
y 
x 
H 1 
H 2 
E = Zseq Js 
 
y 
x 
n 1 
n 2 
E 1 , H 1 
E 2 , H 2 
J1 
J2 
t 
  
 
Fig. 1. Transverse section of multistrip structure in a 
layered media 
             a)                                         b) 
Fig. 2. Equivalence between tangential boundary 
conditions for a strip of thickness t (a) and infinitely 
thin equivalent (b). 
20 mm 
20 mm 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Fig. 3. Propagation constant β and attenuation coefficient α 
versus  number of  basis functions for an YBCO (The London 
penetration depth λL = 0.33 µm, the normal conductivity σn = 
3.0⋅106 (Ohm⋅m)-1) microstrip line on LaAlO3 substrate (f = 
10 GHz). The line with w = 180 µm and t =0.4 µm was placed 
on the substrate with h1Low = 0.508 mm, ε = 23.5 and 
tanδ =  6.0⋅108. 
Fig. 4. The layout of two S-shape YBCO resonators 
on LaAlO3 substrate. The copper (σ = 7⋅107 
(Ohm⋅m)-1) box dimensions are L = 25 mm and  
1
Uph = 7mm. 
          a)         b) 
Fig. 5. The simulated and measured Vendik et al. (6) transmission characteristics of two coupled S-shape resonators 
(Fig. 4). The dielectric substrate of 0.515 mm thickness was used with parameters ε = 23.3,  tanδ =  2.0⋅10-6, the YBCO 
film of 0.45 µm thickness was used with characteristics  λL = 0.144 µm, σn = 5.25⋅106 (Ohm⋅m)-1 (a) and  λL = 0.282 
µm, σn = 12.4⋅106 (Ohm⋅m)-1(b). 
 
REFERENCES 
 
1. I. B. Vendik, O. G. Vendik, A. N. Deleniv, V. V. Kondratiev, M. N. Goubina, and D. V. Kholodniak, “Development of 
CAD tool for a design of microwave planar HTS filters”, 2000, IEEE Trans. Microwave Theory Tech., vol. 48, pp. 1247-
1255. 
2. I. Vendik, “Phenomenological model of the microwave surface impedance of high-Tc superconducting films”, 2000, 
Supercond. Sci. Technol., vol. 13, pp. 974-982. 
3. T. Itoh, “Spectral domain immitance approach for dispersion characteristics of generalized printed transmission lines”, 
1980, IEEE Trans. Microwave Theory Tech., vol. 28, pp. 733-736. 
4. S. Van den Berghe, F. Olyslager, and D. De Zutter, “Accurate modeling of thin conducting layers in FDTD”, 1998, IEEE  
Microwave and Guided Wave Letters, vol. 8, pp. 75-77. 
5. G. Cano, F. Medina, and M. Horno, On the efficient Implementation of SDA for boxed strip-like and slot-like structures”, 
1998, IEEE Trans. Microwave Theory Tech., vol. 46, pp. 1801-1806. 
6. I. B. Vendik, A. N. Deleniv, V. O. Sherman, A. A. Svishchev, V. V. Kondratiev, D. V. Kholodniak, A. V. Lapshin, 
P. N. Yudin, B.-C. Min, Y. H. Choi, and B. Oh, "Narrowband Y-Ba-Cu-O filter with quasi-elliptic characteristic", 2001, 
Appl. Supercond., vol. 11, pp. 477-480. 
811,92 
811,94 
811,96 
811,98 
812,00 
812,02 
812,04 
812,06 
40 
42 
44 
46 
48
50 
52 
0 2 4 6 8 10 12 14 
n 
Chebyshev 
Legendre 
-80
-70
-60
-50
-40
-30
-20
-10
1,77 1,78 1,79 1,8 1,81
S2
1,
 (d
B
)
Frequency, (GHz)
simulation
experiment
T=20K
-80
-70
-60
-50
-40
-30
-20
-10
1,77 1,78 1,79 1,8 1,81
T=77K
simulation
experiment
S2
1,
 (d
B
)
Frequency, (GHz)


Sda full-wave analysis of boxed multistrip lines of finite thickness embedded in a layered lossy medium

SDA FULL-WAVE ANALYSIS OF BOXED MULTISTRIP LINES OF FINITE THICKNESS EMBEDDED IN A LAYERED LOSSY MEDIUM  A.Deleniv*, M. Gashinova**, and I. B. Vendik**  * Chalmers University of Technology, Kemivägen 9, SE-412 96 Göteborg, Sweden, Phone: +4631-772 3605, Fax: +4631-164513, Email:anatoli@ep.chalmers.se **Electrotechical University, St.Petersburg,  197376, Russia, Phone/Fax: +7(812)3460867, Email: mwlab@eltech.ru  ABSTRACT  Spectral-domain approach (SDA) is applied to modeling propagation characteristics of multiconductor structure on superconducting strip lines with signal strips and ground planes of finite thickness in lossy layered media. Equivalent surface impedance is used to transform the superconducting strip of finite thickness into an infinitely thin strip. A novel procedure of a calculation of the equivalent surface impedance is proposed. Two sets of basis functions (Chebyshev and Legendre polynomials) are compared by their ability to accurate modeling of microwave losses in the structure. The model applicability is illustrated by numerical results compared with experimental and simulated data.  THEORETICAL ANALYSIS  A CAD tool for design of filters based on sections of multicoupled planar transmission lines Vendik  et al. (1) is developed and improved. An efficient full-wave method is derived to calculate complex propagation parameters of the boxed multi-conductor planar lines embedded in a layered isotropic media as shown in Fig. 1. The method can be applied to the analysis of superconducting and normal metal lines.  The dielectric layers are characterized by the thickness hiUp ,hjLow and the complex dielectric permittivity   ( ) ( ) ωσεε jiji)j(i )Low(Up j−= .     (1)   For a superconducting layer the dielectric constant is determined as   ωσε )j(isc)j(i)Low(Up j−= ,     (2)  where σsc  is the complex conductivity of a superconductor, Vendik (2).  The goal of the analysis is to define a vector of the complex propagation constants ( ) ( ) ( )ωαωβω jk z −=  in the wave propagation direction z and the characteristic impedance for each mode propagating in the structure.  With the commonly used spectral domain notation, the linear relationships between fields and currents on the conductor interface can be written as  ( )( )( )( )=nxnzxxxzzxzznxeqsnzeqsJJZZZZJZJZαααα~~~~~~~~,    (3) where eqsZ  is the equivalent surface impedance, ijZ~ is the component of Green’s dyad, jiJ ,~ is the component of the surface current density, αn is the Fourier variable. We use the analytical procedure for evaluation of the spectral Green’s dyad based on so-called immitance approach Itoh (3). An extension to the lossy case is made by use of complex permittivity ( )jiε  for i-th (j-th) dielectric layer and the complex propagation constant zk . Contribution of the superconducting strip conductivity to the loss factor can be accurately estimated with a modified set of boundary conditions applied to an infinitely thin strip. The equivalent surface impedance is used to transform a superconducting strip of finite thickness into an infinitely thin strip (Fig. 2). The real superconducting strip is characterized by the tangential electric and magnetic field components on the upper and lower interfaces related by impedance matrix of a thin conducting layer Van den Berghe et al.(4):  −×=212221121121 ~~1001~~HHZZZZnEE r     (4) )csc()cot(21122211tjZZZtjZZZscscscscββ==−== ,    (5)  where Zsc is the surface impedance and βsc is the complex propagation constant of the superconductor. We  define the electric field in the plane 0≤y≤t of the original strip as  21~)sinh())(sinh(~)sinh()sinh()(~ EttyEtyyEscscscscββββ −−=     (6)  and use the following approximate ratio of the magnetic fields on the upper and lower interfaces for lossless dielectric media  1121~~LowUpHHεε= .      (7)  Keeping the equality between the power stored and dissipated within the volume of a real superconducting strip and the equivalent infinitely thin one, the following equation for the eqivalent surface impedance of infinitely thin strip was obtained  ( )( )( )( )( )( )( ) dykttyktyZhscscscscscscUpLowsceqs202112sinhsinh1sinhsinh11 ∫ −−+= βσββσβεεσ,  (8)  where   ,1 111211UpLowZZkεε−=  1211112 ZZkUpLow−=εε .    (9)  In performing solution for planar transmission lines by SDA, we consider two different sets of basis functions ϕ1 and ϕ2 to expand the unknown electric currents of the system (3). The set ϕ1 is the first type of Chebyshev polynomials Tp with Maxwell weighting function. The set ϕ2 is based on orthogonal Legendre polynomials Pp and is suitable for modeling finite edge current. In our codes, we use these sets to model both the axial current density and the derivative of the transversal one  ).()('    ),()( xbxJxaxJk pkpkpxk pkpkpz ∑∑∑∑ == ϕϕ     (10)  The system of linear equations (3) has to be arranged to a homogeneous one. We formulate two different systems to be solved using the sets ϕ1 and ϕ2:  −−−−−=)('~)(~~~~~)(~)(~2nxnzneqsxxnxznzxeqszznxnzJJZZjZjZZZVEααααααα     (11)  for ϕ1 and −−−−−=Π)('~)(~~~~~)(~)(~232nxnzneqsxxnxznzxneqszznxnzJJZZjZjZZZVEαααααααα ,    (12)  for ϕ2, with ∫−=xLzx dxxExV ')'()( ,     (13) tL    x0yw1Dw2 wnm1 mns1εUpIεLowI hLowIhLowIhUpIhUpI∫ ∫− −=ΠxLxL zzdxdxxExE '')'()( .    (14)     The both systems can be presented in general form as [ ] ( )[ ][ ].)(~ ,,~)(~ nznn JkFE αωαα =      (15)  Solving (15) by Galerkin’s method leads to a homogeneous system of linear equations  [ ][ ] 0),( =XkA zω      (16)  with [X] being unknown expansion coefficients of the current density. The entries of [A(ω,kz)] are numerical series , whose general term in (16) is   [ ] [ ]*~),,(~)(~*)(~)(~1,,,,jqznnnkpnkpnnkpjkqp kFA ϕωααϕαϕαϕ µνµν ∑∑=  (17)  where µ,ν stand for x- or z- components of spectral dyad, k, j=1..n are the strip numbers, p,q are the numbers of basis function. In order to reduce numerical efforts for solving the dispersion equation, the effective speed-up algorithm in the case of using the basis function set ϕ1 was implemented Gano et al. (5), while for the set ϕ2 the large argument approximation for Legendre polynomial basis functions was used.  NUMERICAL RESULTS.  The convergence behavior of Galerkin's method for both sets of basis function was studied. The propagation constants versus number of applied basis function are shown in Fig. 3 for superconducting (YBCO) microstrip line. Evidently, the Chebyshev polynomial-based set demonstrates quite poor convergence (typically, more than 15 basis functions have to be applied to provide 1% accuracy). It makes this basis set computationally inefficient. The basis set based on the Legendre polynomials shows better behavior (6-7 basis functions are enough to obtain the result close to converged value). All data presented below are computed using 7 basis functions of this set.  A comparison between results of calculation of the HTS multistrip structure characteristics by the method developed with the experimental characteristics was done for the multistrip structure of two coupled S-shape YBCO resonators on LaAlO3 substrates (Fig. 4). Simulated and measured transmission coefficients at two different temperatures (T = 20 K and 77 K) Vendik et al. (6) are presented in Fig. 5. A good agreement gives a reliable confirmation of correctness of the method developed.  CONCLUSIONS  A computationally efficient full-wave method of analysis of HTS multistrip multilayered structure is proposed and discussed. The method allows taking into account real characteristics of the dielectric layers and conducting strips of a finite thickness. Numerous examples of simulation of different structures revealed a high efficiency and reliability of the method. As result, the complex propagation constants and related modal impedances of the multiconductor structure containing superconducting (normal metal) strip and lossy multilayer dielectric can be calculated. The method can be effectively used for a simulation and design of planar filter structures including tunable devices using ferroelectric layer.   ~ ~ ~ ~ y x H 1 H 2 E = Zseq Js  y x n 1 n 2 E 1 , H 1 E 2 , H 2 J1 J2 t    Fig. 1. Transverse section of multistrip structure in a layered media              a)                                         b) Fig. 2. Equivalence between tangential boundary conditions for a strip of thickness t (a) and infinitely thin equivalent (b). 20 mm 20 mm                   Fig. 3. Propagation constant β and attenuation coefficient α versus  number of  basis functions for an YBCO (The London penetration depth λL = 0.33 µm, the normal conductivity σn = 3.0⋅106 (Ohm⋅m)-1) microstrip line on LaAlO3 substrate (f = 10 GHz). The line with w = 180 µm and t =0.4 µm was placed on the substrate with h1Low = 0.508 mm, ε = 23.5 and tanδ =  6.0⋅108. Fig. 4. The layout of two S-shape YBCO resonators on LaAlO3 substrate. The copper (σ = 7⋅107 (Ohm⋅m)-1) box dimensions are L = 25 mm and  1Uph = 7mm.           a)         b) Fig. 5. The simulated and measured Vendik et al. (6) transmission characteristics of two coupled S-shape resonators (Fig. 4). The dielectric substrate of 0.515 mm thickness was used with parameters ε = 23.3,  tanδ =  2.0⋅10-6, the YBCO film of 0.45 µm thickness was used with characteristics  λL = 0.144 µm, σn = 5.25⋅106 (Ohm⋅m)-1 (a) and  λL = 0.282 µm, σn = 12.4⋅106 (Ohm⋅m)-1(b).  REFERENCES  1. I. B. Vendik, O. G. Vendik, A. N. Deleniv, V. V. Kondratiev, M. N. Goubina, and D. V. Kholodniak, “Development of CAD tool for a design of microwave planar HTS filters”, 2000, IEEE Trans. Microwave Theory Tech., vol. 48, pp. 1247-1255. 2. I. Vendik, “Phenomenological model of the microwave surface impedance of high-Tc superconducting films”, 2000, Supercond. Sci. Technol., vol. 13, pp. 974-982. 3. T. Itoh, “Spectral domain immitance approach for dispersion characteristics of generalized printed transmission lines”, 1980, IEEE Trans. Microwave Theory Tech., vol. 28, pp. 733-736. 4. S. Van den Berghe, F. Olyslager, and D. De Zutter, “Accurate modeling of thin conducting layers in FDTD”, 1998, IEEE  Microwave and Guided Wave Letters, vol. 8, pp. 75-77. 5. G. Cano, F. Medina, and M. Horno, On the efficient Implementation of SDA for boxed strip-like and slot-like structures”, 1998, IEEE Trans. Microwave Theory Tech., vol. 46, pp. 1801-1806. 6. I. B. Vendik, A. N. Deleniv, V. O. Sherman, A. A. Svishchev, V. V. Kondratiev, D. V. Kholodniak, A. V. Lapshin, P. N. Yudin, B.-C. Min, Y. H. Choi, and B. Oh, "Narrowband Y-Ba-Cu-O filter with quasi-elliptic characteristic", 2001, Appl. Supercond., vol. 11, pp. 477-480. 811,92 811,94 811,96 811,98 812,00 812,02 812,04 812,06 40 42 44 46 4850 52 0 2 4 6 8 10 12 14 n Chebyshev Legendre -80-70-60-50-40-30-20-101,77 1,78 1,79 1,8 1,81S21, (dB)Frequency, (GHz)simulationexperimentT=20K-80-70-60-50-40-30-20-101,77 1,78 1,79 1,8 1,81T=77KsimulationexperimentS21, (dB)Frequency, (GHz)

AMS Acta - Alm@DL - Università di Bologna

SDA FULL-WAVE ANALYSIS OF BOXED MULTISTRIP LINES OF FINITE THICKNESS EMBEDDED IN 
A LAYERED LOSSY MEDIUM 
 
A.Deleniv*, M. Gashinova**, and I. B. Vendik** 
 
* Chalmers University of Technology, Kemivägen 9, SE-412 96 Göteborg, Sweden, 
Phone: +4631-772 3605, Fax: +4631-164513, Email:anatoli@ep.chalmers.se 
**Electrotechical University, St.Petersburg,  197376, Russia, Phone/Fax: +7(812)3460867, Email: mwlab@eltech.ru 
 
ABSTRACT 
 
Spectral-domain approach (SDA) is applied to modeling propagation characteristics of multiconductor structure on 
superconducting strip lines with signal strips and ground planes of finite thickness in lossy layered media. Equivalent 
surface impedance is used to transform the superconducting strip of finite thickness into an infinitely thin strip. A novel 
procedure of a calculation of the equivalent surface impedance is proposed. Two sets of basis functions (Chebyshev and 
Legendre polynomials) are compared by their ability to accurate modeling of microwave losses in the structure. The 
model applicability is illustrated by numerical results compared with experimental and simulated data. 
 
THEORETICAL ANALYSIS 
 
A CAD tool for design of filters based on sections of multicoupled planar transmission lines Vendik  et al. (1) is 
developed and improved. An efficient full-wave method is derived to calculate complex propagation parameters of the 
boxed multi-conductor planar lines embedded in a layered isotropic media as shown in Fig. 1. The method can be 
applied to the analysis of superconducting and normal metal lines. 
 
The dielectric layers are characterized by the thickness hiUp ,hjLow and the complex dielectric permittivity 
 
 ( ) ( ) ωσεε jiji)j(i )Low(Up j−= .     (1) 
  
For a superconducting layer the dielectric constant is determined as  
 
ωσε )j(isc
)j(i
)Low(Up j−= ,     (2) 
 
where σsc  is the complex conductivity of a superconductor, Vendik (2). 
 
The goal of the analysis is to define a vector of the complex propagation constants ( ) ( ) ( )ωαωβω jk z −=  in the wave 
propagation direction z and the characteristic impedance for each mode propagating in the structure. 
 
With the commonly used spectral domain notation, the linear relationships between fields and currents on the 
conductor interface can be written as 
 
( )
( )
( )
( )




=


nx
nz
xxxz
zxzz
nx
eq
s
nz
eq
s
J
J
ZZ
ZZ
JZ
JZ
α
α
α
α
~
~
~~
~~
~
~
,    (3) 
where eqsZ  is the equivalent surface impedance, ijZ
~
 is the component of Green’s dyad, jiJ ,
~
 is the component of the 
surface current density, αn is the Fourier variable. We use the analytical procedure for evaluation of the spectral Green’s 
dyad based on so-called immitance approach Itoh (3). An extension to the lossy case is made by use of complex 
permittivity ( )jiε  for i-th (j-th) dielectric layer and the complex propagation constant zk . Contribution of the 
superconducting strip conductivity to the loss factor can be accurately estimated with a modified set of boundary 
conditions applied to an infinitely thin strip. The equivalent surface impedance is used to transform a superconducting 
strip of finite thickness into an infinitely thin strip (Fig. 2). The real superconducting strip is characterized by the 
tangential electric and magnetic field components on the upper and lower interfaces related by impedance matrix of a 
thin conducting layer Van den Berghe et al.(4): 
 








−
×=


2
1
2221
1211
2
1 ~
~
10
01
~
~
H
H
ZZ
ZZ
n
E
E r     (4) 
)csc(
)cot(
2112
2211
tjZZZ
tjZZZ
scsc
scsc
β
β
==
−== ,    (5) 
 
where Zsc is the surface impedance and βsc is the complex propagation constant of the superconductor. 
We  define the electric field in the plane 0≤y≤t of the original strip as 
 
21
~
)sinh(
))(sinh(~
)sinh(
)sinh()(~ E
t
tyE
t
yyE
sc
sc
sc
sc
β
β
β
β −
−=     (6) 
 
and use the following approximate ratio of the magnetic fields on the upper and lower interfaces for lossless dielectric 
media 
 
1
1
2
1
~
~
Low
Up
H
H
ε
ε
= .      (7) 
 
Keeping the equality between the power stored and dissipated within the volume of a real superconducting strip and the 
equivalent infinitely thin one, the following equation for the eqivalent surface impedance of infinitely thin strip was 
obtained 
 
( )( )
( )
( )
( )( )
( ) dykt
tyk
t
yZ
h
sc
scsc
sc
scsc
UpLowsc
eq
s
2
0
211
2
sinh
sinh1
sinh
sinh
1
1 ∫ −−
+
= β
σβ
β
σβ
εεσ
,  (8) 
 
where  
 
,1 1
1
1211
Up
LowZZk
ε
ε
−=  121
1
112 ZZk
Up
Low
−=
ε
ε .    (9) 
 
In performing solution for planar transmission lines by SDA, we consider two different sets of basis functions ϕ1 and 
ϕ2 to expand the unknown electric currents of the system (3). The set ϕ1 is the first type of Chebyshev polynomials Tp 
with Maxwell weighting function. The set ϕ2 is based on orthogonal Legendre polynomials Pp and is suitable for 
modeling finite edge current. In our codes, we use these sets to model both the axial current density and the derivative 
of the transversal one 
 
).()('    ),()( xbxJxaxJ
k p
k
p
k
px
k p
k
p
k
pz ∑∑∑∑ == ϕϕ     (10) 
 
The system of linear equations (3) has to be arranged to a homogeneous one. We formulate two different systems to be 
solved using the sets ϕ1 and ϕ2: 
 













−
−
−
−
−
=


)('~
)(~
~~
~
~
)(~
)(~
2
nx
nz
n
eq
sxx
n
xz
n
zxeq
szz
nx
nz
J
J
ZZ
j
Z
j
Z
ZZ
V
E
α
α
αα
α
α
α     (11) 
 
for ϕ1 and 













−
−
−
−
−
=

Π
)('~
)(~
~~
~~
)(~
)(~
2
32
nx
nz
n
eq
sxx
n
xz
n
zx
n
eq
szz
nx
nz
J
J
ZZ
j
Z
j
ZZZ
V
E
α
α
αα
αα
α
α ,    (12) 
 
for ϕ2, with 
∫
−
=
x
L
zx dxxExV ')'()( ,     (13) 
tL    x0
y
w1
D
w2 wn
m1 mn
s1
εUp
I
εLow
I hLow
I
hLow
I
hUp
I
hUp
I
∫ ∫
− −
=Π
x
L
x
L zz
dxdxxExE '')'()( .    (14) 
    
The both systems can be presented in general form as [ ] ( )[ ][ ].)(~ ,,~)(~ nznn JkFE αωαα =      (15) 
 
Solving (15) by Galerkin’s method leads to a homogeneous system of linear equations 
 
[ ][ ] 0),( =XkA zω      (16) 
 
with [X] being unknown expansion coefficients of the current density. The entries of [A(ω,kz)] are numerical series , 
whose general term in (16) is  
 
[ ] [ ]*~),,(~)(~*)(~)(~
1,,
,,
j
qznn
n
k
p
n
k
p
n
n
k
p
jk
qp kFA ϕωααϕ
αϕαϕ µν
µ
ν ∑∑=
  (17) 
 
where µ,ν stand for x- or z- components of spectral dyad, k, j=1..n are the strip numbers, p,q are the numbers of basis 
function. 
In order to reduce numerical efforts for solving the dispersion equation, the effective speed-up algorithm in the case of 
using the basis function set ϕ1 was implemented Gano et al. (5), while for the set ϕ2 the large argument approximation 
for Legendre polynomial basis functions was used. 
 
NUMERICAL RESULTS. 
 
The convergence behavior of Galerkin's method for both sets of basis function was studied. The propagation constants 
versus number of applied basis function are shown in Fig. 3 for superconducting (YBCO) microstrip line. Evidently, the 
Chebyshev polynomial-based set demonstrates quite poor convergence (typically, more than 15 basis functions have to 
be applied to provide 1% accuracy). It makes this basis set computationally inefficient. The basis set based on the 
Legendre polynomials shows better behavior (6-7 basis functions are enough to obtain the result close to converged 
value). All data presented below are computed using 7 basis functions of this set. 
 
A comparison between results of calculation of the HTS multistrip structure characteristics by the method developed 
with the experimental characteristics was done for the multistrip structure of two coupled S-shape YBCO resonators on 
LaAlO3 substrates (Fig. 4). Simulated and measured transmission coefficients at two different temperatures (T = 20 K 
and 77 K) Vendik et al. (6) are presented in Fig. 5. A good agreement gives a reliable confirmation of correctness of the 
method developed. 
 
CONCLUSIONS 
 
A computationally efficient full-wave method of analysis of HTS multistrip multilayered structure is proposed and 
discussed. The method allows taking into account real characteristics of the dielectric layers and conducting strips of a 
finite thickness. Numerous examples of simulation of different structures revealed a high efficiency and reliability of 
the method. As result, the complex propagation constants and related modal impedances of the multiconductor structure 
containing superconducting (normal metal) strip and lossy multilayer dielectric can be calculated. The method can be 
effectively used for a simulation and design of planar filter structures including tunable devices using ferroelectric layer. 
 
 
~ 
~ 
~ ~ 
y 
x 
H 1 
H 2 
E = Zseq Js 
 
y 
x 
n 1 
n 2 
E 1 , H 1 
E 2 , H 2 
J1 
J2 
t 
  
 
Fig. 1. Transverse section of multistrip structure in a 
layered media 
             a)                                         b) 
Fig. 2. Equivalence between tangential boundary 
conditions for a strip of thickness t (a) and infinitely 
thin equivalent (b). 
20 mm 
20 mm 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Fig. 3. Propagation constant β and attenuation coefficient α 
versus  number of  basis functions for an YBCO (The London 
penetration depth λL = 0.33 µm, the normal conductivity σn = 
3.0⋅106 (Ohm⋅m)-1) microstrip line on LaAlO3 substrate (f = 
10 GHz). The line with w = 180 µm and t =0.4 µm was placed 
on the substrate with h1Low = 0.508 mm, ε = 23.5 and 
tanδ =  6.0⋅108. 
Fig. 4. The layout of two S-shape YBCO resonators 
on LaAlO3 substrate. The copper (σ = 7⋅107 
(Ohm⋅m)-1) box dimensions are L = 25 mm and  
1
Uph = 7mm. 
          a)         b) 
Fig. 5. The simulated and measured Vendik et al. (6) transmission characteristics of two coupled S-shape resonators 
(Fig. 4). The dielectric substrate of 0.515 mm thickness was used with parameters ε = 23.3,  tanδ =  2.0⋅10-6, the YBCO 
film of 0.45 µm thickness was used with characteristics  λL = 0.144 µm, σn = 5.25⋅106 (Ohm⋅m)-1 (a) and  λL = 0.282 
µm, σn = 12.4⋅106 (Ohm⋅m)-1(b). 
 
REFERENCES 
 
1. I. B. Vendik, O. G. Vendik, A. N. Deleniv, V. V. Kondratiev, M. N. Goubina, and D. V. Kholodniak, “Development of 
CAD tool for a design of microwave planar HTS filters”, 2000, IEEE Trans. Microwave Theory Tech., vol. 48, pp. 1247-
1255. 
2. I. Vendik, “Phenomenological model of the microwave surface impedance of high-Tc superconducting films”, 2000, 
Supercond. Sci. Technol., vol. 13, pp. 974-982. 
3. T. Itoh, “Spectral domain immitance approach for dispersion characteristics of generalized printed transmission lines”, 
1980, IEEE Trans. Microwave Theory Tech., vol. 28, pp. 733-736. 
4. S. Van den Berghe, F. Olyslager, and D. De Zutter, “Accurate modeling of thin conducting layers in FDTD”, 1998, IEEE  
Microwave and Guided Wave Letters, vol. 8, pp. 75-77. 
5. G. Cano, F. Medina, and M. Horno, On the efficient Implementation of SDA for boxed strip-like and slot-like structures”, 
1998, IEEE Trans. Microwave Theory Tech., vol. 46, pp. 1801-1806. 
6. I. B. Vendik, A. N. Deleniv, V. O. Sherman, A. A. Svishchev, V. V. Kondratiev, D. V. Kholodniak, A. V. Lapshin, 
P. N. Yudin, B.-C. Min, Y. H. Choi, and B. Oh, "Narrowband Y-Ba-Cu-O filter with quasi-elliptic characteristic", 2001, 
Appl. Supercond., vol. 11, pp. 477-480. 
811,92 
811,94 
811,96 
811,98 
812,00 
812,02 
812,04 
812,06 
40 
42 
44 
46 
48
50 
52 
0 2 4 6 8 10 12 14 
n 
Chebyshev 
Legendre 
-80
-70
-60
-50
-40
-30
-20
-10
1,77 1,78 1,79 1,8 1,81
S2
1,
 (d
B
)
Frequency, (GHz)
simulation
experiment
T=20K
-80
-70
-60
-50
-40
-30
-20
-10
1,77 1,78 1,79 1,8 1,81
T=77K
simulation
experiment
S2
1,
 (d
B
)
Frequency, (GHz)


http://amsacta.unibo.it/413/1/Eug_p134.pdf

Sda full-wave analysis of boxed multistrip lines of finite thickness embedded in a layered lossy medium

Abstract

Similar works

Full text

Available Versions

University of Birmingham Research Portal

AMS Acta

University of Birmingham Research Portal

AMS Acta

AMS Acta - Alm@DL - Università di Bologna