Antennas & Wave Propagation Guide
Antennas & Wave Propagation Guide
WAVE PROPAGATION
2[IRUG8QLYHUVLW\3UHVV
Contents
Preface
Symbols xv
ix
2[IRUG8QLYHUVLW\3UHVV
x Contents
2[IRUG8QLYHUVLW\3UHVV
Contents xi
2[IRUG8QLYHUVLW\3UHVV
xii Contents
2[IRUG8QLYHUVLW\3UHVV
Contents xiii
References 397
Index 399
2[IRUG8QLYHUVLW\3UHVV
CHAPTER 1
Electromagnetic Radiation
Introduction
Most of us are familiar with cellular phones. In cellular communication
systems, there is a two-way wireless transmission between the cellular phone
handset and the base station tower. The cell phone converts the audio sig-
nals into electrical form using a microphone. This information is imposed on
a high frequency carrier signal by the process of modulation. The modulated
carrier is radiated into free space as an electromagnetic wave which is picked
up by the base station tower. Similarly, the signals transmitted by the tower
are received by the handset, thus establishing a two way communication.
This is one of the typical examples of a wireless communication system
which uses free space as a medium to transfer information from the trans-
mitter to the receiver. A key component of a wireless link is the antenna
which eciently couples electromagnetic energy from the transmitter to
free space and from free space to the receiver. An antenna is generally a
bidirectional device, i.e., the power through the antenna can ow in both the
directions, hence it works as a transmitting as well as a receiving antenna.
Transmission lines are used to transfer electromagnetic energy from one
point to another within a circuit and this mode of energy transfer is generally
known as guided wave propagation. An antenna acts as an interface between
the radiated electromagnetic waves and the guided waves. It can be thought
of as a mode transformer which transforms a guided-wave eld distribution
into a radiated-wave eld distribution. Since the wave impedances of the
guided and the radiated waves may be dierent, the antenna can also be
thought of as an impedance transformer. A proper design of this part is
necessary for the ecient coupling of the energy from the circuit to the free
space and vice versa.
One of the important properties of an antenna is its ability to transmit
power in a preferred direction. The angular distribution of the transmitted
2[IRUG8QLYHUVLW\3UHVV
2 Chapter 1 Electromagnetic Radiation
2[IRUG8QLYHUVLW\3UHVV
Introduction 3
2[IRUG8QLYHUVLW\3UHVV
4 Chapter 1 Electromagnetic Radiation
Visible light
Microwave bands
Millimeter wave
Band
designation
L S C X Ku K Ka
1 GHz2 GHz 2 GHz4 GHz 4 GHz8 GHz 8 GHz12.4 GHz 12.4 GHz18 GHz 18 GHz27 GHz 27 GHz40 GHz 40 GHz300 GHz
Fig. 1.3 Radio wave spectrum along with the band designations and typical ap-
plications.
the atmosphere is not a straight line. The air molecules also get ionized
by solar radiation and cosmic rays. The layer of ionized particles in the
atmosphere, known as the ionosphere, reects high frequency (3 MHz to
30 MHz) waves. A multi-hop communication link is established by repeated
reections of the electromagnetic waves between the ionosphere and the
surface of the earth. This is the mode of propagation of shortwave radio
signals over several thousand kilometres.
Both the radiation properties of the antennas and the propagation condi-
tions play a very important role in establishing a successful communication
link. This book addresses both these issues in some detail. It is assumed that
the students have some basic knowledge of electromagnetic theory. However,
in the following section some of the basic concepts of electromagnetic theory
used in the analysis of antennas are presented for easy reference as well as
for introducing the notation used in the book.
2[IRUG8QLYHUVLW\3UHVV
1.1 Review of Electromagnetic Theory 5
EXAMPLE 1.1
2[IRUG8QLYHUVLW\3UHVV
6 Chapter 1 Electromagnetic Radiation
I = 2.24ej1.1071 A
EXAMPLE 1.2
i(t) = Re{Iejt }
The eld vectors that vary with space, and are sinusoidal functions of
time, can also be represented by phasors. For example, an electric eld vector
E(x, y, z, t), a function of space (x, y, z) having a sinusoidal variation with
time, can be written as
E(x, y, z, t) = Re E(x, y, z)ejt (1.5)
2[IRUG8QLYHUVLW\3UHVV
1.1 Review of Electromagnetic Theory 7
Using phasor notation, Maxwells equations can be written for the elds
and sources that are sinusoidally varying with time as1
E = jH (1.6)
H = jE + J (1.7)
D= (1.8)
B=0 (1.9)
The rst two curl equations are the mathematical representations of Fara-
days and Amperes laws, respectively. The divergence equation [Eqn (1.8)]
represents Gausss law. Since magnetic monopoles do not exist in nature,
we have zero divergence for B [Eqn (1.9)].
The current density, J, consists of two components. One is due to the
impressed or actual sources and the other is the current induced due to the
applied electric eld, which is equal to E, where is the conductivity of
the medium (unit: siemens per metre, S/m). In antenna problems, we are
mostly working with elds radiated into free space with = 0. Therefore, in
the analyses that follow, unless explicitly specied, J represents impressed-
source current density.
In an isotropic and homogeneous medium, the electric ux density, D,
and the electric eld intensity, E, are related by
D = E (1.10)
where is the electric permittivity (unit: farad per metre, F/m) of the
medium. 0 is the permittivity of free space (0 = 8.854 1012 F/m) and
the ratio, /0 = r is known as the relative permittivity of the medium. It is
1See Cheng 2002, Hayt et al. 2001, Jordan et al. 2004, and Ramo et al. 2004.
2[IRUG8QLYHUVLW\3UHVV
8 Chapter 1 Electromagnetic Radiation
2[IRUG8QLYHUVLW\3UHVV
1.1 Review of Electromagnetic Theory 9
H = B = A (1.12)
or
1
H= A (1.13)
Substituting this into the equation E = jH, Maxwells rst equa-
tion is also incorporated
E = j(H) = j( A) (1.14)
or
(E + jA) = 0 (1.15)
(E + jA) = V (1.16)
E = (V + jA) (1.17)
2[IRUG8QLYHUVLW\3UHVV
10 Chapter 1 Electromagnetic Radiation
Equations (1.13) and (1.17) relate the H and E elds to the potential
functions A and V . Now, to satisfy Maxwells second equation, H =
jE + J, substitute the expression for the E and H elds in terms of the
potential functions [Eqns (1.13) and (1.17)]
1
( A) = j(V + jA) + J (1.18)
which is valid for a homogeneous medium. Expanding the left hand side
using the vector identity
A = ( A) 2 A (1.19)
we have
2 A + 2 A = J + ( A + jV ) (1.20)
So far we have satised three of Maxwells four equations. Note that only
the curl of A is dened so far. Since the curl and divergence are two
independent parts of any vector eld, we can now dene the divergence
of A. We dene A so as to relate A and V as well as simplify Eqn (1.20)
by eliminating the second term on the right hand side of the equation. We
relate A and V by the equation
A = jV (1.21)
This relationship is known as the Lorentz condition. With this the magnetic
vector potential, A, satises the vector wave equation
2 A + k 2 A = J (1.22)
where
k = (1.23)
is the propagation constant (unit: radian per metre, rad/m) in the medium.
Now, to satisfy Maxwells fourth equation, D = , we substitute
E = (V + jA) in this equation to get
(V jA) = (1.24)
or
2 V + j( A) = (1.25)
2[IRUG8QLYHUVLW\3UHVV
1.1 Review of Electromagnetic Theory 11
Eliminating A from this equation using the Lorentz condition [Eqn (1.21)]
2 V + k 2 V = (1.26)
Thus, both A and V must satisfy the wave equation, the source function
being the current density for the magnetic vector potential, A, and the
charge density for the electric scalar potential function V .
The right hand side of this equation is zero everywhere except at the source
itself. Therefore, in the source-free region, V satises the homogeneous wave
equation
1 V
2
r2 + k2 V = 0 (1.29)
r r r
ejkr
V (r) = V 0 (1.30)
r
travelling in the +r-direction. V 0 is the complex amplitude of the scalar
1(x, y, z):
rectangular co-ordinates; (r, , ): spherical co-ordinates.
2 SeeAppendix E for details on the coordinate systems and vector operations in dierent coordi-
nate systems.
2[IRUG8QLYHUVLW\3UHVV
12 Chapter 1 Electromagnetic Radiation
EXAMPLE 1.3
Show that
ejkr
V (r) = V 0
r
are solutions of
1 V
2
r2 + k2 V = 0
r r r
+ ejkr
V (r) = V 0
r
Substituting into the left hand side (LHS) of the given equation
1 ejkr
+ + ejkr
LHS = 2 r2 V0 + k2 V 0
r r r r r
1
+ ejkr ejkr + ejkr
=V0 2 r2 2 jk + k2 V 0
r r r r r
+ 1
jkr jkr
2 +e
jkr
=V0 e jkre + k V 0
r2 r r
+ 1 jkr jkr 2 jkr
2 +e
jkr
=V0 jke jke k re + k V 0
r2 r
=0
ejkr
V (r) = V 0
r
2[IRUG8QLYHUVLW\3UHVV
1.1 Review of Electromagnetic Theory 13
Substituting into the left hand side of the given dierential equation
1 ejkr ejkr
LHS = 2 r2 V0 + k2 V 0
r r r r r
jkr
1 1 e
r2 2 ejkr + jkrejkr + k 2 V 0
=V0 2
r r r r
1 jkr jkr 2 jkr
2 e
jkr
=V0 jke + jke rk e + k V 0
r2 r
=0
The instantaneous value of the scalar potential V(r, t) for the wave prop-
agating in the +r-direction can be written as
ej(tkr)
+
V(r, t) = Re V 0 (1.31)
r
+
v is the phase angle of V 0 . The equation for the constant phase spherical
wave front is
v + t kr = const (1.32)
The velocity of the wave is the rate at which the constant phase front moves
with time. Dierentiating the expression for the constant phase front surface
with respect to time, we get
dr
j jk =0 (1.33)
dt
+
This follows from the fact that V 0 and, hence, the phase v , is independent
of time, i.e., dv /dt = 0. Therefore, the velocity (v, unit: metre per second,
2[IRUG8QLYHUVLW\3UHVV
14 Chapter 1 Electromagnetic Radiation
dr
v= = (1.34)
dt k
Substituting the value of the propagation constant from Eqn (1.23), the
wave velocity is
1
v= = (1.35)
The velocity of the wave in free space is equal to 3 108 m/s. The distance
between two points that are separated in phase by 2 radians is known as
the wavelength (, unit: metre, m) of the wave. Consider two points r1 and
r2 on the wave with corresponding phases
1 = v + t kr1
2 = v + t kr2
such that
2 1 = k(r1 r2 ) = k = 2 (1.36)
2
k= (1.37)
The velocity can be written in terms of the frequency and the wavelength
of the wave
2f
v= = = f (1.38)
k 2/
EXAMPLE 1.4
50
E(x, y, z, t) = a cos(4 106 t 0.063r) V/m
r
2[IRUG8QLYHUVLW\3UHVV
1.1 Review of Electromagnetic Theory 15
50 j0.063r
E = a e V/m
r
Substituting this in Maxwells equation, Eqn (1.6), and expressing the curl
in spherical coordinates
a ra r sin a
r
1
jH = E = 2 /r / /
r sin
0 rE 0
50
E = a (j0.063)ej0.063
r
Therefore, the magnetic eld is given by
1 0.063 50 0.063r
H= E = a e
j r
0.2 j0.063r
H = a e A/m
r
2[IRUG8QLYHUVLW\3UHVV
16 Chapter 1 Electromagnetic Radiation
0.2
H = a cos(4 106 t 0.063r) A/m
r
Source
q (x', y', z' )
R
P (x, y, z)
Field
r'
point
r
y
o
x
Fig. 1.4 Position vectors of source and eld points
2[IRUG8QLYHUVLW\3UHVV
1.1 Review of Electromagnetic Theory 17
v
V (r, , ) = (1.40)
4R
(x , y , z )v ejkR
V (r, , ) = (1.41)
4 R
where V is the volume over which (x , y , z ) exists, or the source volume.
The instantaneous value of the scalar potential V(r, , , t) is obtained by
1 ejkR+jt
V(r, , , t) = Re V (r, , )ejt = Re (x , y , z ) dv
4 R
V
(1.43)
Using the relation v = /k, this reduces to
1 ej(t R
v
)
V(r, , , t) = Re (x , y , z ) dv (1.44)
4 R
V
It is clear from this expression that the potential at time t is due to the
charge that existed at an earlier time R/v. Or the eect of any change
in the source has travelled with a velocity v to the observation point at a
distance R from the source. Therefore, V is also known as the retarded scalar
potential.
2[IRUG8QLYHUVLW\3UHVV
18 Chapter 1 Electromagnetic Radiation
In Section 1.1.1 it is shown that both, electric scalar potential, V and mag-
netic vector potential, A, satisfy the wave equation with the source terms
being / and J, respectively. Therefore, a similar heuristic argument can
be used to derive the relationship between the current density distribution
J(x , y , z ) and the vector potential A(r, , ). Starting from the expression
for the magnetic vector potential for a static current density we introduce
the delay time R/v to obtain the retarded vector potential expression
for the time-varying current density distribution J. The vector potential at
any time t is related to the current density distribution at time (t R/v).
Further, the vector A has the same direction as the current density J.
The relationship between the current density J(x , y , z ) and the vector
potential A(r, , ) is given by simply multiplying the static relationship
with the ejkR term. Thus, the retarded vector potential is given by
ejkR
A(r, , ) = J(x , y , z ) dv (1.45)
4 R
V
ejkR
A(r, , ) = Js (x , y , z ) ds (1.46)
4 R
S
For a line current I (in A), the integral reduces to a line integral
ejkR
A(r, , ) = I(x , y , z ) dl (1.47)
4 R
C
1
H= A (1.48)
2[IRUG8QLYHUVLW\3UHVV
1.2 Hertzian Dipole 19
1
E= H (1.49)
j
x = y = z = 0 (1.51)
R= (x x )2 + (y y )2 + (z z )2 = x2 + y 2 + z 2 = r (1.52)
2[IRUG8QLYHUVLW\3UHVV
20 Chapter 1 Electromagnetic Radiation
az Az
z
a r Ar
a A
P (x, y, z)
I0 dl
y
rs
in
X
Fig. 1.5 Components of the vector potential on the surface of a
sphere of radius r , due to a z -directed current element kept at the
origin
ejkr
A(r, , ) = az I0 dl = az Az (1.53)
4 r
Note that the vector potential has the same vector direction as the current
element. In this case, the az -directed current element produces only the
Az -component of the vector potential.
The H and E elds of a Hertzian dipole are computed using the rela-
tionships given by Eqns (1.48) and (1.49), respectively. The E and H elds
are generally computed in spherical coordinates for the following reasons
(a) the (ejkr/r) term indicates that the elds consist of outgoing spheri-
cal waves which are simple to represent mathematically in spherical coor-
dinates and (b) the spherical coordinate system allows easy visualization
of the behaviour of the elds as a function of direction and simplies the
mathematical representation of the radiated elds.
2[IRUG8QLYHUVLW\3UHVV
1.2 Hertzian Dipole 21
From Fig. 1.5 we can relate the z-component of the vector potential, Az ,
to the components in spherical coordinates Ar , A , and A as
Ar = Az cos
A = Az sin (1.54)
A = 0
Substituting the components of A from Eqn (1.54) into Eqn (1.55), we get
ar ra r sin a
1
A= 2 /r / / (1.56)
r sin
Az cos rAz sin 0
1
A = a Az sin (jkr + 1) (1.58)
r
Substituting the result in Eqn (1.48) and simplifying, we get the expres-
sions for the components of the H eld of a Hertzian dipole in spherical
coordinates as
Hr = 0 (1.59)
H = 0 (1.60)
I0 dl sin ejkr 1
H = jk 1+ (1.61)
4 r jkr
2[IRUG8QLYHUVLW\3UHVV
22 Chapter 1 Electromagnetic Radiation
The electric eld can be obtained from Maxwells curl equation. Substituting
the expression for H in Eqn (1.49) we get
a ra r sin a
r
1 1 1
E= H= /r / / (1.62)
j
j r2 sin
0 0 r sin H
After performing the indicated derivative operations, Eqn (1.63) can be sim-
plied to give the electric eld components of a Hertzian dipole in spherical
coordinates as
I0 dl cos ejkr 1
Er = 1+ (1.64)
2r r jkr
kI0 dl sin ejkr 1 1
E = j 1+ (1.65)
4 r jkr (kr)2
E = 0 (1.66)
EXAMPLE 1.5
k 2/ 1
= =
(2f ) (f )
The velocity of the wave, v, is related to the frequency, f , and the wave-
length, , by v = f , and v is related to the permittivity, , and permeability,
, of the medium by v = 1/ . Substituting these in the above equation
and simplifying
k 1
= = =
v
2[IRUG8QLYHUVLW\3UHVV
1.2 Hertzian Dipole 23
The impedance of the medium is related to and by = / and,
therefore
k
=
6
10
4
10
2
10
Amplitude
0
10
1/(kr)
2 1/(kr)2
10
1/(kr)3
4
10
6
10 2 1 0 1 2
10 10 10 10 10
kr
Fig. 1.6 Dependence of 1/(kr), 1/(kr)2 , and 1/(kr)3 on (kr)
2[IRUG8QLYHUVLW\3UHVV
24 Chapter 1 Electromagnetic Radiation
distances from the Hertzian dipole, only the terms containing 1/r are re-
tained in the electric and magnetic eld expressions. The electric and the
magnetic elds of a z-directed Hertzian dipole for r are given by
EXAMPLE 1.6
Calculate and compare the r and components of the electric eld inten-
sities at x = 100 m, y = 100 m, and z = 100 m produced by a Hertzian
dipole of length dl = 1 m kept at the origin, oriented along the z-axis,
excited by a current of i(t) = 1 cos(10 106 t) A, and radiating into free
space.
Solution: The frequency of excitation is = 10 106 rad/s, and therefore
f = 5 106 Hz. The dipole is radiating in free space with parameters =
4 107 H/m and = 8.854 1012 F/m. Therefore, the impedance of
the medium is
4 107
= = = 376.73
8.854 1012
The distance r between the eld point and the dipole (which is at the
origin) is
r= x2 + y 2 + z 2 = 1002 + 1002 + 1002 = 173.2 m
2[IRUG8QLYHUVLW\3UHVV
1.2 Hertzian Dipole 25
EXAMPLE 1.7
x = r sin cos
y = r sin sin
z = r cos
2[IRUG8QLYHUVLW\3UHVV
26 Chapter 1 Electromagnetic Radiation
At any point P (r, , ), ar , a , and a denote the unit vectors along the
r, , and directions, respectively. The unit vector along the r-direction is
given by
r
ar =
|r|
where |r| = (r sin cos )2 + (r sin sin )2 + (r cos )2 = r and hence
r/
a = = ax sin + ay cos
|r/|
This transformation from rectangular to spherical coordinates can be ex-
pressed in a matrix form as
ar sin cos sin sin cos ax
a = cos cos cos sin sin ay
a sin cos 0 az
A = ax + ay + az = ar Ar + a A + a A
2[IRUG8QLYHUVLW\3UHVV
1.2 Hertzian Dipole 27
2[IRUG8QLYHUVLW\3UHVV
28 Chapter 1 Electromagnetic Radiation
EXAMPLE 1.8
Derive the expressions for the elds of a current element I0 dl kept at the
origin, oriented along the x-axis, and radiating into free space.
Solution: Since the current element is oriented along the x-direction, the
magnetic vector potential has only the x-component. Following the proce-
dure given in Section 1.2, we can write the magnetic vector potential as
0 ejkr
A = ax I0 dl
4 r
The unit vector ax can be written in terms of the unit vectors of the spherical
coordinates (see Example 1.7)
0 ejkr
A= I0 dl (ar sin cos + a cos cos a sin )
4 r
The magnetic eld is given by
a ra r sin a
r
1 1
H= A= 2 /r / /
r sin
Ar rA r sin A
2[IRUG8QLYHUVLW\3UHVV
1.2 Hertzian Dipole 29
So far we have learnt how to compute the elds due to a current distri-
bution using the vector potential approach. Every antenna can be looked
at as a current distribution producing electric and magnetic elds in the
surrounding space and, therefore, we have learnt the basics of computing the
elds of an antenna. In the following chapter we will learn about properties
of antennas and introduce the various terms associated with antennas.
2[IRUG8QLYHUVLW\3UHVV
30 Chapter 1 Electromagnetic Radiation
Exercises
1.1 Prove that the spherical coordinate system 1.12 Show that at large distances from a radi-
is orthogonal. ating Hertzian dipole (r ), the elec-
1.2 Show that for any twice dierentiable tric and magnetic elds satisfy Maxwells
scalar function, , = 0. equations.
1.3 Show that for any twice dierentiable vec- 1.13 A z -directed Hertzian dipole placed at
tor function A, A = 0. the origin has length dl = 1 m and is ex-
cited by a sinusoidal current of amplitude
1.4 Prove the vector identity A =
( A) 2 A. I0 = 10 A and frequency 1 MHz. If the
dipole is radiating into free space, calcu-
1.5 In a source-free region show that the E late the distance in the x-y plane from
and H elds satisfy 2 E + k 2 E = 0 and the antenna beyond which the magnitude
2 H + k 2 H = 0, respectively. of the electric eld strength is less than
1.6 Show that V (r) = V0 ejkr/r , where V0 is 1 103 V/m.
a complex constant and k is a real number, Answer: 6279 m
represents a wave travelling in the positive 1.14 Derive an expression for the elds of a
r-direction. Hertzian dipole of length dl carrying a
1.7 Plot the equiphase surfaces of the electric current of I0 which is located at the ori-
eld of an EM wave given by (a) E = gin of the coordinate system and oriented
a E0 ejkr/r and (b) E = ay E0 ejkx , along the y -axis.
where E0 is a complex constant. 1.15 Find the strength of the z -component of
1.8 Derive Eqns (1.59)(1.61) from Eqn (1.57). the electric eld at (0, 100 m, 0) produced
1.9 Derive Eqns (1.64)(1.66) from Eqn (1.63). by a z -directed Hertzian dipole of length
dl = 0.5 m, placed at the origin, carrying
1.10 For the Hertzian dipole considered in Sec-
a current of i(t) = 2 cos(6 106 t) A,
tion 1.2, compute the electric and magnetic
and radiating into free space. If the dipole
elds in the rectangular coordinate system
is oriented along the x-axis, what will be
directly by taking the curl of az Az . Now
strength of the x-component of the elec-
convert the elds into spherical coordinates
tric eld at the same point?
and compare them with the results given in
Answer: Ez = 0.01856 99.3 V/m;
Eqns (1.59)(1.61) and Eqns (1.64)(1.66).
Ex = 0.01856 99.3 V/m
1.11 Show that = k , where the symbols
have their usual meaning.
2[IRUG8QLYHUVLW\3UHVV