Lecture 2: Wave Equations, Green’s Function, Huygens’ Principle

and Fresnel Zones

2.1. Wave Equation Presentation

Physically, EM-wave propagation phenomena can be described by use of both the
scalar and vector wave equation presentations. In the case of an isotropic
inhomogeneous medium one can present Maxwell’s equations in the following form,
using their time harmonic presentations (1.7):
   1 ( r )   E( r )   2  (r )E( r )  ij(r )     1 (r )M (r )
(2.1)
1 1
   (r )  H(r )    (r ) H( r )  iM ( r )     ( r ) j(r )
2

Because most problems of wave propagation above the terrain, including built-up
environments, reduce to propagation in a homogeneous, source-free isotropic medium,
this system can be easily simplified from system (1.31) by taking into account the
relations (1.17a)-(1.17c) with  (r )   ,  (r )   , (r)   , that is,
    E(r )   2 E(r )  0
(2.2)
    H(r )   2 H(r )  0

Because both equations are symmetric, one can use one of them, namely that for E ,
and by introducing the vector relation     E  (  E)   2 E and taking into
account that   E  0 , finally obtain
 2 E( r )  k 2 E( r )  0 ,

(2.3)
where k 2   2  . It can be shown that all other electromagnetic vectors satisfy as
well the same wave equation as (2.1) (see, particularly, the equations (1.34) for
Hertzian vectors in a source-free medium).
In special cases of a homogeneous, source-free, isotropic medium, the three
dimensional wave equation reduces to a set of scalar wave equation. This is because in
Cartesian coordinates, E(r )  E x x 0  E y y 0  E z z 0 , where x 0 , y 0 , z 0 are unit
vectors in the directions of the x, y, z coordinates, respectively. Hence, the equation
(2.3) consists of three scalar equations such as
 2  ( r )  k 2  (r )  0 , (2.4)

1
where (r ) can be either E x , E y , or E z . This statement is not true in cylindrical or
spherical coordinate systems. The problems of independent solution of each scalar
wave equation, such as (2.4) is the subject next lecture.

2.2. Boundary Conditions

Equations (2.4) describe all the propagation phenomena within an infinite

inhomogeneous isotropic medium. But if we consider two inhomogeneous finite or
semi-finite regions, we need to introduce boundary conditions at the interface between
these two regions in order to solve one of the two equations (2.2). In this case the
procedure to solve the vector wave equation is as follows.
As a first step this equation is solved separately for each region. Then, in the
second step, by patching the solution together via boundary conditions, we obtain the
solution for two neighboring regions. It can be easily shown that the boundary
conditions follow from one of the two vector wave equations (2.2). To do so, we
integrate the first equation of (2.2) within a small region in the interface of the two
inhomogeneous semi-finite or finite regions, as presented in Fig. 2.1.

Fig. 2.1

Then using Stokes’s theorem for the surface integral of a curl and using the same
integration over surface S for both equations (2.2), we finally obtain after
straightforward derivations and taking the limit   0 (see Fig. 2.1) respectively for
the magnetic-field component
n  H1  n  H 2  j S (2.5)
and for electric-field component

2
 n  E1  n  E 2  M S (2.6)
where M S and j S is a magnetic and electric current sheet at the interface,
respectively. Equation (2.5) states that the discontinuity in the tangential component of
the magnetic field is proportional to the electric current sheet j S . This is the first
boundary condition for solving any one vector electromagnetic equation from (2.6).
Equation (2.6) states that the discontinuity in the tangential component of the electric
field is proportional to the magnetic current sheet M S . This is the second boundary
condition for (2.1).
Both boundary conditions (2.5) and (2.6) can be simplified for the case of wave
propagation above a flat intersection. In the case considered, the first boundary
condition (2.5) for an isotropic non-magnetized (  =1) source-free ( M S  0 , j S  0 )
sub-soil medium reduces to
H1n  H 2 n
and
H1  H 2 

(2.7)
Both conditions are valid in the case of finite conductivity of each medium, which is
satisfied within the air-ground surface. The first condition in (2.7) states, that the
normal components of the magnetic field of the wave are continuous at the interface of
air-conductor intersection surface. The second condition in (2.7) states that the
tangential component of magnetic field is also continuous at the interface of the air-
ground surface.
As for the second boundary condition (2.6), it also can be simplified for the
interface of the air-conductor intersection as
n  E1  n  E 2
or (2.8)
E 1  E 2 

Condition (2.8) states that the tangential components of the electric field of an EM-
wave are continuous at the interface of the air-conductor intersection surface. One may
notice that conditions (2.5) and (2.6) are more general than those described by (2.7)
and (2.8), and satisfy various kinds of isotropic inhomogeneous media that consist of
both electric and magnetic sources.
Example 1: Two homogeneous linear isotropic magnetic materials have an
interface at x=0, as shown in Fig 1E.

3
 Fig. 1E.

At the interface, there is a surface current with density j s  20u y  A / m  . The relative

permeability,  r1  2 , and the magnetic field strength, H1  15u x  10u y  25u z , are
given in the region x<0, and  r 2  5 in the region x>0.
Find: The magnetic field strength H 2 in the region x>0.
Solution
1) The magnetic field can be separated according to (2.7) on two components
with independent boundary conditions. Thus, for the normal component, we
write
 r1 2
 r1  0 H 1n   r 2  0 H 2n  H 2 x  H 1x   H 1x  6 A / m
 r2 5
2) For the tangential component, we write
H 2  H 1  j s  H 2 z  H 1z  j y  25  20  45 A / m

3) Since there is no change in the y-component of the magnetic field, we write

H 2  6u x  10u y  45u z A / m

2.3. Wave Equation Solutions

To analysis different kinds of wave equation solutions, let us consider free space, as a
simple example of an infinite isotropic homogeneous source-free medium, the vector
wave equation of which can be presented in very simple form (2.2)-(2.3) for one of the
component of the EM field, or by use of the Hertzian vector (r ) :

4
  2 (r )  k 2 (r )  0 (2.9)
where, once more, the wavenumber k   (in many practical cases of terrain
propagation  =1 with great accuracy, and one can rewrite it as k  ).

2.3.1. Plane Waves in Homogeneous Medium

For plane waves in a Cartesian coordinate system each of the equation (2.3) or (2.9)
can be rewritten in scalar form (2.4) for any Cartesian component of vectors E( r ) ,
H( r ) or ( r ) . Usually, in the literature there is other form of presentation of
equation (2.4) by introducing instead wavenumber k, the phase velocity

 c
v ph  v   . In this case (2.4) can be rewritten as:
k 

2
 2 (r )  (r )  0 (2.10)
v2
Wave equations (2.4) or (2.10) have the well-known solution [1-3]
( r )  exp ik  r  (2.11)
The waves that satisfy scalar equation (2.10) and are determined by expression (2.11)
are called plane waves. Wave vector k denotes here the direction of propagation of
the plane wave in free space (see Fig. 1.1, Lecture 1). If one considers the plane wave
that propagates in any direction, say along the x-axis, then the fundamental solution of
(2.10) is
( x)  A exp ikx  B exp ikx (2.12)
This solution describes the waves, propagating in the positive direction (with the sign
“+” in the exponent) and in the negative direction (with the sign “-” in the exponent)

c
respectively along the x-axis with phase velocity v ph  
which equals c in an

ideal free space.

But one can note that EM fields have a vector character and satisfy Maxwell’s
equations (1.7) or wave equations such as (2.3), (2.4). Thus one can find the field
vectors in the following form
E(r )  e E E 0 exp ik  r 
or (2.13)
H (r )  e H H 0 exp ik  r 

5
where e E and e H are the constant unit vectors, i.e., | e E | | e H |  1 ; E 0 and H 0
are the complex amplitudes, which are constant in space and time. From conditions in
free space without sources
E  0 and H  0
it follows that
eE k  0 and eH k  0 (2.14)
which denote that E and H are perpendicular to direction of wave propagation k .
Moreover, because in free space the first Maxwell equation (1.1a) or (1.7a) reduces to
  E  iB  0, B  H  H (2.15)
we can finally obtain from (2.43) for a plane wave (2.41) in free space (with   1 ):
i (k  e E ) E 0  ke H B0  exp ik  r   0 (2.16)
Equations (2.16) have solutions:
k  eE
eH  , B0   E0 (2.17)
k
Hence, vectors e E , e H and k form the system of orthogonal vectors, where vectors
E and B oscillate in phase and their ratio is constant (see Fig. 2.3). The wave, which
is described by relations (2.13) and (2.17), is a transverse wave propagating in the k -
direction.

2.3.2. Poynting Vector and Poynting Theorem

This theorem is the simple law of EM-wave energy conservation. It is known from
electrostatics and magnetostatics, that the work of the electric field to move a single
charge q is equal q  v  E , where v is the vector of the charge velocity. The same
work of the magnetic field for this charge is equal to zero, because the magnetic field
direction is perpendicular to the velocity vector [1-3]. For a continuous distribution of
charges and currents in a medium, the total work of the EM field in the volume V in
unit time is equal [1-3]:

 j  Edv
V

(2.18a)
This expression determines the velocity of the decrease in the field energy within the
volume V.
Let us now obtain the law of energy conservation, using Maxwell’s equations
(1.1). We shall substitute the current density j in (2.10a) using (1.1b)

6
  D 
 j  Edv   E    H  E  t dv
V V

(2.18b)
Taking into account the vector equality [4]
   E  H  H    E  E    H

one can easily rewrite (2.18) as:

 D B 
 j  Edv      (E  H)  E 
V V
t
H
t 
dv

(2.19)
If we now present the density of total field energy according to [1-3] as
W   E  D  B  H ,
(2.20)
then (2.11) will be rewritten in the following form:
 W 
  j  Edv    t    ( E  H) dv
V V 

(2.21)
Because expression (2.21) is written for any volume V, it can be presented in
differential form:
W
   S = -j  E
t
(2.22) Equation (2.22) is the equation of EM-field energy conservation, or the
equation of continuity. It can be easily shown that the vector S  E  H in the brackets
on the right hand side of (2.22) has the dimension of watt/ m 2 , which is that of power
density. From (2.22) is clear that it may be associated with the direction of power flow.
The vector that determines the power flow of EM field is called the Poynting vector.
Equation (2.21) is the integral Poynting theorem and equation (2.22) is its vector
presentation.
Using the time harmonic presentation of Maxwell’s equations, one can convert
(2.13) to the time harmonic form. In fact, if we now introduce instead of the derivation

 1
t
the term i, and present the operation of averaging < E  H >, as Re E  H * ,
2
 
taking into account Gauss’s theorem for the term   (E  H * ) , we finally obtain from
(2.21) the Poynting theorem presented in time harmonic form:

7
  ds  ( E  H )  dv(H  B  E  D * )   dvE  j
* *

S V V

(2.23)
or accounting for (2.20), we can present the total energy of the electromagnetic field
as
E 2  H 2
Wtotal  W E  W H   (2.24)
2 2
Then, the Eq. (2.22) can be rewritten as

  E 2 H 2 
S  E  H  ds   t V 2  2 dv  V E dv
2

                 
total power flow storage e / m energy heating
leaving the region

(2.25)

Example 2: The field vectors in free space are given by the following

4 uz E
expressions: E  10 cos(t  z )u x (V/m) and H  Z0
( A / m) . The frequency
3

f=500 MHz, Z 0  120 ()

Find: Poynting vector.
Solution
1) In phasor notations, the fields are expressed as
4  4 
4  j   z j   z
j   z 3  3 
E  10e 3 
ux and H  10 e uy 
e
uy
120 12
2) The Poynting vector is
1 10
S av (r )  Re E(r )  H * (r )  u z  0.133u z (W / m 2 )
2 2  12
Example 3: In Fig. 2E a capacitor of radius a and of width b, which is
connected with the source of voltage V(t) via the corresponding electrical circuit, is
given.

8
 Fig. 2E.

Find: the power flowing via the area of the radial edge of a capacitor using Poynting
theorem, if the density of the total magnetic field equals zero. Show that this power
equals the rate of total energy changes.

Solution
According to conditions of the problem, the total energy equals only the energy of the
electric field, i.e.,
 E 2 
Wtotal  WE  WH  W E  
2
 E 2 
 V  


 2   a b
2
 
   

H 2
According to Poynting's theorem (2.25), since wH   0 and heating is not
2
accounting, we finally get:
d  E 2  
  a 2 b    a 2 b  E
dWtotal dWE dE
    (E1)
dt dt dt  2   dt

On the other hand, according to 2nd Maxwell equation (see Lecture 1):

 B   D   E 
   0 dl    j  ds 
t     ds
t 
L S S

or

2aH   a 2  dE
dt
 H 
a dE
2 dt

Then, the Poynting's vector power density equals:

a dE dE
PS   E  H   S   E   (2ab)   (a 2 b) E (E2)
2 dt dt
Comparing (E1) and (E2), we obtain well-known result:

9
 dWtotal
PS  
dt
(E3)

2.3.3. Wave Polarization

The vector of electric field in the plane wave as described by formula (2.13) is directed
along unit vector e E . To obtain the more general case of wave polarization we need
an additional linear polarized wave independent of the first one. It can be easily shown
that two linear independent solutions, which satisfy wave equation (2.12) or (2.13) can
be presented in the following form
E 1  e 1 E1 exp ik  r 
(2.26)
E 2  e 2 E 2 exp ik  r 

The magnetic field components of the EM wave satisfy, according to (2.13) in free
space (  =1), the following relations:
k Ej
Bj   , Bj Hj, j=1, 2 (2.27)
k
Here amplitudes E1 ( H1 ) and E 2 ( H 2 ) are the complex values, which enable us to
introduce the phase difference between the two components of the EM wave. Thus the
common solution for the plane EM-wave propagated along vector k , can be
presented as a linear combination of E1 and E 2 :
E(r )   e 1 E1  e 2 E 2  exp ik  r  (2.28)
If E1 and E 2 have the same phase, then solution (2.20) describes the linear
polarized wave with polarization vector directed to the e 1 -axis at angle
E 
  tan 1  2  (2.29a)
 E1 

and with amplitude

 
1/ 2
E  E12  E 22 (2.29b)
as presented in Fig. 2.2.
If E1 and E 2 have different phases, then the EM wave (2.28) is elliptically


polarized. If E1  E 2 and phase difference equal , then the elliptically polarized
2
wave becomes a circularly polarized wave.

10
 Fig. 2.2

In this case solution (2.28) can be rewritten as

E( r )  E 0  e 1  e 2  exp ik  r  (2.30)
The sign “+” corresponds to anticlockwise rotation (sometimes called the wave with
left-hand circular polarization). The sign “-” corresponds the waves with right-hand
circular polarization (see Fig. 2.3).

Fig. 2.3

Then two waves with circular polarization can be considered as a basic system for
describing the common case of polarized waves. Let us introduce in the common case
the orthogonal complex unit vectors:
1
e   e1  ie 2  (2.31)
2
Then, the common presentation of a polarized wave (2.28) by the use of linearly
polarized waves, now, by the use of two circularly polarized waves (2.30) and
expression (2.31), can be rewritten as
E(r )   e  E   e  E   exp ik  r  (2.32)

11
where E  and E  are the complex amplitudes of two circularly polarized waves with
opposide directions of rotation. If their modules are different, but their phases are
equal, then expression (2.32) describes, as above, an elliptically polarized wave with
main elliptical axes directed along e 1 and e 2 . The ratio of these semi-axes equals of

E
(1-q)/(1+q), where q  . If the complex amplitudes have different phases, so that
E
E
 q  exp i (2.33)
E

then the ellipses’ axes for E - vector are rotated by angle . In Fig. 2.4 the
2
common case of an elliptical polarized EM wave is presented.

Fig. 2.4

At each spatial point the vector E (the same applies to the vectors H or B )
describes ellipses, as shown in Fig. 2.4. For the case q  1 , we once more return to
the case of linearly polarized wave.

2.3.4. Cylindrical and Spherical Waves in Free Space

As mentioned in the literature [1-3, 5-7], to obtain the common vector presentation of
cylindrical and spherical waves is a very complicated problem, which can be reduced
for the case of an isotropic homogeneous source-free medium (with properties which
limit to those in free space, but with   1 ) to the simple scalar form, as was done
above for the plane wave in the Cartesian coordinate system.
The scalar wave equation in the cylindrical coordinate system  , , z can be written
as:

12
 1   1  2 
   2   k 2   (r)  0 (2.34)
      2
z 2


The above partial differential equation can be solved by separation of variables, and its
one-dimensional (along the z-axis) solution can be presented in the following form [4]
 (r )  n () exp in  ik z z (2.35)
where n is an integer since the wave field has to be 2 periodic in  . Then, by
substituting (2.35) in (2.34), we reduce it to an ordinary differential equation with the
full derivative (d/dx, x is variable) presentation [4]:
 1 d d n2 
   2  k r2  n ()  0 (2.36)
  d d  

where k r   k 2  k z2 
1/ 2
. One may notice that (2.28) is the Bessel equation [10] with
two lineary independent solutions. Its general solution is a linear superposition of the
following four spectral functions [4]: Bessel function, J n ( k r ); Neumann function,

N n ( k r ); Hankel functions of first order, H n(1) ( k r ) , and second order,

H n( 2 ) ( k r ), respectively. Since only two of these four special functions are
independent, they are lineary related each other, that is,

J n (k   ) 
1
2
 
H n(1) ( k   )  H n( 2) ( k   ) ,

N n (k   ) 
1
2i

H n(1) (k   )  H n( 2) (k   ) , 

(2.37a)
H n(1) (k   )  J n ( k   )  iN n ( k   ),
or
H n( 2) ( k   )  J n ( k   )  iN n ( k   ).

(2.37b)
Full information about the properties of these special functions can be obtained from
the special books [4, 10], which describe all mathematical functions. For our purposes
the exponential approximation of those functions is very important.
For example, the general representation of a cylindrical wave (2.27) can be
reduced in a very simple exponential form:

13
 2  n   n 
 (r ) ~ exp i  i  exp i  ()  ik z z  ik r   (2.38)
k r   2 4   

Here  is the arc length in the  direction, and n /  can be thought of as the
component of vector k , if one compares the cylindrical wave presentation (2.38) with
that for a plane wave (2.11). Consequently, expression (2.30) looks like a plane wave
propagating mainly in the direction k  k z zˆ  k r ̂ , when    .
We now consider a spherical wave presentation in free space. In the spherical
coordinate system  r, ,  , the scalar wave equation is [4]:
 1  2  1   1 2 
 2 r  2 sin   2  k 2   (r )  0 (2.39)
 r r r r sin    r sin  
2 2


Following [3], we present the solution of this equation in the form:

 (r )   (r , ) exp im (2.40)
The general equation (2.31) can be further simplified by the separation of variables by
letting
(r , )  bn ( kr ) Pnm (cos )
(2.41)
where Pnm (cos ) is the associate Legendre polynomial satisfying the equation

 1 d d  m2  m
 sin   n ( n  1)    Pn (cos )  0 (2.42)
 sin  d d  sin 2
 

Analogously, bn ( kr ) satisfies the equation

1 d 2 d n(n  1) 
 2 r  k2   bn ( kr )  0 (2.43)
 r dr dr r2 
Equation (2.43) is just the spherical Bessel equation, and b n ( kr ) is either the

spherical Bessel function, jn ( kr ) , spherical Neumann function, n n ( kr ) , or the

spherical Hankel functions, h (1) ( 2)

n ( kr ) and h n ( kr ) [4, 10].

As is shown in [4, 10], the spherical special functions can be approximated by

exp ikr
the spherical functions proportional to ~ . If so, one can represent spherical
r
wave, as a plane one when    .
In an unbounded homogeneous medium using any kind of source, it is a
difficult problem to obtain a strict solution of the wave equation, which describes EM

14
wave propagation in such a medium. Usually, to obtain a particular solution of a wave
equation one can assume that the source is a point with respect to volume metric
dimensions around this source. In the literature [1-3, 6-9] to determine the criterion of
a point source requires that the linear dimensions of the source, l must be smaller


than the wavelength in the considered medium, that is, l << . In this case the so-

called Green’s function, as a solution of the wave equation for a point source, can be
introduced. Moreover, if any real antenna can be represented as a general real source
by a linear superposition of point sources, one can obtain a general solution for the
wave equation with such a source by using the solution of the wave equation for a
point source, in the other words, by use of Green’s function as a point source function.
This result is also connected with the topic of linearity of the wave equation in the
considered medium. Below, we will examine the boundary-value problems both for
scalar and vector waves by employing the Green’s function presentation.

2.4. Green’s Function Presentation

Scalar presentation of Green’s function. First of all we examine the Green’s function
presentation for the scalar wave equation. In this case let us construct the solution of a
scalar wave equation in any volume V of free space having an arbitrary source s( r ).
Such a solution can be written in the same form, as (2.10), but with a source in its
right-hand side:
 2 ( r)  k 2 ( r)  s( r) (2.44)
First, we will introduce the same equation for the Gereen’s function, but with a point
source in its right-hand side:
 2 G ( r, r ' )  k 2 G ( r , r' )    ( r  r' )

(2.45)
The given functions G( r, r' ) and ( r) can be easily found from the principle of
linear superposition, since G( r, r' ) , as was mentioned above, is the solution of (2.44)
with a point source in the right-hand side. In fact, one can notice that an arbitrary
source s( r ) is just
s( r )   dr ' s(r ' )(r  r ' )
(2.46)
which actually a linear superposition of point sources in mathematical terms.
Consequently, the solution of (2.45) is just

15
 (r )    dr ' G ( r , r ' ) s(r ' )
V

(2.47)
which is the linear superposition of the solution of (2.37).
To find the solution of (2.45) for free space, or more correctly, for an
unbounded, homogeneous medium, one can solve it in spherical coordinates with the
origin at vector r' . In this case (2.45) reduces to
 2 G ( r)  k 2 G ( r)    ( r) (2.48)
But due to the spherical symmetry of a point source, G( r) must also be spherically
symmetric. Then, for r  0, the homogeneous, spherically symmetric solution of
(2.48) is given by
exp ikr exp ikr
G ( r)  A B (2.49)
r r
Since sources are absent at infinity, a physically correct solution of (2.48) can be
presented as
exp ikr
G ( r)  A (2.49a)
r
To determine the unknown constant A, we substitute (2.41a) into (2.39) and integrate
it over a small volume about the origin to yield
A exp ikr A exp ikr
 dv       dv  k  1.
2
(2.50)
V
r V
r

Note that the second integral in (2.50) vanishes when V  0 , because dv=4r 2 
dr.
The first integral in (2.50) can be converted into a surface integral using
Gauss’s theorem to obtain
 2 d
exp ikr 
4 r A   1
 dr r  r 0
or
1
A
4
(2.51)
As was mentioned above, the solution of (2.45) must depend only on | r - r' | .
Therefore, Green’s function must be presented, as a solution of (2.45), in the
following form:

16
 exp ik | r - r' |
G ( r, r' )  G ( r - r ' ) 
| r - r' |

(2.52)
Moreover, it can be seen that G ( r, r' )  G ( r' , r) from reciprocity irrespective of the
shape of volume V [3]. This fact and formula (2.52) imply that Green’s function is
translationally invariant for unbounded, homogeneous media. Consequently, a general
solution of inhomogeneous equation (2.44) by using (2.47) can be finally presented as:
exp ik | r  r '|
(r )    dr ' s( r ' )
V
|r - r' |

(2.53)

2.5. Huygens’ Principle – Mathematical Description

Huygens’ principle, which is named by the Danish researcher Christian Huygens, shows

how a wave field on the surface S determines the wave field off the surface S (Fig.

2.5a) or, inversely, inside the area bounded by the surface S (Fig. 2.5b). In the other

words, each point on the surface S can be interpreted as a source of a spherical wave ,

which can be observed at any point A, either in the outside space with volume V, if a

source O is inside it (Fig. 2.5a), or inside the bounded area S, if a source O is outside

the surface S (Fig. 2.5b).

Above we indicated that any scalar wave in an unbounded source-free

homogeneous isotropic medium can be described by the homogeneous equation (2.10)

(or (2.36) without the right-hand side).

If we now multiply (2.10) by G(r , r ' ) and (2.45) by (r ) , subtracting the resulting

equations and integrating over a volume V containing vector r' , we have

 dr G(r , r ' ) 

(r )  (r ) 2 G (r , r ' )  (r ' )
2

(2.54)

Taking into account the following relations [4]

G 2    2 G     G  G 

17
we can rewrite (2.54) by use of the so-called Green’s theorem. Green’s theorem or, as

it is sometimes called, the second Green formula [2] states the equivalence of volume

integral (2.54) with the surface one, i.e.,

 dr G(r, r' )

2
  dsG(r, r' ) 
 (r )   (r ) 2 G (r , r ' ) 
n
 (r )
G ( r , r ' ) 
n 
V S

(2.55)

It is well known [4], that for any two scalar functions f and g there is some relation

g
between them: f g  n  f . Taking into account this relation and Green’s theorem,
n

based on Gauss’ divergence law, one can rewrite the left-hand side of (2.46) as:

 (r ' )   dsn   G(r, r' )(r)  (r)G(r, r' )

S

(2.56)

Equation (2.56) is the mathematical expression of the statement that once (r ) and

n  ( r) is known on surface S, then ( r' ) away from S can be found.

18
 Fig. 2.5

From general equation (2.56) one can obtain two different situations at the bounded

surface S. In fact, if G( r , r ' ) satisfies equation (2.10) or (2.44) (without its right-

hand side) with the boundary conditions n  G( r, r' )  0 for r S , then expression

(2.56) becomes

( r ' )   dsG (r , r ' )n  (r ).

S

(2.57a)

19
On the other hand, if now G( r, r' ) has only to satisfy equation (2.10) or (2.44)

(without its right-hand side) for both r and r' in volume V, and no boundary

condition has yet been imposed on G( r, r' ) , then, in the case of G( r, r' )  0 for

r S , (2.48) becomes

 (r ' )    ds (r )n  G (r , r ' ).

S

(2.57b)

Equations (2.56), (2.57a), and (2.57b) are various forms of Huygens’ principle

depending on the definition of Green’s function G( r, r' ) on the bounded surface S.

For example, equations (2.57a) and (2.57b) state that only n  ( r) or ( r) need

be known, respectively, on the surface S in order to determine wave function (r' ) .

2.6. Huygens’ Principle – Geometrical Description

To show the well-known fact that propagation of rays can be explained by a concept of

electromagnetic wave propagation, let us now formulate Huygens’ principle, presented

above, for free space without obstacles or discontinuities, i.e., for an unbounded

homogeneous medium. This case was an early description of what actually happens

with wave energy when it simply travels in free space in straight-ahead manner. Here,

in simple terms, the principle suggests that the energy from each point propagates in all

forward directions to form many elementary spherical wavefronts, which Huygens

called wavelets. The envelope of these wavelets forms the new wavefront. In other

words, each point on a wavefront acts as the source of secondary elementary spherical

waves described by Green’s function G( r, r' ) . These waves combine to produce a

new wavefront in the direction of wave propagation. With great accuracy each

wavefront can be represented by the plane which is normal to wave vector k (see Fig.

2.6, line AA’, as a starting wave position).

20
 Spherical elementary waves originate from every point on AA’ to form a new

wavefront BB’ which is drawn tangential to all elementary waves with equal radii. As

can be seen from the illustration of Huygens’ principle presented in Fig. 2.6a, the

secondary waves originating from points along AA’ do not have a uniform amplitude in

all directions.

In fact, if  represents the angle between the direction to any point C on the
elementary sphere (see Fig. 2.6b) and the normal to the wavefront (or parallel to k ),

then the amplitude of the secondary wave in a given direction is proportional to

(1  cos  ) . If so, the amplitude in the k -direction is proportional to ~ (1  cos 0)  2

In any other direction the amplitude is less than 2. In particular, the amplitude

of any elementary wave in the backward direction is ~ (1  cos  )  0 , that is, the

waves not propagate backward. The waves propagate forward along straight lines

normal to their wavefronts. Moreover, the consideration of elementary waves

originating from all points on AA’ leads to the expressions for the field at any point on

BB’ in the same form, the solution of which shows that the field at any point on BB’ is

exactly the same as that at the nearest point on AA’.

The phase difference between the oscillations at these neighbouring points of

lines AA’ and BB’ depends on the distance between them, d, and is therefore

proportional to ~ kd  2d /  . If d   , all points at the fronts AA’ and BB’ oscillate

in phase; if d   / 2 , all points oscillate in anti-phase etc.

21
Fig. 2.6a: Geometry of Huygens principle in free space.

Fig. 2.6b

22
Hence, from Huygens’ principle in the particular case of unbounded free space follows

the phenomenon of straight line wave propagation (that is, along the traces where is

the minimum of loss of energy, according to Fermat's principle), as light rays in optics.

Moreover, this also follows from Fermat principle, according to which rays propagate

in free space along the lines where the integral of energy is minimal, that is, along the

straight lines. This concept will be used for the case of optical rays' reflection and

diffraction from intersections by use the same vector and scalar wave equations

described in Lecture 2.

2.7. Fresnel-Zone Concept for Free Space

The existence of Fresnel zones also follows from Huygens’ principle not only in

obstructive conditions for both terminals, transmitter and receiver, when any obstacles

are placed around them and diffraction phenomenon is predominant.

In the case of free unbounded space, let us once more return to the integral

form (2.55) of EM-field presentation at any point r between the observer and the

source of radiation, using the Green’s function source presentation. In the case of free

space, instead of the virtual sources at any boundary limit of the space volume V,

within which the real source exists, we will introduce virtual sources to describe the

EM field in each point of the wavefront in space along the wave propagation path. A

main result, which follows from Huygens’ principle is that at any point of a source-free

unbounded medium, the total field is a superposition of elementary spherical waves,

which are radiated by virtual sources in space and reach to the observation point along

straight paths.

Thus, let us consider that the radiation source is placed in free space at point A

and the receiver is at point B, as shown in Fig. 2.7. We also consider an imaginary

23
plane with area S normal to the line-of-sight path at any point between A and B, which

passes across the point O at the line AB (see Fig. 2.7).

Now, if we “work” with infinite volume V, Green’s theorem (2.55) can be

represented for any vector of the EM wave, namely, the Hertzian vector, as

(R ' ) exp ik | R  R '|

(R )   n | R  R '|
ds

(2.58)

where | R  R '|  r is the distance from any point in the imaginary plane S and

observer at point B.

Fig. 2.7

If the initial radiation source can be assumed to be a point source with Green’s

e ikr1 r1 r
function G ~ , then for any point in the plane S, because   0 , we finally
r1 n r1

obtain:

1 1  
 r0 exp ik r1  r1
'
  ds
2    r1
(R ) ~   ik  (2.59)
 r1 r1r1'

All the mentioned distances are presented schematically in Fig. 2.7.

24
 We will talk about the so-called wave zone or far zone between the plane S and

two terminals A and B, i.e., when

kr1  1 , kr  1

(2.60)

In this case, the first term in brackets within integral (2.59) is less than the second one,

and can be neglected. Moreover, in the process of integration the variables r and r1

are changed. Therefore, because the inequality (2.60) is valid for line-of-sight

propagation links in free space, relatively small changes of variable r cause fast

 
oscillations of the product ~ exp ik  r1  r1  . On the other hand, this fact leads to fast
'

changes of sign both for real and imaginary parts of integrand in (2.59). At the same

time, other products in the integrand of (2.59) are changed more weakly with relatively

small deviations of r and r1 . In this case, the known method of stationary phase is

usually used to derive such an integral, containing both slow and fast terms inside the

integrand [4, 7, 10]. We will not present here all the complicated analyses and

derivations of the integral in (2.59), only its final form for the observed point B, as

   
( B) ~ 
ik
2 r0 r0'
 
exp ik r0  r0'    dxdy expi k2  r1  r1  x
'
2

 0 0 
(2.61)
 k  1 1 
exp i   '  y 2 
 2  r0 r0  

For each integral in the two-dimensional integral of (2.61) one can use the well known

integral presentation [4, 10]



 exp i  d 
2
i



Therefore,

25
 exp ikd 
( B) ~ 
ik
2 r0 r0'  
exp ik r0  r0'  i
k  1 1

d (2.62)
  
2  r0 r0' 

where d  r0  r0' is the distance between the source (point A) and the observer (point

B). Hence, as is shown from expression (2.62), if at the plane S the source A creates a

e ikr1
field ~ , then at the observed point B, the virtual dipoles, uniformly distributed at
r1

e ikd
S, will create a field ~ , the same as for the direct wave from A to B. This is the
d

main content of Huygens’ principle. Additional analysis of integral (2.51), extended

above S for the farthest wave zones

ik r0 exp ik r1  r1 
'

2   r1
( B) ~  ds
r1r1'

(2.63)

shows that the plane S can be split into the concentric circles (hoops) of arbitrary

radius. It is apparent that any wave, which has propagated from A to B via point Ci ,

i  1, 2,... , on any of these hoops has traversed a longer path than AOB (namely, A Ci

B>AOB at Fig. 2.7). While passing from one hoop to another, the real and imaginary

parts of the integrand in (2.63) change its sign. The boundaries of these hoops are

determined by the condition:

  
k r1  r1'  r0  r0'    n 2 , n  1, 2,... (2.64)

The physical meaning of these hoops for wave propagation is that if the virtual sources

of the elementary waves lie within the first hoop, they send to observer radiation with

the same phase for each elementary wave. Sources from two neighboring hoops send

respective radiation, which extinguish each other. Some elementary analysis of integral

26
(2.63) shows that a no-vanishing result exists only from the first, central hoop. The

hoops are usually called the Fresnel zones [1-3, 5-9]. Let us derive the width of these

zones b . For the first Fresnel zone, assuming x 2  y 2  b 2 , and using expressions

(2.63) and (2.64) for n=1, we have:

k  1 1 2 

  ' x  y 
2  r0 r0 
2

2
(2.65)

and

 1 1 r0 r0'
b1  x 2
y 2
    
k  r0 r0'  (r0  r0' )
~ R (2.66)

where R is the minimal range from each r0 and r0' .

Hence, the width of the first Fresnel zone is larger than wavelength, i.e.,

b1 ~ R   . For greater zones (with number n>1), for which the width b is

smaller than the distance to each zone center (simply, the radius of each circle b  bn ),

one can easily obtain after differentiation of (2.66) the following equation:

 1 1 
kbb  '  
 r0 r0  2

(2.67) and

 R b12
b ~ ~  b1 (2.68)
2k b b

Here R is, once more, the minimal range from each r0 and r0' .

Thus the width of the hoops with n>1 decreases with an increase of radius of

each zone bn . At the same time it can be shown that the radius of each Fresnel zone

of any specific number of the family of zones can be expressed in terms of zone

numbers n and the distance between both points A and B and the imaginary plane S as

[1-3, 9]

27
 nr0 r0'
b  bn  (2.69)
(r0  r0' )

from which, introducing in (2.69) n=1, one can immediately obtain expression (2.58)

the radius of the first Fresnel zone. As follows from (2.68) and (2.69) the width of the

Fresnel zones b decreases with increasing zone number n. At the same time, the

area of these zones is not dependent on zone number n, that is,


2 bb ~ R
2

(2.70)

It is clear that the radii of the individual hoops depend on the location of the imaginary

plane with respect to points A and B. The radii are largest midway between points A

and B and become smaller as the points are approached. Moreover, as follows from

(2.64), the family of hoops have a specific property: the path length from A and B via


each circle is n longer than the direct path AOB, i.e., for n=1 (first zone) ACB-
2

 
AOB= , so the excess path length for the innermost circle is . Other zones will
2 2


have an excess proportional to with parameter of proportionality n=2, 3, 4, ....
2


The loci of the points for which the excess A Ci B - AOB= n defines a family
2

of ellipsoids, the radii of which are described by expression (2.69). But in free space

without any obstacles, as we showed above mathematically, only the first ellipsoid is

valid and presents the first Fresnel zone which passes through both points, transmitter

(T) and receiver (R), as illustrated in Fig. 2.8.

28
 Fig. 2.8

This is why, despite the fact that in free space reflection and diffraction phenomena are

not observed, and no effect of interference between neighboring zones exists to

describe the loss-less phenomenon of wave propagation, the concept of Fresnel zones

is also used. This approach allows us to obtain the first hoop’s width in line-of-sight

propagation conditions and then to estimate through formula (2.58), by use of the

“working” frequency for the respective radio link, the range R of wave propagation in

conditions of direct visibility between any receiver and transmitter.

Bibliography
[1] Stratton, J. A., Electromagnetic Theory, New York: McGraw-Hill, 1941.
[2] Jackson, J. D., Classical Electrodynamics, New York: John Wiley & Sons, 1962.
[3] Chew, W. C., Waves and Fields in Inhomogeneous Media, New York: IEEE Press,
1995.
[4] Abramowitz, M., and I. A. Stegun, Handbook of Mathematical Functions, New
York: Dover Publications, 1965.
[5] Elliott, R. S., Electromanetics: History, Theory, and Applications, New York:
IEEE Press, 1993.
[6] Felsen, L. B., and N. Marcuvitz, Radiation and Scattering of Waves, New York:
IEEE Press, 1994.
[7] Al’pert, Ia. L., V. L. Ginzburg, and E. L. Feinberg, Radio Wave Propagation,
Moscow: State Printing House for Technical -Theoretical Literature, 1953.
[8] Kong, J. A., Electromagnetic Wave Theory, New York: John Wiley & Sons, 1986.

29
[9] Blaunstein, N., and Ch. Christodoulou, Radio Propagation and Adaptive Antennas
for Wireless Communication Links: Terrestrial, Atmospheric and Ionospheric, New
Jersey: Wireless InterScience, 2007.
[10] Dudley, D. G., Mathematical Foundations for Electromagnetic Theory, New
York: IEEE Press, 1994.
[11] Optical Fiber Sensors: Principles and Components. Ed. by J. Dakin and B.

Culshaw, Arthech House, Boston-London, 1988.

[12] Palais, J. C., Optical Communications, in Handbook: Engineering

Electromagnetics Applications, Ed. by R. Bansal, New York: Taylor and Frances,
2006.
[13] Kopeika, N. S., A System Engineering Approach to Imaging, Washington: SPIE
Optical Engineering Press, 1998.

30