# Analysis of a Substrate Integrated Waveguide Hybrid Ring (Rat-Race) Coupler

Authors: Ramin Dehdasht-Heydari; Keyvan Forooraghi; Mohammad Naser-Moghadasi
Journal: International Journal of Antennas and Propagation
 (2014)

---

## Abstract

This paper presents an efficient analysis of a substrate integrated waveguide (SIW) single-layer hybrid ring coupler (rat-race) for millimeter-wave and microwave applications. The scattered field from each circular cylinder is expanded by cylindrical eigenfunctions with unknown coefficients that have been solved by electric and magnetic tangential boundary on each metallic via. The coupler S-matrix is calculated by using mode matching that uses the cylindrical vector expansion analysis to minimize the computational time and provides more physical insight. To achieve higher bandwidth, the radiuses of the coupler under analysis have been optimized in Matlab code by invasive weed optimization (IWO) method, and the results have been verified by CST package. The return loss and the isolation are less than −15 dB, and −18 dB, respectively. The insertion loss is divided equally-
3
±
0.2 dB, with 0
±
5 and 180
±
10 degrees in output ports over the operating frequency bandwidth and the agreement of phase differences in output ports has been examined objectively by feature selective validation (FSV) technique.


---

## Body

## 1. Introduction

The concept of substrate integrated waveguides was proposed [1] to overcome the drawbacks of conventional rectangular waveguide and has drawn much attention, for it can offer a promising solution to compact, low-cost, efficient, and hybrid integration of different kinds of conventional microstrip structures [2] in the same substrate at microwave and millimeter wave frequencies. In developing SIW base devices, many passive components, such as transitions, filters, oscillators, antennas, circulators, and couplers, have been studied in microwave and millimeter-wave systems [3–10].The analysis of substrate integrated waveguide (SIW) structures has been carried out in many methods. For instance, the method of moments accompanied with the Floquet mode expansions [11], the finite-difference frequency domain (FDFD) method [12], the boundary integral-resonant expansion (BI-RME) method [13] and the method of lines [14] have been applied to determine the propagation characteristics of the substrate integrated waveguides (SIWs). Moreover, the design of a real device mostly depends on the finite-element method and the method of moment, which can be time consuming and memory demanding for large structures.Recently, several papers have made use of a cylindrical vector wave (CVW) expansion in order to study the simplified 2D case [15–17] and the full 3D case [18–24]. These mode expansions allow an efficient full-wave analysis of SIW structures with metallic and dielectric posts. Furthermore, there is a need for using efficient methods of mode matching and CVW to analyze some specific structures such as couplers. The design of the SIW rat race coupler basically has been presented in [25] with transition from microstrip to SIW rat-race arms. However, in [25], the output phases results (i.e. phase differences of 0 and 180 degrees) of the coupler were not mentioned at all. In addition, the dimensions of the coupler should be optimized to achieve the more relative bandwidth. Because the analysis of the coupler has not been examined, its optimization has not been carried out properly in this case.In this paper, the dyadic Green’s function in a cylindrical medium can be represented in terms of cylindrical vector wave eigenfunctions. So, at first we find the vector wave function for each PEC via in the SIW rat-race coupler structure which is excited by coaxial cables and shown in Figure1. Because of using four coaxial cables instead of microstrip line transitions and a closed form of arms (λ
/
4 distance from the coaxial cables), our structure in this paper is of low loss and profile.(a) Top view of the proposed rat-race coupler under analysis. (b) One via of rat-race coupler and its related waves.
(a)
(b)When these vector wave expansions are known, the dyadic Green’s function for our rat-race coupler can be expressed in terms of these vector wave functions. By imposing boundary conditions on each via and transferring the obtained equations to matrix, we find the unknown coefficient of CVW. In the end, we find admittance matrix from the coefficients and conveniently we can reach scattering parameters of the proposed rat-race coupler. The mode-matching method has been implemented in Matlab software to compare the results with CST package. Also, to achieve better characteristics of rat-race coupler, we use an optimization with IWO and the optimized results are compared. It can be observed that the relative bandwidth has been increased from 9% to 17% over 11.4–13.7 GHz. In addition, the agreement of the output phase results has been shown clearly with feature selective validation (FSV).


## 2. Analysis of Rat-Race Coupler with Cylindrical Vector Wave (CVW) Expansion Formulations

As the structure in Figure1(a) is a homogeneous medium, the electric and magnetic fields fulfill the following equations:(1a)
∇
×
∇
×
E
(
r
)
-
k
2
E
(
r
)
=
0
,
(1b)
∇
×
∇
×
H
(
r
)
-
k
2
H
(
r
)
=
0
,where k is the wave number in dielectric substrate and k
=
ω
μ
0
ε
0
ε
r.Also, in a homogenous medium,E
(
r
) and H
(
r
) are derivable from scalar potential ϑ
(
r
). Now, we can define two vector wave functions P
(
r
) and Q
(
r
) as the equivalence of an electric or a magnetic field. Notice that the above equation is valid for Q
(
r
):
(2)
∇
·
P
(
r
)
=
0
,
∇
×
P
(
r
)
≠
0
.For consideration of the relation betweenϑ
(
r
) and P
(
r
) and Q
(
r
), let us assume the following equation for scalar potential:
(3)
(
∇
2
+
k
2
)
ϑ
(
r
)
=
0
,
and the relation between ϑ
(
r
), P
(
r
), and Q
(
r
) is as follows:
(4)
P
(
r
)
=
∇
×
(
ϑ
(
r
)
a
n
)
,
(5)
Q
(
r
)
=
1
k
∇
×
P
(
r
)
,
where a
n is unit normal vector on the surface.In a cylindrical coordinate we can express in each via in region 1 (outgoing region) this harmonic scalar potentialϑ
n
(
r
):
(6)
ϑ
n
(
r
)
=
H
n
(
2
)
(
k
m
ρ
)
e
j
(
k
z
z
-
n
φ
)
,
and H
n
(
2
) indicates the Hankel function of the second kind, k
m
=
k
2
-
k
z
2, k
z
=
m
π
/
h, and h is the height of substrate in Figure 1(a). In (6), m and n are the modes in vertical direction (z-direction) and azimuthal direction (φ-direction), respectively. Also, in this paper we consider the dependence of the SIW structure on the variation to z, which is expressed as an exponential term.By using (6) in (4) and (5) we obtain P
(
r
) and Q
(
r
) for Figure 1(b):(7a)
P
n
(
r
)
=
(
-
1
ρ
j
n
e
j
(
k
z
z
-
n
φ
)
H
n
(
2
)
(
k
m
ρ
)
)
a
ρ
+
(
-
k
m
ρ
H
n
′
(
2
)
(
k
m
ρ
)
e
j
(
k
z
z
-
n
φ
)
)
a
φ
(7b)
Q
n
(
r
)
=
(
j
k
m
k
z
ρ
H
n
′
(
2
)
(
k
m
ρ
)
e
j
(
k
z
z
-
n
φ
)
)
a
ρ
+
(
1
ρ
n
k
z
e
j
(
k
z
z
-
n
φ
)
H
n
(
2
)
(
k
m
ρ
)
)
a
φ
+
(
k
m
e
j
(
k
z
z
-
n
φ
)
(
-
2
H
n
′
(
2
)
(
k
m
ρ
)
m
m
m
m
m
m
m
m
m
-
ρ
k
m
H
n
′′
(
2
)
(
k
m
ρ
)
)
)
a
z
,where
(8)
H
n
′
(
2
)
(
k
m
ρ
)
=
∂
∂
(
H
n
(
2
)
(
k
m
ρ
)
)
,
H
n
′′
(
2
)
(
k
m
ρ
)
=
∂
∂
(
H
n
′
(
2
)
(
k
m
ρ
)
)
.So, arbitrary electric and magnetic fieldsE
(
r
) and H
(
r
) in region 1 can be expressed as (9a)
H
total
(
r
)
=
H
1
+
H
2
,
(9b)
E
total
(
r
)
=
E
1
+
E
2
,where E
2 and H
2 are scattering electrical and magnetic waves from the via and E
1 and H
1 are incident waves from a magnetic source (M
s
′) in the region 1. To calculate H
1, we need this equation:
(10)
H
1
=
-
j
ω
ε
∭
G
m
2
·
M
s
′
d
V
′
,
where G
m
2 is the magnetic Green function of second kind [26, 27] which for parallel plate waveguide by using residue theorem in cylindrical coordination explicit form is given by
(11)
G
m
2
(
r
,
r
′
)
=
-
ρ
^
ρ
^
k
0
2
δ
(
r
-
r
′
)
-
1
j
4
π
∑
m
=
0
N
z
∑
n
=
0
N
φ
(
2
-
δ
m
0
)
1
k
m
2
h
×
[
P
n
(
k
m
,
k
z
,
ρ
)
P
n
′
(
k
m
,
k
z
,
ρ
′
)
m
m
+
Q
n
(
k
m
,
k
z
,
ρ
)
Q
n
′
(
k
m
,
k
z
,
ρ
′
)
]
,
(12)
δ
m
0
=
{
0
if
m
≠
0
1
if
m
=
0
,
where P
n and Q
n for ρ
>
ρ
′ are similar as outgoing wave in region 1 for (7a) and (7b). For ρ
<
ρ
′, the equations are replaced with (J
n
(
k
m
ρ
)) instead of (H
n
(
k
m
ρ
)). Also, P
n
′ and Q
n
′ can be achieved by interchanging ρ
<
ρ
′ with ρ
>
ρ
′ and vice versa. The m and n in (11) are eigenvalues modes in vertical and azimuthal directions, respectively.After obtainingH
1 from (10) the electrical field E
1 is obtained from
(13)
E
=
1
j
ω
ϵ
∇
×
H
.From (1a), (1b), and (3), we can write TE and TM modes for scattering waves (H
2, E
2) in region 1 as CVW expansions:
(14)
H
2
(
r
)
=
∑
l
=
1
P
via
∑
m
=
1
N
z
∑
n
=
1
N
φ
(
C
l
,
m
,
n
TM
P
n
+
D
l
,
m
,
n
TE
Q
n
)
.
By using (13),
(15)
E
2
(
r
)
=
∑
l
=
1
P
via
∑
m
=
1
N
z
∑
n
=
1
N
φ
(
C
l
,
m
,
n
TM
Q
n
+
D
l
,
m
,
n
TE
P
n
)
.To determine unknown coefficients(
C
l
,
m
,
n
TM
,
D
l
,
m
,
n
TE
), the boundary conditions for the surface of each metallic should be imposed:
(16)
n
^
×
(
E
1
+
E
2
)
=
0
in
surface
of
each
PEC
post
with
radius
a
i
.After imposing boundary conditions that span over the metal posts of sidewalls, arms, inner circle, inner conductor of coaxial cables, and four matching impedance metallic holes which are demonstrated in Figure1(a), we can obtain a matrix (for TE and TM modes separately) for implantation in Matlab code. By using (16) the unknown coefficient vector F
l
,
m
,
n
TM
/
TE
=
[
C
l
,
m
,
n
TM
,
D
l
,
m
,
n
TE
] can be expressed as
(17)
[
F
1
,
m
,
0
TM
/
TE
⋮
F
l
,
m
,
n
TM
/
TE
]
=
[
U
]
-
1
[
S
]
,
where U is interaction matrix between the vias and S is excitation vector.In (17), for propagation modes, we use the direct inverse command in Matlab to obtain unknown coefficients, and, for nonpropagation modes, we use iteration methods because of the singularity.The admittance formulation of four coaxial cables in Figure1(a) is
(18)
Y
(
p
i
,
p
j
)
=
-
∫
H
p
i
·
M
p
j
′
(
ρ
)
d
S
′
|
V
|
2
,
where M
′
(
ρ
) is surface current and has the same distribution as the E field in a coaxial cable:
(19)
E
′
(
ρ
)
=
V
0
2
ρ
L
n
(
b
/
a
)
a
ρ
,
and with Schelkunoff’s equivalence principle:
(20)
M
′
(
ρ
)
=
-
2
n
^
×
E
(
ρ
)
.With excitation of the coupler with coaxial cables, the scattering parameters of the coupler conveniently can be obtained from (18).


## 3. Optimization of the Rat-Race Coupler with IWO Algorithm

Invasive weed optimization (IWO) algorithm has been first used by Mehrabian and Lucas [28] in dynamic and control systems theory. The main parameters involved in improving the bandwidth of a SIW hybrid ring coupler with IWO are R
1, R
2, and R
3. IWO parameters which are used for optimization are shown in Table 1. In our optimization, we consider a goal function (goal: Max (S
11) ≤ −15 dB), and when we get to our goal after passing the necessary iterations, the other output results of the coupler will be satisfied. S
11 variations of the hybrid ring coupler in interest band versus number of iterations are shown in Figure 2. As shown in this figure, the optimized value of S
11 is −15 dB after 175 iterations number which corresponds to the R
1 = 4.5 mm, R
2 = 6.5 mm, and R
3 = 15.5 mm over the 11.4–13.7 GHz bandwidth.Table 1
IWO parameters values for Figure1(a).
Symbol
Quantity
Value
No
Number of initial population
5
itmax
Maximum number of iterations
1000
dim
Problem dimension
3
P
max
⁡
Maximum number of plant population
5
S
max
⁡
Maximum number of seeds
25
S
min
⁡
Minimum number of seeds
0
n
Nonlinear modulation index
3
σ
initial
Initial value of standard deviation
2
σ
final
Final value of standard deviation
0.01
R
1
ini
Initial search area
2 mm to 20 mm
R
2
ini
4 mm to 30 mm
R
3
ini
10 mm to 45 mmFigure 2
S
11 of the rat-race coupler versus IWO iteration number.


## 4. Feature Selective Validation (FSV)

Feature selective validation (FSV) is a very useful technique to compare the agreement of two plots usually to decide that the two plots have excellent, very good, good, fair, poor, or very poor agreement [29]. By using FSV, the comparison between results is to be measured objectively, eliminating the element of subjectivity from the decision making process.The FSV technique is based on the decomposition of the results to be compared into two component measures and then the recombination of the results to provide a global goodness of fit measure. The components used are three figures of merit of the comparison of two data sets [30–33]:(1)
the amplitude difference measure (ADM) that compares the amplitudes and trends of the two data sets,
(2)
the feature difference measure (FDM) that compares the rapidly changing features (as a function of the independent variable),
(3)
the global difference measure (GDM) that is obtained with combination of the ADM and FDM.Accordingly, the comparison of two data sets can be ranked. This ranking is useful for making a selection between multiple comparisons. So, two quality factors for each figure of merit (ADM, FDM, and GDM) are regarded (i.e., the GRADE and the SPREAD).The GRADE provides a numerical indication of the quality of the comparison and the SPREAD gives a numerical indication of the level of confidence that can be placed on this evaluation. In GRADE, the smaller it is, the better the comparison of the results is, and, in SPREAD, the higher the reliability is, the smaller it is. Furthermore, GRADE and SPREAD can be computed for each figure of merit and reported on a GRADE/SPREAD chart.When we use FSV tool [34], there are time domain and frequency domain analyses. A time domain analysis is performed on two sets of data that are generic amplitude values.In the case of frequency domain, there are three options to analyze:(1)
magnitude option that the input data are magnitude values,
(2)
phase option that the input data are phase values,
(3)
combined option in which the 1D FSV combined analysis performs two FSV analyses (i.e., magnitude and phase) and combines them in a single result.


## 5. Result and Discussion

The formulations in Section2 and IWO algorithm have been implemented in Matlab code. The S-parameters from analysis, after optimization of the structure in Section 3, are shown in Figure 3. The S
11 (return loss) and S
14 (isolation) over 11.4–13.7 GHz are less than −15 dB and −18 dB in the order mentioned and insertion loss is divided equally −3 ± 0.2 dB in ports 2 and 3. The analysis results show excellent agreement with the simulation. Figure 4 shows relative phases for the rat-race coupler. The phase difference in Figure 4(a) is 0 ± 5° and in Figure 4(b) is 180 ± 10° over the interest band. Also, there is good agreement between the two proposed methods. Obviously, the relative bandwidth after optimization from Figures 3 and 4 is approximately 17%.Figure 3
TheS-parameters versus frequency from analysis and CST.Phase responses of the SIW rat-race coupler: (a) in-phase response and (b) out-of-phase response.
(a)
(b)As an application of the FSV technique which has been explained in Section4, the agreement of the output phase results in Figure 4 has been performed with 1D FSV tool 2.0.4 [32] in frequency domain. According to [34], excellent, very good, good, fair, poor, and very poor agreements fall in the ranges of [00.1), [0.10.2), [0.20.4), [0.40.8), [0.81.6), and [1.6infinity), respectively. From Figure 5, it is observed that the FSV tool results show good and fair agreements for in-phase responses in the most range of the bandwidth of interest. Also, based on Figure 6, very good and excellent agreements have been obtained for out-of-phase responses in almost all of the operating frequency range.Agreement cases for in-phase responses with 1D FSV tool: (a) ADMi and FDMi, (b) ADMc and FDMc, (c) GDMi and GDMc, and (d) GRADE/SPREAD chart.
(a)
(b)
(c)
(d)Agreement cases for out-of-phase responses with 1D FSV tool: (a) ADMi and FDMi, (b) ADMc and FDMc, (c) GDMi and GDMc, and (d) GRADE/SPREAD chart.
(a)
(b)
(c)
(d)To prove the efficiency of our analysis, simulation time for the considered case and comparing it to CST software are reported in Table2. From this table, the run time ratio of the analysis related to CST is 1 : 10.9.Table 2
Simulation time on a Core i5 with 4 GB RAM (azimuthally mode numbers(
n
)
=
7; vertically mode numbers (
m
)
=
2).
Structure type in Figure1
Number of metallic vias
CST CPU time
Analysis in this paper
Time domain
50 frequency points
Hybrid ring (rat-race) coupler
330
1200 s
110 s


## 6. Conclusion

In this paper, the dyadic Green’s function and vector wave functions for cylindrical coordination have been used for a SIW rat-race coupler. Our output results with the mode matching method after optimization with IWO and comparing it with CST using FSV technique and CPU time show accuracy and efficiency of our theory, respectively. Throughout 11.4–13.7 GHz, the structure has 17% relative bandwidth. By using the proposed analysis and optimizing the dimensions of the coupler, the bandwidth can shift in the interested band.



---

*Source: 803515-2014-04-28.xml*