05877153b961d8 - Basics of Electromagnetics & Transmission Lines - Jagadeeshwar Reddy - ch-1
05877153b961d8 - Basics of Electromagnetics & Transmission Lines - Jagadeeshwar Reddy - ch-1
1
2 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
KQ1Q2
F
R2
where K is proportional constant
In SI, a unit for Q1 and Q2 is coulombs(C), for R meters(m) and for F newtons(N).
1
K
4 0
where 0 = permittivity of free space (or) vacuum
= 8.854 × 10–12 farads/meter
9
= 10 farads/m
36
36 = 9 × 109 m/farads
K
4 10 9
Q1Q2
F …..(1.2.1)
2
4 0 R
Assume that the point charges Q1 and Q2 are located at (x1, y1, z1) and (x2, y2, z2) with
the position vectors r1 and r2 respectively. Let the force on Q2 due to Q1 be F12 which
can be written as
Q1Q2
F12 aR …..(1.2.2)
4 0 R 2 12
where aR12 is unit vector along the vector R12 . Graphical representation of the vectors
in rectangular coordinate system is shown in Fig.1.1
Where a x is the unit vector along X-axis and a y is the
i.e., R12 r2 r1
Fig. 1.1 Graphical
where r1 x1ax y1a y z1az representation of the
vectors
r2 x2 ax y2 a y z2 a z
STATIC ELECTRIC FIELDS 3
Q1Q2 R12
Now F12
4 0 R 2 R12
R12
aR12
R12
Q1Q2 R12
R12 R
4 0 R 2 R
Q1Q2 R12
4 0 R 3
Q1Q2 r2 r1
…..(1.2.3)
4 0 r2 r1 3
Q N
r rK
4 0
QK r rK
3
…..(1.2.4)
K 1
QQ Q
i.e., E …..(1.3.1)
Q 4 0 R 2
4 0 R 2
Consider a point charge Q with position vector r , then the electric field intensity E at
some point with position vector r1 due to point charge Q is
4 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Q
E aR …..(1.3.2)
4 0 R 2
where aR is the unit vector along R . Graphical representation of vector is shown in
Fig.1.2
From Fig. 1.2, R r1 r
Q R
E
4 0 R R
2
Q R
4 0 R R
2
If we have more than one point charge i.e., Q1, Q2,… QN with the position vectors
r1 , r2 ,....rN respectively. Then the electric field intensity E at some point with position
vector r can be written as
Q1 r r1 Q2 r r2 Q r rN
E .... N
4 0 r r1 3
4 0 r r2 3
4 0 r rN 3
1 N
r rK
4 0
QK r rK
3
…..(1.3.3)
K 1
Problem 1.1
Point charges 1 mC and –2 mC are located at (3, 2, –1) and (–1, –1, 4) respectively.
Calculate the electric force on a 10 nC charge located at (0, 3, 1) and the electric field
intensity at that point.
Solution
We know
Q 2
r rK
F
4 0
Q
K 1
K
r rK
3
=
10 109
1 103
3ax a y 2az
2 103
ax 4a y 3az
4 0
3 3
9 1 4 1 16 9
STATIC ELECTRIC FIELDS 5
3 2 1 8 2 6
= 90 103 ax ay az
52.38 132.57 52.38 132.57 52.38 132.57
= 90 103 0.0723 ax 0.0413 a y 0.0834 az
F
Also we know E
Q
0.0065 0.037 0.0075
= a
9 x
a
9 y
az
10 10 10 10 10 109
= 650 ax 370 a y 750 az kV/m.
Problem 1.2
Point charges 5 nC and –2 nC are located at 2ax 4az and –3ax 5az respectively.
(a) Determine the force on a 1 nC point charge located at ax 3a y 7 az . (b) Find the
electric field E at ax – 3a y 7az .
Solution
(a) We know
Q 2
r rK
F
4 0
QK r rK
3
K 1
= 10 9 10 10 5
9 9 9
ax 3a y 3az
2 4ax 3a y 2az
3 3
1 9 9 16 9 4
5 8 15 6 15 4
= 9 109 ax ay az
82.81 156.169 82.81 156.169 82.81 156.169
= 9 109 ax 0.112 a y 0.143 az 0.155
F
(b) E , here Q = 1 nC
Q
E 1.008 ax 1.287 a y 1.395 az V / m
Problem 1.3
Point charges Q1 and Q2 are respectively located at (4, 0, –3) and (2, 0, 1). If Q2 = 4 nC,
Find Q1 such that (a) The E at (5, 0, 6) has no Z-component. (b) The force on a test
charge at (5, 0, 6) has no X-component.
Solution
Q 2
r rK
We have F
4 0
QK r rK
3
K 1
9
F 1 Q1 (5,0,6) (4,0, 3) 4 10 (5,0,6) (2,0,1)
(a) E
Q 4 0
3 3
1 81 9 25
Given E has no Z – component, considering only Z components on both sides
1 Q1 9 4 109 5
0
4 0 82 3 3
34
Q1 9 4 109 5
3 3
82 34
3
20 41
Q1 nC 8.3nC
9 17
(b) Given the force on test charge has no X-component
Q Q1 4 109 3
0
4 0 82 3 3
34
Q1 4 109 3
=
3 3
82 34
3
41
Q1 12 nC 44.95nC
17
STATIC ELECTRIC FIELDS 7
Problem 1.4
Two point charges of equal mass ‘m’, charge ‘Q’ are suspended at a common point by
two threads of negligible mass and length ‘l’. Show that at equilibrium the inclination
angle ‘’ of each thread to the vertical is given by Q2 = 16 0 mgl2 sin2 tan , (or)
tan 3 Q2
,
1 tan 2 16 0 mgl 2
if ‘’ is very small
Q2
Show that 3
16 0 mgl 2
Solution:
When two charges are suspended from a common point with threads of length ‘l’, we
can represent graphically as sown in Fig.1.3, where T is the tension in thread ‘mg’ is the
weight of charge towards ground due to gravitational force and F is force on charge at
‘A’(B) due to charge at ‘B’(A). T cos is the vertical component of ‘T’ which is
upwards and T sin is the horizontal component of ‘T’ which is opposite to F . To form
equilibrium either at ‘A’ or ‘B’
T cos = mg …..(1.3.4)
T sin = F …..(1.3.5)
1.3.4 T sin F
1.3.5 T cos mg
8 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
F
tan
mg
Q2
where F
4 0 r 2
r/2
sin
From Fig.1.3 l
r 2l sin
Q2
tan
4mg 0 r 2
Q2
4mg 0 4l 2 sin 2
Q2
tan
16mgl 2 0 sin 2
Q2
sin 2 tan …..(1.3.6)
16mgl 2 0
Problem 1.5
Two small identical conducting spheres have charges of 2 × 10–9 and – 0.5 × 10 –9 C
respectively. (a) When they are placed 4 cm apart what is the force between them? (b) If
they are brought into contact and then separated by 4 cm. What is the force between
them?
Solution
(a) We know
Q1 Q2
F
4 0 R 2
1
9 109
4 0
2 109 0.5 109 9 109
F
16 104
= – 5.625 N
(b) When they are brought into contact, charges will be added and again when they are
separated charge will be distributed equally
Q1 = 0.758 × 10 –9 C
Q2 = 0.75 × 10 –9 C
F 3.164 N
Problem 1.6
If the charges in the above problem are separated with the same distance in a kerosene
(r = 2), then find (a) and (b) as in the previous problem.
Solution
5.625
(a) Fk N
2
= – 2.8125 N
3.164
(b) Fk = 1.582 N
2
10 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Problem 1.7
Three equal +Ve charges of 4 × 10 –9 C each are located at 3 corners of a square, side
20 cm. Determine the magnitude and direction of the electric field at the vacant corner
point of the square.
Solution
Fig. 1.4
E1 = Electric field intensity at Q4 due to Q1
Q1
=
4 0 R 2
= 900 V/m
E2 = 450 V/m
E3 = 900 V/m
The electric field intensity at vacant point is
E E2 E1 cos 45o E3 cos 45o
900 900
450
2 2
450 900 2
= 1722.792206 V/m
and the dot product of one unit vector with the other
one gives ‘0’.
i.e., ax a y = 0 a y az = 0 az ax = 0
i.e., ax a y az a y az ax az ax a y
Note: If the given problem is of circular symmetry, then it would be better to use
cylindrical coordinates rather than Cartesian coordinates.
Where is the radial distance from origin, is the
azimuthal angle from X-axis to the radial distance and Z
is same as in Cartesian coordinate system. The
cylindrical coordinate system is shown in Fig. 1.6.
Where a , a and az are unit vectors along radial
axis, azimuthal angle and z-direction respectively.
The dot product of any unit vector with itself gives
‘1’.
Fig. 1.6 Cylindrical coordinate
i.e. a a 1 a a 1 az az 1 system
The dot product of any unit vector with the other unit vector gives ‘0’
i.e. a a 0 a . az 0 az a 0
The cross product of any unit vector with the other unit vector, which is next to the
first one in anticlockwise direction, results last unit vector in the anticlockwise direction.
i.e., a a az a az a az a a
Consider a general vector A with components A, A, Az along the three axes, then it
can be represented as
A A a A a Az az
X 2 Y 2 2 X 2 Y 2
a y a sin a cos
The unit vector a z of Cartesian coordinate system and cylindrical Fig. 1.8
coordinate system is same az = az
We know that in Cartesian co-ordinate system
A Ax ax Ay a y Az az .
A Ax cos Ay sin a Ay cos Ax sin a Az az Fig. 1.9
A A a A a Az az
where
A = Ax cos + Ay sin
A = – Ax sin + Ay cos
A z = Az
i.e., in matrix form
A cos sin 0 Ax
A sin cos 0
Ay
Az 0 0 1 Az
Ax ax a ax a ax az A
or
A a a a y a
a y az A
y y
Az az a az a az az Az
Note:
1. Differential displacement or elemental length is
dl d a d a dz az
2. Differential or elemental normal area is dS d dz a
d dz a
d d az
3. Differential or elemental volume is dv = d d dz
The dot product of any unit vector with the other unit
vector gives ‘0’
Fig. 1.10 Spherical
i.e., ar a 0 a a 0 a ar 0 coordinate system
The cross product of unit vectors is: ar a a ,
a a ar , a ar a
y
From Fig.1.11 tan 1
x
x2 y 2
and tan
z z
x2 y 2
tan 1
z
relation between unit vectors of Cartesian and spherical co-ordinate systems is as follows:
ax = sin cos ar + cos cos a – sin a
a y = sin sin ar + cos sin a + cos a
az = cos ar – sin a
Note:
1. Differential displacement or elemental length is
dl dr ar rd a r sin d a
16 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Q= s ds
S
Let us find the electric field at point (x, y, z) due to line charge distribution along
Z- axis. We know electric field intensity due to line charge distribution as
l dl
E aR
l 4 0 R 2
R
where aR , dl = d z ,
R
18 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Since the charge distribution has cylindrical symmetry, we use cylindrical coordinate
system to obtain Electric field intensity.
From the Fig. 1.13
R a z z a z
L dz R
E
l 4 0 R 2 R
L dz R dz a z z az
L
4 0 R
4 0 2 z z 2 3/ 2
3
l l
L dz a z z az
l
4 0 2 z z 2 3/2
From the Fig. 1.13
z z
tan z z tan
2
cos R sec 2 z z ' sec
R
z = OT – z – z = OT – tan
d z = 0 – sec2 d
L sec d a tan az
2
4 0 l
E
3 sec3
L
2
sec2 d sec a cos az sin
4 0 3 sec3
1
L 2
4 0
(cos a sin az ) d
1
L
sin a 2 cos a z 2
4 0 1 1
STATIC ELECTRIC FIELDS 19
L
sin 2 sin 1 a cos 2 cos 1 az
4 0
which is electric field at point (x, y, z) due to line charge distribution from ‘A’ to ‘B’
along Z-axis. If ‘A’ is tending to – then 1 becomes /2 and ‘B’ is tending to then
2 becomes –/2.
L
E sin sin a cos cos az
4 0 2 2 2 2
2a L
4 0
L
E a
2 0 …..(1.5.4)
which is the electric field at point (x, y, z) due to infinite line charge distribution along
Z-axis.
Assume that the elemental surfaces are located on the sheet at ‘1’ and ‘2’ .
Then the elemental charge dQ on elemental surface ds is dQ = s ds.
The elemental electric field at point (0, 0, h) due to the elemental surface ds is
20 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
dQ
dE aR
4 0 R 2
where
R
dQ s ds and aR
R
Since the surface is infinite it has circular symmetry, hence we can use cylindrical
coordinate system to obtain electric field intensity.
Here ds lies on and axises, Hence ds d d
From Fig.1.14
a R haz
R haz a
dQ R
dE
4 0 R 3
dQ a haz
4 0 2 h 2 3/2
Since the sheet is symmetry with respect to origin on XY plane, for every electric field
due to elemental surface (for example elemental surface located at ‘1’) there will be an
equal and opposite electric field due to the elemental surface on the other side(for
example elemental surface located at ‘2’) in the direction of ‘ ’ (radial length), so finally
when we add up the electric fields due to all the elemental surfaces on the sheet the
electric field in the ‘’ direction will get cancelled. We will have only the electric field
perpendicular to the sheet i.e., along Z-direction.
Q haz
By integrating the above equation, E
4 0 2 h 2
3/2
Where Q s d d
2
haz 1
E s d d
0 0 2
h
2 3/ 2 4 0
STATIC ELECTRIC FIELDS 21
s 2 h
E d d az
4 0 0
3/ 2
0 h
2 2
s h
1
3/ 2
2 h2
2 d 2 az
4 0 0
2
2 3
s h 1
h 2 2 1
az
3
2 0 2 1
2 0
2 1/2
s h 1 h a
E
2 0 2 1 / 2
z
s
E az …..(1.5.5)
2 0
If we observe the above equation, the electric field is independent of the height ‘h’ i.e.,
the point can be considered at anywhere on the Z-axis.
The above equation can be generalized as
s
E an …..(1.5.6)
2 0
Where an is the unit vector which is perpendicular to the sheet.
Consider a parallel plate capacitor of equal and opposite charge on each plate, the
electric field due to these parallel plates can be written as
s s s
E an an an …..(1.5.7)
2 0 2 0 0
Problem 1.8
A circular ring of radius ‘a’ carries a uniform charge L C/m and is placed on the XY
plane with axis the same as the Z-axis.
L ah
(a) Show that E 0,0, h a .
3/ 2 z
2 0 h 2 a 2
Solution
(a) Here dl = a d
dQ = L dl
= L ad
dQ
dE ar
4 0 R 2
Fig. 1.15
a R
a r R ; r2 3
R R R
dQ a a ha z
dE
4 0 a 2 h 2 3/2
dQ = L a d
Q L a d
when we add up electric fields, the electric field in direction gets cancelled.
dQ haz
E
4 0 a 2 h 2
3/2
STATIC ELECTRIC FIELDS 23
L a d ha z
4 0 a 2 h 2 3 / 2
2
a ha z L ah
L
4 0 a h 2 3 / 2
2 d 2 a
0
2
h2
3/ 2
az
0
dE
(b) 0
dh
a 3
2 3/2
h2 .1 h (a 2 h 2 )1/2 2h
La 2
az 0
2 0
3
a h2
2
Q
E az
4 0 h 2
Problem 1.9
Derive an expression for the electric field strength due to a circular ring of radius ‘a’ and
uniform charge density L C/m. Obtain the value of height ‘h’ along Z-axis at which the
net electric field becomes zero. Assume the ring to be placed in X-Y plane.
Solution
Derivation is as in Problem. 1.8.
L ah
E az
3/2
2 0 a 2 h 2
La
E 3/2
az
a2
2 0 h 2 2 1
h
From the above equation we can say that for h = , the net electric field becomes zero.
Problem 1.10
A circular ring of radius ‘a’ carries uniform charge L C/m and is in XY-plane. Find the
Electric field at point (0, 0, 2) along its axis.
Solution
Replacing ‘h’ in problem.1.8 with ‘2’ and solving, we get
La2
E a
3/2 z
2 0 a 2 4
Due to symmetry, the electric field in ‘’ direction will be zero. Finally total electric
field will be in Z-direction.
v dv
Ez E az cos
v 4 0 R
2
Ez
v 4 0 R 2
By applying cosine rule in the Fig.1.16
(r ) 2 z 2 R 2 2 zR cos
r z 2 R 2
2
cos
2 zR
Similarly
R 2 z 2 ( r ) 2 2 zr cos
z 2 r R2
2
cos …..(1.5.8)
2 zr
On differentiating equation (1.5.8), we get
2 R
sin d dR
2 zr
R
sin d dR
zr
26 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
v a3 v 4 3
Ez a
0 z 2
3 4 0 z 3 2
Q
E az …..(1.5.9)
4 0 z 2
The electric field due to a sphere of radius ‘a’ with volume charge density v is similar
to the electric field due to a point charge which is placed at origin.
Problem 1.11
A circular disk of radius ‘a’ is uniformly charged with s C/m2. If the disk lies on the
Z = 0 plane with it’s axis along the Z-axis
h
(a) Show that at point (0, 0, h), E s 1 a
2 0 2 1/ 2
z
2
h a
(b) From this derive the E due to an infinite sheet of charge on the Z = 0 plane.
(c) If a << h, Show that E is similar to the field due to a point charge.
STATIC ELECTRIC FIELDS 27
Solution
Fig. 1.17
dQ
(a) dE ar
4 0 R 2
dQ = s ds; ds = d. d,
= s d d
a R haz
R haz a
s d d haz a
E
4 0
32
s h2 2
2 a
s h
E az d d
4 0 0
3/2
0 h
2 2
a
s 1
3/2
az 2 h h 2 2 d 2
4 0 0
2
28 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
a
2 3
s h 1
az
h 2 2 1
3
2 0 2 1
2 0
4 0
s h
a z 2 h 2 a 2
1/ 2
h2
1/2
s h az 1 1
2 0
h 2
a2 h
s h a
E 1
2 0 2 1/ 2
z
2
h a
(b) a ;
s
E az
2 0
(c) when a<< h, the volume charge density becomes a point charge located at origin,
Q
E az
4 0 h 2
Problem 1.12
The finite sheet 0 < x < 1, 0 < y < 1 on the Z = 0 plane has a charge density
s= xy (x2 + y2 + 25)3/2 nC/m2.
Find
(a) the total charge on the sheet
(b) the electric field at (0, 0, 5)
(c) the force experienced by a – 1 nC charge located at (0, 0, 5)
Solution
(a) dQ = s ds
Q= s ds
s
STATIC ELECTRIC FIELDS 29
1 1
3/2
= xy x 2 y 2 25 n dx dy
x 0 y 0
1 1
x
2 3/2 1
= n x y 2 25 d ( y 2 ) dx
x 0 y 0
2
1
5/2 1
= n x x 2 y 2 25 2 1 dx
x 0
0 5 2
1
n 2
1 d ( x2 )
5/2 5/2
5 x0
= x 26 x 2 25
2
7/2 1
n 2
1
7/2
= x 26 x 2 25
5 0 7
n
=
27 7/2 2 26 7/2 257/ 2
35
n
= 102275.868136 179240.733942 78125
35
Q = 33.15 nC
(b) Electric field at (0, 0, 5)
s ds
dE aR ; on Z-plane point is (x, y, 0)
4 0 R 2
R 0,0,5 x, y,0 xa x ya y 5a z
aR R xa x ya y 5az
2 3 3
R R x 2 y 2 25
ds R
E s
s 4 0 R
3
3/2
1 1 xy x 2 y 2 25 109 xax ya y 5az
3
dx dy
x 0 y 0
4 0
x 2 y 2 25
30 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
1 1
1
4 0 x 2 yax xy 2 a y 5 xyaz dx dy 109
x 0 y 0
1 1 1
1
1
y2 y3 y2
4 0 x 2 ax x a y 5 x az dx 109
2 0 3 0 2 0
x 0
1
1 x2 x 5
4 0
2
ax a y xaz dx 109
3 2
x 0
1 x 3 1 x2
1
5 x2
1
ax a y az 109
4 0 6 0 6 0 2 2 0
1 1 1 5 9
6 ax 6 a y 4 az 10
4 0
1 1 5
9 109 ax a y az 109
6 6 4
1.5 ax 1.5 a y 11.25 az V/m
(c) F = qE
Problem 1.13
A square plane described by –2 < x < 2, –2 < y < 2, z = 0 carries a charge density
12|y| mC/m2. Find the total charge on the plate and the electric field intensity at (0, 0, 10)
Solution
dQ = s ds
Q= s ds
s
2 2
= 12 | y | 103 dxdy
x 2 y 2
STATIC ELECTRIC FIELDS 31
2 0 2
= 103 12 y dy 12 y dy dx
x 2 y 2 y 0
0 2
y2 2
y2
= 10 12 12 dx
3
x 2 2 2 2 0
2
= 103 12 2 12 2 dx
x 2
2
= 48 103 dx 48 103 4 192 mC
x 2
s ds
dE aR ; R 0, 0, 10 x, y,0 xax ya y 10az
4 0 R 2
s ds R
dE
4 0 R 3
s ds R
E
s
4 0 R 3
12 | y | 103 xax ya y 10az
2 2
3
dx dy
x 2 y 2
4 0
x 2 y 2 100
0
2
xyax y 2 a y 10 yaz 2 xya y 2 a 10 ya
dy dx
6 x y z
9 10 12 dy
3/2 3/2
x 2 y 2 2 2
x y 100 y 0 2 2
x y 100
2 2 2
3/2 3/2
108 106 x 2 yax x 2 y 2 100
y 0
dy 10 2 yaz x 2 y 2 100 dy dx
x 2 y 0
2 2 2
az x y 100
3/2 3/2
108 106 ax x y 100
2 2
x d ( y 2 ) 10 2 2
d ( y 2 ) dx
x 2 y 0 y 0
2 2 2
x 2 y 2 100
2 1/2 1/2
x y 100
2
a 10
a dx
108 10 x
6
1/ 2 x 1/ 2 z
x 2
0 0
2
2 x x 2 104
a 20 x 2 104
a dx
1/2 1/2 1/2 1/2
108 106
2 x x 2 100
x
x 2 100
z
x 2
1/2 1/2
x x 2 104 & x x 2 100 are odd functions
x 104 x 100
2 1/2 2 1/ 2
and & are even functions
a
0 if f is odd
a
f ( x)dx 2 f ( x)dx if f is even
a
0
2
1 1
E = –20 × 108 × 10 × 2
6
az dx
2
x 0 x 2 104 x 2 102
2
x 1 x
40 108 106 sinh 1 sinh 10 az
104 0
2 1 1
40 108 106 sinh 1 sinh 5 az
104
E = 16.46 az MV/m.
STATIC ELECTRIC FIELDS 33
A cos dL A.dL
L L
= A. ds …..(1.6.2.1)
s
The formulae for D can be obtained by multiplying the formulae of E with 0.
A ds
div A A lim S
…..(1.7.1)
v 0 v
Where is the del operator or gradient operator. can be operated on a vector or
scalar. It has got different meanings when it is operating on a vector and scalar. If it is
operating on a scalar V then it can be written as V which is called as scalar gradient. If
it is operating on a vector A with dot product then it is . A and it is called as divergence
of vector A and If it is operating on a vector A with cross product then it is A and it
is called as curl of vector A .
Physically divergence can be interpreted as the measure of how much field diverges or
emanates from a point. Let us consider the Fig.1.21(a) in which field is reaching to the
point. Divergence at that point is –Ve or it is also called as convergence. In Fig.1.21(b)
the field is going away from the point, therefore divergence is +Ve. In Fig.1.21(c) some
of the flux lines or field lines are reaching to the point and same number of field lines are
leaving from the point hence the divergence is zero.
To determine .A let us consider the volume in
Cartesian co-ordinate systems as shown in the Fig.1.22. In
Cartesian co-ordinate system, the vector A with it’s unit
vectors and components along X, Y, Z is
A Ax ax Ay a y Az az
Ay Ay
Ay y x z . So the total flux on Y-axis is Ayxz + xyz – Ayxz
y y
Ay
= xyz
y
Similarly on X and Z-axises also.
A Ay Az
The entire flux in all the directions is x xy z . We know
x y z
A ds
s
a a y a z
A x .Ax a x Ay a y Az a z
x y z
Statement
This theorem states that the outward flux flows through a closed surface is same as the
volume integral of divergence of a vector.
A ds A dv ….. (1.7.3)
s v
Proof:
Consider a vector A Ax ax Ay a y Az az .
Assume dv = dx dy dz
consider the volume integral
Ax Ay Az
A dv x
y
dx dy dz
z
v v
Ax dsx
and Az dsz
s s
38 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Adv
( Ax dsx Ay ds y Az dsz )
v s
= Ax ax Ay a y Az az dsx ax ds y a y dsz az
A ds
s s
Hence proved
Formulae for Gradient
in Cartesian co-ordinate system
V V V
V ax ay az …..(1.7.4)
x y z
in cylindrical co-ordinate system
V 1 V V
V a a az …..(1.7.5)
z
in spherical co-ordinate system
V 1 V 1 V
V ar a a …..(1.7.6)
r r r sin
Problem 1.15
Find the gradient of the following scalar fields
(a) V = e–z sin 2x cos hy
(b) U = 2z cos2
(c) W = 10r sin2 cos
Solution
(a) Since given V is in x and y, consider gradient in Cartesian co-ordinate system
V V V
V ax ay az
x y z
r r sin
Formulae for Divergence of a Vector
in Cartesian co-ordinate system
Ax Ay Az
A
x y z …..(1.7.7)
in cylindrical co-ordinate system
1 A 1 A Az
A …..(1.7.8)
z
in spherical co-ordinate system
A
2
1 r Ar
1 sin A
1 A
…..(1.7.9)
r 2
r rsin rsin
Problem 1.16
Determine the divergence of the following vector fields.
(a) P x 2 yz ax x3 zy a y xy 2 z 3 az
(b) Q sin a 2 z a zcos az
1
(c) T cos ar r sin cos a cos a
r2
(d) N r 3 sin ar sin 2 cos 2 a cos r 2 a
Solution
(a) Given P x 2 yz ax x3 zy a y xy 2 z 3 az
40 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
1 Q 1 Q Qz
Q
z
1 1
2 sin 0 cos
= 2sin + cos
1
(c) Given T cos ar r sin cos a cos a
r2
2
1 r Tr
T 2
1 sin T
1 T
r r r sin r sin
1 1 1
2 0
0 r 2sin cos cos
r r sin r sin
= 2 cos cos
(d) Given N r 3 sin ar sin 2 cos 2 a cos r 2 a
N
2
1 r Nr
1 sin N
1 N
r 2
r r sin r sin
1 4 1 1 sin 3 2 1
5r sin sin cos 0
r 2
r sin 2 3 r sin
1 sin 3
5r 2 sin cos 2 cos 2
2r 6r sin
The charge enclosed within a volume or closed surface whose volume charge density
v is
Q v dv
v
D ds D dv …..(1.8.1b)
s v
By comparing the volume integrals in equations (1.8.1a) and (1.8.1b) we can write as
v D …..(1.8.2)
which is the Maxwell’s first equation for electrostatics (time in-varying fields)
Consider unsymmetrical distribution as shown in
Fig. 1.23a. The flux flowing through the closed surface
shown in Fig. 1.23a is = 5 – 2 = 3 nC. The charge
enclosed by the surface is Q = 3 nC.
Consider an empty closed surface as shown in Fig. 1.23 Closed surface
Fig. 1.23b. Flux flowing through the closed surface shown
in Fig. 1.23b is = 0 and hence charge enclosed by the closed surface is zero.
Conclusion
Gauss law holds good even if the charge distribution is unsymmetrical as shown in
Figs.1.23a & b. But to find either E or D , the charge distribution must be symmetrical. It
can be rectangular symmetry or cylindrical symmetry or spherical symmetry.
42 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
If the continuous charge distribution depends on either ‘x’ or ‘y’ or ‘z’, then the
distribution will have rectangular symmetry. So to find either E or D , we can use
rectangular co-ordinates.
If the continuous charge distribution depends only on and is independent of and z
then the distribution will have cylindrical symmetry. So, to find either E or D , we can use
cylindrical co-ordinates.
If the continuous charge distribution depends on ‘r’ and is independent of and then
the symmetry it will have is spherical. So to find either E or D , we can use spherical
co-ordinates.
d s r 2sin d d ar
Flux flowing through the sphere is
D.ds
s
2
Dr ar .r 2 sin d d ar
0 0
2
Dr r 2 sin d d
0 0
Fig. 1.24 Gaussian surface
2
about a point charge
r cos 0 d
2
Dr
0
2 Dr r 2 2
4 r 2 Dr
STATIC ELECTRIC FIELDS 43
Here D the electric flux density is perpendicular to the surface of the cylinder i.e., it
will be in ‘’ direction in cylindrical co-ordinate systems.
44 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
D D a
l 2
D a . d dza
z 0 0
l 2
D dz d
z 0 0
D 2 l …..(1.8.3)
L
D a and
2
D L
E a …..(1.8.5)
0 2 0
Which is similar to the formula derived by using Coulomb’s law.
Problem 1.17
Given D z cos 2 az C/m2. Calculate the charge density at (1, /4, 3) and the total
charge enclosed by the cylinder of radius 1m with –2 < z < 2 m.
STATIC ELECTRIC FIELDS 45
Solution
We know
v D
in cylindrical co-ordinate system the divergence can be written as
1 D 1 D Dz
v
z
Dz
v since D has only Z- component
z
v cos 2
v 1, , 3 1 cos2
1
C/m3
4 4 2
1 2 2
Qenc cos 2 d d dz
0 0 z 2
1 2
2 cos 2 4 d d
0 0
1
1
2 sin 4 d
1
4 2
2 2
0
1
1
3 4
4 2 d 4
3 0 3
C
0
Problem 1.18
If D 2 y 2 z a x 4 xy a y x a z C/m2. Find
Solution
According to Maxwell’s I equation
v D
1 1 1
4 x dx dy dz
x 0 y 0 z 0
1 1
4 x 1 dx dy
x 0 y 0
1
4 x 1 dx
x 0
1
x2 4
4 2C
2 0 2
Problem 1.19
Given the electric flux density D 0.3r 2 ar nC/m2, in free space. Find
(a) E at point (2, 25o, 90o)
(b) the total charge within the sphere r = 3
(c) the total electric flux leaving the sphere r = 4
Solution
(a) Given D 0.3r 2 ar nC/m2
D 0.3r 2 ar
E
0 8.854 1012
STATIC ELECTRIC FIELDS 47
0.3 4
E 2 , 25o ,90o 8.854 1012 ar 1.355 1011 ar 109 135.5 ar V/m
(b) we know v D
2
1 r Dr
1 sin D
1 D 1
= 2 0.3 4 r 3 n = 1.2r n
r 2
r r sin r sin r
= r2 sin d d dr
3 2
Q 1.2r n r 2 sin d d dr
r 0 0 0
3
n 1.2 r
3
sin 2 d dr
r 0 0
r cos 0 dr
2.4 n 3
r 0
3
r4
2.4 n 2 305.4 nC
4 0
4 2
1.2 nr sin d d dr
3
(c) Q
r 0 0 0
Dz ds ds
top bottom
Assume that the area of the elemental surface as A, then
Dz A A
2 ADz
Charge enclosed by the rectangular box is
Qenc s ds
Qenc s ds
Qenc s A
STATIC ELECTRIC FIELDS 49
Qenc v dv
v
2 r
v r sin d d dr
2
Fig. 1.27 Gaussian surface for
0 0 r 0 uniformly charged sphere
4
= v r 3
3
The flux flowing through the spherical surface
D.ds
s
As the flux density is normal to the surface it will have components only in ‘r’
direction.
50 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
2
= Dr r 2 sin d d
0 0
= Dr 4 r2
According to Gauss’s law charge enclosed = flux flowing through the surface
i.e., Qenc =
4
v r3 = Dr 4 r2
3
v
Dr r
3
v
D r ar
3
D
and E v r ar …..(1.8.7)
0 3 0
Case II (r > a)
To find the electric flux density out side the sphere of radius ‘a’, consider a sphere of
radius ‘r’, which is treated as Gaussian surface as shown in Fig.1.28.
Charge enclosed by the sphere of radius ‘r’ is
Qenc v dv
v
Qenc v dv
v
2 a
v r sin d d dr
2
Fig. 1.28 Gaussian surface
0 0 r 0 for uniformly charged sphere
4
= v a3
3
Flux flowing through the surface
2
Dr r 2sin d d
0 0
= Dr 4 r2
Qenc = according to Gauss’s law
4
v a3 = Dr 4 r2
3
STATIC ELECTRIC FIELDS 51
v
Dr a3
3r 2
v a3
D ar
3r 2
v a3
and E ar
3r 2 0
v a 3 4
E ar
3r 2 0 4
Q
E ar
4 0 r 2
…..(1.8.8)
which is similar to the formula derived by using Coulomb’s law.
Problem 1.20
A charge distribution with spherical symmetry has density
r
0 , 0 r R
v R
0, rR
Determine E everywhere
Solution:
Replace ‘a’ with ‘R’ in Fig.1.27, Then
Case I: Inside the sphere of radius ‘R’
The charge enclosed by the sphere of radius ‘r’ is Qenc v dv
v
r
Qenc 0 dv
v
R
2 r
r3
0 R sin d d dr
0 0 r 0
2 r
0
d sin d r dr
3
R 0 0 r 0
52 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
4 r 4 0
4R
0
Qenc r4
R
The flux flowing through the spherical surface
D.ds
s
As the flux density is normal to the surface it will have components only in ‘r’
direction.
2
= Dr r 2 sin d d
0 0
= Dr 4 r2
According to Gauss’s law charge enclosed = flux flowing through the surface i.e.,
Qenc =
0
r 4 Dr 4 r 2
R
0
Dr r2
4R
0
D r 2 ar
4R
D 0 2
and E r ar
0 4 R 0
Case II: Outside the sphere of radius ‘R’
Charge enclosed by the sphere of radius ‘r’ is
Qenc v dv
v
r
Qenc 0 dv
v
R
2 R
0
r sin d d dr
3
R 0 0 r 0
= 0 R 3
Flux flowing through the surface
STATIC ELECTRIC FIELDS 53
2
Dr r 2 sin d d
0 0
= Dr 4 r2
Qenc = according to Gauss’s law
0 R3 = Dr 4 r2
0 R 3
Dr
4r 2
R3
D 0 2 ar
4r
R3
and E 02 ar
4r 0
Problem 1.21
A sphere of radius ‘a’ is filled with a uniform charge density of v C/m3. Determine the
electric field inside and outside the sphere.
Solution
The answer is as derived in section 1.8.5 case-I(inside the sphere) and case-II(outside the
sphere).
Problem 1.22
A charge distribution in free space has v = 2r nC/m3 for 0 < r < 10 m and ‘0’ otherwise.
Determine E at r = 2 m and r = 12 m
Solution
Replace ‘a’ with ’10 m’ in Fig.1.27, Then
At r = 2 m
The charge enclosed by the sphere of radius ‘2m’ is Qenc v dv
v
Qenc 2rndv
v
2 2
r sin d d dr
3
2n
0 0 r 0
= 32 nC
54 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
As the flux density is normal to the surface it will have components only in ‘r’
direction.
2
= Dr r 2 sin d d
0 0
= Dr 16
According to Gauss’s law charge enclosed = flux flowing through the surface i.e.,
Qenc =
32 n = Dr 16
Dr 2n
D 2n ar and
D
E 226ar V m
0
At r = 12 m
Charge enclosed by the sphere of radius ’12 m’ is
Qenc v dv
v
Qenc 2rndv
v
2 10
r sin d d dr
3
2n
0 0 r 0
= 20 C
Flux flowing through the surface
2
Dr r 2 sin d d
0 0
= Dr 4 122
Qenc = according to Gauss’s law
STATIC ELECTRIC FIELDS 55
20 = Dr 4 122
Dr 0.0347
D 0.0347 ar and
E 3.92 ar k V m
B
W Q E dL
A
B
W
E dL
Q A
which is work done per unit charge and it is also called potential difference VAB.
Q
We know that electric field intensity E due to a point charge is ar
4 0 r 2
Where rA and rB are position vectors of point A and point B from origin
Q 1 1
VAB
4 0 rB rA
Q Q
VB V A …..(1.9.1)
4 0 rB 4 0 rA
Where VB and VA are absolute potentials at point B and A respectively. From the
above equation VAB is the potential at B with reference to the potential at A.
If A is at then VA = 0.
The above equation can be generalized for a potential (V) at any point having distance
‘r’ as
Q
V (Here Q is located at origin) …..(1.9.2)
4 0 r
If the point charge is placed at a distance r , then the electric potential at point ‘r’ can
be written as
Q
V …..(1.9.3)
4 0 r r
If we have ‘n’ number of point charges Q1, Q2, …,Qn with position vectors r1, r2, ….,
rn respectively, then the potential at ‘r’ is
Q1 Q2 Qn
V .... …..(1.9.4)
4 0 r r1 4 0 r r2 4 0 r rn
STATIC ELECTRIC FIELDS 57
For line charge distribution with charge density L, in the above equation Q can be
replaced by L dL .
For surface charge distribution with charge density s, in equation (1.9.4), Q can be
replaced by s ds .
Similarly Q can be replaced by v dv , For volume charge distribution with charge
density v .
Problem 1.23
Two point charges – 4 C and 5 C are located at (2, –1, 3) and (0, 4, –2) respectively.
Find the potential at (1, 0, 1). Assuming ‘0’ potential at infinity.
Solution
Q1 Q2
V
4 0 r r1 4 0 r r2
4 106 5 106
V
4 0 (1,0,1) (2, 1,3) 4 0 (1,0,1) (0, 4, 2)
Simplifying, we get
V = –5.872 kV
Problem: 1.24
A point charge 3 C is located at the origin in addition to the two charges of previous
problem. Find the potential at (–1, 5, 2). Assuming V() = 0.
Solution:
r – r1 = 1 25 4 = 5.478
r – r2 = 9 36 1 = 6.782
r – r3 = 16 1 1 = 4.243
3 103 4 103 5 103
V 9
5.478 6.782 4.243
= 10.23 kV
Problem 1.25
A point charge of 5 nC is located at the origin if V = 2 V at (0, 6, –8) find
58 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Solution
Q 1 1
(a) VA V
4 0 rA r
rA = (–3, 2, 6) – (0, 0, 0) = 32 22 62 7
rB = (1, 5, 7) – (0, 0, 0) = 1 52 7 2 = 75
5 109 1 1
VB 2 9
10 75 10
4
36
VB = 2.696 V.
(c) VAB = VB – VA = –1.233 V
*Problem 1.26
A point of 5 nC is located at (–3, 4, 0), while line y = 1, z = 1 carries uniform charge
2 nC/m.
(a) If V = 0 V at O(0, 0, 0), find V at A(5, 0, 1).
(b) If V= 100 V at B(1, 2, 1), find V at C(–2, 5, 3).
(c) If V = – 5 V at O, find VBC.
Solution:
Let the potential at any point be
STATIC ELECTRIC FIELDS 59
V = VQ +VL
Where VQ is potential due to point charge
Q
i.e., VQ
4 0 r
by neglecting constant of integration
and VL is potential due to line charge distribution,
for infinite line, we have
L
E a
2 0
L
VL E .dl a .d a
2 0
L
VL ln
2 0
by neglecting constant of integration.
Here ρ is the perpendicular distance from the line y = 1, z = 1(which is parallel to the
x-axis) to the field point.
Let the field point be (x, y, z), then
( x, y, z ) ( x, 1, 1) ( y 1) 2 ( z 1) 2
L Q
V ln
2 0 4 0 r
by neglecting constant of integration.
(a) O (0, 0, 0) (0,1,1) 2
A (5,0,1) (5,1,1) 1
rO (0,0,0) (3, 4,0) 5
L Q 1 1
VO VA ln O
2 0 A 4 0 rO rA
60 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
2 109 2 5 109 1 1
0 VA ln
109 1 109 5 9
2 4
36 36
1 1
VA 36ln 2 45
5 9
VA 36ln 2 4 8.477 V
(b) B (1, 2,1) (1,1,1) 1
C (2,5,3) (2,1,1) 20
circulation = A.dL . This circulation depends upon the area chosen in the non
conservative field. Let the area be S. Then the ratio of A.dL to S can be considered
as one unit. As the field is normal to this unit we can write the above expression as
A.dL a .
n
S
In general will be from point to point. This can be denoted by taking Limit S0
which gives curl of vector A i.e.,
A lim
A dL a …..(1.10.1)
n
S 0 S
The curl of vector A gives circulation that exists on the chosen closed surface.
As A or curl of a vector A is a vector. It can be
represented with three components in a rectangular
co-ordinate system i.e., [curl A ]1, [curl A ]2, [curl A ]3 along
X, Y & Z axises with ax , a y & az as unit vectors respectively.
[curl A ] = lim
A dL a
x
yz 0 yz
Let the components of vector A at P be (Ax, Ay, Az)
Q
The line integral from P to Q is A dL .
P
At ‘R’ we have moved by a distance z. As RS line is parallel to Y-axis, consider the
Ay
Y-component at ‘R’ as Ay z
z
S
Ay
A dL Ay z y
R
z
P
At ‘S’ to find A dL consider the ‘z’ component at ‘S’ which is Az
S
P
A dL Az z
S
Q R S P
A dL P A dL Q A dL R A dL S A dL
Az Ay
= Ay y Az z y z Ay y z y Az z
y z
A Ay
= z z y
y z
A Ay
z z y
y
z
A Ay
[Curl A ]1 = ax = z ax
zy
y z
Similarly we can also construct equations for [Curl A ]2 by considering the elemental
surface on ZX plane which is perpendicular to Y-axis
A A
[Curl A ]2 = a y x z
z x
and for [Curl A ]3 we have to consider the elemental surface on XY plane which is
perpendicular to Z-axis
Ay Ax
[Curl A ]3 = az
x y
Curl A = [Curl A ]1+ [Curl A ]2+ [Curl A ]3
A Ay Ax Az Ay Ax
A ax z ay az x y
y z z x
Which can be written in matrix form as
Cartesian co-ordinate system:
ax ay az
A
x y z
Ax Ay Az
Cylindrical co-ordinate system:
a a az
1
A
z
A A Az
Spherical co-ordinate system:
ar ra r sin a
1
A 2
r sin r
Ar rA r sin A
64 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Stoke’s Theorem
Stoke’s theorem gives the relation between the line integral and surface integral as
A ds A dL …..(1.10.2)
S L
where A is the field vector. According to above equation finding curl of a vector at
every point in a chosen surface and adding all those values will be equal to the contour
line integral of the boundary of the chosen surface.
Proof:
Let us consider a rotational field and choose a surface on it as shown in Fig.1.33.
We know that
A lim
A dL a
n
s 0 S
The above equation can be written as
A ds A dL
S L
A ds1 A dL A dL A dL A dL
S A B C D
D E F A
A ds2 A dL A dL A dL A dL
S A D E F
A ds A ds1 A ds2
S S S
B C D A A E F A
= A dL A dL A dL A dL A dL A dL A dL A dL
A B C D D D E F
B C D E F A
= A dL A dL A dL A dL A dL A dL
A B C D E F
From the above equation by finding curl of a vector A at all the points in a chosen
surface and adding up all the values will be equal to the contour line integral of the
chosen boundary surface. Adding up all the curls is nothing but integrating the curl of a
vector w.r.t. chosen surface.
A ds A dL
S L
…..(1.11.1)
A
A L
The total work done in moving a point charge from A to B and back to A is ‘0’.
From equation (1.11.1) we can say that the electrostatic fields are conservative fields
or irrotational fields.
According to Stoke’s theorem E ds
E dL
S L
66 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
E ds 0 or E 0 …..(1.11.2)
S
dv Ex dx E y dy Ez dz …..(1.11.3)
E V
Which is the relation between E and V.
Problem 1.27
10
Given the potential V sin cos
r2
(a) Find the electric flux density D at (2 , / 2 ,0)
(b) Calculate the work done in moving a 10 mC charge from point A(1, 30o, 120o) to
B(4, 90o, 60o)
STATIC ELECTRIC FIELDS 67
Solution
(a) We have
E V
Since V is given in spherical co-ordinate system, consider V in spherical co-
ordinate system
v 1 v 1 v
E ar a a
r r r sin
1 10cos cos 1 10sin sin
= 10 2r 3 sin cos ar
r r2
a
r sin r2
a
1 10cos cos 1 10sin sin
= 10 2r 3 sin cos ar
r r2
a
r sin r2
a
20sin cos 10cos cos 10sin
= 3
ar 3
a a
r r r3
10
=
r3
2sin cos ar cos cos a sin a
D E 0
8.825 1011
= 2sin cos ar cos cos a sin a
r3
8.825 1011
= 2.1.1ar 0 0
r3
D (2,
, 0) = 22.1 a pC/m2
r
2
B
(b) Work done = Q E dL Q VAB
A
= Q (VB – VA)
10 1
VB 1. 0.3125 V
16 2
10 1
VA 0.5 5 0.5 2.5 V
1 2
VB – VA = 2.8125 V
W = 10–3 × 10 × (VB – VA) = 28.125 mJ
68 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Problem 1.28
Given that E 3 x 2 y ax xa y k V m . Find the work done in moving a –2 C charge
from (0, 5, 0) to (2, –1, 0) by taking the path
(a) (0, 5, 0) (2, 5, 0) (2, –1, 0)
(b) y = 5 – 3x
Solution
(a) Line equation for (0, 5, 0) to (2, 5, 0) is
x x1 y y1 z z1
x1 x2 y1 y2 z1 z2
x0 y 5 z 0
02 55 00
y=5 z=0
dy = 0 dz = 0
2 ,5 ,0
3x y ax xa y dx ax dya y dzaz
2
W1 QK
0 , 5, 0
2 , 5 , 0
3x y dx xdy
2
QK
0 , 5, 0
2
3x
3 2
2 10 5 dx 0
0
x3 2
2 10 3 5 2
3
3 0
= 36 mJ
Line equation for (2, 5, 0) to (2, –1, 0)
Z=0 dz = 0
x2 y 5 z 0
2 2 5 1 0 0
x=2 dx = 0
STATIC ELECTRIC FIELDS 69
2 , 1, 0
3x y dx xdy
2
W2 QK
2 , 5, 0
1
W2 QK 2dy 2QK 1 5 24 mJ
5
W = W1 + W2 = 12 mJ
(b) Line equation for (0, 5, 0) to (2, 5, 0) is y = 5 – 3x
dy = –3dx
2 , 1, 0
3x y dx xdy
2
W QK
0 , 5, 0
2
W 2 103 3 x 2 5 3 x dx 3x dx =12 mJ
0
Let us find potential at P ( r, , ) due to electric dipole. We know the potential at ‘P’
Q Q
due to a point charge +Q is VQ and potential at ‘P’ due to –Q is VQ
4 0 r1 4 0 r2
Potential at ‘P’ due to electric dipole is
Q Q
V
4 0 r1 4 0 r2
Q 1 1
4 0 r1 r2
x d
from Fig.1.34 cos x cos
d 2 2
r1 r x
d
r1 cos r
2
y d
cos y cos
d /2 2
r r2 y
d
r2 r cos
2
Q 1 1
V
4 0 r d cos r d cos
2 2
Q d cos
4 0 2 d
2
r cos
2
if r >>d
Q d cos
V …..(1.12.1)
4 0 r 2
STATIC ELECTRIC FIELDS 71
1
V due to dipole
r2
d .ar = d cos and here define electric dipole moment p Q d whose unit is C-m.
V
Q d ar p ar
1
p.
r
…..(1.12.2)
4 0 r 2
4 0 r 2
4 0 r 2
r
If the electric dipole center is other than origin, let it be at r then the above equation
can be generalized as
1 r r' p.( r r ')
V p. …..(1.12.3)
4 0 r r '
2
r r ' 4 0 r r ' 3
The electric field due to dipole with center at the origin can be obtained as
E V
Since V in equation (1.12.1) is in terms of r and consider V in spherical co-ordinate
system, then
v 1 v
E ar a
r r
1 1 Qd
E 2 r 3 cos ar sin 2 a
rr 4 0
Qd 3 sin
2r cos ar 3 a
4 0 r
Qd
2 cos ar sin a
4 0 r 3
p
2cos ar sin a …..(1.12.4)
4 0 r 3
Problem 1.29
An electric dipole located at the origin in free space has a moment
p 3ax 2a y az n Cm
Solution
(a) We have
1 r r'
V p.
4 0 r r '
2
r r'
r ' = (0, 0, 0)
r r' 4 9 16 29
(3ax 2a y az ) (2ax 3a y 4az )
V 9 109 109
29 29
9 4
= 0.235 V
29 3/ 2
(b) r = 2.5 = 30o = 40o
x = r sin cos = 0.958
y = r sin sin = 0.8035
z = r cos = 2.165
upon simplifying we get
V = 1.97 V
If the surface is having continuous charge distribution then the above equation becomes
1
2 L
WE LV dL for line charge distribution …..(1.12.9)
1
2 S
WE sV dS for surface charge distribution …..(1.12.10)
1
2 v
WE V V dv for volume charge distribution …..(1.12.11)
A V AV A V
i.e., D V DV D V
1
from (1.12.12) WE
2 v
DV D V dv
1 1
WE
2v DV dv D V dv
2v
According to divergence theorem, first integral can be written as
1 1
2 S
WE DV . dS D V dv
2v
1 1
For point charges the potential V , E 2
r r
STATIC ELECTRIC FIELDS 75
1 1
For dipoles the potential V 2
,E 3
r r
Surface ds r2
1
If we consider the point charges the product of V and E and product of
r3
1
DV and dS . For very large surface the first integral will become zero.
r
1
2 v
WE D V dv
1 1
2v D.( E )dv D.E dv
2v
D 0 E
1
2 v
Energy = WE 0 E E dv Joules …..(1.12.13)
dWE 1
The energy density J/m3 is 0 E 2 wE J/m 3 …..(1.12.14)
dv 2
Problem 1.30
Three point charges –1 nC, 4 nC and 3 nC are located at (0, 0, 0), (0, 0, 1) and (1, 0, 0)
respectively. Find the energy in the system.
Solution
WE = W1 + W2 + W3
= 0 + Q2V21+ Q3(V31+V32)
Q1 Q3 Q1 Q2
Q2 .
4 0 r2 r1 4 0 r3 r1 r3 r2
1 QQ
Q1Q2 Q1Q3 2 3
4 0 2
1 12 18
9
4 3 .10
10 2
4 .
36
76 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
12
9 7 nJ 13.37 nJ
2
Problem 1.31
Point charges Q1 = 1 nC, Q2 = – 2 nC, Q3 = 3 nC and Q4 = – 4 nC are positioned one at a
time and in that order at (0, 0, 0), (1, 0, 0), (0, 0, –1) and (0, 0, 1) respectively. Calculate
the energy in the system after each charge is positioned.
Solution
Energy after Q1 is positioned is W1 = 0
Q1 2 1 1018
W2 = Q2V21 Q2 . = –18 nJ
4 0 r2 r1 109
4 . (1,0,0) (0,0,0)
36
Energy after Q2 is positioned W = W1+ W2 = –18 nJ
2
Energy after Q3 is positioned
W3 = W2 +Q3(V32+V31)
dQ
i
dt
Current Density
The current i flowing through a surface s is denoted as Jn = i/S A/m2.
i = Jn S
If current density Jn is perpendicular to the surface S
i = Jn S
If Jn is not perpendicular to S, then i J S
The total current flowing through surface is I J ds .
S
Based on how the current I is produced, the current densities are classified in to
(i) convection current density (ii) conduction current density and (iii) displacement
current density.
Convection Current Density
Conductors are not involved for flowing current in case of convection current. Hence it
will not satisfy ohm’s law. The current flowing through an insulating material like liquid
or vacuum is convection current. A beam of electrons through a vacuum tube is an
example of convection current.
Consider a filament which is having volume charge density v as shown in Fig.1.37
We know that
Fig.1.37 Current in a filament
Q = v V = v s l
Dividing with t
Q l
v S
t t
Q l
I v SU y I and Uy
t t
I
the current density J vU y
S
78 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
eE
Conduction current density J ne
m
J ne 2 E …..(1.13.4a)
m
J E
where
= conductivity of the conductor = ne 2 …..(1.13.4b)
m
STATIC ELECTRIC FIELDS 79
Problem 1.32
1
If J 2cos ar sin a A/m2. Calculate the current passing through
r3
(a) Hemispherical shell of radius 20 cm.
(b) A spherical shell of radius 20 cm.
Solution
I J .ds
2 2
1
r 0 0
= sin 2 d d
2 /2
1 cos 2
= d
r 0 2 0
1 2
= 1 1 2 10 31.4 A
2r 0.2
(b) = 0 to 2, = 0 to and r = 0.2 m for spherical shell
2
1
I sin 2 d d
r 0 0
2
1 cos 2
r 0 2 0
d
2
1
1 1 d 0 A
2r 0
80 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
Problem 1.33
For the current density J 10 z sin 2 a A/m2. Find the current through the cylindrical
surface = 2, 1< z < 5 m.
Solution
Since it is cylinder ds d dz a
We have
I J ds
5 2
10 z sin d dz
2
z 1 0
5
10 z (1 cos )
z 1
= 754 A
*Problem 1.34
In a cylindrical conductor of radius 2 mm, the current density varies with distance from
the axis according to J 103 e 400 r A/m2. Find the total current I.
Solution
Since it is cylinder ds d dz a
Here r = ρ = 0.02 m,
J 103 e 400 a A m 2
2 z
10 e
3 400
I d dz
0 z 0
I 2 z103 e 400
I 4 ze0.8 z5.65 A
STATIC ELECTRIC FIELDS 81
Problem 1.35
1
If the current density J (cos ar sin a ) A m 2 , find the current passing through a
r2
sphere of radius 1.0 m.
Solution
We know the total current I J .ds
s
I sin 2 d
0
cos 2
0 A
2 0
Consider a volume which has dielectric material with volume charge density v. Then
the total volume charge density T v pv D
v pv 0 E
v 0 E pv
v 0 E P pv P
v 0 E P
v D
where D 0 E P …..(1.14.2)
The electric flux density D in free space is 0 E i.e., P = 0 in free space.
From the above equation we can say that D is getting increased by P in dielectric
materials.
STATIC ELECTRIC FIELDS 83
0 1 X E E
0 r E
E
where = 0.r r 1 XE
0
Where is the permittivity of dielectric material and 0 is the permittivity of free
space and r is the dielectric constant or relative permittivity. The dielectric constant
r can be defined as the ratio of to 0 .
If electric field strength is more such that it pulls the electrons from the outer shells of
dielectric molecules, then the dielectric material becomes conducting material and we can
say dielectric material has been broken.
Dielectric strength can be defined as the maximum electric field with which dielectric
material can tolerate or withstand.
Isotropic dielectrics: If does not change with the direction then we can say the
dielectric as isotropic dielectric.
Similarly conducting materials are classified as
If ‘’ is independent of E then the conducting material is linear conducting material.
If ‘’ is independent of direction then the conducting material is isotropic conductor.
If ‘’ does not change from point to point then the conducting material is
homogeneous conductor.
dQ dQ d
v dv
dt
can be written as
dt dt
= v dv …..(1.16.1b)
V
t
equations(1.16.1a) = (1.16.1b)
Jdv t v dv
v v
v
J …..(1.16.1c)
t
which is the continuity current equation.
The left side of the equation is the divergence of the Electric Current Density J .
This is a measure of whether current is flowing into a volume (i.e., the divergence of J is
positive if more current leaves the volume than enters).
Recall that current is the flow of electric charge. So if the divergence of J is positive,
then more charge is exiting than entering the specified volume. If charge is exiting, then
STATIC ELECTRIC FIELDS 85
the amount of charge within the volume must be decreasing. This is exactly what the right
side is a measure of how much electric charge is accumulating or leaving in a volume.
Hence, the continuity equation is about continuity - if there is a net electric current is
flowing out of a region, then the charge in that region must be decreasing. If there is more
electric current flowing into a given volume than exiting, then the amount of electric
charge must be increasing.
E v
v
E …..(1.16.2)
Consider the conduction current equation (point form of ohm’s law)
J E …..(1.16.3)
v
From (1.16.2) E
v
J from (1.16.3)
v
v from continuity equation
t
v
t
v
on integrating
ln v t ln v0
v t
e e
t / /
…..(1.16.4)
v0
v0 = initial volume change density
Tr
86 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
we know E V
V v …..(1.17.1a)
which is the Poisson’s equation for in-homogeneous medium.
For charge free medium v = 0
V 0 …..(1.17.1b)
which is the Laplace’s equation for in-homogeneous charge free medium.
For homogeneous medium since is constant
2 v
V …..(1.17.2)
which is the Poisson’s equation for homogeneous medium.
STATIC ELECTRIC FIELDS 87
2V 0 …..(1.17.3)
which is the Laplace’s equation for homogeneous charge free medium.
We know
V V V
V ax ay az
x y z
2V 2V 2V
V 2V 0 …..(1.17.4)
x 2 y 2 z 2
Which is Laplace’s equation in rectangular co-ordinate system.
where 2 is Laplacian operator
In cylindrical co-ordinate system is
1 V 1 2V 2V
2V 0 …..(1.17.5)
2 2 z 2
In spherical co-ordinate system
1 V 1 V 1 2V
2V 2 r 2 2 sin 2 2 0 …..(1.17.6)
r r r r sin r sin 2
Problem 1.36
Write Laplace’s equation in rectangular co-ordinates for two parallel planes of infinite
extent in the X and Y directions and separated by a distance ‘d’ in the Z-direction.
Determine the potential distribution and electric field strength in the region between the
planes.
Solution
2V = 0
2V 2V 2V
0
x 2 y 2 z 2
since the potential is constant in X and Y directions
V V 2V 2V
0
x y x 2 y 2
2V Fig. 1.40
0
z 2
88 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
V
A
z
V = Az + B
At Z = 0 V = V1
V1 = 0 +B
At Z = d V = V2
V2 = Ad +B
V2 = Ad+V1
V2 V1
A
d
V2 V1
The Potential distribution is V z V1
d
The Electric field strength is
V V V V V
E V = az 2 1 az 1 2 az
z d d
The parameter of the capacitor i.e., ‘capacitance’ is Fig. 1.41 Two conductors
connected to V
defined as the ratio of charge on one of the conductors to
the potential difference between two conductors.
STATIC ELECTRIC FIELDS 89
Q
D ds E ds
C S
S
…..(1.18.1)
V E dL E dL
1.18.1 Parallel Plate Capacitor
Consider two conductors whose area as ‘A’ and are separated by a distance ‘d’ as shown
in Fig.1.42.
We know the electric field intensity E between
s
parallel plate capacitors in free space as E an
0
s
But from the Fig. 1.41 E (ax )
E will be in opposite direction of x-axis
Q Fig. 1.42 Parallel plate
Q = s. A s
A capacitor
where dL dx ax
d
Q
V ax .dx ax
0
A
d
Q Qd
V dx
0
A A
Q A
C …..(1.18.2)
V d
Energy stored in the parallel plate capacitor is
90 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
1
2 v
WE E E dv
1
WE
2v
s a x s a x dv
1 s2
2 v 2
WE dv
1 s2 1 s2 s2 Ad
2 V
WE dv A d
2 2
Q
Replace s by
A
1 Q 2 Ad 1 Q 2 1
WE VQ …..(1.18.3)
2 A2 2 C 2
*Problem 1.37
Calculate the capacitance of a parallel plate capacitor with a dielectric, mica filled
between plates. r of mica is 6. The plates of the capacitor are square in shape with
0.254 cm side. Separation between the two plates is 0.254 cm.
Solution
A
We have C
d
Here 0r 8.854 1012 6
Q D ds where D D a and ds d dz a
2 L
Q D d dz 2 D L 2 E L
0 z 0
Fig. 1.43 Co-axial capacitor
Q Q
i.e., E E a
2 L 2 L
To find the capacitance of co-axial capacitor. We need to find the potential difference
between the two cylinders.
a
V E.dl where dl d a
b
a
V E .d a
b
a
Q
d
b
2 L
Q b
V ln
2 L a
Q 2 L
C …..(1.18.4)
V ln b
a
Which is the expression for Coaxial capacitance.
D.d s
Q
S
where D D r a r ;
ds r 2 sin d d ar
2
Q d r 2 sin Dr d Fig. 1.44 Spherical capacitor
0 0
Q
Dr
4 r 2
Q
Er
4 r 2
Q
E ar
4 r 2
To find the capacitance of spherical capacitor. We need to find the potential difference
between the two spheres.
a
V E dl
b
where dl dr ar
a
Q 1
4 b r 2
V dr
Q 1 1
4 a b
Q 4
C …..(1.18.5)
V 1 1
a b
Which is the expression for Spherical capacitance.
94 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
A ds A dL
S L
The left side of the equation is the divergence of the Electric Current Density J .
This is a measure of whether current is flowing into a volume (i.e., the divergence
of J is positive if more current leaves the volume than enters).
Recall that current is the flow of electric charge. So if the divergence of J is
positive, then more charge is exiting than entering the specified volume. If charge
is exiting, then the amount of charge within the volume must be decreasing. This is
exactly what the right side is a measure of - how much electric charge is
accumulating or leaving in a volume. Hence, the continuity equation is about
continuity - if there is a net electric current is flowing out of a region, then the
charge in that region must be decreasing. If there is more electric current flowing
into a given volume than exiting, then the amount of electric charge must be
increasing.
18. Write the expression for energy density in electrostatic field.
1
Ans. wE E2
2
19. Write down the expression for capacitance between two parallel plates.
A
Ans. C
d
20. What is meant by displacement current?
Ans. Displacement current is the current flowing through the capacitor.
3. The intensity of the field due to a line charge pL1 at a distance r1 = 1 cm away from
it is E1 = 1 V/m. What is the intensity E2 of the field of the line charge pL2 = 4 at a
distance r2 = 2 cm from it?
(a) E2 = 1 V/m (b) E2 = 4 V/m
(c) E2 = 2 V/m (d) E2 = ½ V/m
4. Charge Q is uniformly distributed in a sphere of radius a1. How is the charge
density going to change if this same charge is now occupying a sphere of radius
a2 = a1/4 ?
(a) It will increase 4 times (b) It will increase 64 times
(c) It will increase 16 times (d) It will increase 2 times
5. A line charge pL = 5 × 10–3 C/m is located at (x, y) = (0, 0), and is along the z-axis.
Calculate the surface charge density ps (ps > 0) and the location xp (xp > 0) of an
infinite planar charge distributed on the plane at x = xp, so that the total field at the
point P (0. 5 × 10–3, 0) m, is zero.
(a) s 1/(2 ) C/m 2, x p 5 103 m (b) s 1/(2 ) C/m 2 , x p
6. The volume charge density associated with the electric displacement vector in
spherical coordinates sin sin ar cos sin a cos a is
(a) 0 (b) 1
(c) Not compatible (d) sin
7. The divergence theorem
(a) Relates a line integral to a surface integral
(b) Holds for specific vector fields only
(c) Works only for open surfaces
(d) Relates a surface integral to a volume integral
8. The flux of a vector quantity crossing a closed surface
(a) is always zero
(b) is related to the quantity’s component normal to the surface
(c) is related to the quantity’s component tangential to the surface
(d) is not related in any way to the divergence of that vector quantity
98 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
9. The flux produced by a given set of fixed charges enclosed in a given closed region
is
(a) Dependent on the surface shape of the region, but not the volume
(b) Dependent on the total volume of the region, but not the surface shape
(c) Dependent on the ratio of volume to surface area of the region
(d) Not dependent on any of these as long as the charges are inside the region
10. Consider charges placed inside a closed hemisphere. Consider the flux due to these
charges through the curved regions (Flux A) and through the flat region (Flux B)
(a) Flux A = Flux B
(b) Flux A > Flux B
(c) Flux A < Flux B
(d) Not enough information to decide the relation between Flux A and Flux B
11. An electron (qe = 1.602 × 10–19 C) leaves the cathode of a cathode ray tube (CRT)
and travels in a uniform electrostatic field toward the anode, which is at a potential
Va = 500 V with respect to the cathode. What is the work W done by the
electrostatic field involved in moving the electron from the cathode to the anode?
(a) W 5kJ (b) W 8 1019 J
(c) W 8 1017 J (d) W 5J
12. In the previous question, what is the electric field strength E E if the distance
between the cathode and the anode is 10 cm?
(a) E 5 V/m (b) E 500 V/m
(c) E 50 V/m (d) E 5 kV/m
13. The electrostatic potential due to a point charge Q1 at a distance r1 = 1 cm away
from it is V1 = 1 V. What is the potential V2 of a charge Q2 = 4Q1 at a distance
r2 = 2 cm from it?
(a) V2 0.5 V (b) V2 1 V
(c) V2 4 V (d) V2 2 V
14. The electrostatic potential due to a dipole p1 = p1a2 at a distance r1 = 1 cm away
from it along the z-axis, is V1 = 1 V. What is the potential V2 of a dipole
p2 = 4p1aZ at a distance r2 = 2 cm from it along the z-axis?
(a) V2 0.5 V (b) V2 1 V
(c) V2 4 V (d) V2 2 V
STATIC ELECTRIC FIELDS 99
2 103
15. The electrostatic potential V x V. Where x is measured in meters and
0
0 is the permittivity of vacuum, exists in a region of space (vacuum) in the shape
of a parallelogram of size 10 × 10 × 1 cm. What is the electrostatic energy
WE stored in this region?
(a) WE 2 1010 J (b) WE 1 109 J
(c) WE 4 1010 J (d) WE 3 109 J
16. Which statement is not true?
(a) The static electric field in a conductor is zero
(b) The conductor surface is equipotential
(c) Zero tangential electric field on the surface of a conductor leads to zero
potential difference between points on the surface
(d) The normally directed electrical field on the surface of a conductor is zero
17. The “skin” effect results in
(a) Current flowing in the entire volume as frequency increases
(b) Current flowing only near the surface as frequency increases
(c) Current flowing only near the surface as frequency decreases
(d) Current flowing near the surface at any frequency
18. As frequency increases, skin effect results in
(a) Decreased resistance
(b) Increased resistance
(c) No change in resistance
(d) Increase or decrease depending on material properties.
19. In a parallel-plate capacitor, the charge on the plates is C. What is the electric flux
density magnitude D, if the area of each plate is A = 10–4 m2. Assume uniform field
distribution.
(a) D 105 C/m 2 (b) D 105 /0 C/m 2
(c) D 105 0 C/m 2 (d) D 1013 C/m 2
20. For the capacitor in Previous question, find the voltage between its plates, provided
its capacitance is C = 10 pF.
(a) V 885V (b) V 0 V
(c) V 100 V (d) V 105 V
100 BASICS OF ELECTROMAGNETICS AND TRANSMISSION LINES
21. The capacitor in above Q no. 19 and 20 uses dielectric of permittivity 0 . The
maximum allowable field intensity (dielectric strength) of this dielectric is
Eds = 3 MV/m. (If E > Eds, the material breaks down.) What is the maximum
voltage Vmax, up to which the capacitor can operate safely (its breakdown voltage)?
(a) Vmax 885V (b) Vmax 1000 V
(c) Vmax 3 106 V (d) Vmax 265V
22. A coaxial capacitor whose cross-section is shown in the figure below has a central
conductor of radius r1 and an outer conductor of radius r3. The region between the
two conductors consists of two regions: (i) the region r1 < p < r2 has a relative
permittivity of r1 2 and (ii) the region r2 < p < r3 has a relative permittivity of
r 2 1 . The radius r2 is such that r2 r1 e2 and r3 r2 e where e 2.71 .
Answers
Exercise Questions
1. State the Coulomb’s law in SI units and indicate the parameters used in the equations
with the aid of a diagram.
2. State Gauss’s law. Using divergence theorem and Gauss’s law, relate the density D
to the volume charge density v .
3. Explain the following terms:
(a) Homogeneous and isotropic medium and
(b) Line, surface and volume charge distributions.
4. State and Prove Gauss’s law. List the limitations of Gauss’s law.
5. Express Gauss’s law in both integral and differential forms. Discuss the salient
features of Gauss’s law.
6. Derive Poisson’s and Laplace’s equations starting from Gauss’s law.
7. Using Gauss’s law derive expressions for electric field intensity and electric flux
density due to an infinite sheet of conductor of charge density C/m.
8. Find the force on a charge of –100 mC located at P(2, 0, 5) in free space due to
another charge 300 C located at Q(1, 2, 3).
9. Find the force on a 100 C charge at(0, 0, 3) m, if four like charges of 20 C are
located on X and Y axes at 4 m.
10. Derive an expression for the electric field intensity due to a finite length line charge
along the Z-axis at an arbitrary point Q(x, y, z).
11. A point charge of 15 nC is situated at the origin and another point charge of –12 nC
is located at the point (3, 3, 3) m. Find E and V at the point(0, –3, –3).
12. Obtain the expressions for the field and the potential due to a small Electric dipole
oriented along Z-axis.
13. Define conductivity of a material. Explain the equation of continuity for time
varying fields.
14. As an example of the solution of Laplace’s equation, derive an expression for
capacitance of a parallel plate capacitors.
15. In a certain region J 3r 2 cos ar r 2 sin a A m , find the current crossing the
surface defined by 30o ,0 2 ,0 r 2 m .