Electrodynamics and Insulating Materials B
Electrodynamics and Insulating Materials B
1 Introduction 1
2 Wave classification 5
2.1 Reduced Maxwell’s equations - A recap . . . . . . . . . . . 5
2.2 Application . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.2.1 TEM or Plane waves . . . . . . . . . . . . . . . . . . 8
2.2.2 TE (H) waves . . . . . . . . . . . . . . . . . . . . . . 11
2.2.3 TM (E) waves . . . . . . . . . . . . . . . . . . . . . 13
4 Waveguides 89
4.1 General wave characteristics along uniform guiding structures 89
4.1.1 TEM modes . . . . . . . . . . . . . . . . . . . . . . . 91
4.1.2 TM modes or E modes . . . . . . . . . . . . . . . . . 92
4.1.3 TE modes or H modes . . . . . . . . . . . . . . . . . 95
4.2 Wave velocities and dispersion in guiding structures . . . . 97
4.3 Group velocity . . . . . . . . . . . . . . . . . . . . . . . . . 99
This forms the second part of the course "Electrodynamics and Insulating
Materials" in which we will treat topics in
• Coaxial line
• Microstrip
• Coplanar line
Waveguides include
Fig. 1.1 shows a section of the electromagnetic spectrum and the desig-
nated application. Table 1.1 presents some traditional applications of the
low-frequency end of the electromagnetic spectrum.
Note that Channels 54 − 69 (734M Hz − 862 M Hz) have now been
switched off for analogue TV (European Union). The frequency band from
734 MHz−790 MHz is now allocated to digital terrestrial TV (DTT)/terrestrial
1. Introduction 3
790MHz 791MHz
FIGURE 1.2. Old and new frequency plan for the 800 M Hz band (European
Union)
digital video broadcasting (DVB-T) applications. The band from 791 MHz−
862 M Hz has been allocated for international mobile telecommunications.
An important point to note is that the 4G long term evolution (LTE) ap-
plication is now collocated with DTT/DVB-T applications and presents
potential interference problems. This new frequency plan is shown in Fig.
1.2.
An alternative designation of some of the microwave bands are presented
in Table 1.2.
4 1. Introduction
2
Wave classification
→ −
− →
H = H t + zHz exp (−γz)
6 2. Wave classification
where we drop the factor exp (jωt) and use the phasor formulation for
convenience.
Now consider Faraday’s induction law in a source-free region of space.
We can write
−
→ − → −→ − → −→
∇ × E = ∇ t − zγ × E t + zEz = −jωµ H t + zHz
or
−
→ −
→ −
→ −
→ −
→
∇ t × E t + ∇ t × zEz − γ
z × E t − γ
z × zEz = −jωµ H t − jωµ
z Hz
We see that this separates the fields into transverse and axial components
as,
−
→ −
→
∇ t × E t = −jωµ z Hz , axial field components
and
−
→ →
− −
→
z × ∇ t Ez + γ E t = jωµ H t , transverse field components
or
−
→ −
→ −
→ −
→ −
→
∇ t × H t + ∇ t × zHz − γ
z × H t − γ
z × zHz = jωǫ E t + jωǫ
z Ez
Thus,
−
→ −
→
∇ t × H t = jωǫ
z Ez , axial field components
and −
→ →
− −
→
z × ∇ t Hz + γ H t = −jωǫ E t , transverse components
or
−
→ − → −
→ −
→
∇ t · H t + ∇ t · zHz − γ
z · H t − γ
z · zHz = 0
Since
−
→ −
→
∇ t · zHz = z · H t = 0
then
−
→ − → −
→ − →
∇ t · H t − γ
z · zHz = 0 ⇒ ∇ t · H t = γHz
In a region of space containing free charges, we can write,
−
→ − → −
→ −→
∇ · D = ǫ ∇ t − γ z · E t + zEz = ρ
2.1 Reduced Maxwell’s equations - A recap 7
or
−
→ − → −
→ −
→ ρ
∇ t · E t + ∇ t · zEz − γ
z · E t − γ
z · zEz =
ǫ
or
→ −
− → ρ −
→ − → ρ
z · zEz = ⇒ ∇ t · E t = γEz +
∇ t · E t − γ
ǫ ǫ
When the region contains no free charges, we have
−
→ − →
∇ t · E t = γEz
Thus −
→ −
→
∇t × −E t = −jωµ z Hz
→ →
− −
→
z × ∇ t Ez + γ E t = jωµ H t
−
→ −
→
∇t × −H t = jωǫ
z Ez
→ −
→ −
→
z × ∇ t Hz + γ H t = −jωǫ E t
−
→ − →
∇ t · H t = γHz
−
→ − →
∇ t · E t = γEz
We can express the transverse fields as functions of the axial field com-
ponents alone. In the transverse curl equations, we write
−→ → −
− → →
− −
→
z × z × ∇ t Ez + γ E t = ∇ t Ez + γ E t = −jωµ z × Ht
and −
→ →
− −
→
−jωµ z × ∇ t Hz + γ H t = (−jωµ) (−jωǫ) E t
or →
− ω2 µǫ −
→ jωµ −
→
−jωµ z × H t = − Et + z × ∇ t Hz
γ γ
Thus, subtracting, we obtain
−
→ → ω 2 µǫ −
− → jωµ −
→
∇ t Ez + γ E t + Et − z × ∇ t Hz = 0
γ γ
or
ω2 µǫ −
→ −
→ jωµ −
→
γ+ E t = − ∇ t Ez + z × ∇ t Hz
γ γ
Hence,
−
→ γ− → jωµ −
→
E t = − 2 ∇ t Ez + 2 z × ∇ t Hz
kc kc
where
kc2 = γ 2 + ω 2 µǫ = γ 2 + k2
Similarly,
−
→ → −
− → −
→ −
→
z × z × ∇ t Hz + γ H t = ∇ t Hz + γ H t = jωǫ
z × Et
8 2. Wave classification
and −
→ →
− −
→
jωǫ z × ∇ t Ez + γ E t = (jωǫ) (jωµ) H t
or
−
→ ω 2 µǫ −
→ jωǫ −
→
jωǫ
z × Et = − z × H t − z × ∇ t Ez
γ γ
Thus, subtracting,
−
→ ω2 µǫ −→ jωǫ −
→
∇ t Hz + γ + Ht + z × ∇ t Ez = 0
γ γ
or
−
→ γ− → jωǫ −
→
H t = − 2 ∇ t Hz − 2 z × ∇ t Ez
kc kc
2.2 Application
2.2.1 TEM or Plane waves
−
→ −
→
For TEM waves, E and H lie in a plane transverse to the axis along which
the wave propagates. Thus, Ez = Hz = 0. Hence
−
→ −
→
∇t × E t = 0
−
→ −
→
z × E t = j ωµ
γ Ht
−
→ −
→
∇t × H t = 0
−
→ −
→
z × H t = −j ωǫγ Et
−
→ − →
∇t · H t = 0
−
→ − →
∇t · E t = 0
−
→ −
→ −
→ −
→
The fact that ∇ t × E t = 0, ∇ t × H t = 0 implies that the line integral
−
→ −
→
of either E t or H t along a closed path in the xy-plane is zero equivalent
to the case of static electromagnetic fields. This has a physical basis since
there is neither an axial magnetic nor electric field threading the contour.
We can thus define the transverse field components as gradients of scalar
potentials i.e.,
−
→ −
→ −
→ −
→
E t = − ∇ t Φe , H t = − ∇ t Φm
where Φe and Φm are analytical scalar potential functions.
Consider the transverse divergence of the transverse field components.
We obtain
−
→ − → −
→ − →
∇ t · E t = − ∇ t · ∇ t Φe = ∇2t Φe = 0
and
−
→ − → −
→ − →
∇ t · H t = − ∇ t · ∇ t Φm = ∇2t Φm = 0
2.2 Application 9
which are the two-dimensional Laplace’s equations for the transverse elec-
tric and magnetic fields in a sourceless region of space.
These fields must satisfy the Helmholtz equation, i.e.,
−
→ −
→
∇2 E t + k2 E t = 0
and
−
→ −
→
∇2 H t + k2 H t = 0
where k2 = ω2 µǫ. Thus, since
−→ −
→
∇2 = ∇ t − γ T EM z · ∇ t − γ T EM z = ∇2t + γ 2T EM
or
→ 2
−
∇ t ∇t Φe + k2 + γ 2T EM Φe = 0
Thus, for TEM waves,
√
γ T EM = jβ = ±jk = ±jω µǫ
−
→
We can then obtain H t from
−
→ jωµ − →
z × E t = Ht
γ T EM
or
−
→ −
→
−
→ γ T EMT EM −
→ z × E t z × ∇ t Φe
Ht = z × E t = =−
jωµ ηT EM ηT EM
where
jωµ µ
η T EM = = =η
γ ǫ
η is the intrinsic wave impedance of the unbounded medium. For a TEM
wave propagating in the +z-direction, we hence have,
−
→ −
→
E t = − ∇ t Φe
−
→ −
→
−
→ z × E t z × ∇ t Φe
Ht = =−
η η
ηT EM is called the wave impedance of the medium and is used in situations
where the wave is confined. For the TEM waves, we see that the wave
impedance is equal to the intrinsic impedance of the unbounded medium.
10 2. Wave classification
or
∂ 2 −
→ −
→
− ∇ t Φm = − (jωµ) (jωǫ) ∇ t Φm
∂z 2
or
∂2 ∂2
2
Φm = −ω 2 µǫΦm ⇒ 2 I (z) = −ω 2 µǫI (z) = −ω2 L′ C ′ I (z)
∂z ∂z
we can write
2
−
→ ∂Φm ω µǫ −
→ jωµ kc2 −→
∇t = 2
− 1 ∇ t Hz = − 2 jωǫ + ∇ t Hz
∂z kc kc jωµ
or
∂Φm k2 jωµ
= − jωǫ + c Hz
∂z jωµ kc2
We note that the quantity jωµHz /kc2 has dimensions of voltage which
we designate V (z) while Φm has those of current which we designate I (z)
and therefore express the equations as
∂
V (z) = −jωµI (z) = −Z ′ I (z)
∂z
and
∂Iz kc2
=− + jωǫ V (z) = −Y ′ V (z)
∂z jωµ
where
jωµ
V (z) = Hz
kc2
I (z) = Φm
Z ′ = jωµ
and
kc2
Y′ = + jωǫ
jωµ
are the per unit length series impedance and shunt admittance of a trans-
mission line respectively. The characteristic impedance of the line is hence,
Z′ µ 1 η
Z0,T E = = =
Y′ ǫ
1 − (kc /k)2 1 − (kc /k)2
and for these waves, γ T E = 0 at k = kc in which case the wave does not
propagate. kc is hence called the cut-off wavenumber of a TE mode.
The equivalent circuit for TE wave propagation is shown in Fig. 6.7.
2.2 Application 13
jωµ
k2c
jωε
jωµ
and
−
→ jωǫ −
→ jωǫ −
→
Ht = z × E t = − 2 z × ∇ t Ez
γT M kc
We note that
jωǫ k ǫ k
=j =j
γT M γT M µ γT M η
so that
−
→ k −
→ 1 −
→
Ht = j z × E t = j z × E t
γT M η ηT M
where
γT M
ηT M = η
k
is the wave impedance of the E-waves.
−
→
Because E t is curl-free, we can express it as the gradient of an analytical
scalar function Φe as
−
→ −
→ γ −
→ ∂ 1−→
E t = − ∇ t Φe = − T2M ∇ t Ez = ∇ E
t z
kc ∂z kc2
and from
−
→ −
→ −
→ jωǫ −
→ ω2 µǫ −
→
∇ t Ez +γ T M E t = −jωµ z × − 2 z × ∇ t Ez = 2 ∇ t Ez
z × H t = −jωµ
kc kc
we can write
−
→ ∂Φe ω 2 µǫ −
→ −
→
−γ T M ∇ t Φe = = 2 ∇ t Ez − ∇ t Ez
∂z kc
or
2
∂Φe ω µǫ 1 jωǫ
= 2
− 1 Ez = − jωµ + Ez
∂z kc jωǫ kc2
We note that the quantity jωǫEz represents an axial displacement current
density while 1/kc2 has dimensions of area. Thus jωǫEz /kc2 is a current in
the axial direction which we denote by an equivalent current I (z). Hence
∂ jωǫ ∂
Ez = I (z) = −jωǫV (z)
∂z kc2 ∂z
and
∂ kc2
V (z) = − jωµ + I (z)
∂z jωǫ
These have the same form as the differential equations for a lossless trans-
mission line with series impedance per unit length given by
kc2
Z ′ = jωµ +
jωǫ
and shunt admittance per unit length
Y ′ = jωǫ
Thus, we can represent the wave propagation characteristics of the TM
waves as shown in Fig. 6.6.We note that TE and TM waves are duals of
each other since
η T E ηT M = η2
as long as γ is the same for both waves.
For energy flow along the z-axis, power is defined in terms of the trans-
verse field components alone. Thus, for TE waves
1 −
→ −
→
Pave = Re E t × H ∗t · z dx dz
2
S
1 → −
− →
= Re ηT E z × H t × H ∗t · z dx dz
2
S
ηT E −
→ − → 1 −
→ − →
= H t · H ∗t dx dz = E t · E ∗t dx dz
2 2ηT E
S S
2.2 Application 15
k2
c
jωµ jωε
jωε
−
→ −
→
• The coaxial transmission line in which the E and H fields are com-
pletely confined in the dielectric region
FIGURE 3.2. A two-conductor line showing the electric and magnetic fields
V = V0 exp (−jkz)
−
→
and since the line integral of E t is independent of the path chosen between
−
→
the conductors, then E t is the gradient of the scalar potential.
−
→
The line integral of H t around one conductor, say S2 is given by
−
→ −
→
Ht · d l = Js dl = I0
C C
−
→ −
→
where n × H t = J s due to the assumption of a perfect conductor with
n
the outward unit normal vector on the conductor surface. Since n and
−
→ −
→
H t are coplanar in the transverse plane, J s is axial. Far away from the
−
→ − → → −
− →
conductor surfaces, ∇ × H = 0 but the line integral H t · d l = 0 because
C
−
→
of J s . Associated with the magnetic field is a unique current wave
I = I0 exp (−jkz)
−
→ −
→
Since Φ is independent of frequency, both E t and H t are as well.
We now consider the coaxial line of the transverse cross-section shown in
Fig. 3.3. In cylindrical coordinates, the transverse Laplace’s equation takes
the form
1 ∂ ∂Φ 1 ∂Φ
r + 2 =0
r ∂r ∂r r ∂φ
and if Φ has angular symmetry, ∂Φ/∂φ = 0 and
1 ∂ ∂Φ d dΦ
r = r =0
r ∂r ∂r dr dr
20 3. Theory and applications of multi-conductor transmission lines
Hence,
dΦ C1
=
dr r
or
Φ = C1 ln r + C2
Since Φ = V0 @ r = a and Φ = 0 @ r = b, then C2 = −C1 ln b and
C1 = V0 / ln (a/b) so that
ln (r/b)
Φ= V0
ln (a/b)
and
−
→ dΦ −r V0
Et = −
r exp (−jkz) = exp (−jkz)
dr r ln (a/b)
r V0
= exp (−jkz)
r ln (b/a)
and
−
→ 1 −
→
φ V0
H t = r × E t = exp (−jkz)
η r η ln (b/a)
Since the potential difference between the outer and inner conductors is
V0 , then the unique voltage wave associated with the electric field is
V = V0 exp (−jkz)
and the current density on the inner conductor is
−
→ −
→ −
→ z V0
Js=n
× H t = r × H t = exp (−jkz)
a η ln (b/a)
The magnitude of the current on the conductor is thus,
2π
1 V0 V0
|I0 | = a dφ = 2π
a η ln (b/a) η ln (b/a)
0
3.2 The TEM wave along a parallel-plate transmission line 21
where dS = a dφ. The unique current wave associated with the magnetic
field is thus
I = I0 exp (−jkz)
The rate of energy flow or power along the line is from Poynting theorem
b 2π
1 −
→ −
→
P = Re E t × H ∗t · rr dr dφ
2
a 0
2 b 2π
1 V0 dφ dr V02
= =π
2η ln (b/a) r η ln (b/a)
a 0
V0 η
Z0 = = ln (b/a)
I0 2π
h y z
x
and −
→
E µ
η = − −
→ = ǫ
H
−
→ − → −
→ E0
−y × H = J s,h ⇒ J s,h = zHx = − z exp (−jβz)
η
Thus, the surface current and free surface charge densities on the plates
vary sinusoidally with z.
The electric and magnetic field intensities in the dielectric satisfy Maxwell’s
equations, i.e.
−
→ − → −
→ dEy −
→ − → −
→ dHx
∇ × E = −jωµ H ⇒ = jωµHx , ∇ × H = jωǫ E ⇒ = jωǫEy
dz dz
−
→ − →
Thus, integrating ∇ × E over y from y = 0 to y = h, we obtain
h h
d
Ey dy = jωµ Hx dy
dz
0 0
3.2 The TEM wave along a parallel-plate transmission line 23
or
d h
− V (z) = jωµJs,h h = jω µ (Js,h W ) = jωL′ I (z)
dz W
where
h
V (z) = − Ey dy = −hEy (z)
0
I (z) = W Js,h
W x
d
Hx dx = jωǫ Ey dx
dz
0 0
or
d W
− I (z) = −jωǫEy (z) W = jω ǫ (−Ey (z) h) = jωC ′ V (z)
dz h
where
W
C′ = ǫ
h
is the capacitance per unit length of the parallel-plate transmission line
structure.
These equations can be combined into a single second-order differential
equation in either V (z) or I (z) as
d2
V (z) = −ω2 L′ C ′ V (z)
dz 2
or
d2
I (z) = −ω2 L′ C ′ I (z)
dz 2
which have general solutions of the form
where √ √
β = ω L′ C ′ = ω µǫ
Since the line is infinitely long, then
V (z) = V + exp (−jβz)
I (z) = I + exp (−jβz)
We define a line impedance at a point z along the line as
∆ V (z) V+ L′
Z (z) = = + =
I (z) I C′
If only the −z-propagating wave exists on the line, we have
V (z) = V − exp (jβz)
I (z) = I − exp (jβz)
Thus,
d −
− V exp (jβz) = −jβV − exp (jβz) = jωL′ I − exp (jβz)
dz
or
V (z) V− ωL′ L′
Z (z) = = − =− =−
I (z) I β C′
We note that the line impedance is independent of the position along the
line for the case of voltage and current waves propagating in only one
direction. The line impedance for this special case is called the characteristic
impedance Z0 of the line and depends only on the line’s parameters and
geometry. Using the values for L′ and C ′ of a parallel-plate transmission
line, we obtain,
h µ h
Z0 = =η , Ω
W ǫ W
This is the impedance that will be seen at any location on the line looking
towards an infinitely long (reflection-free) line and is known as the charac-
teristic impedance of the line. The wave propagates along the line with a
phase velocity
ω 1 1
vp = = √ = √ , ms−1
β LC′ ′ µǫ
and −1/2
µ/ǫ′ ǫ′′
Z0 = = η 1 − j
1 − jǫ′′ /ǫ′ ǫ′
When the losses are small such that ǫ′′ /ǫ′ << 1, the propagation constant
can be expanded in a binomial series as
ǫ′′ ǫ′′
γ = α + jβ = jk 1 − j ′ ∼ jk + k ′
ǫ 2ǫ
Hence
ǫ′′
α=k
2ǫ′
and
β=k
and the characteristic impedance is given by
ǫ′′
Z0 ∼ η 1 + j ′
2ǫ
where
β = k Z0 = η, ǫ′′ ∼ 0
We saw previously that if ǫ′′ is small but different from zero, it can be
expressed as an equivalent conductivity σd = ωǫ′′ which results in a shunt
−
→ −
→
current J = σd E t between the conductors. The power loss per unit length
of the line is hence,
1 −
→ − → ωǫ′′ → −
− →
Pdiel = J · J ∗ dS = E t · E ∗t dS
2σd 2
S S
and
∂P
− = Pdiel = 2αPin (0) exp (−2αz) = 2αP (z)
∂z
which states that the power loss at any plane z is directly proportional
to the total power present at this plane. From Poynting’s theorem, the
time-average power propagated along the line is given by,
1 −
→ −
→
P (z) = Re E t × H ∗t · z dS
2
S
1 −
→ →
− 1 −
→ − →
= Re E t × z × E ∗t · z dS = Re E t · E ∗t dS
2η 2η
S S
Thus,
Pdiel σdη ωǫ′′ η kǫ′′
αdiel = = = =
2P (z) 2 2 2
the same as before.
where
2
δs =
ωµc σ c
is the skin depth.
The tangential electric field at the conductor surface must therefore be
−
→ −
→ −
→
Zm J s where J s = n × H is the surface current density. Thus, in the case
of a finite σ, the axial electric field is no longer zero and the wave then no
longer truly TEM. This axial field gives rise to a component of Poynting’s
vector directed into the conductor which accounts for the conductor power
−
→
loss. If |Zm | is much smaller than Z0 , Zm J s is also small relative to the
transverse components so that the field solution is very nearly TEM. This
mode of wave propagation is referred to as quasi-TEM.
In applying the perturbation theory, we let
−
→ −
→
Js=n
× Ht
−
→
where n is the unit outward normal to the conductor surface and H t the
lossless transverse field. Power loss per unit length of line due to the surface
impedance is then,
1 −
→ − → Rm →
− →
−
Pc = Re Zm J s · J ∗s dl = n × H ∗t dl
× Ht · n
2 2
C1 +C2 C1 +C2
where
µ
Z0 = η =
ǫ
in the lossless case.
When both the conductors and dielectric are imperfect, the total atten-
uation coefficient is the sum of the individual attenuations provided these
are separately small. Thus, in the coaxial line, with permittivity ǫ = ǫ′ −jǫ′′
and conductors of conductivity σc , the pertinent quantities in the loss-free
case are replaced by
Z0 ∼ η, β ∼ k
Thus
−
→ r V0 −
→
φ V0
Et = exp (−jkz) H t = exp (−jkz)
r ln (b/a) r Z0 ln (b/a)
and
b 2π
1 −
→ −
→ V02
P (z) = Re E t × H ∗t · rr dr dφ = π
2 Z0 ln (b/a)
a 0
b 2π
ωǫ′′ −
→ −
→ ωǫ′′ V02
Pdiel = E t × H ∗t · rr dr dφ = π
2 ln (b/a)
a 0
2 2π
2
Rm Y V0 1 1 V0 1 1
Pc = + dφ = πRm +
2 ln (b/a) a b Z0 ln (b/a) a b
0
Thus,
ǫ
R′ C ′ =
σd
or the conductance per unit length of line,
1 σd ′
G′ = = C
R′ ǫ
For the parallel-plate line, we thus have
W
G′ = σd
h
If the plates have a large but finite conductivity σ c , the ohmic losses in
the conductor require a non-zero tangential electric field at the conductor
surface such that
−
→ 1
P ave = Re [ z Ez × x
Hx ]
2
is non-zero and is equal to the average power per unit area dissipated in
each conducting plate.
Consider the plate at y = h having a surface current density Js,h = Hx .
We can thus define a surface impedance on the conductor as
∆ Etan
Zm = , Ω
Js
so that, for the plate at y = h,
Ez Ez
Zm = = = ηc
Js,h Hx
30 3. Theory and applications of multi-conductor transmission lines
The ratio of the axial to transverse components of the electric field is given
by
Ez |ηc Hx | ǫ ωµc ǫ ωǫ
= = |η | = =
Ey |ηHx | c
µ µσ c σc
so that for Cu in air with σCu ∼ 5, 8 × 107 Sm−1 , ǫ ∼ ǫ0 at f = 3 GHz,
W
t grounded
h
ε h t W
ε
h
grounded grounded
microstrip line triplate line
W W
ε h ε h
slotline
coplanar strip line
a
W
s
ε h
3.3.7 Stripline
Finds common use in microwave filters but has been largely replaced by the
coplanar and microstrip lines. It is advantageous for fabricating couplers
but much more difficult to fabricate. The two-conductor structure and ho-
mogeneous dielectric allows TEM propagation which is the desired mode.
30π b
Z0 = √
ǫr Wef f + 0, 441b
Wef f W W
= , > 0, 35
b b b
2
W W
= 0, 35 − , < 0, 35
b b
which assumes strips of zero thickness. Thus, the design formulas can be
expressed as
W √
= x, Z0 ǫr < 120
b √
= 0, 85 − 0, 6 − x, Z0 ǫr > 120
where
30π
x= √ − 0, 441
Z0 ǫr
The attenuation due to dielectric loss is the same as for other TEM lines.
The attenuation due to conductors may be obtained by perturbation theory
or Wheeler’s incremental inductance rule
2, 7 × 10−3 √
αc = Rm ǫr Z0 A, Z0 ǫr < 120
30π (b − t)
0, 16 √
= Rm B, Z0 ǫr > 120
Z0 b
where
2W 1 b+t 2b − t
A=1+ + ln
b−t π b−t t
and
b 0, 414t 1 4πW
B = 1+ 0, 5 + + ln
0, 5W + 0, 7t W 2π t
where t is the strip thickness. The cross-section of the stripline is illustrated
in Fig. 3.12.
Example: Find the width for a 50 Ω copper stripline conductor with b =
0, 32 and ǫr = 2, 20. If the dielectric has a loss tangent of 0, 001 and the
operating frequency is 10 GHz, calculate the attenuation in dB/λ. Assume
conductor thickness of t = 0, 01 mm.
3.3 Important practical planar transmission lines 37
√ √
Z0 ǫr = 74, 2 < 120 ⇒ x = 30π/ (Z0 ǫr ) − 0, 441 = 0, 83 cm
W = bx = 0, 266 cm
2πf √
k= ǫr = 310, 6 m−1
c
k tan δ
αdiel = = 0, 155 Np/m
2
2, 7 × 10−3 Rm ǫr Zc A
Rm = 0, 026 Ω ⇒ αcond = = 0, 122 N p/m, A = 4, 74
30π (b − t)
α = αdiel + αcond = 0, 277 N p/m ⇒ αdB = −20 log10 eα = 2, 41 dB/m
c
λ = √ = 2, 02 cm ⇒ αdB /λ = (2, 41 dB/m) (0, 0202 m) = 0, 049 dB/λ
f ǫr
• Low loss
• Good machinability
• Good thermal conduction since heavy metal heat sinks are undesir-
able at microwave frequencies as they interact undesirably with the
electromagnetic fields.
38 3. Theory and applications of multi-conductor transmission lines
unit length of line. For the moment, we satisfy ourselves with the following
expressions for the transmission line parameters
µ −
→ − →
L′ = H · H ∗ dS
II ∗
S
derived by equating
1 ∗ ′
II L = Um
4
ǫ′ −
→ − →
C′ = E · E ∗ dS
VV∗
S
derived by equating
1
V V ∗ C ′ = Ue
4
′ Rm −
→ − →
R = ∗ H · H ∗ dl
II
S1 +S2
and
ωǫ′′ −
→ − →
G′ = E · E ∗ dS
VV∗
S
∂
v (z, t) − v (z + ∆z, t) = i (z, t) R′ ∆z + L′ ∆z i (z, t)
∂t
40 3. Theory and applications of multi-conductor transmission lines
L’∆z R’∆z
i(z,t) i(z+∆z,t)
∂
v (z + ∆z, t) ∼ v (z, t) + ∆z v (z, t)
∂z
Hence
∂ ∂
−∆z v (z, t) = i (z, t) R′ ∆z + L′ ∆z i (z, t)
∂z ∂t
or
∂ ∂
v (z, t) = −i (z, t) R′ − L′ i (z, t)
∂z ∂t
Similarly
∂
i (z, t) − i (z + ∆z, t) = v (z + ∆z, t) G′ ∆z + C ′ ∆z v (z + ∆z, t)
∂t
or
∂ ∂ ∂ ∂
−∆z i (z, t) = v (z, t) + ∆z v (z, t) G′ ∆z + C ′ ∆z v (z, t) + ∆z v (z, t)
∂ ∂z ∂t ∂z
∂
= v (z, t) G′ ∆z + C ′ ∆z v (z, t)
∂t
since ∆z is small. Thus
∂ ∂
i (z, t) = −v (z, t) G′ − C ′ v (z, t)
∂z ∂t
We differentiate the first equation w.r.t. z and the second one w.r.t. t to
obtain
∂2 ′ ∂ ′ ∂
2
v (z, t) = −R i (z, t) − L i (z, t)
∂z 2 ∂z ∂z ∂t
∂2 ∂ ∂2
i (z, t) = −G′ v (z, t) − C ′ 2 v (z, t)
∂t ∂z ∂t ∂t
3.4 A distributed circuit analysis 41
or
∂2 ′ ′ ′ ∂ ′ ′ ∂ ∂2
v (z, t) = −R −v (z, t) G − C v (z, t) −L −G v (z, t) − C v (z, t)
∂z 2 ∂t ∂t ∂t2
or
∂2 ∂ ∂2
2
v (z, t) = R′ G′ v (z, t) + (R′ C ′ + L′ G′ ) v (z, t) + L′ C ′ 2 v (z, t)
∂z ∂t ∂t
or
∂2 ′ ′ ′ ′ ′ ′ ∂ ′ ′ ∂
2
v (z, t) − R G v (z, t) − (R C + L G ) v (z, t) − L C v (z, t) = 0
∂z 2 ∂t ∂t2
For harmonic time-dependent voltage and current, we can write
2
γ − R′ G′ − jω (R′ C ′ + L′ G′ ) + ω2 L′ C ′ V (z) = 0
where
V (z) = V + exp (−γz) + V − exp (γz)
is the line voltage in phasor form at z down the line and
γ = jω (R′ C ′ + L′ G′ ) − ω 2 L′ C ′ + R′ G′
′
√ 1 R′ G′ 1 R′ G
′
= jω L C 1 + ′ + ′ − 2
jω L ′ C ω L′ C′
In the forward going and reflected waves, I + , V + and I − , V − are the inci-
dent and reflected quantities respectively.
The characteristic impedance of the line is given by
V+ V− R′ + jωL′ R′ + jωL′ L′ 1 + (R′ /L′ ) /jω
Zc = + = − − = = =
I I γ G′ + jωC ′ C ′ 1 + (G′ /C ′ ) /jω
or
′
L′ 1 R G′
Zc ∼ 1+ − ′
C′ j2ω L′ C
In the loss-free case √
γ = jω L′ C ′
and
L′
Zc = Z0 =
C′
d
− Hy (z) = (ωǫ′′ + jωǫ′ ) Ex (z)
dz
Comparing with
d
− V (z) = (R′ + jωL′ ) I (z)
dz
and
d
− I (z) = (G′ + jωC ′ ) V (z)
dz
Then the equations governing Ex (z), Hy (z) and V (z), I (z), are similar
in form. We combine to obtain
d2
Ex (z) = γ 2 Ex (z)
dz 2
d2
Hy (z) = γ 2 Hy (z)
dz 2
with
γ = α + jβ = (ωµ′′ + jωµ′ ) (ωǫ′′ + jωǫ′ )
and the intrinsic impedance of the lossy medium for a wave propagating in
the +z-direction is
µ′′ + jµ′
η=
ǫ′′ + jǫ′
Hence
1 ′ C′ L′ 1 R′ L′ √
α∼ R + G′ = ′
+ G Z0 , Z0 = , β = ω L′ C ′
2 L′ C′ 2 Z0 C′
ω 1
=√ vp ∼
β L′ C ′
−1/2 ′
′
L′ R′ G′ L 1 R G′
Zc = 1+ 1+ ∼ 1+ − ′
C′ jωL′ jωC ′ C′ j2ω L′ C
or
L′
Zc ∼ = Z0
C′
since the frequency dependent imaginary part is negligible when R′ <<
ωL′ and G′ << ωC ′ .
(c) Distortionless line (R′ /L′ = G′ /C ′ )
′ ′ ′
RC C
′
γ = α + jβ = (R + jωL ) ′ + jωC = ′ (R′ + jωL′ )
L′ L′
′ C′ √
α=R ′
, β = ω L′ C ′
L
ω 1
vp ∼ = √
β L′ C ′
R′ + jωL′ L′
Zc = ′ ′ ′ ′
= = Z0
R C /L + jωC C′
(a) Determine the resistance, inductance and conductance per unit length
of the line
(b) Find the velocity of wave propagation and the velocity factor
(c) Determine the percentage to which the amplitude of a travelling voltage
wave decreases in 1 km and 5 km
(a) For a distortionless line,
R′ G′
′
= ′
L C
Thus
L′
Z0 = = 50 Ω
C′
3.5 Propagation characteristics on an infinitely long transmission line 45
Thus
R′ = αZ0 = 0, 057 Ωm−1
L′ = C ′ Z02 = 0, 25 µHm−1
R′ C ′
G′ = = 22, 8 µSm−1
L′
(b)
1
vp = √ = 2 × 108 ms−1
L′ C ′
and
∆ v 2 × 108 ms−1
velocity factor = = = 0, 7
c 3 × 108 ms−1
(c)
V2
= exp (−αz) = 31, 7% @ 1 km
V1
= 0, 32% @ 5 km
V+ V+
I (z) = exp {− (α + jβ) z} = exp (−αz) exp (−jβz)
Z0 Z0
The time-average power propagated down the line across the plane at z is
given by
1 1 V +2
P (z) = Re [V (z) I ∗ (z)] = Re [Z0 ] exp (−2αz)
2 2 |Z0 |2
46 3. Theory and applications of multi-conductor transmission lines
(a) For a lossy transmission line, the time-average power loss per unit
length is
1 2 V02 ′
Pl (z) = I (z) R′ + |V (z)|2 G′ = 2 R + G′ |Z0 |2 exp (−2αz)
2 2 |Z0 |
Thus
1 ′
α= R + G′ |Z0 |2 , N pm−1
2Z0
(b) For a low-loss line, Z0 ∼ L′ /C ′
1 R′ ′ 1 ′ C′ ′ L′
α∼ + G Z0 = R ′
+G
2 Z0 2 L C′
and for a distortionless line, Z0 = L′ /C ′ , R′ /L′ = G′ /C ′ so that
′ C′
α=R
L′
Zg
Ii IL
Vg Vi Zi ZL
VL
γ,Z0
-l -z 0
The reflected waves V − (z) exp (γz) and I − (z) exp (γz) vanish when the
line is terminated in its characteristic impedance. Consider a load ZL con-
nected to a sinusoidal generator Vg of source impedance Zg through a
finite-length line of length l as depicted in Fig. 3.15. For convenience, we
assume the load to be at z = 0 so that the source is at z = −l. Thus, at
the load, we can write
V+ V−
V (0) = VL = V + + V − , I (0) = IL = I + + I − = −
Z0 Z0
where the characteristic impedance of the line Z0 is defined as
V+ V−
Z0 = +
=− −
I I
By Ohm’s law,
VL = IL ZL
Thus, we can eliminate the wave amplitudes by writing
VL = IL ZL = V + + V −
IL Z0 = V + − V −
from which we see that
IL IL
V+ = (ZL + Z0 ) , V − = (ZL − Z0 )
2 2
IL IL
I (z) = (ZL + Z0 ) exp (−γz) − (ZL − Z0 ) exp (γz)
2Z0 2Z0
IL IL
= Z0 cosh γz − ZL sinh γz
Z0 Z0
48 3. Theory and applications of multi-conductor transmission lines
Zg Ii
Vg Vi Zi
Xi,oc
inductive
λ/2 3λ/4 λ
l
0 λ/4 5λ/4
capacitive
ZL + jZ0 tan βl
Zi = Z0
Z0 + jZL tan βl
Z0
Zi,oc = = −jZ0 cot βl
j tan βl
Thus, depending on the sign of cot βl, Zi,oc can have either inductive or
capacitive characteristics. If βl is small, tan βl ∼ βl and
Z0 1 L′ 1
Zi,oc ∼ = √ ′
=
jβl jω L′ C ′ l C jωC ′ l
Xi,sc
inductive
λ/4 λ 5λ/4
l
0 λ/2 3λ/4
capacitive
A quarter-wavelength line
If l is such that
λ
l = (2n + 1)
4
then
2π λ π
βl = (2n + 1) = (2n + 1)
λ 4 2
and
π
tan βl = tan (2n + 1) = ±∞
2
In this case,
R20
Zi =
ZL
52 3. Theory and applications of multi-conductor transmission lines
A half-wavelength line
For l = nλ/2, n = 1, 2, . . .,
2π nλ
βl = = nπ
λ 2
and
tan βl = 0
so that
Zi = ZL
We can determine Z0 and γ of a given line by measuring the input
impedance of a section of the line with open- and short-circuit terminations.
For ZL → ∞, Zi,oc = Z0 coth γl and for ZL = 0, Zi,sc = Z0 tanh γl. Thus
Z0 = Zi,oc Zi,sc
and
1 Zi,sc
γ = tanh−1
l Zi,oc
These expressions hold regardless of whether the line is lossy or not.
Example: The open- and short-circuit impedances measured at the input
terminals of a lossless transmission line of length 1, 5 m less than λ/4 are
−j54, 6 Ω and j103 Ω resp.
(a) Determine Z0 and γ of the line
(b) Without changing the operating frequency, determine the input im-
pedance of a short-circuited section of this line that is twice as long
as the one given.
(c) How long should the short-circuited line be in order for it to appear as
an open circuit at the input terminals at this frequency?
(a)
Z0 = Zi,oc Zi,sc = −j54, 6 × j103 = 75 Ω
and
1 Zi,sc 1 j103
γ = tanh−1 = tanh−1 = j0, 628 m−1
l Zi,oc 1, 5 −j54, 6
(b)
γl = j0, 628 × 2 × 1, 5 m = j1, 884 rad
Zi,sc = 75 tanh (j1, 884) Ω = j75 tan 1, 884 = −j231 Ω
Thus
λ
1, 5 m < < 3, 0 m
4
3.6 Propagation characteristics on finite-length transmission lines 53
and −2
|Zi,sc |2 nπ 2
∆f 2
= 1+
|Zi,sc,max |2 2αl 2f0
Therefore, the half-power frequencies f1 and f2 , where f2 > f1 are where
this ratio is 1/2. i.e.,
nπ ∆f β ∆f
= = 1, n odd
2αl 2f0 2α f0
The Q of the parallel-resonant circuit is hence,
f0 β
Q= =
∆f 2α
and assuming low loss,
′
ωL′ R G′ 1 ′ L′ ′
Q= = + = R + G
R′ + G′ L′ /C ′ ωL′ ωC ′ ωL′ C′
For a well insulated line, G′ L′ /C ′ << R′ . Hence
ωL′
Q∼
R′
Example: The measured attenuation of an air-dielectric coaxial line at
400 M Hz is 0, 01 dBm−1 . Determine the Q and half-power bandwidth of
a quarter-wavelength section of this line terminated in a short-circuit.
At f = 4 × 108 Hz, λ = 0, 75 m, β = 8, 38 radm−1
0, 01
α =, 01 dBm−1 = N pm−1
8, 69
Thus
β 8, 38 × 8, 69
Q= = = 3641
2α 2 × 0, 01
f0
∆f = = 0, 11 × 106 Hz
Q
where
ZL − Z0
ΓL = = |ΓL | exp (jθ Γ )
ZL + Z0
is the voltage reflection coefficient of the load. Similarly, for the current
wave,
IL IL
I (z) = (ZL + Z0 ) exp (−γz) − (ZL − Z0 ) exp (γz)
2Z0 2Z0
IL ZL − Z0
= (ZL + Z0 ) exp (−γz) 1 − exp (2γz)
2Z0 ZL + Z0
IL
= (ZL + Z0 ) exp (−γz) {1 − ΓL exp (2γz)}
2Z0
Note that the current reflection coefficient defined as
∆ I− V−
ΓI = = −
I+ V+
is equal to the negative of the voltage reflection coefficient. Thus
IL
V (−z) = (ZL + Z0 ) exp (γz) {1 + ΓL exp (−2γz)}
2
and
IL
I (−z) = (ZL + Z0 ) exp (γz) {1 − ΓL exp (−2γz)}
2Z0
For a lossless transmission line, γ = jβ. Hence
IL
V (−z) = (ZL + Z0 ) exp (jβz) [1 + |ΓL | exp {−j (2βz − θ Γ )}]
2
IL
I (−z) = (ZL + Z0 ) exp (jβz) [1 − |ΓL | exp {−j (2βz − θΓ )}]
2Z0
The line current and voltage phasors are obtained as
which are standing wave patterns with minima and maxima located at fixed
points along the line. The voltage standing wave ratio (VSWR) is defined
as,
∆ |Vmax | 1 + |ΓL |
S= =
|Vmin | 1 − |ΓL |
|Vmax | and |Imin | occur together at z = −zM such that
θ Γ − 2βzM = 2nπ, n = 0, 1, 2, . . .
θΓ − 2βzm = (2n + 1) π, n = 0, 1, 2, . . .
(a) Show how the value of the terminating resistance on a lossless line of
known characteristic impedance Z0 can be determined by measuring
V SW R.
(b) What is the impedance of the line looking towards the load at a distance
λ/4 where λ is the operating wavelength?
Z0
|Vmax | = VL , |Vmin | = VL
RL
RL
|Imin | = IL , |Imax | = IL
Z0
Thus
|Vmax | |Imax | RL
= =S=
|Vmin | |Imin | Z0
or
RL = SZ0
If RL < Z0 , θ Γ = π and |Vmin |, |Imax | occur at βz = 0 and |Vmax |,
|Imin | at βz = −π/2.
Z0
|Vmin | = VL , |Vmax | = VL
RL
RL
|Imax | = IL , |Imin | = IL
Z0
3.6 Propagation characteristics on finite-length transmission lines 57
Thus
|Vmax | |Imax | Z0
= =S=
|Vmin | |Imin | RL
or
Z0
RL =
S
(b) λ is determined from twice the distance between to neighbouring voltage
(or current) maxima or minima. At z = −λ/4, βz = −π/2, cos βz =
0, sin βz = −1. Thus
λ λ VL
V − = jIL Z0 , I − =j
4 4 Z0
λ V (−λ/4) Z2
Z − = Zi = = 0
4 I (−λ/4) RL
θ Γ = π − 2βzm , n = 0
58 3. Theory and applications of multi-conductor transmission lines
Z0 ZL
z’m z’=0
lm
Z0 Rm
z’m z’=0
λ/2
(iii) Determine ZL as
1 + |ΓL | exp (jθΓ )
ZL = RL + jXL = Z0
1 − |ΓL | exp (jθΓ )
θΓ − 2βzM = −2nπ, n = 0
λ = 2 × 0, 2 m = 0, 4 m
2π
β= = 5π rads−1
λ
S−1
|ΓL | = = 0, 5
S+1
3
θ Γ = π − 2βzm = π rad
2
ΓL = |ΓL | exp (jθΓ ) = −j0, 5
(b)
1 − j0, 5
ZL = 50 = 30 − j40 Ω
1 + j0, 5
(c)
Rm + j50 tan βlm
30 − j40 = 50
50 + jRm tan βlm
Equate real and imaginary parts or recall
λ Z0
|zm | + lm = , Rm =
2 S
Hence
λ
lm = − |zm | = 0, 15 m
2
60 3. Theory and applications of multi-conductor transmission lines
50
Rm = = 16, 7 Ω
3
Note that
′ λ
lm = lm − = 0, 05 m
4
and
′
Rm = Z0 S = 150 Ω
are also valid when we use the locations of the maxima instead.
Z0 Vg 1 + ΓL exp (2γz)
V (z) = exp {−γ (z + l)}
Z0 + Zg 1 − Γg ΓL exp (−2γl)
and
Vg 1 − ΓL exp (2γz)
I (z) = exp {−γ (z + l)}
Z0 + Zg 1 − Γg ΓL exp (−2γl)
Example: A 100 M Hz sinusoidal generator with Vg = 10∠0◦ and source
resistance 50 Ω is connected to a lossless 50 Ω air line 3, 6 m long and
terminated in a 25 + j25 Ω load. Determine
Thus
2π 2πf 2π
β= = = radm−1 ⇒ βl = 2, 4π rad
λ c 3
ZL − Z0
ΓL = = 0, 447∠116, 6◦ , Γg = 0
ZL + Z0
At −z from the load, we have
Z0 Vg
V (z) = (1 + ΓL exp (−j2βz)) exp {−jβ (l − z)}
Z0 + Zg
Thus, measured from the generator, this will be the point z − l and
we therefore have
Z0 Vg
V (z − l) = [1 + ΓL exp {j2β (z − l)}] exp {jβ (z − l)}
Z0 + Zg
where z ≥ 0 and measured from the generator.
(b)
Z0 Vg
Vi = V (0) = {1 + ΓL }
Z0 + Zg
Z0 Vg
= {1 + 0, 447 exp (−j4, 8π + j0, 6478π)} exp (−j2, 4π)
Z0 + Zg
Z0 Vg
= {1 + 0, 447∠ − 27, 4◦ } exp (−j2, 4π)
Z0 + Zg
= 5 (1 + 0, 3969 − j0, 2057) exp (−j2, 4π) V
= 7, 06∠ − 80, 4◦ V
and
Z0 Vg
VL = V (l) = [1 + ΓL ]
Z0 + Zg
= 5 {1 + 0, 447∠116, 6◦ } V
= 5 {1 − 0, 2001 + j0, 3997}
= 4, 47∠26, 6◦ V
(c)
1 + |Γ|
S= = 2, 62
1 − |Γ|
62 3. Theory and applications of multi-conductor transmission lines
(d)
2
1 VL
Pave = RL = 0, 2 W
2 ZL
If ZL = 50 Ω, Γ = 0 and |VL | = |Vi | = Vg /2. Maximum power would
thus be delivered to the load and is called the available power. This
has the value
V2
Pave,max = L = 0, 25 W
2RL
The reflected power is hence
2
|Γ| Pave,max = 0, 05 W
Zg t=0
VDC Z0 Z0
z1
z=-l z=0
and
I (z1 ) = I1+ + I1− + I2+ + I2− + . . .
The ultimate voltage across the load as time progresses will be
VL = V1+ 1 + ΓL + Γg ΓL + Γg Γ2L + . . .
= V1+ 1 + Γg ΓL + Γ2g Γ2L + . . . + ΓL 1 + Γg ΓL + Γ2g Γ2L + . . .
or
1 + ΓL
VL = V1+
1 − Γg ΓL
64 3. Theory and applications of multi-conductor transmission lines
For RL = 3Z0 , Rg = 2Z0 say, Γg = 1/3, ΓL = 1/2 and V1+ = VDC /3. Thus
VL → 95 V1+ = 35 VDC as t → ∞. Similarly
vg (t) = V0 {u (t) − u (t − T0 )}
The transient response of the line is the superposition of the two step
functions.
Example: A rectangular pulse of an amplitude 15 V and duration 1 µs is
applied through a series resistance of 25 Ω to the inputs of a 50 Ω lossless
coaxial line. If the line is 400 m long and is terminated in a short-circuit
at the far end, determine the voltage response midway down the line as a
function of time up to 8 µs. The dielectric constant is 2, 25.
ΓL = −1, Γg = −1/3. We have
vg (t) = 15 u (t) − u t − 10−6
c
v = √ = 2, 0 × 108 ms−1
ǫr
l
T = = 2 µs
v
Z0
V1+ = 15 V = 10 V
Z0 + RG
Fig. 3.23 illustrates
the solution
process. The values indicated in brackets
apply to the −15u t − 10−6 function applied 1 µs later.
3.7 Transients on transmission lines 65
t5 P5
Γ2g Γ2L V+1
4T
P4
t4 ΓgΓ2L V+1
3T
t3 P3
ΓgΓLV+1
2T
P2
t2 ΓLV+1
T
V+1
P1
t1 slope=1/v
z1 z
0 l
ΓgΓ2L V+1
t4 l
Z0 V0
V1+ = − V0 = −
Z0 + R 2
At t = l/v, this wave arrives at the open-circuit having lowered the voltage
everywhere on the line to V0 /2. A reflected wave V1− is sent back to the
sending end where
V0
V1− = ΓV1+ = −
2
and arrives here at t = 2l/v and lowers the line voltage to zero.
The current through the resistor
IR = −I1
3.7 Transients on transmission lines 67
9 +1/9
(+1/3)
8
-1/3
+1/3 6
5 (-1/3)
(+1)
4
-1
+1 2
1 (-1)
z
0 200 400
Z0 t=0 t=0
IR I+1
R VR Z0 LINE
t=0 t=0
I1
IR
R LINE R LINE
z=0 z=0
(c) VOLTAGES OPPOSE HENCE NO TRANSIENTS (d) INTERESTING PART
where
V1+ V0
I1 = I1+ = =−
Z0 2Z0
in the time interval 0 ≤ t ≤ l/v. At t = l/v, I1+ reaches the open circuit
and the reflected current must change sign to make the total current zero
as required by the open circuit. Thus
V0
I1− = −I1+ =
2Z0
V0
V1+ =
2
3.7 Transients on transmission lines 69
Z0 t=0 Z0 t=0
iL(t) iL(t)
iL(t) vL(l,t)
t t
T T
V1- (l,t)
v(z,t)
V+1
2V+1
t
T
V+1
-V+1
z
z1 l
Z0 t=0
iC(t)
V0 Z0 C vC(t)
z=-l z=0
vC(l,t) iC(t)
t t
T T
V1- (l,t)
v(z,t)
V+1
2V+1
t
T
-V+1 z
z1 l
(rL − 1) + jxL
Γr + jΓi =
(1 + rL ) + jxL
or equivalently
(1 + Γr ) + jΓi
rL + jxL =
(1 − Γr ) − jΓi
From which we see that
1 − Γ2r − Γ2i 2Γi
rL = 2 , xL =
(1 − Γr ) + Γi 2 (1 − Γr )2 + Γ2i
For a given value of rL , the locus of the reflection coefficient describes a
circle specified through the equation
2
2
rL 1
Γr − + Γ2i =
1 + rL 1 + rL
with the centre located in the complex Γ-plane at
rL
,0
1 + rL
and a radius
1
1 + rL
xL forms another family of circles
2
2
1 1
(Γr − 1)2 + Γi − =
xL xL
3.8 The Smith-chart 73
x=1.0
r=0
r=0.5 r=1.0 r=2.0
x=-0.5
x=-1.0 x=-2.0
unit circle
FIGURE 3.29. Simplified Smith chart
of radius
1
xL
and centred at
1
1,
xL
The Smith-chart is thus a graph of r- and x-circles in the complex Γ-plane
for|Γ| ≤ 1. Fig. 3.29 presents an example of the impedance Smith-chart.
The fractional distance from the center of the chart to a given point
relative to the radius |Γ| = 1 of the circle is a measure of the value of |Γ|
at that point. θΓ is the angle the line from the centre passing through the
given point makes with the positive real axis. Thus,
• A |Γ| circle intersects the real axis at two points: at PM towards the
open-circuit point at (r = 1, x = 0), and Pm towards the short-circuit
point at (r = 0, x = 0)
• PM and Pm represent a case of pure resistive loading on the line with
PM representing RL > Z0 and Pm representing resistive loading with
RL < Z0 . Since
RL
S= = rL
Z0
74 3. Theory and applications of multi-conductor transmission lines
For a lossless line therefore, the impedance seen looking toward the load
at a point z < 0 from the load is given by
where like before, we locate the load at the origin of the z-axis. Using
normalised impedance, we can write
where
ϕ = θΓ − 2βz
Thus keeping |Γ| constant and subtracting (rotating clockwise) from θΓ
through 2βz = 4πz/λ locates the point |Γ| exp (jϕ), the reflection coeffi-
cient of the input impedance a distance z from the load zL .
Example: Use the Smith-chart to determine the input impedance of a
section of a 50 Ω lossless line that is 0, 1λ long and terminated in a short
circuit.
We have zL = 0, Z0 = 50 Ω, −z = 0, 1λ. On the Smith-chart, mark Psc
at the intersection of r = 0 and x = 0. Move along the circumference of
the chart ( |Γ| = 1) a distance 0, 1λ towards the generator in a clockwise
direction to point P1 as shown on Fig. 3.30. At P1 read off r = 0 and
x ∼ 0, 725 or zi = j0, 725 ⇒ Zi = zi Z0 = j36, 3 Ω. Confirm this result
using
Zi = jZ0 tan βl, (Ω)
Given
−l = 0, 434λ, Z0 = 100 Ω, ZL = 260 + j180 Ω
Thus
3.8 The Smith-chart 75
P1
Psc
(a) We have
ZL
zl = = 2, 60 + j1, 8
Z0
zl − 1, 0
|Γ| = = 0, 60
zl + 1, 0
We plot on the Smith chart, the point P2 = (2.6, 1.8) and draw a
circle of radius OP2 = |Γ| centred at the origin of the chart ( (0, 0)).
We then draw the line OP2 and extend it to the point P2′ on the unit-
circle and read off the length of line toward the generator as 0, 220.
Phase angle of reflection coefficient is hence
(b) The |Γ| = 0, 60 circle intersects the positive real axis OPoc at r = S =
4 so that the voltage standing-wave ratio is 4
76 3. Theory and applications of multi-conductor transmission lines
(c) For the impedance, advance P2′ at 0, 220 through l = 0, 434λ towards
the generator to point P3′ , first to 0, 50 equivalent to 0, 0 and a further
0, 154
0, 50 − 0, 220 = 0, 280
and then
0, 280 + 0, 154 = 0, 434
since a complete revolution around the Smith chart corresponds to a
length of λ/2.
Join the origin with P3′ which then intersects the |Γ| = 0, 6 circle at
point P3 and read off the value of r and x at P3
zi = 0, 69 + j1, 2
Zi = Z0 zi = 69 + j120 Ω
(d) Going from P2 to P3 , the |Γ| = 0, 6 circle intersects the positive real
axis OPoc at PM where the voltage is maximum. Hence, maximum
voltage occurs at
(0, 25 − 0, 22) λ = 0, 03λ
from the load.
(a) For a lossy line of a sufficient length l such that 2αl cannot be ignored,
the expression for input impedance must be modified to
P’3
P3
P2 P’2
toward generator
Psc Poc
O PM
toward load
on the chart as point P1 as shown in Fig. 3.32 and draw a straight line
from the centre of the circle through P1 to P1′ on the circumference
of the chart. Now measure OP1 /OP1′ = exp (−2αl) = 0, 89. Hence
determine α as
1 1 1
α= ln = ln 1, 124 = 0, 029 N pm−1
2l 0, 89 4
Record the length of the arc Psc P1′ in wavelengths toward the generator
(i.e. clockwise) which gives the number of wavelengths l/λ = 0, 20 and
l
2βl = 4π = 4π × 0, 20 = 0, 8π ⇒ β = 0, 20π
λ
0, 364 + 0, 20 − 0, 5 = 0, 064
and join P3′ to the centre of the chart by a straight line. This line
intersects the |Γ|-circle at P3 . Mark on OP3 a point Pi such that
78 3. Theory and applications of multi-conductor transmission lines
P’3
P1 P’1
P3
Pi
P2
P’2
λ/4
64Ω
B
Z’01
Z0
Z’02
B’ 25Ω
λ/4
(b) Under matched conditions, there are no standing waves on the main
line. On the matching lines however, we have
′
RL1 − Z01 1 + |Γ1 |
Γ1 = ′ = −0, 11 ⇒ S1 = 1 − |Γ | = 1, 25
RL1 + Z01 1
′
RL2 − Z02 1 + |Γ2 |
Γ2 = ′ = −0, 33 ⇒ S2 = 1 − |Γ | = 1, 99
RL2 + Z02 2
yB + ys = 1
80 3. Theory and applications of multi-conductor transmission lines
d
B
Z0 yB
yi yL
ys
Z0
B’
where
yB = Z0 YB , ys = Z0 Ys
with YB and Ys the admittances of the load section and short-circuited
stub respectively.
The short-circuited stub is purely susceptive so that the normalised ad-
mittance is satisfied only if
yB = 1 + jbB , ys = −jbB
Thus, we seek d such that yB looking towards the load at the junction has
unity real part and find l for the stub such that its susceptance cancels that
of the load. Analytically, we can write
(rL + jxL ) + jt
zB =
1 + j (rL + jxL ) t
and
2
rL t − (1 − xL t) (xL + t)
bB = 2 + (x + t)2
rL L
yB = 1 + jbB , ys = −jbB
Thus
2
2xL t rL − rL − x2L
gB = 1 ⇒ t2 − + =0
rL − 1 rL − 1
or
1
2 2
t = xL ± rL (1 − rL ) + xL , rL = 1
rL − 1
xL
= − , rL = 1
2
and
d 1
= tan−1 t, t ≥ 0
λ 2π
1
= π + tan−1 t , t < 0
2π
l 1 1
= tan−1 , bB ≥ 0
λ 2π bB
1 −1 1
= π + tan , bB < 0
2π bB
Double-stub matching
Consider Fig. 3.35. d0 is chosen arbitrarily e.g. λ/16, λ/8, 3λ/16, 3λ/8,
etc. lA and lB are adjusted to match ZL to the main line. Total admittance
looking into the terminals B − B ′ looking towards the load must be equal
to the conductance of the main line, i.e.
1
Yi = YB + YSB = Y0 =
Z0
or
1 = yB + ySB
where ySB is purely imaginary. This hence requires that
d0
B A
yA
Z0
yB ZL
yi
ySA
ySB
Z0
B’ Z0
A’
lB lA
and that on S1 is then −Q′ , where C1 and C2 are the peripheries of con-
ductors S1 and S2 . Since the potential difference between the conductors
3.9 ♣Transmission line parameters 83
− −
→ → −
→
since H t = η−1 E t = η−1 n
· E t at the surface S2 and the tangential
electric and normal component of the magnetic fields are separately equal
to zero when the conductor is perfect. The characteristic impedance of the
line is then
V0 V0 ǫ
Z0 = = ′ ǫη = ′ η
I0 Q C
Thus −
→ 2
H t dl
Rm C1 +C2
R′ = 2
2 2
−
→
H t dl
C2
or
2
− Xm R′ I0 R ′ I0
Xm → 2
Um = H t dl = =
4ω 4ω Rm ω 2
S1 +S2
Thus
ǫ′ Q′ 2πǫ′
C′ = =
ǫ V0 ln (b/a)
86 3. Theory and applications of multi-conductor transmission lines
−
→ − →
φ V0
H = H t exp (−jβz) = exp (−jβz)
r η ln (b/a)
and the current
2π
−
→ 2πV0
I0 = H t · φa dφ =
η ln (b/a)
0
Hence
V0 ln (b/a)
Z0 = =η
I0 2π
and the flux linking the center conductor is
b b
−
→ V0 dr µ
Ψ = µ H t · φdr =µ = V0
η ln (b/a) r η
a a
Thus
Ψ
′ µηV0 b µ b
L = = ln = ln
I0 2πηV0 a 2π a
and
ωǫ′′ ′ ωǫ′′ 2πǫ′ 2πωǫ′′
G′ = C = =
ǫ′ ǫ′ ln (b/a) ln (b/a)
and the power loss in the conductors
1 ′ 2 Rm πV02 1 1
R I0 = Pcond = 2 +
2 η ln (b/a) a b
so that
Rm 1 1
R′ = +
2π a b
and
R′
Lm =
ω
Hence
′ µ b 1 1 1
L + Lm = ln + +
2π a 2πωσδ s a b
If D/a >> 1,
D D
cosh−1 ∼ ln
a a
The inductance per unit length is
µ D
L′ = cosh−1 , Hm−1
π a
and
πσd
G′ = , Sm−1
cosh−1 D/a
Assuming a surface current density Js flows in a very thin surface layer of
the conductor, the current in each wire is given by
I = 2πaJs
and
1 Rm
Pc = I 2 , W m−1
2 2πa
and
′ Rm 1 πfµc
R =2 − , Ωm−1
2πa πa σc
Note that Js is not uniform over the circumference of the wires in the
strict sense as the proximity of the two wires will tend to distort current
distribution.
88 3. Theory and applications of multi-conductor transmission lines
4
Waveguides
y
x z
or
−
−
→ →
∇2t E t + γ 2 + k2 E t + z ∇2t Ez + γ 2 + k2 Ez = 0
or
−
→ −
→
∇2t E t + kc2 E t = 0, kc2 = γ 2 + k2
and
∇2t Ez + kc2 Ez = 0, kc2 = γ 2 + k2
For simple, media, we have already seen that the reduced Maxwell’s curl
equations can be expressed as,
−
→ −
→ −
→ −
→
∇ t × E t = −jωµ
z Hz , ∇ t × H t = jωǫ
z Ez
and
−
→ →
− −
→ −
→ →
− −
→
z × ∇ t Ez + γ E t = jωµ H t , z × ∇ t Hz + γ H t = −jωǫ E t
Also multiplying the curl equation of the magnetic field by −jωµ, we obtain
−
→ −
→ −
→
−jωµγ
z × H t = (−jωµ) (−jωǫ) E t + jωµ
z × ∇ t Hz
or
−
→ jωµ −
→ γ−→
E t = 2 z × ∇ t Hz − 2 ∇ t Ez
kc kc
kc is a characteristic wavenumber of the guide. In a similar way we can
express the transverse component of the magnetic field intensity as
−
→ jωǫ −
→ γ−→
H t = − 2 z × ∇ t Ez − 2 ∇ t Hz
kc kc
These equations show that the transverse field components can be deter-
mined from the axial components alone.
Depending on the values of Hz and Ez , we have
which is the same as the phase velocity in the unbounded medium filling
the guide. The intrinsic impedance or wave impedance for TEM waves is
given by
jωµ γ µ
η T EM = = T EM = =η
γ T EM jωǫ ǫ
just like for the unbounded medium. Thus,
−
→
−
→ jωǫ −
→ z × E t
Ht = z × E t =
γ T EM ηT EM
Since magnetic field lines always close on themselves, TEM waves exist in
−
→ −
→
a waveguide only if the field lines of B , hence of H form closed loops in
the transverse plane.
92 4. Waveguides
From the general Ampere’s circuital law, the line integral of the magnetic
field intensity in any closed contour in a transverse plane must be equal
to the sum of the axial conduction and displacement currents enclosed by
the path of integration. In the absence of an inner conductor in a guide,
no axial current can exist inside a waveguide. Hence no closed magnetic
loops will exist in any transverse plane. Therefore, hollow waveguides of
any shape do not support the TEM mode of propagation.
and since
−
→ k2 −→
∇ t Ez = − c E t
γT M
then
−
→
−
→ jωǫ −
→ z × E t
Ht = z × E t =
γT M ηT M
where
γT M
ηT M =
jωǫ
is the wave impedance for the TM wave.
Since we have
kc2 = γ 2T M + ω 2 µǫ ⇒ γ T M = kc2 − ω 2 µǫ
We define kc2 = ω 2c µǫ or
kc kc c
fc = √ =
2π µǫ 2π
4.1 General wave characteristics along uniform guiding structures 93
where
1
c= √
µǫ
is the phase velocity of a plane wave in the unbounded medium. At an
operating frequency f = fc , γ T M = 0 and the wave does not propagate in
the waveguide. fc is hence called the cut-off frequency of the guide. The
value of fc for a given waveguide mode depends on the eigenvalue of that
mode, i.e. the particular value of kc for the mode. Thus,
2
fc
γ T M = jk 1 −
f
The phase velocity of the mode propagating in the guide is then given
by
ω c λg
vp = = = c>c
β λ
1 − (fc /f )2
The phase velocity in the guide is hence always greater than c in
the unbounded medium and frequency-dependent as well through the
ratio λg /λ. Waveguides are therefore dispersive.
94 4. Waveguides
ηTM/η
evanescent
region
f/fc
1
and with
−
→ k2 −→
∇ t Hz = − c H t
γT E
we have
−
→ jωµ −
→ jωµ −
→ −
→
E t = 2 z × ∇ t Hz = − z × H t = −η T E z × H t
kc γT E
where
jωµ
ηT E =
γT E
is the wave impedance of the TE waves.
From
2
fc
γ T E = kc2 − ω 2 µǫ = jk 1 −
f
we note,
is purely resistive.
96 4. Waveguides
ηTE/η
f/fc
1
and
jωµ
ηT E =
k (fc /f)2 − 1
which is purely reactive. The variation of the ηT E normalised to η
with the frequency normalised to fc is shown in Fig. 4.3.
Example:
(b) Repeat for a frequency equal to one half the cut-off frequency
(c) What are the wave impedance and guide wavelength for the TEM mode?
Thus
ηT M = 0, 866η < η, λT M = 1, 155λ > λ
ηT E = 1, 155η > η, λT E = 1, 155λ > λ
(b) At f = 0, 5fc no propagating mode exists and the guide wavelength has
no physical meaning.
2
kc f 0, 276
ηT M = 1− = −j kc
jωǫ fc fǫ
2 −1
jωµ f 3, 63fc µ
ηT E = 1− = −j
kc fc kc
propagating TE
and TM modes
P TEM mode
ωc
α/kc
f/fc
∞
F (ω) = f (t) exp (−jωt) dt
−∞
and
∞
1
f (t) = F (ω) exp (jωt) dω
2π
−∞
We regard the guide in which this signal propagates as a linear system with
an impulse response H (ω) ⇋ h (t). At the receiving end (load side) of the
guide, we can express the output signal as
Suppose that β (ω) can be expanded in a Taylor’s series about ω0 , valid for
narrow-band signals with ω max << ω0 . Then,
dβ (ω) 1 d2 β (ω)
β (ω) ∼ β (ω 0 ) + |ω=ω0 (ω − ω0 ) + |ω=ω0 (ω − ω 0 )2 + . . .
dω 2 dω 2
4.3 Group velocity 101
β (ω) ∼ β 0 + β ′0 (ω − ω0 )
where
dβ (ω)
β 0 = β (ω 0 ) , β ′0 = |ω=ω0
dω
Letting Ω = ω − ω 0 , we obtain
ωm
1
so (t) = A Re exp {j (ω 0 t − β 0 z)} F (Ω) exp j t − β ′0 z Ω dΩ
2π
−ωm
'
(
= A Re f t − β ′0 z exp {j (ω 0 t − β 0 z)}
= Af t − β ′0 z cos (ω0 t − β 0 z)
This is a time shifted replica of the original modulated signal envelope and
the velocity of the envelope, i.e., velocity of a constant amplitude point on
the wavefront, is defined as the group velocity,
−1
1 dβ
vg = = |ω=ω0
β ′0 dω
d2
Ez (y) + kc2 Ez (y) = 0
dy 2
Ez (y) = 0 @ y = 0, b
Thus,
nπ
Ez (y) = An sin y, n = 0, 1, 2, . . .
b
104 5. Waveguide modes in the parallel-plate transmission line
x z
where
nπ
kc = , n = 0, 1, 2, . . .
b
Since
−
→ γ− →
E t = − 2 ∇ t Ez
kc
and
−
→ jωǫ −
→
H t = 2 z × ∇ t Ez
kc
then the only other non-zero components are,
jωǫ nπ
Hx (y) = An cos y, n = 0, 1, 2, . . .
kc b
and
γ nπ
Ey (y) = − A cos y, n = 0, 1, 2, . . .
kc b
since
∂Ez
=0
∂x
The propagation constant is expressed as
nπ 2
γT M = − ω 2 µǫ
b
and the cut-off frequencies for the different TM modes are
n nc
fc,T M,n = √ =
2b µǫ 2b
We note that TM0 for n = 0 exhibits no cut-off characteristics and is
therefore TEM as a consequence of the two conductors in the structure
5.1 Waveguide modes in a parallel-plate transmission line 105
which allow an axial current flow. The mode having the lowest cut-off
frequency is called the dominant mode and is the TEM mode in this case.
The first TM waveguide mode in the parallel-plate guide is the TM1 mode
and starts to propagate when the operating frequency is f > c/2b.
Example:
(a) Write the instantaneous field expressions for the TM 1 mode in a parallel-
plate waveguide.
(b) Sketch the electric and magnetic field lines in the yz-plane.
(a) We have
π π
Ez (y, z, t) = Re A1 sin y exp (−jωt − jβz) = A1 sin y cos (ωt − βz)
b b
b π
Ey (y, z, t) = β A1 cos y sin (ωt − βz)
π b
b π
Hx (y, z, t) = ωǫ A1 cos y sin (ωt − βz)
π b
where π 2
β = ω 2 µǫ −
b
−
→
(b) In the yz-plane, E has both y- and z-components. At a given time t,
the slope of the electric field lines is given by
dy Ey
=
dz Ez
At t = 0 say,
dy Ey (y, z, 0) b π
= = −β cot y tan βz
dz Ez (y, z, 0) π b
We can re-write this in the form
sin (πy/b) b sin βz
dy = −β dz
cos (πy/b) π cos βz
or
d (cos (πy/b)) d (cos βz)
− = dz
cos (πy/b) cos βz
which on integration yields
π
− ln cos y = ln [cos βz] + const
b
or
π
cos y cos βz = const, 0 ≤ y ≤ b
b
106 5. Waveguide modes in the parallel-plate transmission line
y/b
βz
x π/2 π 3π/2
5.1.2 TE modes
We have Ez , Hz = 0. We can then write
Hz (y, z) = Hz (y) exp (−γz)
5.1 Waveguide modes in a parallel-plate transmission line 107
y/b
βz
x
π/2 π 3π/2
or
π
sin y cos βz = const, 0 ≤ y ≤ b
b
with the constant lying between ∓1. Fig. 5.3 illustrates the field lines.
and
′
Uave = {ue + um } dS
S
and
−
→ 1 ωǫβ 2 nπ
z · P ave = − Re [Ey Hx∗ ] = 2
An cos2 y
2 2kc b
b
−
→ ωǫβb 2
Pz,ave = P ave · zdy = A
4kc2 n
0
−
ǫ → − → ǫ nπ β2 nπ
ue = Re E · E ∗ = A2n sin2 y + 2 cos2 y
4 4 b kc b
and
−
µ → − → µ ω2 ǫ2 nπ
um = Re H · H ∗ = A2n cos2 y
4 4 kc2 b
Hence
b b
ǫ 2 2 nπ β2 2 nπ
ue dy = A sin y + 2 cos y dy
4 n b kc b
0 0
2
2
ǫb 2 β ǫb k
= An 1 + = A2n
8 kc 8 kc
and
b b
µ ω 2 ǫ2 nπ
um dy = A2n cos2 ydy
4 kc2 b
0 0
2
µb 2 2
2 ǫb k
= 2
ω ǫ An = A2n
8kc 8 kc
Thus
2
ωβ ω β fc
ven = 2 = = vp 1 −
k k k f
110 5. Waveguide modes in the parallel-plate transmission line
5.3.2 TM modes
We have,
nπ 2 ǫ′′
γ = − ωµǫ′
1−j ′
b ǫ
σd nπ 2
= j ω 2 µǫ′ 1 − j ′ −
ωǫ b
nπ 2 nπ 2 −1
∼ j ω 2 µǫ′ − 1 − jωµσd ω2 µǫ′ −
b b
or −1
nπ 2 jωµσ d 2 ′ nπ 2
γ∼j ω 2 µǫ′ − 1− ω µǫ −
b 2 b
where we have assumed
nπ 2
ωµσ d << ω2 µǫ′ −
b
5.3 Attenuation in the parallel-plate the guide 111
We have
nπ √
= 2πfc µǫ
b
Hence
nπ 2
2
fc
ω 2 µǫ′ − = ω µǫ′ 1−
b f
Thus
2
σd µ 1 fc
γ = αd + jβ = ′
+ jω µǫ′ 1−
2 ǫ 2 f
1 − (fc /f)
or
σdη
αd = , N pm−1
2
2 1 − (fc /f )
and
2
fc
β = ω µǫ′ 1− , rads−1
f
Thus, for TM modes, αd decreases with increasing frequency.
From the principle of energy conservation,
Pl (z)
αc =
2P (z)
b
1
P (z) = W − Ey Hx∗ dy
2
0
2 b
2
ωǫ′ β bAn nπ bAn
= W cos2 ydy = W ωǫ′ β
2 nπ b 2nπ
0
and hence
2ωǫ′ Rm 2Rm
αc = =
βb
ηb 1 − (fc /f )2
With
πf µc
Rm =
σc
then
2 −1/2
2 πf µc fc
αc = 1−
ηb σc f
2 −1/2
2 πfc µc fc fc
= 1−
ηb σc f f
5.3.3 TE modes
The propagation constant is the same as for the TM modes so that αd has
the same value. For αc , we note that
b
1
P (z) = W Ex Hy∗ dy
2
0
2 b
2
W ωµβ bBn nπ bBn
= sin2 ydy = W ωµβ
2 nπ b 2nπ
0
and
1
Pl (z) = 2W |Js,x |2 Rm = W |Hz (y = 0)|2 Rm = W Bn2 Rm
2
Thus,
2
Pl (z) 2Rm nπ 2 2Rm kc 2Rm (fc /f)2
αc = = = =
P (z) ωµβb b ηb k
ηb 1 − (fc /f)2
6.1.1 TM or E Modes
We assume the guide walls are perfectly conducting and that the enclosed
dielectric has parameters ǫ, µ. For a wave propagating in the +z-direction
with a propagation factor of the from exp (jωt − γz), the axial component
of the electric field in a charge-free region in the dielectric satisfies the
114 6. Rectangular section waveguides
b
y
x z
or
∇2t Ez + kc2 Ez = 0, kc2 = γ 2 + k2
with the boundary conditions that the tangential electric field components
vanish at the conducting walls.
We assume that Ez can be expressed as
Ez = X (x) Y (y) = 0
Hence
1 d2 X (x) 1 d2 Y (y)
+ = −kc2
X (x) dx2 Y (y) dy 2
This means that each of the terms must be separately equal to a constant,
say
1 d2 X (x) 2 1 d2 Y (x)
= −k x , = −ky2
X (x) dx2 Y (x) dx2
such that the separation constants kx and ky satisfy
for the propagating modes. For a given TMmn mode, the cut-off frequency
is given by
c mπ 2 nπ 2 c m 2 n 2
fc,mn = + = +
2π a b 2 a b
116 6. Rectangular section waveguides
or alternatively,
2
λc,mn =
π (m/a)2 + (n/b)2
Neither m nor n can be zero and the dominant TM mode is hence the
TM11 .
Example: Write the instantaneous field expressions for the TM 11 mode
hence sketch the electric and magnetic field lines in a typical xy-plane and
a typical yz-plane in a waveguide with cross-section dimensions a, b.
We have
β π π π
Ex (x, y, z, t) = E11 cos x sin y sin (ωt − βz)
kc2 a a b
β π π π
Ey (x, y, z, t) = 2
E11 sin x cos y sin (ωt − βz)
kc b a b
π π
Ez (x, y, z, t) = E11 sin x sin y cos (ωt − βz)
a b
ωǫ π π π
Hx (x, y, z, t) = − 2 E11 sin x cos y sin (ωt − βz)
kc b a b
ωǫ π π π
Hy (x, y, z, t) = 2 E11 cos x sin y sin (ωt − βz)
kc a a b
where
mπ 2 nπ 2
β = k2 − kc2 = ω2 µǫ − +
a b
In a typical xy-plane, the slopes of the electric and magnetic field lines
are given by
dy a π π
= tan x cot y
dx E b a b
dy b π π
= − cot x tan y
dx H a a b
and we note that
dy dy
= −1
dx E dx H
−
→ −
→
i.e. the E and H lines are orthogonal everywhere in the xy-plane and the
−
→ −
→
E lines are normal and the H lines parallel to the conducting guide walls.
In a typical yz plane, say at x = a/2 or sin (πx/a) = 1 and cos (πx/a) =
0, we have
dy β π π
= 2 cot x tan (ωt − βz)
dz E kc b a
−
→
H has only an x-component. Fig. ?? illustrates the field lines of the TM 11
mode.
6.1 Rectangular section waveguides 117
y/b y/b
x/a βz
π/2 π 3π/2 2π
x/a
6.1.2 TE or H modes
Here, we have
∇2t Hz + γ 2 + k2 Hz = 0
with boundary conditions
∂Hz
= 0, @ x = 0, a ⇒ Ey = 0 @ x = 0, a
∂x
and
∂Hz
= 0, @ y = 0, b ⇒ Ex = 0 @ y = 0, b
∂y
−
→
since tangential components of E on the walls must vanish.
A solution of the homogeneous Helmholtz equation is hence readily seen
to be
mπ nπ
Hz = Hmn cos x cos y, m = 0, 1, 2, . . . , n = 0, 1, 2, . . .
a b
The transverse field components are from
−
→ γ− → γ ∂ γ ∂
H t = − 2 ∇ t Hz = − 2 x
Hz − 2 y Hz
kc kc ∂x kc ∂y
and
−
→ jωµ −
→ jωµ ∂ jωµ ∂
E t = 2 z × ∇ t Hz = − 2 x
Hz + 2 y Hz
kc kc ∂y kc ∂x
Hence,
jωµ nπ mπ nπ
Ex = 2
Hmn cos x sin y
kc b a b
jωµ mπ mπ nπ
Ey = − 2
Hmn sin x cos y
kc a a b
γ mπ mπ nπ
Hx (x, y) = 2 Hmn sin x cos y
kc a a b
γ nπ mπ nπ
Hy (x, y) = 2 Hmn cos x sin y
kc b a b
118 6. Rectangular section waveguides
Ex (x, y, z, t) = 0
ωµ π π
Ey (x, y, t) = 2 H10 sin x sin (ωt − βz)
kc a a
β π π
Hx (x, y, t) = − 2
H10 sin x sin (ωt − βz)
kc a a
Hy (x, y, t) = 0
where π 2
β= ω 2 µǫ −
a
Only Ey , Hx , Hz are non-zero with Ey and Hx independent of the y-
dimension. At x = a/2 say, sin (πx/a) = 1 and cos (πx/a) = 0 so that
only Ey and Hx varying sinusoidally with z are non-zero. Hence, at t = 0,
the electric and magnetic field lines in the xy-plane are as shown in Fig.
6.3.
In the xz-plane, all the field components are non-zero and the slope of
−
→ −
→
the H and E lines are governed by
dx β π π
= 2 tan x tan βz
dz H kc a a
Thus, the field lines are independent of y. The surface current density on
−
→
the guide walls is related to H by
−
→ −
→
Js=n
×H
6.1 Rectangular section waveguides 119
y/b y/b
x/a βz
π/2 π 3π/2
x/a
βz
π/2 π 3π/2 2π
FIGURE 6.3. The TE10 electric and magnetic field lines in a rectangular
waveguide
−
→
where n is the outward normal from the inner surfaces of the walls and H
the field intensity at the walls.
Hence, at t = 0,
−
→
J s (x = 0) = −yHz (0, y, z, 0) = −
y H10 cos βz
−
→ −
→
J s (x = a) = yHz (a, y, z, 0) = −
y H10 cos βz = J s (x = 0)
−
→
J s (y = 0) = x
Hz (x, 0, z, 0) − zHx (x, 0, z, 0)
π β π π
H10 cos x cos βz − z 2 H10 sin x sin βz
= x
a kc a a
and
−
→ −
→
J s (y = b) = − J s (y = 0)
Fig. 6.4 shows the surface currents on the guide walls for the TE 10 in a
rectangular waveguide.
Example: Standard air-filled waveguides have been designed for the radar
bands. One type designated the WG-16 is suitable for X-band applications
and has dimensions, a = 2, 29 cm, b = 1, 02 cm. If it is desired that a
WG-16 guide operate only in the dominant TE 10 mode and the operating
frequency be at least 25% above the cut-off frequency of this mode but no
higher than 95% of the next higher cut-off frequency, what is the allowable
operating frequency range?
120 6. Rectangular section waveguides
βz
π/2 π 3π/2 2π
FIGURE 6.4. Surface current distribution for the TE10 mode in a rectangular
waveguide
c 30 GHzcm
fc,10 = = = 6, 55 GHz
2a 2 × 2, 29 cm
30 GHzcm
fc,20 = = 13, 10 GHz
2, 29 cm
Thus, the range is
1, 25fc,10 ≤ f ≤ 0, 95fc,20
or
8, 19 GHz ≤ f ≤ 12, 45 GHz
Now
a b a b
2 mπ nπ mπ nπ
sin x cos2 y dx dy = cos2 x sin2 y dx dy
a b a b
0 0 0 0
ab/4, m=0, n=0
= {ab/2, m or n=0
6.1 Rectangular section waveguides 121
Thus
2
2 ab β mn mπ 2 nπ 2
Pmn = |Hmn | 2
ηT Emn +
ǫ0m ǫ0n kc,mn a b
2
2 ab β mn
= |Hmn | ηT Emn
2ǫ0m ǫ0n kc,mn
ǫ0m = 1, m = 0
= 2, m > 0
and
ǫ0n = 1, n = 0
= 2, n > 0
P = Pmn + Plp
since
a
mπ lπ
sin x sin x dx = 0, l = m
a a
0
For low-loss waveguides, the power orthogonality will still be valid except
for the case when two or more degenerate modes are present. The presence
of low losses may also result in strong coupling between the degenerate
modes.
For the conduction losses, we note that only Ey , Hx and Hz are non-zero
for the TE10 mode so that
b a
1
P (z) = − Ey Hx∗ dxdy
2
0 0
a 2 b a
1 π
= ωµβ H10 sin2 xdxdy
2 π a
0 0
a 2
= ωµβ H10 ab
2π
The power loss Pl (z) is computed by noting that we have four walls to
consider with
−
→ −
→
J s (x = 0) = J s (x = a) = −
yHz (x = 0) = −
y H10
and
−
→ −
→ π βa π
J s (y = 0) = − J s (y = b) = x
H10 cos x − z H10 sin x
a π a
Thus
Pl (z) = 2Pl (z) |x=0 + 2Pl (z) |y=0
where
b 2
1 −
→ b 2
Pl (z) |x=0 = J s (x = 0) Rs dy = H10 Rm
2 2
0
a 2 − 2
1 −→ →
Pl (z) |y=0 = J s,x (y = 0) + J s,x (y = 0) Rm dx
2
0
2
a βa 2
= 1+ H10 Rm
4 π
Thus,
2
a βa 2
Pl (z) = b+ 1+ H10 Rm
2 π
2
a f 2
= b+ 1+ H10 Rm
2 fc
6.1 Rectangular section waveguides 123
and
1 + (2b/a) (fc /f )
αc,T E10 = Rm
ηb 1 − (fc /f)2
2
1 πµc 2b f
= 1+
ηb σc 1 − (fc /f )2 π fc
with fc,11 = 16, 10 GHz and fc,10 = 6, 55 GHz. The curves for αc,T M11
and αc,T E10 show a broad minimum in the operating range f > fc with
αc,T E10 < αc,T M11 .
Theoretical formulae for attenuation are in good agreement with mea-
sured values for frequencies up to ∼ 5 GHz. At higher frequencies however,
measured values are generally higher depending on the surface smoothness
of the guide walls. When the roughness is on the order of the skin depth at
the operating frequency, the effective surface area is increased thus, increas-
ing the losses. This may explain the higher measured attenuation values.
Example: A TE 10 wave at 10 GHz propagates in a brass rectangular
waveguide with inner dimensions a = 1, 5 cm, b = 0, 6 cm. If the dielectric
is polyethylene with ǫ = 2, 25, µr = 1 and loss tangent 4 × 10−4 , determine
(a)
2
ω fc
β= 1− = 74, 5π = 234 radm−1
v f
λ
λg = = 0, 0268 m
1 − (fc /f )2
(b)
v
vp = = 2, 68 × 108 ms−1
2
1 − (fc /f )
1 µ
ηT E10 = = 337, 4 Ω
ǫ
1 − (fc /f)2
π 2 2 π 2
2πf
β 10 = k2 − = − = 158, 05 m−1
a c a
6.1 Rectangular section waveguides 125
b d jBC b jBL
d
a a
(a) CAPACITIVE IRIS (b) INDUCTIVE IRIS
d d
b b
a a
(c) PROTRUDING POST (d) TUNING POST OR SCREW
σ Cu = 5, 8 × 107 S/m
ωµ
Rm = = 0, 026 Ω
2σCu
Rm
α= 3 µ = 26π 2 + a3 k2 = 0, 0125 N p/m
a bβk ǫ
attenuation
−20 log10 exp (−αl) = 0, 11 dB
λg πd
bL = BL ηT E10 = − cot2
a 2π
Another discontinuity that provides shunt susceptance is a conducting
post protruding into the waveguide on the broad face as shown in Fig.
6.5(c). If d is small, the susceptance is capacitive. When d becomes appre-
ciable compared to b, considerable current flows along the post causing an
inductive effect. Resonance occurs near d = 3b/4 and for d longer than this,
the susceptance becomes inductive.
In practice, the post is a screw installed in an axial slit in the centre of
the broad face where H is minimum and does not therefore disturb the
field pattern significantly. By adjusting d and the position along the guide,
impedance matching similar to stub matching or tuning can be provided.
Example: Measurements on a WG-10 S-band waveguide ( a = 7, 21 cm,
b = 3, 40 cm) feeding a horn antenna show a VSWR of 2, 0 at the operating
frequency of 3 GHz and existence of an electric field maximum 12 cm from
the neck of the horn. Find the location and dimensions of a symmetrical
inductive iris necessary to achieve a perfect match. Assume the guide to be
lossless.
c
fc,10 = = 2, 08 GHz
2a
λ c
λg = = = 13, 9 cm
1 − (fc /f )2 f − fc2
2
12 cm
= 0, 863λg
13, 9 cm
RL
rL = =S
Z0
Hence
YL 1
gL = = = 0, 5
Y0 S
Proceed as in the single-stub tuning scheme to obtain d = 4, 72 cm located
at (0, 863 − 0, 348) λg = 7, 16 cm.
6.2 Power capacity of guiding structures 127
A higher Pmax requires the use of a larger coaxial line while keeping the
ratio b/a fixed to ensure the same characteristic impedance. However a
and b determine the highest frequency at which the line can be operated
mode-free. Thus, for a given maximum operating frequency fmax ,
2
2
0, 025 cEd Ed
Pmax = = 5, 8 × 1012
µ0 /ǫ0 fmax fmax
ab E02 ab Ed2
Pmax = =
4 Z 4 Z
where Z is the wave impedance of the guide. Many rectangular waveguides
are constructed for single-mode operation with a ∼ 2b to avoid propagation
of the TE20 , which is the next higher mode. Thus, we must have
c
a<
fmax
128 6. Rectangular section waveguides
so that
2
2
0, 11 cEd 13 Ed
Pmax = = 2, 6 × 10
µ0 /ǫ0 fmax fmax
Arcing and voltage breakdown are high frequency processes so that Pmax
represents the peak value. Hence, in constructing a guiding structure, pro-
vide at least a safety factor of 2 so that practical Pmax values are approxi-
mately half the theoretical ones. In case the reflection coefficient becomes
unity, the maximum guide voltage will be twice the value in a reflection-
free case. Hence Pmax allowable is reduced by a factor of 4. Power handling
capacity can be improved by using pressurised air, inert gases or using a
dielectric with a higher dielectric strength.
kc2
Z ′ = jωµ +
jωǫ
Y ′ = jωǫ
k2c
jωµ jωε
jωε
kc2 = ω2c µǫ
where ωc is the radian cutoff frequency from wave theory. The characteristic
impedance for the TM equivalent circuit is given by
Z′ jωµ + kc2 /jωǫ µ ω 2
c
Z0 = = = 1 − = ZT M
Y′ jωǫ ǫ ω
130 6. Rectangular section waveguides
jωµ
k2c
jωε
jωµ
Z0 Z0
Z0 Z0
Z01, Z02,
γ1 γ2
FIGURE 6.8. Example field problems where transmission line analogy proves
useful
132 6. Rectangular section waveguides
ε ε0
a x
t
Partially filled waveguide
the boundary conditions, the current associated with the dominant mode is
discontinuous by the amount contributed by the higher order modes. This
is what is accounted for by the shunt reactance connected at the disconti-
nuity.
y y
ky,ZTM
b b
ε ,µ Zin,l
y
Zin,g
x x
a Transverse-resonance
equivalent circuit
and
Zin,l (y) = jηT M tan ky y
Thus
jηT M {tan ky (b − y) + tan ky y} = 0
or
sin ky b
jηT M =0
cos ky (b − y) cos ky y
6.3 Some useful ideas for waveguide analysis 135
y y
b b
ky,a,Z0,a
ε0
t t
ε ky,d,Z0,d
x x
Partially filled waveguide Transverse-resonance
equivalent circuit
Hence
nπ
kc = ky = , n = 0, 1, 2, . . .
b
For a partially loaded rectangular waveguide, the situation for the TE0n
modes can be represented as shown in Fig. 6.11. In the region 0 < y < t,
we have β = ky,d and
k k0
Z0,d = η = η0
ky,d ky,d
where
√ µ0
k0 = ω µ0 ǫ0 , η 0 =
ǫ0
k0
Z0,a = η0
ky,a
Thus
ky,a {tan ky,d t + ky,d tan ky,a (b − t)} = 0
and together with the phase matching condition of the tangential field
components
β 2 = k2 − ky,d
2
= k02 − ky,a
2
Ey −ηT E,mn Hx − jβ
bk2
mn n
sin mπx nπx
a cos b exp (−jβ mn z)
c,mn
k0
ηT E,mn β mn Z0
β mn
ηT M,mn k0 Z0
2
m 2
kc,mn π a + nb
λc,mn √ 2ab 2 2
(mb) +(na)
β mn 2 2
k0 − kc,mn
2Rm m2 b3 + n2 a3
αT M =
2 2
2
bZ0 1 − kc,mn /k0 m b a + n2 a3
2
where
ωµ0
Rm =
2σ
and
ǫ0m = 1, m = 0
= 2, m > 2
z
d
Consider the rectangular cavity of Fig. 6.12. We choose the z-axis as the
reference direction of propagation but because of the conducting termina-
tion at z = 0, d, only a standing wave pattern exists in the cavity. We use
three indices m, n, p to denote the TM or TE standing wave pattern in the
cavity.
1 mπ mπ nπ pπ
Ey = − Emn cos x sin y sin z
kc2 a a b d
where mπ 2 nπ 2
kc2 = +
a b
For the TMmnp modes, resonance frequency of the cavity is given by
1 mπ 2 nπ 2 pπ 2
ω mnp = + + = 2πfmnp
µǫ a b d
maximum Hx inside the front or back faces as shown in Fig. 6.12. The
best position for the probe or loop is affected by the impedance matching
requirements of the microwave circuit of which the cavity is part of. A
common method of coupling energy out of the cavity is by introducing a
hole or iris at an appropriate location in the cavity wall. The field at the
hole must have a component favourable for exciting the desired mode in
the cavity.
Example: Determine the dominant modes and their frequencies in an air-
filled rectangular cavity resonator for a > b > d, a > d > b, a = b = d,
where a, b, and d are the x-, y-, z-directed dimensions of the cavity resp.
With the z-axis chosen as the reference direction of propagation, neither
m nor n can be zero for TM modes but TE modes allow either m or n
to be zero but not simultaneously but p cannot be zero. Thus, the lowest
order modes are TM 110 , TE 011 and TE 101 with TM 110 having the lowest
resonance frequency for a > b > d given by
c 1 1
f110 = +
2 a2 b2
TE 101 will have the lowest resonance frequency for a > d > b given by
c 1 1
f101 = +
2 a2 d2
while all three modes will have the same resonance frequency for a = b = d
given by
c
f110 = √
a 2
U
Q=ω
Pl
140 6. Rectangular section waveguides
and
πf101 µ0 abd a2 + d2
Q=
Rm {2b (a3 + d3 ) + ad (a2 + d2 )}
Example: What should be the size of a hollow cubic cavity made of copper
in order for it to have a dominant mode resonance frequency of 10 GHz?
Proceed to determine the Q of the resonator at this frequency.
30 GHzcm
f101 = √ = 10 GHz ⇒ a = 21, 2 mm
2a cm
πf101 µ0 a a
Q= = πf101 µ0 σ Cu = 10700, σCu = 5, 8 × 107 Sm−1
3Rs 3
142 6. Rectangular section waveguides
7
Circular section waveguides
and
1 ∂ ∂ ∂2
r Hz + r2 2 Hz + kc2 Hz = 0
r ∂r ∂r ∂φ
where
kc2 = γ 2 + ω 2 µǫ
7.1.1 TM modes
Consider the conducting hollow pipe of inner radius a shown in Fig. 7.1.
We propose that Ez (r, φ) can be expressed as a product of analytic one-
dimensional functions as,
Ez (r, φ) = R (r) Φ (φ)
Thus
r d dR 1 d2 Φ
r + kc2 r2 = −
R (r) dr dr Φ (φ) dφ2
Let ν = 0 be a constant such that
1 d2 Φ
= −ν 2
Φ (φ) dφ2
called the angular differential equation and
r d dR
r + kc2 r2 = ν 2
R (r) dr dr
or
d2 R dR 2 2
r2 2
+r + kc r − ν 2 R (r) = 0
dr dr
7.1 Circular section waveguides 145
or
∞
/ ∞
/
Cl (p + l)2 − ν 2 xp+l + Cl xp+l+2 = 0
l=0 l=0
or
Cl−2
Cl = − , l = 2, 3, 4, . . .
(p + l)2 − ν 2
We note that with C1 = 0, Cl = 0 for l odd. We have three possible
solutions depending on the nature of p1 = ν and p2 = −ν,
146 7. Circular section waveguides
If ν is not an integer, Rν (x) and R−ν (x) represent two linearly independent
solutions of the second-order Bessel’s equation. A general solution can the
be expressed as
R (x) = A1 Jν (x) + A2 J−ν (x)
7.1 Circular section waveguides 147
Replacing n by −n yields,
∞
/ (−1)m x 2m−n
J−n (x) = , n = 0, 1, 2, . . .
m=0
m!Γ (m − n + 1) 2
∞
/ (−1)m x 2m−n
= , n = 0, 1, 2, . . .
m=n
m!Γ (m − n + 1) 2
since
1
= 0, m = 0, 1, 2, . . . n − 1
Γ (m − n + 1)
Thus, putting k = m − n, we obtain
∞
/ (−1)k+n x 2k+n
J−n (x) = , n = 0, 1, 2, . . .
(k + n)!Γ (k + 1) 2
k=0
∞
/ (−1)k x 2k+n
= (−1)n , n = 0, 1, 2, . . .
k!Γ (k + n + 1) 2
k
= (−1)n Jn (x)
and now, Jn (x) and J−n (x) are linearly dependent so that we have only
one solution of the second-order differential equation. To generate a second
solution, we define a function
∆ Jv (x) cos vπ − J−v (x)
Yv (x) =
sin vπ
called the Bessel function of the second kind and order v.
Yv (x) is linearly independent of Jv (x) when ν is not an integer. For
integer ν, Yn (x) reduces to the indeterminate ratio 0/0. We thus define
Yn (x) as a limit as
Yn (x) = lim Yν (x) , n = 0, 1, 2, . . .
ν→n
148 7. Circular section waveguides
and use L’Hospital’s rule to evaluate this limit since, for two functions f (x)
and g (x) for which
f (x) 0
lim =
x→x0 g (x) 0
then
f (x) f ′ (x)
lim = lim ′
x→x0 g (x) x→x0 g (x)
For n = 0, we have
x 2/
(−1) x 2k
∞ k
2 1 1 1
Y0 (x) = J0 (x) ln + γ − 1 + + + . . . +
π 2 π
k=1
(k!)2 2 2 3 k
where
n
/ 1
γ = lim − ln n = 0, 577215 . . .
n→∞ k
k=1
is Euler’s constant.
For small x,
2
ln x, x → 0+
Y0 (x) →
π
The asymptotic behaviour near x = 0 is given by
ν
Γ (ν + 1) 2
Yν (x) ∼ − , ν > 0, x → 0+
π x
and unbounded at the origin. Thus, for problems involving the origin, Yν (x)
is not a valid solution and must be excluded. For arbitrary ν, we can write
Fig. 7.2 illustrates the Bessel functions of the first kind for the first three
integer orders.
The choice of the value of C0 arose from earlier studies of Bessel func-
tions in which
Jn (x) was taken as the coefficient of tn in the expansion of
exp x2 t − 1t , i.e.
n
1 x n
/∞
x 1 1
exp t− = t−
2 t n=0
n! 2 t
/ 1 x n
∞
n
1 1
= n
C0 tn −n C1 tn−1 + . . . + (−1)n n Cn
n=0
n! 2 t t
x 1
The function exp 2 t− t is known as the generating function of the
first order.
7.1 Circular section waveguides 149
Jn(kca)
1,0
J0(kca)
J1(kca)
J2(kca)
kc a
0
FIGURE 7.2. Bessel functions of the first kind for the first few orders
Some problems are simplified when Jν (kc r) and Yν (kc r) are linearly
combined to form the Hankel functions of the first and second kind and
order ν defined as
and
Hν(2) (kc r) = Jν (kc r) − jYν (kc r)
respectively. For kc r → ∞,
2 π νπ
Hν(1) (kc r) → exp j kc r − −
πkc r 4 2
and
2 π νπ
Hν(2) (kc r) → exp −j kc r − −
πkc r 4 2
The large-argument approximations of the Hankel functions are useful in
(1)
the study of radiation and scattering problems. Note that Hν (kc r) for
large kc r is a spherical wave propagating inward toward the origin while
(2)
Hν (kc r) is one propagating outward from the origin.
When kc2 is negative, the solutions take the form
Iν (ξr) = j −ν Jν (jξr)
150 7. Circular section waveguides
and
π
Kν (ξr) = j ν+1 Hν(1) (jξr)
2
where jξ = kc . These are the modified Bessel functions of the first and
second kind and order ν. For ξr → ∞, these have the asymptotic approxi-
mations
1
Iν (ξr) → exp (ξr)
2πξr
and
π
Kν (ξr) → exp (−ξr)
2πξr
Kν (ξr) for large r therefore exhibits an exponential decay which is a char-
acteristic of evanescent waves. This function is useful in the study of surface
waves in rod dielectric waveguides.
For the TM modes, we have
where
Ez (r, φ) = R (r) Φ (φ)
Since all field components are periodic in φ with period 2π, the only ad-
missible solution is either sin nφ or cos nφ or a linear combination of these
in the expression of Φ (φ). We may chose either since this only translates
to a shift in the choice of the φ = 0 reference point.
Let Ez (r, φ) = An Jn (kc r) cos nφ. Then
−
→ φ=−γ − →
E t,T M = rEr + φE ∇ t Ez
kc2
where
−
→ ∂ ∂
φ
∇ t Ez = r + Ez
∂r r ∂φ
Thus
jβ
Er = − An Jn′ (kc r) cos nφ
kc
jβn
Eφ = 2 An Jn (kc r) sin nφ
kc r
jωǫn
Hr = − An Jn (kc r) sin nφ
kc2 r
jωǫ
Hφ = − An Jn′ (kc r) cos nφ
kc
where
d
Jn′ (kc r) = Jn (kc r)
dr
7.1 Circular section waveguides 151
TABLE 7.1. The first few zeros of the Bessel’s function of the first kind and order
n
n p1n p2n p3n
0 2, 405 5, 520 8, 654
1 3, 832 7, 016 10, 174
2 5, 135 8, 417 11, 620
7.1.2 TE modes
Hz (r, φ, z) = Hz (r, φ) exp (−γz)
where
∇2t Hz + kc2 Hz = 0
152 7. Circular section waveguides
FIGURE 7.3. TM01 electric and magnetic field lines in a hollow circular waveguide
with
∂
Jn (kc r) |r=a = Jn′ (kc a) = 0
∂r
Let the roots of this equation be p′mn . Thus
p′mn
kc,mn =
a
and the first few roots are presented in Table 7.2
Now
d
J0 (x) = −J1 (x)
dx
thus p′0n = p1n and T E0n and T M1n are therefore degenerate.
7.1 Circular section waveguides 153
TABLE 7.2. The first few zeros of the derivative of Bessel functions of the first
kind and order n
m p′m1 p′m2 p′m3
0 3, 832 7, 016 10, 174
1 1, 841 5, 331 8, 536
2 3, 054 6, 706 9, 970
We see that the dielectric loss dominates. The total attenuation is then
α = αc + αd = 0, 547 N p/m
and loss
−20 log10 [exp (−αl)] = 2, 38 dB
7.1 Circular section waveguides 155
TM (E) Modes
TM (E) Modes
Hz 0
Hr −Eφ /ηT M,mn
Hφ Er /η
T M,mn
−jβ z
Ez Jn pmn a
r
e mn (cos nφ or sin nφ)
jβ mn pmn ′ pmn r
−jβ mn z
Er − ak2 Jn a e (cos nφ or sin nφ)
c,mn
jβ mn pmn pmn r
−jβ z
Eφ − rk2 nJn a e mn (cos nφ or sin nφ)
c,mn
β mn
ηT M,mn k0 η 0
pmn
kc,mn a
2πa
λc,mn pmn
k0 β mn π 2 ′2 ′
power 4
2η0 kc,mn ǫom pmn Jm (pmn )
α √ 2Rm2
aZ0 1−kc,mn /k02
where
ωµ0 µ0
Rm = , η0 =
2σ ǫ0
and
ǫ0m = 1, m = 0
= 2, m > 2
156 7. Circular section waveguides
cavity of radius a and length d, the energy stored in the cavity for the
TM010 mode is given by
a
ǫ0 2 ǫ0 C02 2, 405
U = 2Ue = |Ez | dV = (2πd) J02 r rdr
2 2 a
V 0
2
a 2
= πǫ0 dC02 J (2, 405)
2 1
The time-average power loss due to dissipation in the wall arising from a
finite conductivity is given by
1 1
Pave = |Js |2 Rm = |Hφ |2 Rm
2 2
−
→
There are radial currents J r on the face plates and uniform axial surface
−
→
currents J z on the inside of the cylindrical wall. Thus,
a
Rm
Pl = 2 |Jr |2 (2πr) dr + (2πad) |Jz |2
2
0
a
2 2
= πRm 2 |Hφ | rdr + ad |Hφ (r = a)|
0
a
πRm 2 2 2, 405 2
= C 0 2 J1 r rdr + adJ1 (2, 405)
η20 a
0
or
πRm 2
Pl = C (a + d) J12 (2, 405)
η20 0
so that
η0 2, 405
Q=
Rs 2 (1 + a/d)
Note that
r2 n2
Jn2 (kc r) rdr = ′
Jn (kc r) + 1 − 2
Jn (kc r)
2 (kc r)2
1
J1′ (kc r) = J0 (kc r) − J1 (kc r) , J0 (kc a) = 0
kc r
With
πf µc
Rm =
σ
160 7. Circular section waveguides
2, 405 0, 115
fT M,010 = √ = GHz
2πa µ0 ǫ0 a
Example: A hollow circular cylindrical cavity resonator is constructed of
copper such that its length d equals its diameter 2a. Determine a and d
for a resonance frequency of 10 GHz for the TM 010 mode hence estimate
its Q.
We have
2, 405 0, 115
fT M,010 = √ = GHz = 10 GHz
2πa µ0 ǫ0 a
Hence
a = 1, 15 cm, d = 2a = 2, 30 cm
πfµ0
Rm = = 2, 61 × 10−2 Ω
σ
so that
η0 2, 405
Q= = 11580
Rs 2 (1 + a/d)
Comparing with rectangular and circular resonators of comparable size at
the same frequency, we note
Thus, these have approximately the same volume but the total surface area
of the rectangular cavity is 8, 2% larger leading to a lower Q. The Q of the
circular cavity is ∼ 8, 2% higher than that of the rectangular cavity.
8
Dielectric waveguides
8.1.1 TM modes
We have
d2
Ez (y) + kc2 Ez (y) = 0
dy 2
where
kc = γ 2 + ω 2 µǫ
with γ = jβ inside the slab. Hence, we can write
d
Ez (y) = Eo sin ky y + Ee cos ky y, |y| ≤
2
where Eo and Ee are amplitudes of the odd and even field components
resp. and
ky = ω2 µd ǫd − β 2
162 8. Dielectric waveguides
d/2 µ0,ε0
z
0
µd,εd
-d/2 µ0,ε0
Odd TM modes
(i) In the dielectric region,
Ez (y) = Eo sin ky y
jβ
Ey = − Eo cos ky y
ky
jωǫd
Hx (y) = Eo cos ky y
ky
(ii) In the region y > d/2,
d d
Ez (y) = Eo sin ky exp −α y −
2 2
8.1 Planar dielectric slab in free space 163
jβ d d
Ey = − Eo sin ky exp −α y −
α 2 2
jωǫ0 d d
Hx (y) = Eo sin ky exp −α y −
α 2 2
α ǫ0 ky d
= tan
ky ǫd 2
and
α2 + ky2 = ω 2 {µd ǫd − µ0 ǫ0 }
or
α= ω 2 {µd ǫd − µ0 ǫ0 } − ky2
Ez α
Zm = − =j
Hx ωǫ0
Even TM modes
We have
d
Ez (y) = Ee cos ky y, |y| ≤
2
and invoking the continuity conditions at the dielectric surface, we obtain
α ǫ0 d
= − cot ky
ky ǫd 2
The phase constant of the propagating TM even modes lies between the
intrinsic phase constant of free space and that of the dielectric, i.e.
√ √
ω µ0 ǫ0 < β < ω µd ǫd
√
As β → ω µ0 ǫ0 , α → 0 and the wave is no longer confined within the
dielectric. The limiting frequencies under this condition are known as the
cut-off frequencies of the dielectric waveguide. At cutoff, we have for odd
TM modes,
ω co d
tan µd ǫd − µ0 ǫ0 = 0
2
πfco d µd ǫd − µ0 ǫ0 = (n − 1) π, n = 1, 2, 3, . . .
n−1
fco = √
d µd ǫd − µ0 ǫ0
and
ω ce d
cot µd ǫd − µ0 ǫ0 = 0
2
1
πfce d µd ǫd − µ0 ǫ0 = n − π, n = 1, 2, 3, . . .
2
n − 1/2
fce = √
d µd ǫd − µ0 ǫ0
for even TM modes. Note that fc,o = 0 for n = 1. Hence, the dominant
mode on a slab of any thickness is an odd TM mode.
8.1.2 TE modes
We have
d2
Hz (y) + kc2 Hz (y) = 0
dy 2
with a solution
d
Hz = Ho sin ky y + He cos ky y, |y| ≤
2
8.1 Planar dielectric slab in free space 165
and
d d
Hz (y) = CA exp −α y − , y>
2 2
d d
= CB exp α y + , y<−
2 2
Odd TE modes
We have
(i)
jβ
Hz (y) = Ho sin ky y, Hy (y) = − Ho cos ky y
ky
jωµd
Ex (y) = − Ho cos ky y
ky
inside the dielectric
(ii) For y ≥ d/2
ky d d
Hz (y) = Ho sin exp −α y −
2 2
jβ ky d d
Hy (y) = − Ho sin exp −α y −
α 2 2
jωµd ky d d
Ex (y) = − Ho sin exp −α y −
α 2 2
Thus, at y = d/2,
α µ ky d
= 0 tan
ky µd 2
Ex jωµ0
Zm = =−
Hz α
which is capacitive but increases with increasing frequency. Hence, a TE
surface wave can be supported by a capacitive surface.
166 8. Dielectric waveguides
Even TE modes
d
Hz (y) = He cos ky y, |y| ≤
2
and
α µ ky d
= 0 cot
ky µd 2
Therefore, the lowest order odd TE mode just like the lowest order odd
TM mode has no cutoff frequency.
Example: A dielectric slab waveguide with constitutive parameters µd =
µ0 , ǫd = 2, 50ǫ0 is situated in free space. Determine the minimum thickness
of the slab so that a TM or TE wave of the even type propagates along the
guide at 20 GHz.
Lowest even TM and TE modes have the same cut-off frequency
n − 1/2
fc = √
d µd ǫd − µ0 ǫ0
Thus, for n = 1,
c c
fc = ⇒ dmin = = 6, 12 cm
2d ǫd /ǫ0 − 1 2fc ǫd /ǫ0 − 1
Example:
(a) Obtain an approximate expression for the decaying rate of the dominant
TM surface wave outside a very thin dielectric slab waveguide.
(b) Determine the time-average power per unit slab width transmitted along
the guide
(a) fco = 0 for n = 1 for all d. If ky d/2 << 1, tan ky d/2 ∼ ky d/2. Thus,
for
d ǫd
α <<
2 ǫ0
then
1 ǫ0 2 1 ǫ0 2
α= ky d ∼ ω {µd ǫd − µ0 ǫ0 } d
2 ǫd 2 ǫd
(b)
−
→ 1
P ave = Re [−
y Ey × x
Hx ]
2
8.1 Planar dielectric slab in free space 167
µ0,ε0
θi
θi µd,εd
µ0,ε0
Thus
d/2 d/2
−→ ωǫd β 2
Pave = 2 P ave dy = 2 Eo cos2 ky ydy
ky
0 0
ωǫd β 2 1
= E d+ sin ky d , W m−1
4ky2 o ky
√
ky ∼ ω µd ǫd − µ0 ǫ0 , β ∼ ω µ0 ǫ0
−
→ 1
P ave = Re [
z Ez × x
Hx ] = 0
2
Consider the dielectric slab shown in Fig. 8.2. There will be total internal
reflection of a wave incident obliquely on the lower surface of the slab when
ǫd > ǫ0 and θi > θ i,c , where
ǫ0
θi,c = sin−1
ǫd
and
Ez = A2 exp {−ka (y − a)} , y > a
and since
β k −
→
Hx = −j z × ∇ t Ez · x
kc2 βη
where
k k0
= ǫr = ωǫ, 0 < y < a
η η0
k0
= = ωǫ0 , y > a
η0
then at y = a
ǫr ∂Ez 1 ∂Ez
2 =− 2
kd ∂y ka ∂y
i.e.
ǫr 1
A1 cos kd a = A2
kd ka
or dividing the two transcendental equations and multiplying through by
a, we have
kd a tan kd a = ǫr ka a
When a solution for real ka exists, this will be a surface wave since its
field will decay exponentially in the y-direction while propagating in the
z-direction.
kd
η T M,y = η
k
in the dielectric and
ka
ηT M,y = η
jk0 0
in air. The equivalent circuit is then a short-circuited transmission line of
length a and impedance −jka η0 /k0 . At the reference plane, the resonance
170 8. Dielectric waveguides
condition requires that the sum of the impedances seen looking toward the
short circuit and the input to the infinite line vanish. Hence
kd η ka η0
j tan kd a − j =0
k k0
or
kd tan kd a = ka ǫr
since η/k = η 0 / (k0 ǫr ).
The T M0 mode surface-wave does not have a low-frequency cut-off fre-
quency but this is not a characteristic of all surface-wave guiding structures.
However, all will have only a finite number of propagating modes at any
given frequency. When the mode exists, its phase velocity will be less than
the plane wave velocity in air since
k0 k0
vp = c= 2 c<c
β k0 + ka2
For this reason, a surface-wave is also referred to as a slow wave. The field
components are then
Ez = A sin (kd y) e−jβz
jβ
Ey = − A cos (kd y) e−jβz
kd
k0
Hx = −ǫr A cos (kd y) e−jβz
kd η 0
ωǫ
= − A cos (kd y) e−jβz
kd
8.2 ♣The cylindrical optical fibre 171
k0
Hx = − A sin (kd a) e−ka (y−a)−jβz
ka η0
ωǫ0
= − A sin (kd a) e−ka (y−a)−jβz
ka
in air.
Ez,i = C0 J0 (kc r) , r ≤ a
where
kc2 = γ 2 + kd2 = ω2 µ0 ǫd − β 2
The corresponding Hφ,i is given by
jωǫd
Hφ,i = − C0 J0′ (kc r)
kc
Outside the rod, the fields must be evanescent and we may thus choose
Ez = K0 (ξr)
the modified Bessel function of the second kind and order zero. We write
Ez,0 = D0 K0 (ξr) , r ≥ a
where
ξ 2 = β 2 − k02 = β 2 − ω2 µ0 ǫ0
and
jωǫ0
Hφ,0 = D0 K0′ (ξr) , r ≥ a
ξ
Ez and Hφ must be continuous across the boundary at r = a. Hence
C0 J0 (kc a) = D0 K0 (ξr)
172 8. Dielectric waveguides
and
ǫd ǫ0
C0 J0′ (kc r) = D0 K0′ (ξr)
kc ξ
and combining,
J0 (kc a) ǫd ξK0 (ξr)
′ =
J0 (kc r) ǫ0 kc K0′ (ξr)
with
kc2 + ξ 2 = ω2 µ0 {ǫd − ǫ0 }
Commercially available optical fibres are of two types however,
or
∇2 Ez + k02 n2 − β 2 Ez = 0
where
√
n= ǫr
is the refractive index of a non-magnetic dielectric medium.
8.2 ♣The cylindrical optical fibre 173
An intuitively important point is that Ez does not couple to the other two
components as the electromagnetic field propagates even after reflection at
the cylindrical surfaces. Therefore, since Ez couples only to itself, the wave
equation can be written in cylindrical coordinates as
∂2 1 ∂ 1 ∂2
2
Ez + Ez + 2 2 Ez + k02 n2 − β 2 Ez = 0
∂r r ∂r r ∂φ
Thus, as long as
kc = k02 n2 − β 2 > 0
valid solutions are the Bessel functions of the first kind and order m since
the fields must remain finite at the centre of the core.
When
k02 n2 − β 2 < 0
the valid solution will be the modified Bessel function of the second kind,
Km (ξr) since the fields must vanish at infinity, where
ξ 2 = β 2 − k02 m2
and outside
we can write
j β ∂Ez ∂Hz
Eφ = − − ωµ0
k02 n2 − β 2 r ∂φ ∂r
j ωµ0 ∂Hz ∂Ez
Er = − 2 2 +β
k0 n − β 2 r ∂φ ∂r
j ∂Ez β ∂Hz
Hφ = − 2 ωǫ +
k0 n2 − β 2 ∂r r ∂φ
j ∂Hz ωǫ ∂Ez
Hr = − 2 2 β −
k0 n − β 2 ∂r r ∂φ
In the core region
jβ ′ jmωµ0
Er = − 2 A1 Jm (kc r) + A2 Jm (kc r) exp {j (mφ − βz)}
kc βr
jβ jm jωµ0 ′
Eφ = − 2 A1 Jm (kc r) − A2 kJm (kc r) exp {j (mφ − βz)}
kc r β
jβ ′ jmωǫcore
Hr = − 2 A2 kJm (kc r) + A1 Jm (kc r) exp {j (mφ − βz)}
kc βr
jβ jm jmωǫcore ′
Hφ = − 2 A2 Jm (kc r) + A1 kJm (kc r) exp {j (mφ − βz)}
kc r β
and in the cladding
jβ ′ jmωµ0
Er = 2 B1 γKm (ξr) + B2 Kn (ξr) exp {j (mφ − βz)}
ξ βr
jβ jn jωµ0 ′
Eφ = 2 B1 Km (ξr) − B2 γKm (ξr) exp {j (mφ − βz)}
ξ r β
jβ ′ jmωǫcladd
Hr = 2 B2 γKm (ξr) + B1 Km (ξr) exp {j (mφ − βz)}
ξ βr
jβ jm jmωǫcladd ′
Hφ = 2 B2 Km (ξr) + B1 γKm (ξr) exp {j (mφ − βz)}
ξ r β
and on matching the boundary conditions
Jm (kc a) 0 −Km (γa) 0 A1
0 J (k a) 0 −K m (ξa)
m c A2
=0
β ωµ0 ′ βn ωµ0 ′
akc2 Jm (kc a) j kc Jm (kc a) K (ξa)
aξ2 m
j ξ Km (ξa) B1
ωǫcore ′ β ωǫcladd ′ βm B2
−j kc Jm (kc a) ak2 Jm (kc a) −j ξ Km (ξa) aξ2 Km (ξa)
c
8.3 ♣Dielectric-coated conducting wire 175
In the region r > b, the field must decay with increasing r to correspond to a
surface-wave. We recall that for large ka r we can make the approximations
2 π
J0 (jka r) ∼ cos jka r −
jπka r 4
and
2 π
Y0 (jka r) ∼ sin jka r −
jπka r 4
Hence,
2 π
J0 (jka r) + jY0 (jka r) ∼ 2 exp j jka r −
jπka r 4
In (y) = j −n Jn (jy)
and
π
Kn (y) = j n+1 {Jn (jy) + jYn (jy)}
2
We note that Kn (y) decays exponentially with increasing y making it a
suitable solution for Ez in the region r > b for n = 0. In (y) on the other
hand, increases exponentially and is therefore not valid for physical fields.
Thus,
Ez = A3 K0 (ka r) , r > b
Since Ez = 0 at r = a, (electric wall), then it is necessary that
A1 J0 (kd a) + A2 Y0 (kd a) = 0
or
J0 (kd a)
A2 = −A1
Y0 (kd a)
At r = b, Ez is continuous so that
J0 (kd a)
A3 K0 (ka b) = A1 J0 (kd b) − Y0 (kd b)
Y0 (kd a)
8.3 ♣Dielectric-coated conducting wire 177
we can write
kd2 + ka2 = (ǫr − 1) k02
We invoke the continuity condition of the tangential magnetic field at r = b
so that, in terms of Ez , we have
−
→ jβ −
→
E t = − 2 ∇ t Ez
kc
With
−
→ ∂ ∂
φ ∂
∇ t = r + = r
∂r r ∂φ ∂r
we have
−
→ jβ ∂Ez
E t = rEr = −
r 2
kc ∂r
and
−
→
−
→ φ = z × E t , ηT M = β η
H t = φH
ηT M k
Noting that
k ǫr k0
kc = kd , = = ωǫ
η η0
in dielectric and
k k0
kc = jka , = = ωǫ0
η η0
in air, the transverse fields are found to be
jβ
Er = {kd A1 J0′ (kd r) + kd A2 Y0′ (kd r)}
kd2
and
ǫr k0 ωǫ
Hφ = Er = Er
βη0 β
in the region a < r < b and
jβ k0 ωǫ0
Er = ka A3 K0′ (ka r) , Hφ = Er = Er
ka2 βη0 β
k0 ǫr k0
j A3 K0′ (ka b) = −j {A1 J0′ (kd b) + A2 Y0′ (kd b)}
ka η 0 kd η0
178 8. Dielectric waveguides
or
J0 (kd a)
A3 kd K0′ (ka b) ′ ′
= −ǫr ka A1 J0 (kd b) + Y0 (kd b)
Y0 (kd a)
Hence, we can write
kd K0′ (ka b) J ′ (kd b) Y0 (kd a) − Y0′ (kd b) J0 (kd a)
= −ǫr ka 0
K0 (ka b) J0 (kd b) Y0 (kd a) − Y0 (kd b) J0 (kd a)
and
β = 2, 02 rad/cm
For large ka r,
π
K0 (ka r) ∼ exp (−ka r)
2ka r
Thus at r ∼ 2λ0 , the field is down a factor 0, 052 from its value at r = b.
The effective waveguide is thus about 4λ0 . Therefore surface-wave guides
are practical only at high frequencies.
and
−
→
n
· E =0
We can insert these walls safely along planes for which the tangential mag-
netic field is zero without disturbing the field. Such planes correspond to
certain symmetry planes in a given field problem.
We assume that the electric field does not vary in the x-coordinate and is
therefore a function of y and z only. Thus, only Ey , Ez and Hx are present.
This then allows us to place a magnetic wall along any x = const. plane
without disturbing the field. In our present case, we choose the planes
x = ±W as shown in Fig. 8.7. The problem can now be treated as a
waveguide structure with
• Magnetic walls at x = ±W
intermediate value between its value in air and the dielectric. i.e.,
√
k0 = ω µ0 ǫ0 < β < k
We note that β must of necessity be the same in both air and the dielec-
tric regions at the interface since the tangential components of the fields
are continuous across the air/dielectric interface for all z. Thus, let kc = kd
in the dielectric and kc = ka in air. Thus,
or
(ǫr − 1) k02 = kd2 − ka2
In the two regions,
∂ 2 Ez
+ kd2 Ez = 0, 0 ≤ y ≤ a
∂y 2
and
∂ 2 Ez
+ ka2 Ez = 0, a ≤ y ≤ b
∂y2
with boundary conditions
Ez (y) = 0, y = 0, b
Ez = A1 sin kd y, 0 ≤ y ≤ a
= A2 sin {ka (b − y)} , a ≤ y ≤ b
Hence
A1 sin kd a = A2 sin ka d
and continuity of Hx at y = a requires
ǫr 1
A1 cos kd a = − A2 cos ka d
kd ka
Dividing the two equations we obtain
most of the modes will be evanescent since for real β, we must have ka < k0
and kd < k. Thus, for β lying between k0 and k, ka must be imaginary, say
ka = −jka,0 . We assume a corresponding value kd = kd,0 and hence
and
2 2
kd,0 + ka,0 = (ǫr − 1) k02
The problem can be attacked separately at low and high frequencies.
and
tanh ka,0 d ∼ ka,0 d
182 8. Dielectric waveguides
Thus, we have
2
kd,0 a = ǫ2r ka,0 d
or
2 (ǫr − 1) k02 a
ka,0 =
a + ǫr d
and
ǫr b
2 2 √
β = k0 + ka,0 = k0 = k0 ǫef f
a + ǫr d
where ǫef f is an effective dielectric constant.
If the current density Jz on the inner surface of the upper plate is uniform
then the magnetic field between the plates is given by Hx = Jz (electric
wall). Thus, the time-average energy stored in the magnetic field is given
by
b W
µ0 µ 1
Um = Hx2 dx dy = 0 W bJz2 = L′ Iz2
4 2 4
0 −W
where
Iz = 2W Jz
and
µ0 b
L′ =
2W
The capacitance per unit length can be regarded as two series capacitances,
′ ′
one for the air region, Cair and the other for the dielectric, Cdiel , where
′ 2ǫ0 W ′ 2ǫr ǫ0 W
Cair = , Cdiel =
d a
so that −1
a 2ǫr ǫ0 W
C ′ = 2ǫ0 W d+ =
ǫr a + ǫr d
and
ǫr ǫ0 µ0 b ǫr
L′ C ′ = = ǫ0 µ0 b
a + ǫr d a + ǫr d
At low frequencies, we can of course write
sin kd,0 a kd,0 a
A2 = A1 ∼ −jA1
j sinh ka,0 d ka,0 d
In the dielectric, 0 ≤ y ≤ a, we then have
Ez = A1 kd,0 y
jβ b
Ey = − A1 = −jA1
kd,0 (ǫr − 1) d
8.4 ♣Inhomogeneously-filled parallel-plate transmission line 183
jǫr k0 A1 ǫr (ǫr d + a)
Hx = A1 = j
kd,0 η 0 η0 (ǫr − 1) d
and in the air region, a ≤ y ≤ b,
kd,0 a
Ez = A1 (b − y)
d
jβkd,0 a jβ b
Ey = − 2 A1 = −j ǫr A1 = jA1 ǫr
ka,0 kd,0 (ǫr − 1) d
jk0 A1 ǫr (ǫr d + a)
Hx = A1 = j
kd,0 η 0 η0 (ǫr − 1) d
We note therefore Ez → 0 as k0 , kd,0 → 0 while Ey and Hx remain constant,
hence a TEM situation. Thus, letting
b
b
V = − Ey dy = jA1 (a + ǫr d)
(ǫr − 1) d
0
FIGURE 8.8. Graphical solution for the TM modes in a grounded dielectric slab
since kd = k @ β = 0,
√ √
tan (k0 a ǫr ) = −k0 ǫr tan k0 d
Since the current vanishes at the ends of the wire, charges must be deposited
there. We have
d
i (t) = ± Q (t)
dt
where
Q (t) = Re [Q (z) exp (jωt)]
and
I (z) = ±jωQ (z)
or
I (z)
Q (z) = ±
jω
The pair of positive and negative charges separated by a short distance
constitutes an electric dipole of dipole moment
−
→
p = zQ (z) l
and
−
→ 1 −→ −
→ 1 1 ∂
θ ∂
E = ∇×H = r (Hφ sin θ) − (rHφ )
jωǫ jωǫ r sin θ ∂θ r ∂r
188 9. Antenna theory
or
2
3
Il 2 1 1
Er = − ηk cos θ + exp (−jkr)
2π jkr jkr
2
3
Il 2 1 1 1
Eθ = − ηk sin θ + + exp (−jkr)
4π jkr jkr jkr
and
Eφ = 0
Near field
In the region near the Hertzian dipole (near zone),
2πr
kr = << 1
λ
so that,
1 1
exp (−jkr) ∼ 1 + (−jkr) + (−jkr)2 + (−jkr)3 + . . . ∼ 1
2! 3!
Thus, the term
Idl
Hφ ∼ sin θ
4πr2
dominates and
Il η
Er ∼ 2 cos θ
4πr3 jk
Il η
Eθ ∼ sin θ
4πr3 jk
These expressions describe quasi-static fields and we see that Hφ is in phase
with the current near the current element and can be identified with the
induced field of Ampere’s law while the near-field electric field is identified
with the field obtained from an electric dipole. The H- and E-field compo-
nents are 90◦ out of phase and represent no net energy flow according to
Poynting’s theorem. This zone is referred to as the reactive zone.
Far field
In the far field, we have
2πr
kr = >> 1
λ
Hence
Il
Hφ ∼ j sin θ exp (−jkr)
4πr
Il
Eθ ∼ j ηk sin θ exp (−jkr)
4πr
9.1 Radiation from elemental dipoles 189
in free space. The Poynting vector is radial and the time-average radiated
power is hence
1 −
→ →
− k 2 I 2 l2
Pr = Re E θ × H ∗φ = η0 sin2 θ
2 32π2 r2
The total power over a spherical surface at r from the current element is
π π
2
−
→ − →
k2 I 2 l2 3 l
PT = P ·d S = Pr 2πr2 sin θ dθ = η 0 sin θ dθ = 40π I2 2
W
16π2 λ
S 0 0
since
π π π
3 2
sin θ dθ = sin θ sin θ dθ = 1 − cos2 θ sin θ dθ
0 0 0
1 ' 3 (π
= [− cos θ]π0 + cos θ 0
3
4
=
3
The far field zone is also called the Fraunhofer zone.
In the intermediate zone called the Fresnel zone, both radiative and
stored energy exists. At a distance rf z from the current element such that
2
1 1 1 λ λ
− = 0 ⇒ rf z = = ∼
krf z krf z k 2π 6
the two components are equal. rf z therefore marks the boundary between
the near and far zones.
i (t) = I cos ωt
−
→
d l ′ = −
x sin φ + y cos φbdφ
z P(r,θ,π/2)
r
r’
b
ϕ
I O y
φ’ P’ P’’
b
dl’
x
−
→ µm 2
A =φ β (1 + jkr) exp (−jkr) sin θ, m = Iπb2
4πr2
and
−
→ − → −
→
∇ × H = jωǫ E
−
→ −
→ − → −
→
Thus, exchanging H for E , − E for H , µ for ǫ and ǫ for µ in the two
equations leaves the equations unchanged so that the case of an electric
source can be adapted to one with a magnetic source. We can therefore
represent a current loop of radius b and current I by a magnetic dipole
m = πb2 I which has a magnetic field of the same form as the electric field
of an electric dipole p = Ql. From continuity, the electric charge on one
end of the dipole is
I I
Q= ⇒p= l
jω jω
Replacing Il/jω by µm, the magnetic fields for a magnetic dipole can be
found from those of the Hertzian dipole.
−
→ −
→
If E e and H e are the electric and magnetic field intensities of the
−
→ −
→
Hertzian dipole and E m and H m those of the magnetic dipole, then
−
→
−
→ −
→ −
→ Em
E e = η H m, H e = −
η
as long as the electric and magnetic moments are related as
√
Il = jkm, k = ω µǫ
or
l
Im
Eθ = j ηk sin θ exp (−jkr) sin k (l − z) cos (kz cos θ) dz
2πr
0
60πIm
= j exp (−jkr) F (θ)
r
where
l
cos (kl cos θ) − cos kl
F (θ) = k sin θ sin k (l − z) cos (kz cos θ) dz =
sin θ
0
is the E-plane pattern function of the linear dipole antenna and is sym-
metrical about the plane θ = π/2.
Note that we can determine the fields alternatively by recognising that,
since ϕ = θ anywhere along the line, then
l
−
→
N = z I exp (jkr′ cos θ) dz ′
−l
0
= z Im sin k (l + z) exp (jkr′ cos θ) dz ′
−l
l
z Im sin k (l − z) exp (jkr′ cos θ) dz ′
+
0
and with
exp (αx)
exp (αx) sin (bx + c) dx = 2 {a sin (bx + c) − b cos (bx + c)}
a + b2
194 9. Antenna theory
z r’
θ
r
dz
l ϕ
z
Im
l
−
→ 2Im
N = z {cos (kl cos θ) − cos (kl)}
k sin2 θ
and
Nθ = −Nz sin θ
We therefore have
jηIm cos (kl cos θ) − cos (kl)
Eθ = exp (−jkr)
2πr sin θ
and
Eθ
Hφ =
η
as before
and
2Pr
Rr = 2
∼ 73, 1 Ω
Im
The input reactance of the dipole is small and inductive and can be elimi-
nated by making 2l slightly less than λ/2. The half-dipole is a near-resonant
structure for thin wires for which the open-circuit uniform transmission line
approximation is valid. The resonant length and resonant terminal resis-
tance decrease with increasing wire radius to wavelength ratio.
z y
90o x x
z
FIGURE 9.4. The E-and H-plane radiation pattern of an isolated Hertzian dipole
|Eθ |n = |sin θ|
where dΩ = sin θdφ. The directive gain of the antenna is hence defined as
K (θ, φ) K (θ, φ)
GD (θ, φ) = = 4π
Pr /4π
KdΩ
Ω
The maximum directive gain is defined as the directivity of the antenna and
is the ratio of the maximum radiation intensity to the average radiation
intensity and usually denoted by
Kmax Kmax
D = GD,max = = 4π
Kave Pr
This can be expressed as
4π |Emax |2
D=
2ππ
|E (θ, φ)|2 sin θdθdφ
0 0
FIGURE 9.5. The radiation pattern of a half-wave and the Hertzian dipoles
the radiation efficiency is the ratio of the power gain to the directivity,
Gp Pr
ηr = =
D Pi
The amount of power radiated can be given in terms of an equivalent
radiation resistance Rr which is a fictitious resistance that would dissipate
the same amount of power as Pr when the current through this resistance
is the maximum current along the antenna. The radiation resistance can
be obtained from the input impedance of the antenna regarded as a series
frequency-dependent impedance of the form
Zi = R + jX
2Pr
Rr = 2
= 36, 54 Ω
Im
and
Kmax Kmax
D= = = 1, 64
Kave Pr /4π
the same as the λ/2 antenna.
λ/4 λ/4
I
λ/4
I
l
exp (−jkr)
Eθ = η0 Hφ = j30k sin θ I (z) exp (jkz cos θ) dz
r
−l
Let I (0) be the input current at the feed point of the antenna. Then,
exp (−jkr)
Eθ = j30I (0) k le (θ)
r
where
l
sin θ
le (θ) = I (z) exp (jkz cos θ) dz
I (0)
−l
l
1
le = I (z) dz
I (0)
−l
λ/4
λ
le (θ) = sin θ sin k − |z| exp (jkz cos θ) dz
4
−λ/4
9.4 Antenna array theory 203
or π
2 cos 2 cos θ λ λ
le (θ) = ⇒ le = le,max (θ) = <
k sin θ π 2
Note that le (θ) is only meaningful for relatively short antennas with
maximum current at the feed point. For a receiving antenna, the effective
length is defined as the ratio of the open-circuit voltage Voc induced at the
antenna terminals to the effective field intensity |Ei | at the antenna that
induces it,
Voc
le (θ) = −
Ei
This is the same as when the antenna is used for transmitting. If Ei is not
parallel to the dipole, there is a polarisation mismatch and the magnitude
of the open-circuit voltage will be
−
→
|Voc | = le E i
exp (−jkr0 )
E0 = Em F (θ, φ)
r0
exp (−jkr1 )
E1 = Em F (θ, φ)
r1
and the electric field intensity at P is hence
exp (−jkr0 ) exp (−jkr1 )
E = E0 + E1 = Em F (θ, φ) +
r0 r1
204 9. Antenna theory
r0
r1
φ
0 1
x
d
In the far-field, r0 >> d/2 and the difference between r0 and r1 is suffi-
ciently small such that r1 ∼ r0 in the magnitude expression. In the phase
factor, we retain this difference however since k may be large. Hence,
r1 ∼ r0 − d sin θ cos φ
and
exp (−jkr0 )
E = Em F (θ, φ) [1 + exp {j (kd sin θ cos φ + ϕ)}]
r0
exp (−jkr0 ) ϕ
ψ
= Em F (θ, φ) exp j 2 cos
r0 2 2
where
ψ = kd sin θ cos φ + ϕ
and
2Em ψ
|E| = |F (θ, φ)| cos
r0 2
|F (θ, φ)| may be called the element factor and cos ψ2 the space factor of
the array or array factor.
The pattern function of an array of identical elements is described by
the product of the element and array factors (principle of pattern mul-
tiplication). For two parallel z-directed half-wave dipoles, the total field
magnitude is hence,
2Em cos π2 cos θ ψ
|E| = cos
r0 sin θ 2
and since ψ is also a function of θ, the E-plane pattern is not the same
as that of a single half-wave dipole except for ψ = ±π/2. In the H-plane,
θ = π/2 and the pattern is determined by cos ψ2 .
9.4 Antenna array theory 205
Example: Plot the H-plane radiation patterns of two parallel dipoles for
the following cases
(a) d = λ/2, ϕ = 0
Let the dipoles be z-directed and placed along the x-axis. In the H-plane
where θ = π/2, each dipole is omnidirectional and the normalised pattern
function is equal to the normalised array factor |A (φ)|. Thus,
ψ kd
|A (φ)| = cos = cos cos φ + ϕ
2 2
(a) d = λ/2, kd = π, ϕ = 0
π
|A (φ)| = cos cos φ
2
The pattern has a maximum at φ0 = ±π/2, i.e. broadside direction.
Fig. 9.8.
1 2 1
x
λ/2 λ/2
φ
(a) 1 2
1
0o
λ/2
1 1
1
x
1
(b)
(c)
0 φ 1 2 N-1
x
d
where
ψ = kd cos φ + ϕ
Thus
1 1 − exp (jN ψ) 1 sin Nψ/2
|A (φ)| = =
N 1 − exp (jψ) N sin ψ/2
As φ varies from 0 to 2π, ψ varies from kd + ϕ to −kd + ϕ covering a
range of 2kd or 4π/λ and defines the visible range of the radiation pattern.
Thus,
kd cos φ0 + ϕ = 0
or
ϕ
cos φ0 = −
kd
The values of ϕ lead to
ψ
N = ±nπ, n = 1, 2, . . .
2
Note that the locations of the nulls in φ are different for broadside and
endfire arrays due to the different values of ϕ for the null condition.
(3) An angular width of the main beam between the first nulls for large
N given by
ψ 2π
N 01 = ±π ⇒ ψ01 = ±
2 N
where ψ01 denotes the value of ψ at the first nulls.
2∆φ = 2 (φ01 − φ0 )
9.4 Antenna array theory 209
At φ01 , we have
ψ01
cos φ01 = cos (φ0 + ∆φ) =
kd
which for φ0 = π/2, gives
π 2π
cos + ∆φ = − sin ∆φ = −
2 N kd
or
λ λ λ
∆φ = sin−1 ∼ =
Nd Nd l
where l = N d is the length of the array. Thus, for a long uniform broadside
array, the beam width of the main lobe is approximately twice the reciprocal
of the array length in wavelengths.
∆φ2
cos φ01 = cos ∆φ ∼ 1 −
2
Hence
∆φ2 λ 2λ λ
∼ ⇒ ∆φ ∼ = 2
2 Nd Nd l
which is larger than that of the broadside array.
The side lobes of the radiation occur where the numerator |sin Nψ/2| = 1
or
ψ π
N = ± (2m + 1) , m = 1, 2, . . .
2 2
The first side lobes occur when
ψ π
N = ±3 , m = 1
2 2
where we note that
ψ π
N =± , m=0
2 2
does not represent a location of side lobes because these points are still
within the main lobe region.
210 9. Antenna theory
for large N . Thus, for a long uniform array, the first side lobes are 13, 5 dB
below the main lobe. The level of the side lobe can be lowered by tapering
the current distribution in the array elements.
Example: Determine the array factor and plot the normalised radiation
pattern of a linear array of five isotropic elements spaced λ/2 apart and
having excitation amplitude ratio 1 : 2 : 3 : 2 : 1. Compare the level of the
first side-lobe with that of a five-element uniform array.
The normalised array factor is given by
1
|A (φ)| = |1 + 2 exp (jψ) + 3 exp (j2ψ) + 2 exp (j3ψ) + exp (j4ψ)|
9
1
= |exp (j2ψ) {3 + 2 exp (jψ) + exp (j2ψ)}|
9
1
= |{3 + 4 cos ψ + 2 cos 2ψ}|
9
Fig. 9.12 illustrates the radiation pattern of a non-uniform array. This
expression holds for a general ψ = kd cos φ+ϕ. Broadside radiation requires
ϕ = 0, ψ = kd cos φ and for an element spacing d = λ/2, ψ = π cos φ. The
normalised array factor is hence given by
1
|A (φ)| = |{3 + 4 cos (π cos φ) + 2 cos (2π cos φ)}|
9
an = a0 exp (−jn∆ϕ) , n = 0, 1, 2, . . . , N − 1
9.4 Antenna array theory 211
π 0o
Thus
N−1
/
|A (φ)| = an z −n
n=0
212 9. Antenna theory
FIGURE 9.13. Super gain array showing the zero positioning array synthesis
technique
Real φ then corresponds to points on the unit circle in the complex plane
with phase angles between −kd and kd. When the zeros of |A (φ)| lie on
this part of the unit circle, they correspond to the zeros of the pattern or
the cones of silence.
The zeros are uniformly distributed over the entire unit circle for the
broadside array except for the missing one at φ = 0, where the main lobe
builds up. This suggests a technique for array synthesis by positioning zeros
on the unit circle so that they lie close together in the region where the
pattern is to have minimum magnitude and wider apart where it is to
have maximum magnitude. For an element spacing of less than λ/2, real φ
covers only part of the unit circle and since the zeros of the uniform array
are spread over the entire unit circle, some of these are invisible because
they fall on the part of the circle where φ is imaginary.
An arbitrarily high gain antenna can be synthesised by moving the in-
visible zeros to the visible part of the circle but this then creates a very
large lobe in the position vacated by the zeros. While this may seem incon-
sequential, it represents reactive energy which can be a problem especially
with broadside arrays. Uniform arrays with in-phase currents have very
high gain and synthesised arrays with gains higher than these are known
as supergain antennas. Fig. 9.13 shows the zero pattern for an example of
a supergain antenna.
We note that the half-wave dipole can be analysed as an array of Hertzian
dipoles. A common array antenna is the Yagi-Uda array used for TV re-
ceptions in many homes.
9.5 Antennas for reception 213
A r B
1
Z11-Z12 Z22-Z12
I1 I2
V1 Z12 V2 ZL
V2 = Z21 I1 + Z22 I2
A bilateral space is such that reciprocity holds i.e., Z21 = Z12 . In practice,
A and B will be separated by large distances such that Z12 = Z21 ∼ 0 in
the limit as r → ∞. Hence
V1 ∼ Z11 I1 = ZA IA
V2 ∼ Z22 I2 = ZB IB
However, the coupling from the transmitting to the receiving antenna can-
not be neglected as it is by this means that energy arrives at B. We can
modify the equivalent two-port as shown in Fig. 9.16 and apply Thevenin’s
theorem at the load ZL looking toward antenna A, whereby the open-circuit
voltage is then seen to be
Z12
Voc = V1
Z11
and the source impedance as
Pt
Pave = GD,A
4πr2
Thus
Pr Ae,B GD,A
=
Pt 4πr2
so that
RA RB Ae,B GD,A
|Z21 |2 =
πr2
If B transmits instead while A receives, then
RA RB Ae,A GD,B
|Z12 |2 =
πr2
Hence,
GD,A GD,B
=
Ae,A Ae,B
The ratio of the directive gain to the effective aperture of an antenna is
the universal constant determined by the directive gain and effective area
of any antenna.
Example: Determine the effective area, Ae (θ) of a Hertzian dipole of
length 2l << λ used to receive an incident plane electromagnetic wave of
wavelength λ polarised as shown in Fig. 9.17.
Let the amplitude of the electric field intensity at the dipole be Ei . Thus,
Ei2
Pave =
2η0
9.5 Antennas for reception 217
where
2
πl
Rr = 80 , (Ω)
λ
The effective aperture is hence
Pr η 3 λ2
Ae (θ) = = (l)2 sin2 θ = sin2 θ
Pave 4Rr 2 4π
which is independent of the dipole length.
For an elemental electric dipole,
3 2
GD (θ) = sin θ
2
Thus
4π
GD (θ) = Ae (θ)
λ2
and
4π
GD (θ, φ) = Ae (θ, φ)
λ2
The universal constant is hence given by
GD (θ, φ) 4π
= 2
Ae (θ, φ) λ
In a thin linear antenna of effective length le (θ), power available to the
antenna under matched conditions is
2
Voc (−le Ei )2
Pr = =
8Rr 8Rr
Hence
30π 2
Ae (θ) = l (θ)
Rr e
−
→ ωp exp (−jkr)
Es =
θEs = − θ η0 β sin θ
4π r
2 3 exp (−jkr) ǫr − 1
= −θk b Ei sin θ
r ǫr + 2
Ei2
Pi =
2η0
Hence
4
2 ǫr − 1
σbs = 4πb kb , m2
ǫr + 2
At f = 15 GHz, λ = 20 mm, the radius of the rain drop is given by
3/2 mm = 1, 5 mm << λ. Thus
We see that
σ bs
= 0, 177
πb2
where πb2 is the geometrical cross-section area of the drop.
Pt = 0, 90 × 5 × 104 W = 4, 5 × 104 W
−12 1
Pr = 1, 5 × 10 = 1, 58 × 10−12 W
1 − 0, 05
σ bs = 1 m2
σbs A2e Pt
r4 = ⇒ r = 4, 36 × 104 m = 43, 6 km
4πλ2 PL
In communication systems using geostationary satellites, the geosynchro-
nous orbit has a radius of 42300 km. The Earth’s radius is ∼ 6380 km
and satellites are therefore ∼ 36000 km from the Earth’s surface.
Now consider the situation shown in Fig. 9.18, where Ae is the effective
scattering aperture of the object. For a large and flat conducting sheet
oriented for specular reflection, Ae = Ap where Ap is the physical aperture.
For most objects, Ae < Ap since scattering is isotropic. The radar equation
can now be modified as
Pr A2 σ bs
= e2 4
Pt 4πλ r
where
scattered power (assumed isotropic) 4πr2 Psc
σ bs = =
incident power density Pi
is the radar cross-section, Psc the scattered power density and Pi the inci-
dent power.
9.5 Antennas for reception 221
We note that
λ2
A′e =
4π
then
σ bs λ2
A2e = σbs A′e =
4π
For a large perfectly conducting sphere of radius a, σbs = πa2 , the physical
cross-section area. It will be smaller for spheres of a finite conductivity
however.
The noise power per unit bandwidth available at the terminals of a resis-
tance R at an absolute temperature T is given by the Johnson noise power
formula
Pn = kT
T is then referred to as the noise temperature of the resistor. Replacing the
resistor by a lossless antenna, with a radiation resistance Rr = R, the input
222 9. Antenna theory
impedance looking into the terminals remains unchanged but the noise will
be different unless the region from which the antenna receives power is at
the same temperature T . In radio telescopes, the radiated beam is pointed
to a region of the sky whose temperature is close to 0 K and the antenna
temperature is then this value and not the physical antenna temperature.
Let the antenna temperature be designated as TA , then the noise power is
given by
Pn = kTA W/Hz
The radio telescope can thus be regarded as a temperature measuring
device (radiometer). This assumes that the antenna has negligible thermal
losses and little side and back lobes. Reception with the telescope is then
conveniently described in terms of received power flux and the flux density
found as
Pn kTA
S= = W m−2 Hz −1
Ae Ae
In astronomy, the flux unit (f u) = 10−26 W m−2 Hz −1 is used. This ex-
pression gives the flux density for a celestial source of extent Ω that is
small compared with the antenna beam solid angle ΩA and the antenna
aligned with the source. Hence the expression applies to point sources. The
actual flux density for other objects is twice this value since the antenna
will receive only half of the radiated unpolarised power.
For small angular size sources of known magnitude, the source tempera-
ture can be determined from
ΩS
TA = TS
ΩA
where ΩS is the source solid angle, ΩA the antenna beam solid angle and
TS the source temperature.
B
z r Eθ2
Eθ1
θ1
A
r’ h2
h1
C
d
h1
θ2
A’
and
K
Eθ2 = exp (−jkr′ ) sin θ 2
r′
where
Il
K=j ηk
4π
and assuming a perfectly conducting Earth.
−
→ −
→
In general, E θ1 and E θ2 are not parallel but for d >> h1 and h2 ,
θ1 ∼ θ2 = θ so that
K
Eθ ∼
θ F sin θ exp (−jkr)
r
where
F = 1 + exp {−jk (r′ − r)}
We then have
2
2 2 h2 − h1 (h2 − h1 )2
r = d + (h1 − h2 ) = d 1 + ∼d+
d 2d
and
(h2 + h1 )2
r′ = d2 + (h2 + h1 )2 ∼ d +
2d
so that
(h2 + h1 )2 − (h2 − h1 )2 h1 h2
r′ − r ∼ =2
2d d
and
h1 h2
|F | = 1 + exp −j2k
d
224 9. Antenna theory
FIGURE 9.20. Antennas above ground showing reflection into the plane
FIGURE 9.21. Separation of radiation field into a skywave and ground wave for
a finite σ ground plane
π/2
ηI 2 [cos (kl cos θ) − cos (kl)]2
PT = m dθ
4π sin θ
0
When the conductor is not perfect as will be the case when the ground
plane is the Earth’s surface, for which the displacement and conduction
currents are approximately equal at f ∼ 5 MHz, the radiated field splits
into a space or sky wave and a ground or surface wave as shown in Fig.
9.21 for the Hertzian dipole. The space wave varies as 1/r with distance
from the antenna while the surface wave varies as 1/r2 .
θ
I
excitation terminating
source h impedance
z
ground plane
If we ignore the effect of the nearby ground, the far-field will have an electric
field intensity given by
L
jηk exp (−jkr)
Eθ = sin θ I (z) exp (jkz cos θ) dz
4π r
0
L
jηIm exp (−jkr)
= sin θ I (z) exp {−jkz (1 − cos θ)} dz
4π r
0
or
j60Im L
Eθ = exp −jk r + (1 − cos θ) F (θ)
r 2
where
sin θ sin (kL (1 − cos θ))
F (θ) =
1 − cos θ
is the pattern factor of an isolated travelling-wave antenna of length L.
The pattern is not symmetrical about the θ = π/2 plane unlike in resonant
antennas such as the λ/2 dipole. The main lobes tilt in the direction of the
current wave with the tilt angle becoming greater with increasing antenna
length. The side lobes of this antenna are also large.
The effect of ground can be accounted for using the method of images as
shown in Fig. 9.23. This is then an array of two travelling-wave antennas
carrying equal but opposite currents. Using the principle of pattern multi-
plication, the resulting pattern factor is F (θ) |cos ψ/2| of the two-element
array with ϕ = π. Thus,
cos ψ = cos kh sin θ cos φ + π = |sin (kh sin θ cos φ)|
2 2
9.7 Other types of practical antenna systems 227
2h
h ground plane
I end-view
(two-element array)
−
→
φ = N kI η exp (−jkr) j
E = θEθ + φE θs +
φkπb 2
sin θ
4π 0 r
228 9. Antenna theory
z
FIGURE 9.24. Helical antenna, its equivalent visualisation and E-plane radiation
pattern
V1 = I1 Z11 + I2 Z12
V2 = I1 Z21 + I2 Z22
where the Zij are constants dependent on the element length l1 and l2 as
well as the separation d. Since the voltage at the drive point of the parasitic
element is zero, I2 is determined by I1 through induction. We can therefore
write
Z21
I2 = − I1
Z22
and the array factor then given by
d ∼ 0, 16λ
and
2l2
0, 51 < < 0, 52
λ
with 2l1 ∼ λ/2 in this case.
When the intention is to have maximum radiation in the +z-direction,
element 2 is known as a director and has a slightly smaller length than
the active element. The optimum length and separation from the active
element with 2l1 ∼ λ/2 are in this case
d ∼ 0, 11λ
and
2l2
0, 38 < < 0, 48
λ
The Yagi-Uda array in typical use has both a director and reflector and
the three elements complicate the analysis somewhat. However, conclusions
drawn from the two element array remain fairly valid. Practical Yagis use
a longer reflector and shorter lengths resp. than those given for maximum
radiation. The directivity of the array then turns out to be insensitive
to the director and reflector spacing from the active element. The input
impedance of the array is then given by
I2 I3
Zin = Z22 + Z21 + Z32
I1 I2
9.7 Other types of practical antenna systems 231
z
2a
b34
h3
2 3 4 5
1
and due to the current phasing, Zin tends to be rather low. If the driven
element is a folded dipole, Z22 is four times that of a single dipole and
the overall input impedance can be raised. Additional directors improve
antenna characteristics and directivities in the range 10 − 100 are common
depending on the number of directors in the array. No significant improve-
ment is achievable from providing a second reflector, though. Fig. 9.27
illustrates a five-element Yagi-Uda antenna.
Since all the elements in the Yagi-Uda antenna are coupled, the current
distribution on each element depends on the length and spacing of all the
other elements. From experience, little advantage accrues from using more
than one reflector but directivity improves with the number of director
elements. The Yagi-Uda is an endfire array with the main lobe pointing
away from the reflector as shown in Fig. 9.28 for a three element array. A
good Yagi-Uda antenna should have
• High directivity
• Narrow beam width
232 9. Antenna theory
• Low side-lobes
For a Yagi with a = 0, 003369λ, (ln λ/2a = 5), a typical well dimensioned
six-element Yagi with four directors of equal length and equal spacing may
have the parameters shown in Table 9.1. The directivity of the λ/2 dipole
used as reference is of course 1, 64 or 2, 15 dBi referenced to an omnidirec-
tional antenna.
Uniformly spaced directors of equal length do not translate into an opti-
mal design. An optimally designed Yagi with dipole radius a = 0, 003369λ
would have the parameters shown in Table ??.
rn+1
rn
may be a disadvantage. In the strict sense, the spiral must be infinite for
the antenna to be truly frequency-independent.
The log-periodic antenna of Fig. 9.30 is another example of a frequency-
independent radiative structure. The teeth present discontinuities that tend
to localise the region of maximum radiation and cause the current to dimin-
ish rapidly beyond the region. Tooth lengths are determined by the angles
between the lines from the origin and the region of strongest radiation is
where the tooth length is ∼ λ/4.
The spacing between successive edges of the teeth follow the same rule
governing the distance between neighbouring conductors in the equiangular
spiral, i.e.
fn = τ fn+1
or
1
ln fn+1 = ln fn + ln
τ
Wire or tubular log-periodic antennas may be constructed when the
tubes outline the cut-out design of the planar log-periodic antenna. The-
oretically, sheet thickness and wire diameter should increase linearly with
distance from the feed point in accordance to τ and this becomes important
when the bandwidth demand is severe. The planar log-periodic antenna is
also bidirectional i.e. the main lobe is broadside on either side. However,
9.7 Other types of practical antenna systems 235
rn
rn+1 dn
ln
ln+1
folding the two arms into a wedge-like structure will make the antenna uni-
directional. Without the teeth, we have the bowtie antenna with limited
broadband characteristics and is useful in UHF television reception.
The log-periodic dipole array of Fig. 9.31 is frequency-independent as
well. The dipoles are unequally spaced and have lengths given by
ln+1 rn+1
= =τ
ln rn
where τ is a design parameter. Since the element spacings are related to
the distances to an imaginary apex point and apex angle α, i.e.
dn = rn − rn+1 = rn (1 − τ )
then
dn+1
=τ
dn
as well.
The angle α or spacing factor κ given by
dn
κ=
2ln
is another necessary design parameter. The relationship between τ , α and
κ is
α 1 ln ln (1 − τ ) 1−τ
tan = = =
2 2 rn dn 4κ
236 9. Antenna theory
Plane of aperture
Actual
source Aperture
(S)
−
→
Similar to the vector potential A associated with electric sources, we can
−
→
define a vector potential F associated with the magnetic sources and write,
−
→
−
→ µ J s exp (−jkr)
A= dS ′
4π r
S′
and −
→
−
→ ǫ M s exp (−jkr)
F = dS ′
4π r
S′
and
−
→ −
→
L = M s exp (jkr′ cos ϕ) dS ′
S′
′
where r and ϕ have the previously defined meaning. The only components
of the field decreasing no faster than 1/r are
exp (−jkr)
Eθ = ηHφ = −j {ηNθ + Lφ }
2λr
exp (−jkr)
Eφ = −ηHθ = j {−ηNφ + Lθ }
2λr
9.8 Radiating apertures 239
and
Lθ = −E0 dS sin φ cos θ, Lφ = −E0 dS cos φ
and the far-field radiation is then
jE0 dS
Eθ = exp (−jkr) (1 + cos θ) cos φ
2λr
jE0 dS
Eφ = − exp (−jkr) (1 + cos θ) sin φ
2λr
and
E02 (dS)2 2 2
K = (− cos φ cos θ − cos φ) + (sin φ + sin φ cos θ)
8ηλ2
2
E02 (dS) θ
= cos4
2ηλ2 2
where
(xx′ − yy ′ )
R= (x − x′ )2 + (y − y ′ )2 + z 2 ∼ r −
r
and since the difference between r and R is important only in the phase
consideration, then
j exp (−jkr) ′ ′ (xx′ + yy ′ )
E (x, y, z) = E (x , y ) exp jk dx′ dy′
λr r
S′
• Only small angles about the polar axis are considered so that the
paraxial approximation is valid
a/2 b/2
exp (−jkr) xx′ yy ′
E (x, y, z) = j E0 exp jk exp jk dx′ dy ′
λr r r
−a/2−b/2
exp (−jkr) sin (kax/2r) sin (kby/2r)
= j E0 ab
λr (kax/2r) (kby/2r)
exp (−jkr) kax kby
= j E0 ab sinc sinc
λr 2r 2r
In the plane x = 0 the null angles have a similar form with b and y replacing
a and x respectively. The main lobe hence becomes narrower as the aperture
dimensions increase.
The directivity of the rectangular aperture is then
Piso
D=
P
where
|Emax |2 E 2 (ab)2
Piso = 4πr2 = 2π 0 2
2Z0 λ0 η
is the isotropically radiated power. The actual power radiated, neglecting
currents in the surrounding aperture plane is
E02
P = ab
2η
and
4π 4π
D= ab = 2 A
λ2 λ
where A = ab is the aperture area.
If the aperture dimensions are much larger than the wavelength, almost all
the energy of the radiated field will be contained in a small region about
the z-axis (paraxial field). Thus
−
→
Ep = x
Ep
where
j
Ep = Ea (x′ , y′ ) exp (−jkr) dx′ dy ′
λr0
9.8 Radiating apertures 243
a
b
r xx′ + yy ′ ) · (
= r0 − ( x sin θ cos φ + y sin θ sin φ)
= r0 − (x sin θ cos φ + y ′ sin θ sin φ)
′
Thus,
j
Ep = exp (−jkr0 ) F (θ, φ)
λr0
where
F (θ, φ) = Ea (x′ , y′ ) exp {jk sin θ (x′ cos φ + y ′ sin φ)} dx′ dy ′
apert.
Ea (x′ , y ′ ) = f1 (x′ ) f2 (y ′ )
and
a/2
= C1 f1 (x′ ) exp (jkx′ sin θ) dx′
−a/2
where C1 is a constant.
(b) In the yz-plane where φ = π/2,
a/2 b/2
′ ′
Fyz (θ) = f1 (x ) dx f2 (y ′ ) exp (jky ′ sin θ) dy ′
−a/2 −b/2
a/2
= C2 f2 (y ′ ) exp (jky ′ sin θ) dy ′
−a/2
and
Pr = total power radiated
1 2
= |Ea (x′ , y ′ )| dx′ dy ′
2η
apert.
Thus 2
E (x ′ ′
, y ) dx ′
dy ′
a
4π apert.
D= 2
λ
|Ea (x′ , y ′ )|2 dx′ dy′
apert.
9.8 Radiating apertures 245
(a) The pattern function in a principal plane and the half-power beam
width
(b) The location of the first nulls and the level of the first side lobes
a/2
′ ′ sin ψ
Fxz (θ) = b exp (jkx sin θ) dx = ab
ψ
−a/2
where
πa
ψ=sin θ
λ
Fyz (θ) has the same form, where b replaces a. Note that the pattern
function is similar to the array factor of a uniform array for small
ψ.
The half-power points are obtained by setting
sin ψ1/2 1
=√
ψ1/2 2
Hence,
πa
ψ1/2 = sin θ1/2 = 1, 39
λ
or
λ
sin θ1/2 = 0, 0442
a
For a sufficiently large aperture, sin θ 1/2 ∼ θ1/2 and
λ λ
2θ1/2 ∼ 0, 88 rad ∼ 50◦
a a
Hence
sin ψ1 sin 1, 43π
=
ψ 1, 43π = 0, 217
1
When referred to unity at ψ = 0, the first side lobes are 20 log10 (1/0, 217) =
13, 3 dB below the level of maximum radiation.
Mz = Ex = Em cos kz
λ/4
Lz = 2gEm cos kz ′ exp (jkz ′ cos θ) dz ′
−λ/4
cos π2 cos θ
= 4gEm
k sin2 θ
and
j exp (−jkr) cos π2 cos θ
Eφ = −ηHθ = gEm
πr sin θ
−
→ −
→
which is similar to the case of the resonant dipole except that H and E
are interchanged (i.e. polarisations are orthogonal).
The radiated power is given by
π π
πgEm cos2 2 cos θ
Pr = dθ
2πη0 sin θ
0
or
2Rr,dipole
Gr,slot = ∼ 0, 00103
η2
9.8 Radiating apertures 249
η2
Zstrip Zslot =
4
Thus
2πa
exp (−jkr)
E (r, θ, φ) = j E r′ , φ′ exp jkr′ sin θ cos φ − φ′ dr′ dφ′
λr
0 0
where a is the aperture radius. If E r′ , φ′ is independent of φ′ , E (r, θ, φ)
will be also and φ = 0 in the integral. We use the integration result
2π
exp (jx cos ϕ) dϕ = 2πJ0 (x)
0
with J0 the Bessel function of the first kind and order 0. Thus
a
2π exp (−jkr)
E (r, θ) = j E (r′ ) J0 (kr sin θ) dr′
λr
0
since the imaginary parts are odd functions and integrate to zero. We can
thus write
exp (−jkr) 2 k 2 x2 + y 2 w 2
E (x, y, z) = j E0 πw exp −
λr r2
2 2 2
exp (−jkr) k w sin θ
= j E0 πw2 exp −
λr 4
which is tabulated. The angle, θ0 where the radiation is 1/e of the value at
θ = 0, is given by
−1 2 2 λ
θ 0 = sin ∼ =
kw kw πw
x′ = r′ cos φ′ , y ′ = r′ sin φ′
9.8 Radiating apertures 251
and
x′ cos φ + y ′ sin φ = r′ cos φ cos φ′ + sin φ sin φ′ = r′ cos φ − φ′
Thus
b 2π
F (θ, φ) = E0 exp jkr′ sin θ cos φ − φ′ r′ dφ′ dr′
0 0
b
J1 (kb sin θ)
= E0 2πJ0 (kr′ sin θ) r′ dr′ = E0 2πb2
kb sin θ
0
Thus
−
→ 2πb2 J1 (u)
Ep = x
jE0 exp (−jkr0 )
λr0 u
where
2πb
u = kb sin θ = sin θ
λ
The first null occurs at the first zero of J1 (u), u11 say, where
u11 = 3, 832
corresponding to an angle
3, 832λ 3, 832λ λ
θ 1 = sin−1 ∼ = 1, 22
2πb 2πb D
where D = 2b diameter of the aperture. Thus, the width between the first
nulls is
2, 44λ
2θ1 = rad
D
This is wider than that of the rectangular aperture having a = D.
The level of the first side lobes for the circular aperture is 0, 13 ( 17, 7 dB)
down from the maximum which is better than that of the rectangular aper-
ture with a = D. Note that
2π
exp jw cos φ′ dφ′ = 2πJ0 (w)
0
and
wJ0 (w) dw = wJ1 (w)
252 9. Antenna theory
FIGURE 9.39. Parabolic dish for the case with focal length equal to dish depth
x2
y=
4p
where x is the dish radius or half the dish aperture (diameter). The focus
of the dish shown in Fig. 9.39 is then at
D2 f D
f= ⇒ =
16d D 16d
where D is the diameter of the aperture and d the depth of the dish. If
f = d, then f /D = 0, 25.
A dish may be deep or shallow depending on the paraboloid slice envis-
aged during manufacture. Depending on the location of the focal point, we
can have the following situations
f 1
<
D 4
Part of the dish is not illuminated i.e. we have under illumination. This
may shield the feed system from spurious radiation.
f 1
=
D 4
Illumination of the reflector is more uniform.
f
=1
D
We have over illumination and some signal may miss the dish altogether.
The gain of the dish antenna depends on the size of the dish and the
wavelength and can be written as
(πd)2 η
Gd = 10 ln
λ2
where d is the reflector aperture, λ the wavelength and η the efficiency,
usually 0, 6-0, 8. Deep dishes with low f/D ratios have higher efficiencies
and are more shielded from noise. From manufacturability considerations,
0, 25 ≤ f/D < 1. Shallow dishes are easier to manufacture and transport
while it is also difficult to illuminate over 180◦ . Beam width design may be
used to counter the effects of over illumination.
According to the feed system location, we can have prime focus with
the focus at the centre of the dish or offset focus antennas with the focus
off centre. Large dishes use prime focus because of the greater mechanical
stability offered by the structure. Medium to high power antennas operating
above 10 GHz employ the offset focus configuration. This is because, for
small dishes, the shadow cast by the feed element is a problem. The offset
focus antenna is installed with the dish inclined towards the cooler sky
hence reducing antenna temperature and allowing for debris to slide off the
dish easily.
According to the Rayleigh limit encountered in optical mirrors, there
is little gain in making the reflector more accurate than λ/8. Hence we
can get away with quite some hefty deviations in the focal length without
compromising performance. Furthermore, the dish need not be made from
a continuous metal sheet but may be out of a mesh as long as the mesh is
small compared to the wavelength at the operating frequency. This allows
for lighter dish construction.
254 9. Antenna theory
T C
A
In the Cassegrain antenna, the rays reflecting off the parabolic dish to-
ward the dish focus are trapped before they reach the focus by a hyperbolic
sub-reflector. The feed element is now at the back end of the dish (back-fire
configuration). The antenna then has a slimmer profile than the traditional
front end feed system allowing for better shielding of the element from noise.
Due to the extra hyperbolic sub-reflector, the antenna system will be more
expensive however. Since the sub-reflector also blocks signal rays, the main
dish must have an aperture larger than 100λ. This explains why parabolic
antennas are to be found only in the microwave range of operation.
The sub-reflector is installed with one of its focus coinciding with the
focus of the parabola. Note that a hyperbola is described by
y 2 x 2
+ =4
a b
where
a2 + b2 = c2
The hyperbolic sub-reflector is replaced by an elliptical one which then
allows small reflectors, hence higher efficiency at longer wavelengths. The
elliptical reflector is governed by the same equation as the hyperbolic one
with the foci at
y = ± a2 + b2
At resonance,
1
P = I 2R
2
where R is the input resistance at resonance. This implies that for the
folded dipole, R is four times that of the single dipole so that for the folded
half-wave dipole of Fig. 9.42, R ∼ 292 Ω, constant over a wider range of
frequencies than in the single half-wave dipole.
9.9.3 BALUN
A dipole in isolation or one in parallel with a ground plane is balanced with
respect to ground, since the dipole has equal currents in its two arms. To
avoid unbalancing the radiating structure, the dipole should be fed using
a balanced transmission line such as a two-conductor line. For frequencies
above ∼ 200 M Hz, the feed line of preference is a coaxial line which is an
unbalanced line and will therefore lead to distortion of the radiation pat-
tern. The coaxial line must therefore be connected to the antenna through
a structure that adapts the unbalanced to the balanced structure. Such
a device is known as a BALanced-to-UNbalanced (BALUN) impedance
transformer. In TV systems with a 75 Ω input impedance and using a
folded dipole antenna, the balun is usually a 4 : 1 impedance transformer.
The balun may take various forms the most common ones being the
transmission line transformer of Fig. 9.43 and the structure shown in Fig.
9.44. For the balun in Fig. 9.44, the centre conductor is connected in parallel
with the short-circuited outer conductor of a coax line a distance λ/4 away
from the feed point. The short-circuit is thus seen as an open-circuit at the
feed point.
9.9 Feed techniques for wire antennas 257
EM wave
dl
I
2 2 3
I R C
Va E 4
1 1
dl wire loop
(a) Distributed source (b) Localised source
applied field is distributed over the entire circuit through the mechanism
of the incident electromagnetic waves. In this case, the integration path is
the closed curve C enclosing the circuit and we can therefore write,
−
→ −
→
Va = E a · d l
C
When the sources are localised, the integral path is across the source. The
applied voltage is in this case independent of the path of integration as
shown in see Fig. 9.45 where
1 2
−
→ −
→ −
→ −
→
Va = Ea · d l = − Ea · d l
2 1
The voltage drop across the internal impedance is the complex contribu-
tion of the internal resistance and reactance. The internal impedance per
unit length of the circuit Zi′ is related to the tangential surface electric field
intensity Em
Em
Zi′ =
I
where I is the total current through the circuit. The total internal im-
pedance is then obtained as
−→
1 J −
→ 1 − → −
→
Zi I = ·d l = E m · d l = Zi′ dl
I σ I
where the integral is carried out only over the conductive portions of the
circuit path.
If we select the path of integration to include only the contributions from
the magnetic flux external to the conductor and ignore gaps such as those
at the localised source and capacitors, we can use a closed-path formulation
with negligible error as shown in Fig. 9.46. Thus,
− →
∂A − → ∂ − → − →
·d l = A ·d l
∂t ∂t
260 9. Antenna theory
integration path
V +
-
wire loop
This is the contribution from the flux threading the chosen circuit path
and is valid for a circuit of a size much smaller than the signal wavelength.
−
→
In Fig. 9.45, we integrate ∇Φ over the conductive portion of the circuit
path between points 3 and 4 and assume charges due to the capacitance are
negligible at the source. Recalling that the integral of a gradient function
around a closed path is zero, then the voltage across the capacitor is related
to the charge Q on the capacitor as
Q
V =
C
Hence
3
−
→ −
→ Q 1
∇Φ · d l = Φ3 − Φ4 = = I dt
C C
4
where harmonic time dependence has been assumed as well. For a filamen-
tary current,
−
→ µI exp (−jkr) −→
A= dl ′=0
4πr
9.10 Circuits of a size comparable to signal wavelength 261
where r is the distance of the observation point from the current element
under consideration. Recalling that
−
→ − → − →
B =∇×A
and
−
→ ω−→ −→ − → −
→
E = −j 2 ∇ ∇ · A − jω A
k
we can therefore write
−
→ −
→ µI cos kr − j sin kr −
→ − →
E a · d l − jω d l ·d l ′ =0
4π r
Thus, even with a current that is in phase throughout the circuit, a finite
value of kr leads to a real and imaginary contribution. The real part of the
contribution is related to energy radiation and the imaginary part to an
inductive reactance.
We therefore see that, in a circuit of a size comparable to the signal
wavelength, the retardation effect from one part of the circuit to another,
causes a phase shift which leads to induced effects with real and imaginary
parts unlike the purely reactive effects in a circuit that is small relative to
a wavelength.
arms of the rhombus are usually 2-7λ long and the antenna placed 1-2λ
above ground.
A wave propagating in a waveguide operated near cut-off, has a phase
velocity approximately equal to the velocity in the unbounded medium. The
surface waves propagating on a dielectric guide may then excite appreciable
fields in space. This structure is also shown in Fig. 9.49.
FIGURE 9.50. The electromagnetic horn and leaky waveguide radiating struc-
tures
The antenna field can be divided into the near- or Fresnel field and the
far- or Fraunhofer field. The far-field is important in wireless communica-
tion and we will confine ourselves to the study of this field. In the far-field
region, we assume
• The differences in the lengths of the radial vector − →r at different
points on the antenna as observed from an observation point P in
the far-field are insignificant in describing the magnitude of the field
contributions. This implies that the wavefront far away from the ra-
diating structure can be regarded as planar.
• The differences in the lengths of the radial vector are important when
regarding the phase differences between different points on the an-
tenna as observed from the far-field zone.
• All field components decreasing with distance from the radiating
structure faster than 1/r are insignificant in the far-field (Sommerfeld
radiation condition).
Consider the volume V ′ enclosing a point source a short distance r′ away
from the origin of the spherical coordinate system as shown in Fig. 9.52.
The vector potential at P at r >> r′ from the source is given by
−
→
−
→ J exp (−jkR) ′
A =µ dV
4πR
V′
where R is the distance of P from the point source. Applying the cosine
rule,
R = r2 + r′2 − rr′ cos ϕ ∼ r − r′ cos ϕ
where ϕ is the direction angle between −
→
r and −
→r ′ . Hence
−
→ µ exp (−jkr) −
→
A= J exp (jkr′ cos ϕ) dV ′
4πr
V′
266 9. Antenna theory
so that
−
→ µ exp (−jkr) −
→
A= N
4πr
−
→
Generally, N may have components in any direction so that
→ µ exp (−jkr)
−
A= rNr + φ
θNθ + φN
4πr
In spherical coordinates, we have
−
→ −
→ − → r ∂ ∂Aθ θ 1 ∂Ar ∂
B = ∇×A = (Aφ sin θ) − + − (rAφ )
r sin θ ∂θ ∂φ r sin θ ∂φ ∂r
∂
φ ∂Ar
+ (rAθ ) −
r ∂r ∂θ
Field components not decreasing faster than 1/r are thus,
1 ∂ jk
Hθ = − (rAφ ) = Nφ exp (−jkr)
µr ∂r 4πr
and
1 ∂ jk
Hφ = (rAθ ) = − Nθ exp (−jkr)
µr ∂r 4πr
−
→ −
→ −
→ − → −
→
Since E = ∇ ∇ · A − jω A , then only
jωµ
Eθ = − Nθ exp (−jkr)
4πr
9.12 Antenna analysis by field computation 267
and
jωµ
Eφ = − Nφ exp (−jkr)
4πr
decrease at a rate not faster than 1/r. Thus, there are no radial electric
fields.
The time-average Poynting vector is then
1 ' ( η 2 2
Pr = Re Eθ Hφ∗ − Eφ Hθ∗ = 20 2 |Nθ | + |Nφ |
2 8λ r
where
Eθ Eφ
η0 = =−
Hφ Hθ
and the time-average radiated power is
π 2π π 2π
η0
PT = 2
Pr r sin θ dθ dφ = 2 |Nθ |2 + |Nφ |2 sin θ dθ dφ
8λ
0 0 0 0
−
→
The Poynting vector P gives the actual power density at any point. We
define a radiation intensity K independent of the radial distance r from the
source as the power radiated in a given direction per unit solid angle and is
hence the time-average Poynting vector on a unit-radius sphere. Therefore,
we can write
η
K = 02 |Nθ |2 + |Nφ |2
8λ
so that
π 2π
PT = K sin θ dθ dφ
0 0
If all the currents are in one direction which may be taken as the axis of a
−
→ −
→
set of spherical coordinates, then N and hence A has only one component,
say Nz . Thus
Nφ = 0, Nθ = −Nz sin θ
and
η0 2 2
K= 2 |Nz | sin θ
8λ
in this case. With the currents directed azimuthally about an axis however,
the axis could be taken as the axis of a spherical coordinate system. If
symmetry in φ exists,
η
K = 02 |Nφ |2
8λ
and
π
PT = 2π K sin θ dθ
0
268 9. Antenna theory
dSθ
dSφ dSr
rdθ
dθ dr
θ
y
φ
rsinθ
rsinθdφ
x
Note that as shown in Fig. 9.53, the Cartesian coordinates (x, y, z) are
expressed in terms of the spherical coordinates as
x = r cos φ sin θ
y = r sin φ sin θ
z = r cos θ
and
Nφ = −Nx sin φ + Ny cos φ
If (θ, φ) are the angular coordinates of P and θ′ , φ′ those of the source,
then
cos ϕ = cos θ cos θ ′ + sin θ sin θ′ cos φ − φ′
9.13 Waves guided by biconical guide 269
1 ∂ 1 ∂Er
(rEθ ) − + jωµHφ = 0
r ∂r r ∂r
1 ∂
(sin θHφ ) − jωǫEr = 0
r sin θ ∂θ
1 ∂
− (rHφ ) − jωǫEθ = 0
r ∂r
and
Er = 0
η
rEθ = {A exp (−jkr) + B exp (jkr)}
sin θ
1
rHφ = {A exp (−jkr) − B exp (jkr)}
sin θ
where
Eθ
η=
Hφ
and ∂/∂φ = 0 because of the φ-symmetry.
270 9. Antenna theory
0
π−θ 0
π−θ
1
V = Eθ rdθ = −η {A exp (−jkr) + B exp (jkr)} dθ
sin θ
θ0 θ0
θ0
= 2η ln cot {A exp (−jkr) + B exp (jkr)}
2
where l is the length of the transmission line and c the speed of light, the
system is said to be a low frequency circuit and the small perturbations in
the currents then referred to as quasi-static phenomena. The low frequency
circuit, can then be modeled with an RLC circuit. Additional issues arise
in high frequency systems such as propagation delays, reflections and radi-
ation, which are no longer negligible.
EMC
Emission (source) Susceptibility (victim)
EME EMS
CE (conducted) CS (conducted)
RE (radiated) RS (radiated)
• Common-mode interference
• Differential-mode interference
• asymmetrical interference.
C
- I
Common-mode interference which is mainly a high frequency (≥ 1 MHz)
phenomenon, occurs between all lines in a cable and the reference potential.
If the potentials on the conductors of a two-conductor transmission line are
V1 and V2 resp., the common-mode potential is defined as
V1 + V2
VCM = − Vref
2
where Vref is the reference potential. Common-mode currents cannot occur
when the paths between the conductors are ideal but this is not possible
in practice.
D
- I
Differential-mode interference occurs symmetrically between two lines
and is primarily a low-frequency (below several hundreds kHz) phenom-
enon. If the conductors of a two-conductor transmission line have potential
V1 and V2 resp., the differential-mode potential is defined as
VDM = V1 − V2
A I
Asymmetrical interference describes the interference between one line
and a reference potential.
In the case of radiated interference, we can have differential radiation
currents which are a natural consequence of the physical layout of the
circuit. Different PCB tracks (traces) may carry different current densities
with differential-mode currents leading to H-field radiation. Common-mode
radiation currents mainly lead to E-field radiation due to losses through
the cables and PCB tracks.
10.1 Basic information on EMC 275
VDM
VCM V1 V2
• Mode of propagation
• Coupling mechanism
from the common-mode source. In general, this means that the maximum
measured emission must be 30 µV m−1 .
The common-mode currents are related to the electric field as
60ICM
E=
r
where ICM is the common-mode current and r, the distance from the
common-mode source where the field is measured. Emission may be due to
capacitive coupling between the lines or from a line to ground. Fig. 10.3
shows a test setup for the common-mode current measurement. Differential-
mode currents could be measured with the same setup when one line is
wrapped around the transformer once.
Common-mode currents cause unwanted radiation and can be suppressed
by placing a choke or clamp such as ferrite beads, in the current path. Dig-
ital ribbon cables can give significant problems with common-mode cur-
rent radiation and clamps may be the more feasible approach to suppress
this radiation. Note that ferrite beads or clamps will have no effect on
differential-mode radiation.
D
- C
M
The far-field H-field is related to the far-field E-field as
E
H=
η
where the symbols carry their usual meaning. The H-field is measured at
0, 1 m from the differential-mode source and the E-field at 3 m, which forms
the legal length of cables supplied with consumer products for operation
below 16-300 MHz band.
10.1 Basic information on EMC 277
The H-field can be measured using a current loop made out of a coaxial
cable. The outer conductor is removed and a 50 Ω chip resistor connected
between a loop formed by the inner conductor. The H-field induces a volt-
age across the resistor which is measured by the measuring system. The
sensitivity of the probe is high for low frequencies but becomes indepen-
dent of frequency above about 300 M Hz. The reason for this is that while
the voltage increases with frequency, the inductive impedance of the probe
increases in the same direction preventing any further increase in current.
The E-field measured at a distance r = 3 m from the source is related to
the 0, 1 m H-field as
−
→
E = 2, 7 × 10f × H
where f is the frequency in Hz. A measurement setup for the H-field is
shown in Fig. 10.4. Fig. 10.5 shows the general EMI environment.
ICM
IDM
Cpara R Cpara
ICM
IDM
ICM
through the connecting conductors and return through the ground connec-
tion. Since the parasitic capacitive impedances tend toward short-circuits
at high frequencies, the coupling effects in the connecting cables and the
equipment itself will increase with frequency. Hence, common-mode inter-
ference currents become dominant in the M Hz range. Asymmetrical inter-
ference consists of both symmetrical (differential-mode) and common-mode
components.
10.1.6 Cross-talk
When the coupling path between the EMI source and victim is short, we
have cross-talk. Examples are when both the source and victim share a
common impedance which then forms the path for the cross-talk signal as
well as inductive or capacitive coupling through the magnetic or electric
fields resp. Cross-talk must be particularly addressed when designing PCBs
or multi-cable assemblies.
Fig. 10.7 illustrates a network in which the output and input circuits
share a common impedance Zm . The currents in the two circuits flow
through Zm to the reference node which is usually the circuit ground. By
superposition, the voltage established across the load ZL2 by the current
from Vs1 is given by
ZL2 Zm // (ZL2 + Zs2 )
VZL2 = Vs1
ZL2 + Zs2 Zm // (ZL2 + Zs2 ) + ZL1 + Zs1
Zm ZL2
= Vs1
Zm (ZL2 + Zs2 ) + (ZL1 + Zs1 ) (ZL2 + Zs2 )
and when Zm is small, we can write
VZL2 Zm ZL2
=
Vs1 (ZL1 + Zs1 ) (ZL2 + Zs2 )
We therefore see that cross-talk depends on Zm and there will be no cross-
talk if Zm = 0.
The effect of cross-talk on the signal is specified through the distortion
ratio defined as the voltage drop across ZL2 for Zm = 0 to the resulting
drop when Zm = 0. Thus,
[VZL2 ]Zm =0 Zs1 + ZL1 Vs2
=
[VZL2 ]Zm , small Zm Vs1
I1 Zm I2
ZL1 ZL2
• Circuit emission
• Circuit immunity
• Circuit conductors
• Component placement in the circuit
• Group digital components and use the slowest components the appli-
cation can tolerate.
• Place cables and connectors on one side of the PCB and as close to
each other as possible
• Keep the straight parts of the traces as short as possible and use 45◦
bends in the traces whenever required. To illustrate, a trace of length
0, 6 mm resonates and hence radiates strongly at ∼ 5 GHz, since
c
f=
l
where l is the trace length, c is the speed of light and f is the oper-
ating frequency. Reducing the longest trace to 0, 18 mm pushes the
resonance frequency to ∼ 16, 6 GHz.
284 10. Fundamentals of electromagnetic compatibility
3500
2500
voltage in µV
2000
1500
1000
1,705
500
250
0
0,01 0,1 1 10 100
frequency in MHz
65
60
60
56,5
54 54
electric field in dBµ V/m
55
50 49,5
46
45 43,5
40
40
216
35
88
30
30 100 1000
frequency in MHz
FIGURE 10.11. A comparison of the FCC limits for radiated emissions in Class
A and B digital devices (measurements at 3 m)
45
electric field in dBµV/m
40 40 FCC
37
CISPR-22
35 36,5
34
30
29,5 216
88 230
25
30 100 1000
frequency in MHz
FIGURE 10.12. A comparison of the FCC and CISPR-22 regulations for radiated
emissions in Class A devices (measurement at 30 m)
10.2 The EMC guidelines 289
45
43,5
37
35 35,5
33
30
29,5
230
88
216
25
30 100 1000
frequency in MHz
FIGURE 10.13. A comparison of the FCC and CISPR-22 regulations for radiated
emissions in Class B devices (measurement at 10 m)
500 kHz. The rationale behind this extension is the increasing importance
of switching power supplies over their linear analogues occasioned by their
higher efficiency and light weight. The CISPR-22 regulations for conducted
emissions are specified using a quasi-peak detector (QP) and an average
detector (AV). FCC conducted and CISPR-22 and FCC radiated emissions
limits require the use of a quasi-peak detector. Figs. 10.14 and 10.15 illus-
trate these differences for Class A and B devices.
Because of the great importance of EMC issues in military applications
to ensure that missions will not be compromised, the military has developed
military EMC regulations. In addition to regulating both conducted and
radiated emissions, the military regulates susceptibility which is very im-
portant in military applications since military equipment must be immune
to outside interference. The military EMC regulations are more stringent
than the FCC or CISPR regulations and cover a much wider frequency
range with several subdivisions within this range. In addition, the military
may require that EMC requirements be waived for certain applications if
this is deemed necessary to mission success. On the other hand, CISPR
and FCC regulations cannot be waived for commercial products.
80 79,5
75
73
CISPR-22 (QP)
voltage in dBµV
70 69,5
FCC
66
65
60
60 0,15 CISPR-22 (AV)
0,5 1,705
55
0,1 0,45 1 10 30
frequency in MHz
FIGURE 10.14. A comparison of the FCC and CISPR-22 regulations for con-
ducted emissions in Class A devices
70
66
65
voltage in dBµV
60 CISPR-22 (QP)
56
55
50 CISPR-22 (AV)
48 FCC
46
45
0,1 0,15 0,45 1 5 10 30
0,5
frequency in MHz
FIGURE 10.15. A comparison of the FCC and CISPR-22 regulations for con-
ducted emissions in Class B devices
10.2 The EMC guidelines 291
repeatable. For radiated emissions, the FCC specifies that the measurement
of radiated and conducted emissions be performed on the complete system.
All cables interconnecting to peripheral equipment must be connected and
the system must be in a typical configuration and both the cable and the
system be configured in a way that maximises the emissions. For example,
if the unit under test has interior wire harnesses, these harnesses must be
configured in such a way that of all the possible ways the unit can be as-
sembled with those wire harnesses, this way causes the maximum radiated
emissions. This ensures that for a mass production unit, the worst-case
scenario is taken into consideration.
Testing standards set forth by the FCC for radiated emissions testing
are very specific and difficult to automate. Radiated emissions are to be
measured at 10 m for Class A and at 3 m for Class B devices. These
measurements are to be made over a ground plane using a tuned dipole
antenna at an open field test site. In addition, the tests are to be made with
the measurement antenna in both the vertical and horizontal polarisation
positions. During product development, most companies test their products
in a semi-anechoic chamber, a shielded room with radio frequency absorbing
cones on the walls and ceiling. This semi-anechoic chamber shown in Fig.
10.16 simulates an open field test site by minimising reflections from the
walls and eliminates any ambient signals that may be present in an open
field environment.
Companies sometimes simplify the FCC test procedure using a broad-
band antenna such as a log-periodic or discone antenna. Such antennas
are desirable since, unlike a tuned dipole, their length does not need to
be adjusted with each frequency change allowing companies to test their
products using a frequency sweep rather than measuring at each frequency
separately and adjusting the dipole lengths with each measurement. An
additional test requirement for radiated emissions testing is that the band-
width of the measurement receiver must be at least 100 kHz. A wide re-
ceiver bandwidth ensures that the test picks up emissions from broadband
sources such as the arcing at the brushes of a DC motor in addition to the
intended narrowband signals such as clock signals.
The FCC and CISPR test procedures require that the measurement re-
ceiver use a quasi-peak detector while typical spectrum analysers use peak
detectors. The quasi-peak detectors ensure that fast changing, momentary
signals such as randomly occurring spikes will not charge up the detector to
as high a level as periodic signals would. The FCC is concerned primarily
with the more significant and frequent emissions that would cause inter-
ference with radio and wire communications rather than with randomly
occurring one-time signals.
292 10. Fundamentals of electromagnetic compatibility
shielded room
DUT
(3 or 10 m)
scan height
1 - 4 m for
both vertical
& horizontal
polarisations
ground plane
spectrum analyser
or receiver
LISN
IP L1
C1 C1
IP L1
L IN IN
product N to AC
under net
G
test C2
50 Ω 50 Ω C2
R1 R1
VP VN
green wire
L
product N to AC
under net
G
test
green wire
short circuits while the inductors appear as open circuits and the equivalent
circuit is as shown in Fig. 10.19. The currents on the neutral and phase
lines can therefore be isolated and measured at the 50 Ω resistors. The
currents on the phase and neutral lines have no path onto the AC power
net.
LISN
L
product N to AC
under net
G
test C2
50 Ω 50 Ω
R1 R1
VP VN
green wire
11.1 Introduction
A standard is a document established by consensus and approved by a recog-
nised body which provides for common and repeated use, rules, guidelines or
characteristics for activities or their results aimed at achieving an optimum
degree of order in a given context.
Standards should be based on consolidated results of science, technology
and experience with the aim of promoting optimum community benefits.
While voluntary documentary standards are crucial in facilitating and reg-
ulating industry and commerce, this role is often grossly overlooked. Stan-
dards do not dictate how products are designed or manufactured but may
be used to influence safety requirements through incorporation into regu-
lations. It is important to remember that standards derive their legitimacy
from the voluntary and consensual nature of their development process
however.
Most standards provide utility, i.e. when one uses a standard, it is often
because it provides a reliable though not necessarily the only or the best
solution to a given problem. In research, standards are one of the most
important tools for taking new technology to market and in transform-
ing research findings into guidance documents thereby providing a bridge
between research and industry. This is a crucial link for successful com-
mercialisation of research outcomes. By providing this bridge, standards
are invaluable as a tool for promoting innovation and commercialisation
through
298 11. Standards and standardisation
• A process that seeks to take into account views from all parties con-
cerned and reconciliation of any conflicting arguments
11.3.2 Standardisation
This is the activity of establishing a standard with regard to the actual po-
tential problems, provisions for common and repeated use aimed at achiev-
ing the optimum degree of order in a given context. In particular, standard-
isation consists of the process of formulating, issuing and implementing the
standard. Important benefits from standardisation are
Working draft
NWIP received NWIP approved submitted by
from member and registered in proposer, further
body database developed by
experts
3 month ballot by
members of relevant
TC.
Approval if at least 5 3 month review Consensus
members agree to and ballot by among
participate and 50% members of TC experts
of those voting support - may take 2 to
proposal 4 months by
agreement with TC
Resolution
of comments
PAS & TR published received
if at least simple majority
of voting members approve
document. TS published if
at least 66% of voting members
approve document. Same
requirement for formal route
Published as IS
Most organisations such as SMEs, NGOs or consumer groups may find this
untenable.
While physical meetings might provide direct and immediate input to a
project, participation via correspondence and other electronic means can be
more effective since it can be accomplished at a time and place convenient
to the individual.
Is the output of a
new and repeatable
technique for
- identification
- characterisation
- manipulation
- preparation No
Ignore output
- verification, etc.
Or
a modification of an
existing technique or
procedure to allow its
use under previously
untested conditions?
Yes
No
Yes
Yes
Yes
Yes
Yes
to provide a project leader (PL) who will work with experts nominated by
other members to develop the document to the final working draft (FWD)
stage. The country proposing the project is also expected to provide a draft
or outline of the document.
In ISO, approval of a NWIP requires that at least five (5) P members
agree to participate in the work by nominating experts together with at
least a simple majority of those P members voting being in favour of the
development. Abstentions in all ballots are not counted. The approved work
items (AWIs) are then developed by nominated experts acting in their
individual capacity within a project group (PG) under the leadership of
the PL. Project groups covering the same general area of work are usually
grouped together under working groups (WG).
Where a large volume of work is undertaken in an area that does not im-
pact on other parts of the TC’s work and can be managed more effectively
as a separate group, a subcommittee (SC) can be formed to handle the
work. Subcommittees act simultaneously and are essentially independent
of their parent TC. Both the TC and its SC can establish a WG, approve a
NWIP , develop documents to the committee draft (CD) stage etc. Work is
undertaken electronically in general but most PGs have a physical meeting
during the once or twice yearly TC planning meetings.
Once experts reach consensus on the contents of a document, their final
working draft (FWD) is submitted to the membership of the TC or SC
for approval as a CD. This is achieved through a three-month ballot of P
members with recent changes to the ISO/IEC directives permitting two-,
three- or four-month balloting periods depending on the circumstances. In
this ballot, P members are required to provide a national position through
their mirror TCs, i.e. approve, reject or abstain. Where appropriate, this
may be accompanied by comments on the document contents. Comments
on documents under ballot can be of a general, editorial or technical na-
ture and are submitted on a standard comments template. Comments must
be accompanied by a recommendation on how each comment is to be ad-
dressed.
Committee drafts are approved if they garner support of at least two-
thirds of those voting with abstentions not counted. Depending on the type
of document originally proposed, the CD can be published as a PAS, TR
or TS. Only a simple majority of voting P members is required to approve
a PAS or TR. A formal route toward full international standard is via
preparation of a DIS. Whichever route is chosen, it is necessary that any
comments received are satisfactorily resolved before publication or progress
to the DIS stage. This is normally done at a comments resolution meeting
under the leadership of the PL. Once comments have been addressed to
the satisfaction of all commenting members, the document is either sent
for publication or submitted for its DIS ballot. The DIS ballot runs over a
five-month period during which each ISO member has an opportunity to
review and comment on the document. Approval of the document requires
11.3 How standards are developed 307
Requirements for approval of an FDIS are the same as those for a DIS. In
CEN, all members are automatically members of all TCs. However, to take
part in the work of a given TC, a member must appoint a representative
to the TC in order to access relevant committee documents. Unlike in ISO
which operates on the principle of one member, one vote, CEN operates
a weighted voting system giving larger economies a higher percentage of
the votes. Weighted voting applies to NWIP ballots for TS and EN and to
approve ballots of both types of publications.
Inventive step refers to something that is not obvious to those skilled in
the art, i.e. something that another expert would not do as a matter of
course as a consequence of their training.
All standards have the following, specific to each document
(a) A unique number that includes the source of the document and/or its
status as well as year of approval. For example, ISO 9001:2008, EN
13485:2003, EN 1071-3:2001, ISO TS 80004-1:2010, CWA 14243:2003
(b) Title consisting of up to three elements
(i) Introductory element giving the general field to which the docu-
ment belongs. Often, the title of the committee responsible for
the document.
(ii) Main element indicating the principal subject of the document
(iii) Complementary element giving specific aspect of the main ele-
ment addressed in the document
(c) Scope defining the subject of the document and aspects covered. This
should indicate any limits to applicability of the document or partic-
ular parts of it and should reflect the title while the title reflects the
scope. It is important to agree on the scope as early as possible for
all experts to be fully aware of what they are trying to develop.
308 11. Standards and standardisation
(d) Normative references, i.e. a list of other standards that form or parts
of which form an integral part of the document and without reference
to which, it is impossible to comply with the standard established by
the document.
(g) Apparatus
(h) Procedures written in the imperative form and not the narrative
(k) Annexes
11.3.7 Consensus
There are minimum ballot thresholds for a document to be published as a
TS, TR, PAS or IS. In CEN, the minimum requirement for EN or TS is 71%
of the weighted votes of members voting. A TR requires a simple majority of
those members voting. CEN does not have an equivalent to PAS. However,
the overriding objective should always be to reach a consensus amongst all
involved. Remember that by consensus, we mean that there is no sustained
opposition to substantial issues by any important part of the concerned
interests but not necessarily unanimity.
The principle of consensus is vital to standardisation as it imparts legiti-
macy to the process and the documents produced. People and organisations
use standards not through coercion but because they are in general agree-
ment with their content. Even without requiring unanimity, consensus does
require those involved to respond from a position of knowledge not igno-
rance. Members therefore have the responsibility to ensure that they only
support or reject documents when they have sufficient national expertise
to make informed judgement. Members should always abstain when the
necessary expertise does not exist. Note that if a member cannot provide
11.3 How standards are developed 309
• Low viscosity, low density, low volatility, low solvent power but high
flash point
• Paraffins
• Asphatics
• Mined
After extraction from the ground, oils are refined to remove metals and
other impurities. In the refined form, oils can have dielectric strengths be-
tween 400 and 500 kV cm−1 which could go as high as 1000 kV cm−1 with
further refining. The dielectric strength will however fall during use due
to a combination of factors which include aging and contamination. Most
insulating oils have an ǫr ∼ 2, 2 and power factors below 0, 0005 at 50 Hz
and 20 ◦ C.
F
Fluorinated liquids are characterised by high chemical stability allowing
application at temperatures up to 200 ◦ C. They have low loss tangents
of < 0, 0005, high resistivity 1014 -1018 Ωcm and ǫr ∼ 1, 8. Their cooling
properties are better than those of mineral oils and silicones. These liquids
have therefore found applications in some electronic equipment such as
radio transformers in electronic systems. However, fluorinated liquids suffer
degradation through moisture absorption and are therefore to be found
primarily in sealed equipment.
S
Synthetic hydrocarbons have no major advantages over the natural min-
eral oils apart from the relative ease of guaranteeing their chemical compo-
sition. They have found use mainly in smoothing capacitors for DC power
and in high pressure gas-filled cables.
O
Esters are useful mainly at high frequencies owing to their relatively low
dielectric losses and high dielectric constants, properties that are desirable
in many high frequency capacitors. Chemically, esters are manufactured
by reacting acids and alcohols but some occur naturally in plants, e.g.
castor oil. Since esters have a high moisture absorption capacity, they are
to be found mainly in sealed devices. Esters have dielectric constants in the
range 2-3, 5 at low frequencies and power factors at 50 Hz of ∼ 0, 05 which
decreases to about 0, 001 at 1 MHz.
V
Vegetable oils may be classified under
This classification depends on the ease with which the oils form films
resulting from the ease with which they can be oxidised and polymerised
12.1 Liquid insulation 313
when exposed to air. From a dielectric point of view, drying and non-
drying oils are important with linseed oil finding use in the formulation
of insulation varnishes and castor oil as a plasticiser in insulating resin
compositions employed in structural insulation materials and in coating
compositions in electrical equipment.
A! (
)
These are chlorinated aromatic (benzene-based) hydrocarbons. They are
fire-resistant and have been used in transformers and capacitors that re-
quire high dielectric constants. The breakdown voltage is fairly high and
lies between 80 and 180 kV cm−1 and power factors in the range 0, 02-0, 05.
Askarels have resistivities in excess of 10 GΩcm and dielectric constants in
the range 4, 8-5, 3. However, askarels are now recognised as accumulative
persistent ecological contaminants and their use is now largely discontin-
ued.
S
Silicone fluids are polyorganic silicones with a unique combination of
properties. They are clear liquids with an oily consistency and have ǫr ∼
2, 2-2, 8. Their viscosity can vary over a wide range but silicon fluids are
highly stable at high temperatures. Polyorganic silicones have electrical
properties that are largely independent of frequency and temperature mak-
ing them excellent coolants in pulsed radar systems, aircraft and radio
transformers. However, silicones are detrimentally affected by arcing mak-
ing them unsuitable for use in switch-gear.
When choosing a liquid dielectric, several factors such as environmental
factors and space-saving properties must be considered beside their physical
properties. Cost and chemical stability are important factors as well. Hence,
mineral oils remain the dominant liquid insulators in active application to
date.
B T
Mechanically pure liquids contain no moisture or other impurities. They
may however contain dissolved gases which tend to collect around the elec-
trodes in the form of bubbles. The gas bubbles may then become ionised
and elongated from polarisation, leading to breakdown.
L G
B!&
where P is the partial pressure of the solute gas above the solution in atm
of absolute pressure, c is the solute concentration in solution in mol L−1
(molarity) and kH is Henry’s law constant in L atm mol−1 , atm per mole
fraction or P a m3 mol−1 . Thus,
P = kH c
Table 12.1 presents kH values for common solute gases in water at 298 K.
Since there are various formulations of Henry’s law, kH has different di-
mensions in these formulations and some of these are presented in Table
12.1 as well. Note that kH,cp and kH,pc are reciprocals of each other.
where, caq is the concentration of solute in mol L−1 , Lsol is volume of
solution in L, Pgas is the gas partial pressure above solution in atm of
absolute pressure and atm is atmospheres of absolute pressure.
Henry’s law is a limit law and applies only to dilute solutions. The range
of applicable concentrations becomes increasingly narrow with increasing
divergence of the system from ideal behaviour, i.e. roughly as the solute
and solvent become more chemically different. Further, Henry’s law applies
only when the solute and solvent do not chemically react. Therefore, since
CO2 rapidly forms hydrated CO2 and finally carbonic acid (H2 CO3 ) in
water, Henry’s law does not apply in the strict sense.
12.1 Liquid insulation 315
TABLE 12.1. Some forms of Henry’s law and constants for some common gases
in water at 298 K
P c Pgas caq
equation kH,pc = cgas kH,cp = Pgas
aq
kH,pn = naq kH,cc = cgas
aq
Lsol ·atm molgas atm·molsol
dimension molgas Lsol ·atm molgas dimensionless
−3 4
O2 769, 23 1, 3 × 10 4, 259 × 10 3, 180 × 10−2
CO2 29, 41 7, 8 × 10−4 7, 099 × 104 1, 907 × 10−2
H2 1282, 05 3, 4 × 10−2 0, 163 × 104 0, 8317
N2 1639, 34 6, 1 × 10−4 9, 077 × 104 1, 492 × 10−2
He 2702, 7 3, 7 × 10−4 14, 97 × 104 9, 051 × 10−3
Ne 2222, 22 4, 5 × 10−4 12, 30 × 104 1, 101 × 10−2
Ar 714, 28 1, 4 × 10−3 3, 955 × 104 3, 425 × 10−2
CO 1052, 63 9, 5 × 10−4 5, 828 × 104 2, 324 × 10−2
where ρmelt and ρgas are the densities of the solute gas in the melt and gas
phases resp., µmelt,E and µgas,E are the excess chemical potentials of the
solute gas in the melt and gas phases and
1
β=
kT
where k = 1, 38 × 10−23 JK −1 is Boltzmann’s constant.
Raoult’s law is similar to Henry’s law and states that the partial pressure
of a component in a solution is proportional to the concentration of that
component in the solution. Using the molar fraction n = N/NA , where
N is the number of molecules of the component in solution and NA =
6, 02 × 1023 mol−1 is Avogadro’s number, to express the concentration of
the solute in solution, we can express Henry’s law in the form
P = kH n
and
P
lim = kH , Henry’s law
n
n→0
• Resins or polymers
• Composite materials
• Thickness
• Finish
• Porosity
• Tensile strength
• Tearing resistance
Paper insulators are prone to water absorption and are therefore usually
impregnated with mineral or vegetable oil. This treatment lends these insu-
lators an ǫr of ∼ 3, power factor in the range 0, 0017-0, 0025 at densities of
0, 75 - 1, 7 gcm−3 and a conductivity of 10−17 Scm−1 at 25 ◦ C. The impulse
breakdown voltage in paper insulators is approximately 3000 kV cm−1 . The
usual application areas of these insulators are hence
F M
Fibrous materials combine strength and durability and therefore are to
be found in many applications. Fibrous materials include cotton, silk and
synthetic fibres such as jute, viscous rayons, nylon and fibre glass. They are
flexible, cheap, possess high mechanical strengths and are easy to process.
They however have low dielectric strengths and are highly hygroscopic.
However, these materials can have improved electrical properties with an
oil impregnation.
The dielectric constants of the fibrous materials range from 3-8 in a
perfectly dry condition but many of their electrical properties are frequency
dependent. For example, ǫr decreases with increasing frequency and shows
a fair amount of anisotropy. Fibrous insulators are affected by moisture and
their resistivity will decrease with increasing moisture content.
These insulators find such applications as conductor insulation, winding
and coil insulation, backing for mica etc. Inorganic fibres such as asbestos
and fibre glass withstand high temperatures but have low elasticity and
machinability.
I(
, C
, F
B
M
These materials can be classified as waxes, insulating varnishes, filling
and bonding materials. Waxes are soft in nature and have low mechanical
strength but can be used as impregnating materials. Artificial waxes such as
paraffins, ceresine and halowax have excellent insulating properties. How-
ever, natural ones such as beeswax enjoy limited insulation applications.
Insulating varnishes are solutions of resins, bitumen, etc. and are mainly
used as impregnating, coating or adhesive materials. Bonding materials are
used to seal off separately insulated parts while adhesive tapes are to be
found extensively in wiring.
R
12.2 Solid-state insulating materials 319
Resins are mainly natural and synthetic polymers (plastics) and can be
classified on the basis of their behaviour under thermal stress as thermo-
plastic or thermosetting resins. Thermoplastic resins can be formed into
fibres and will retain their properties after melting and solidification on
cooling. However, they have low mechanical strength. Thermosetting resins
on the other hand, lose their properties after a heat treatment and have
poor insulating properties but good mechanical properties. Natural resins
find extensive application and include such materials as shellac, resins,
copal and amber. Shellac is obtained from tropical trees and has good ad-
hesive properties. Synthetic resins include polystyrene, polyethylene, PVC,
acrylic, polyimide and resins derived from cellulose, polyester epoxy resins
etc.
N
S
R
Natural rubber is a polymer of hydrocarbon isoprene and has poor ther-
mal properties and is prone to oxidation in air. It has a low loss tangent,
∼ 0, 01-0, 03 and an ǫr of ∼ 2, 5-5. Synthetic rubber includes butadiene
rubber, butyl rubber, polychloroprene rubber and silicon rubber and find
use in wiring and aircraft cabling.
• Surface conditions
FIGURE 12.2. A solid dielectric capacitor with cavities and its electrical equiva-
lent circuit
Hence
Vc 1 Va
Ec = =
tc 1 − (1 − tm /tc ) /ǫr tc
which is significantly higher than Va /tm and hence will cause partial dis-
charge in the cavity. The accumulated effect is to heat up the material
gradually until there is total breakdown. The life of the material with in-
ternal discharges will therefore depend on the number of electrical stress
application cycles.
C
E B!&