ELECTROSTATICS

Introduction
In the study of electrostatics, the charge, the source of electric field, is assumed to be
fixed in space and remains constant in time.
The aim of this part of study is to determine the force exerted by the charge(s) on other
charge(s), the electric field intensity, the potential distribution and the energy stored by a
capacitor.
There are numerous practical applications of electrostatic field. The static electric field is
used, for example, to accelerate a charged particle. The use of electrostatic field in the
design of an oscilloscope and an ink-jet printer is widely familiar.

1.1

Coulombs law

In an experiment, French army engineer Colonel Charles Coulomb observed that there
exists a force between two charged bodies. He showed that the electric force between two
point charges is
(a) Directly proportional to the product of the charges,
(b) Inversely proportional to the square of the distance between the charges,
(c) Directed along the line joining the charges,
(d) The magnitude of the force depends on the medium and
(e) Repulsive for charges of identical (like) sign, but attractive for charges of
opposite (unlike) sign.
On the basis of the experimental result, Coulomb formulated the force, known as
Coulombs electrostatic force, between two point charges in a boundless free space as
1 Qq
FQ =
a R Newton (N)
(1)
4 0 R 2
where Q and q are two point charges in Coulomb (C), R=|R| is the distance between them
in meter (m) and 0 is the permittivity of the free space (k~bvbi f`bhvMZv) in which the
charges are situated. If FQ is the force exerted on q (a test charge) by Q (source charge),
the unit vector aR will point in the direction from Q to q as shown in Fig. 2(a).

aR

(x', y', z')

Q
r

R
q (x, y, z)
Q(x', y', z')
(a)

z
R=r-r
r

q
(x, y, z)
y

(b)

Fig. 2
Here R=(x- x')ax+(y- y')ay+(z- z')az and
R ( x x )a x + ( y y )a y + ( z z )a z
aR = =
R
( x x ) 2 + ( y y ) 2 + ( z z ) 2
If r and r are respectively, the position vectors of Q and q as shown in Fig. 2(b), (1) can
be expressed as

Qq
1 QqR
1 Qq (r r )
aR =
=
2
3
4 0 R
4 0 R
4 0 r r 3
1

FQ =

(2)

Here r=xax + yay + zaz, r= x'ax + y'ay + z'az and R=r r=(x- x') ax + (y- y') ay + (zr r
z')az and
is the unit vector aR in the direction of the vector R=r r.
r r
Coulombs law can be used over the wide range from 10-13 cm to many km. Even it
correctly describes,
1) The electric force that bind the electrons of an atom to its nucleus,
2) The forces that bind atoms together to form molecules, and
3) The forces that bind atoms or molecules together to form solids or liquids.
However, the coulombs law fails to explain the stability of the nucleus. The nuclear
force which is larger than the electric force described by the coulombs law holds the
proton of the nucleus together.

Superposition principle in Electrostatic

1.2

The force on a test charge q is a linear function of the amount of charge Q of the source.
The principle of superposition can, therefore, be applied in a system consisting of three or
more discrete charges to determine the electric force exerted on any one of the charges by
the others. If Fi be the force exerted on a test charge q by the ith source charge Qi, and
aRi is the unit vector acts in the direction from Qi to q as shown in Fig. 3, then the
resultant force on q due to several point charges Qi (where i=1,2,3.n) is given by:
n

F = Fi =
i =1

Qiq
1
a Ri =

2
4 0 i =1 R i
4 0
1

R1

Q1

Q i qR i
1
=

3
4 0
Ri
i =1

i =1

Q i q(r ri)
r ri

(3)

aR2
q

Q2

aR1

R2

Fig. 3
Here, in (3) Ri is the vector between the ith source charge and q, r is the position vector
of q on which the force is to be determined and ri is the position vector of the ith source
Qi.
In the above discussion, only the point charges are considered, where all the charges are
assumed to be concentrated at a point having zero volume. However, the charges may be
distributed uniformly along a line, or on a surface or in a volume.
When the charges are spread over a finite volume with Volume charge density v (Cm-3),
then it is called volume charge; where the total charge Q contained within a volume v is
given by:
Q = v dv C , where the integral is to be performed over the volume v.
v

If the charge is spread thinly over a surface of a sheet, it is said sheet charge. If the
Surface charge density at a given point on the surface of the sheet is s (Cm-2), the total
charge Q on the whole surface S is then given by:
Q = s ds C , where the integral is to be performed over the surface S.
s

When the charge is distributed along a line with a Line charge density l (Cm-1), it is
called line charge; the total charge Q in a total length L is given by:
Q = l dl C , where the integral is to be performed along the line L.
l

Direct use of equation (1) or (2) is not possible to determine the force exerted by a line or
surface or volume charge distribution as the source charge is not in the point form.
However, the superposition principle helped us to overcome this problem easily. In case
of line, surface or volume charge distribution, we may divide the entire line, surface or
volume into several numbers of small segments. The charge of the small segment of line
or surface or volume will behave like a point charge. The forces for the individual
segments are calculated using equation (1) or (2) and then add them vectorically to
determine the resultant. This is explained in the following in details.
For the line charge distribution shown in Fig. 4(a), if we consider a small segment dl then
l dl may be assumed as a point charge (l is the line charge density in C/m). The force
dF on q due to this small segment of line charge l dl C is therefore,

dF =

q l dl(r r )

N.
3
r r
The total force on q due to the whole line charge is then given by
l (r r )
q
F = dF =
dl N.
4 0 l r r 3
l

4 0

(4)

In Fig. 4(b), if s is uniform surface charge density in C/m2, the charge sds C on a small
surface ds, can be assumed as a point charge. The force dF on q for this point charge is
then given by
1 q s ds(r r )
dF =
.
3
4 0
r r
The total force on q due to the whole surface is therefore given by
s (r r )
q
ds N.
F = dF =
3

r
0
s
s

(5)

In Fig. 4(c), if v is uniform volume charge density in C/m3, the charge vdv C on an
elemental volume dv, can be assumed as a point charge. The force dF for this point
charge on q is then given by
1 q v dv(r r )
dF =
.
3
4 0
r r
The total force on q due to the whole volume is therefore given by
v (r r )
q
F = dF =
dv N.

4 0 v r r 3
v

(6)

In the above equations, r is the position vector of the point P where the test charge q is
assumed to be placed and r is the position vector of the segmented element as shown in
Fig. 4.

Fig. 4

Electric Field (Zwor ) and Electric Field Intensity (Zwor i

cvej)

1.3

If a point charge q (test charge) is brought into the vicinity of a fixed charge Q, the
charge q experiences a force. It is then said that Q has a field around it which can exert a
force on another charge q. Thus, the space surrounding a charge where forces act is
defined as the electric field. And the electric field intensity/electric field strength/electric
field E at some point in the field is defined as the force on a unit positive charge at that
point. To determine E at a point P, we have to place a very small positive charge q at that
point and measure the force F acting on it. Then
F
E = lim . As because q also creates its own electric field and perturb/alter the initial
q 0 q
field, its magnitude must be as small as possible in order to minimize the distortion.
The unit of E is N/C or more usual V/m.
For a single point charge Q
1 Q
1 Q(r r )
E=
a
or
E
=
V/m as obtained from (1) or (2).
R
4 0 R 2
4 0 r r 3
The electric fields E due to several point charges, line, surface, and volume charge
distribution are, respectively, given by the following eqs. (7), (8), (9), and (10):
n

E = Ei =
i =1

E = dE =
l

E = dE =
s

E = dE =

4 0
1

4 0
1
4 0
1

Q i (r ri)

r ri
l (r r )

i =1

r r

s (r r )
r r

v (r r )

(7)

dl

(8)

ds

(9)

dv
(10)
3
r r
The electric field at a given point around a charge represents graphically by the lines
along which the test charge moves when placed in the field. These line are known as
electric field lines (Zwor ejiLv). It points outward from the charge if it is positive,
v

4 0

otherwise it is inward as shown in Fig.5. The density of field lines, however, is

proportional to the magnitude of the electric field.

Q+

Q-

Q+ >> Q-

Fig. 5
The electric field lines have the following properties
The tangent to the lines at any point gives the direction of the E-field at that point.
Lines start on positive charges and finish on negative ones.
The density of lines gives an indication of the field strength at a given point.

1.4

Electric potential or Voltage (Zwor wefe)

The energy required to move a unit positive charge from one point to another, is defined
as the potential difference between the points.
If a positive test charge q is placed in an electric field E, there will be a force F=qE on
the charge q that acts in the direction of E. Assume that under this force, the charge q
moves a differential distance dl. The incremental energy expended by the electric field,
or simply the amount of work done by E is thus given by dWe= qE.dl where the subscript
e signifies that the work is done by the electric field. However, if an external force Fext =
-F= - qE acts on the charge q, it will move against the E field; the differential work done
by the external force is given as dW= - qE.dl. Assuming the test positive charge q is
moved from point b to point a against E, then by the definition of the potential difference
between the points a and b, Vab, is given by
a
a
dW
Vab = Va Vb =
= E dl volts.
(11)
q
b
b
The ve sign in the above equation signifies that the work has to be done by the external
source against the field in moving the charge q from b to a. To determine the potential
experimentally, the test charge q should be considered sufficiently small so that it does
not perturb the original field distribution.
In eq. (11), dl, in various coordinate systems, is given by,
dl=dx ax+ dy ay+ dz az in Cartesian coordinate system,
dl=d a+ d a+ dz az in Cylindrical coordinate system and
dl=dr ar+ rd a+ r sin d a in Spherical coordinate system.
In moving the test charge q in a closed path, the total work done and hence the potential
difference must be zero. In other words,
E dl = 0 ( KVL) which on using the Stocks theorem takes the following form
E = 0 . It states that the static electric field is irrotational or conservative.

1.4.1 Potential due to a point charge

1

Q(r r )

where r= x'ax +
4 0 r r 3
y'ay + z'az and r=xax + yay + zaz are the position vectors of the location of Q and the
point where the field is to be determined.
Now, r- r=(x- x')ax + (y- y')ay + (z- z')az and | r- r|=[(x- x')2+ (y- y')2+ (z- z')2]1/2.
Following eq. (11) with dl=dx ax+ dy ay+ dz az the potential difference between points a
a
W
and b is Vab = Va Vb = ab = E dl
q
b
Q [(x - x )dx + (y - y )dy + (z - z )dz ]
Now E dl =
4 0 (x - x )2 + (y - y )2 + (z - z )2 3/2
The electric field for a point charge Q is given by E =

Vab =

Q
4 0

[(x - x )dx + (y - y)dy + (z - z)dz] =

[(x - x ) + (y - y) + (z - z )
2

2 3/2

Q
4 0

1
1
)

ra r rb r

where ra and rb are respectively, the position vectors of points a and b. If point b is
assumed to be at infinity ie. rb=, then the potential of point a is called the absolute
potential and is given by
Q
Q
Q
Va =
or, more generally V =
.=
. .
(12)
4 0 ra r
4 0 r r
4 0 R
where R is the distance of the observation point from the charge Q. In eq. (12), V is the
absolute or simply potential at a point of position vector r.
Following the above eq. and the superposition theorem, the (absolute) potential at a point,
due to a group of point charges, line charge, surface charge and volume charge are given
by eqs. (13)-(16), respectively:
n
Qi
1
V=
(13)
i =1 4 0 r ri
V=
V=
V=

1
4 0
1
4 0
1
4 0

l dl

r r

(14)

s ds

r r

(15)

v dv

r r

(16)

In the above equations, ri is the position vector of the ith source charge, r is the position
vector of the observation point where the potential is to be determined, and r is the
position vector of the elemental source.
Equation (12) shows that the potential remains unchanged on a circle or a surface of
constant radius. A line (surface) on which the potential is same is known as an
equipotential line (surface). Thus for a point charge, the equipotential surfaces are
concentric spheres and those for a line charge are concentric cylinders as shown in Fig. 6.

Fig. 6.
1.4.2 Path-independence of potential

Fig. 7
Let us consider few concentric circles around a point charge at the origin as shown in Fig.
7. The magnitude of electric field E at each and every point on each circle is constant. As
the electric force lines/field lines of a point charge is a radial one, work is only done for
movement along the line joining the two charges (one of the charges is the source of
electric field and the other is the test charge used to measure the potential). No work is
done, i.e., dW=0 for any displacement of the test charge along the circular path which is
normal to the field lines, because qE.dl=0. In the above figure, therefore, no work is done
along the arc segments AB, CD, EF and GH. Hence dW for path ABCDEFGH is
dW(BC)+ dW (DE)+ dW (FG) = dW (AH). Hence dW is independent of the path taken
in moving between two positions.

1.5

Electric Field as gradient of Potential

In case there are two points which are separated by an small distance dl, then the work
done by an external force in moving a unit positive charge from one point of potential V
to another point of potential V+dV is
dW=(V+dV)-V= - E dl

dV= - E dl. As V is a function of position (say x,

y and z),
V
V
V
dx +
dy +
dz = E dl
x
y
z
V
V
V
(
ax +
ay +
a z ) (dxa x + dya y + dza z ) = E dl
x
y
z
or, V dl = E dl

where =
a x + a y + a z in Cartesian coordinate system,
x
y
z

=
a +
a + a z in Cylindrical coordinate system and

ar +
a +
a in Spherical coordinate system.
r
r
rsin
Therefore E= - V .
=

1.6

(17)

Equipotential surfaces and E-Field lines

For a +ve point charge, the lines of force/field point radially outwards and the
equipotential lines form a series of concentric circles as shown in Fig. 8. It is easy to
prove that at all points the two types of lines (force line and equipotential line) are normal
to each other.

Fig. 8
In a direction tangential along the equipotential surface/line there can be no change in V.
Therefore there can be no component of E tangential to the surface (as E=-
V) and hence
the only component of E must be normal to the surface.
The properties of equipotential surfaces can be summarized as follows:
(i) The electric field lines are perpendicular to the equipotentials and point from higher to
lower potentials.
(ii) By symmetry, the equipotential surfaces produced by a point charge form a family of
concentric spheres, and for constant electric field, a family of planes perpendicular to the
field lines.
(iii) The tangential component of the electric field along the equipotential surface is zero,
otherwise non-vanishing work would be done to move a charge from one point on the surface
to the other.
(iv) No work is required to move a charged particle along an equipotential surface.
Conceptually, a field line is a directed curve in the electric field such that the forward
drawn tangent at any point on the curve has the direction of electric field E at that point.
Therefore, if dl is a small segment on this curve, then dl=kEE where kE is a constant. By
equating the values of kE for different components of dl and E, one finds the following
differential equations in different coordinate systems:
dx dy dz
=
=
in Cartesian coordinate system,
Ex Ey Ez
d d dz
=
=
in Cylindrical coordinate system and
(18)
E
E
Ez
dr rd rsind
=
=
in Spherical coordinate system.
Er
E
E
The field lines can then be obtained by integrating the above differential equations.
Similarly, the equation of equipotential line/surface can be obtained by equating the
equation of potential to a constant, i. e, V=kv.

1.7

Electric dipole

Many physical systems are electrically neutral but still produce an electric field and are
affected when placed in an electric field. This arises because the positive and negative
charges of the molecules are physically separated when they are brought to an electric
field. Such molecules are said to be polarized. This leads to the concept of an electric
dipole which means a pair of equal charges of opposite polarities with very small
separation between them.
Two equal point charges Q and Q when separated by a small distance s, thus results in
an electric dipole or simply dipole. In an ideal dipole s is very small compared to the
distances to any other charges and also to any points where we wish to find the resultant
electric potential or E-field.
az
az
sin
r-d/2
P
ar
Q
r
1
a
d0

cos

r+d/2
d/2=d/2 az
az= cos ar - sin a
p=Qd=Qd az=p az

-Q

Let us find potential V and electric field E at point P due to the dipole. We can determine
E first and then V using (11) with radial distance of point b as infinity. The 2nd method is
to determine V first using (13) and then E using (17). The two methods are described in
the following.
(a) Determination of E first and then V
Using (7), E at P is given by
1 2 Q i (r ri)
Q (r d / 2) (r + d / 2)

E=
=

4 0 i =1 r ri 3
4 0 r d / 2 3 r + d / 2 3

Now, r d / 2

= [(r d / 2) (r d / 2)]

3 / 2

3 / 2

= r2 r d + d2 / 4

3 / 2

3 / 2

r d
rdcos
r -3 1 2
= r -3 1
r
r 2

As the 2nd term in above eq. is small in comparison to 1 for an ideal dipole, using Binomal
expansion we have r d / 2

3
3dcos
r -3 1 +
. Similarly, we have r + d / 2

2r

3dcos
r -3 1
.
2r

Q (r d / 2) (r + d / 2)

E =

4 0 r d / 2 3 r + d / 2 3

Qr -3
4 0

3dcos

3dcos
1 + 2r (r d / 2) 1 2r (r + d / 2)

Q 3dcos
Q
3dcos

ra r da z =
r d =
[3dcosa r da z ]
3
3

4 0 r r
4 0 r r
4 0 r
Q
=
[3dcosa r d(cosa r sina )] = Q 3 [2dcosa r + dsina ]
3
4 0 r
4 0 r
Q

Qd
[2cosa r + sina ] = p 3 [2cosa r + sina ]
3
4 0 r
4 0 r

(19)
Now the potential at P, which is r distance away from the origin of coordinate system, is
P
P
p

V = E dL =
[2cosa r + sina ] (dr a r + rda + r sin da )
3

4 0 r
If we consider the path of integration along the radial path only
p ar
2pcos
pcos
dr =
=
3
2
4 0 r
4 0 r 2
4 0 r
r

V =

(20)
(b) Determination of V first and then E
Using (13), V at P is given by
2

Qi
1
Q
1
1
V=
=

4 0 r d / 2 r + d / 2
i =1 4 0 r ri
Now, r d / 2

= [(r d / 2) (r d / 2)]

1 / 2

1 / 2

= r2 r d + d2 / 4

1 / 2

1 / 2

r d
dcos
r 1 2
= r -1 1
r
r

As the 2nd term in above eq. is small in comparison to 1 for an ideal dipole, using Binomal
-1

expansion we have r d / 2

1
dcos
r -1 1 +
. Similarly, we have r + d / 2

2r

dcos
r -1 1
.
2r

Qr 1 dcos dcos
Q
1
1
V =

1 +
1

=
4 0 r d / 2 r + d / 2 4 0
2r
2r
p ar
Qdcos
pcos
=
=
=
2
2
4 0 r
4 0 r
4 0 r 2
Using (17) we have

pcos

E = V = - a r +
a +
a
r
rsin 4 0 r 2
r
=

2pcos
psin
p
ar +
a =
[2cosa r + sina ]
3
3
4 0 r
4 0 r
4 0 r 3

In contrast to a single point charge the potential falls off as r-2 (not r-1). The product Qd is
defined as the electric dipole moment, symbol p. This is actually a vector p=Qd where d
is defined as pointing from -Q to +Q.
The total field is
p (4 cos 2 + sin 2 )1 / 2 p(1 + 3 cos 2 )1 / 2
2
2
E = E R + E =
=
4 0 r 3
4 0 r 3
In contrast to a single point charge E falls off as r-3 for an electric dipole (not r-2). It is
very important to note that the potential on the bisecting plane for which =900 is zero;
and the electric field acts perpendicularly on this surface.

1.7.1 E-field and equipotential lines for an electric dipole

The equation for an equipotential surface for any charge distribution is obtained by
setting the expression of V to equal a constant. Since r and are only variable in (20), we
have

r = kV cos where kV is a constant.

By plotting r versus for various values of kV we draw the equipotential lines (dashed
lines) in above figure. On the other hand the electric field lines represent the direction of
E in space. We, therefore get from (19)
p
E=
[2cosa r + sina ] = k E (2 cos a r + sin a )
4 0 r 3
Therefore using (18) we have
2 cos d
dr
rd
dr 2 cos d
dr
=

=
2 cos sin
r
sin
r
sin
2
2
ln r + const1 = ln sin + const 2 r = C E sin
where CE is a constant. By plotting r versus for various values of CE we draw the field
lines (solid lines) in above figure.
Example 1.1 A line charge density l(x) is given by l(x)=3x2 Cm-1. Calculate the total
charge contained between the points x=0 and x=L.
Solution:
L

The total charge Q= 1 ( x) dx = 3x 2 dx = x 3

0

L
0

= L3 C

Example 1.2 A surface charge density s (x,y) is given by s (x,y)=3x2+4y2-xy Cm-2.

Calculate the total charge contained within the area bounded by x=0+a, y=0+a.
Solution:
a a

a a

The total charge Q = s ( x, y ) dxdy = (3 x 2 + 4 y 2 xy ) dxdy

0 0

= ( x 3 + 4 xy 2 x 2 y / 2)
0

0 0

x =a

dy = (a 3 + 4ay 2 a 2 y / 2)dy = (a 3 y + 4ay 3 / 3 a 2 y 2 / 4)

x =0

a
0

= (a 4 + 4a 4 / 3 a 4 / 4) = 25a 4 / 12 C
Example 1.3 Charges of +Q, +5Q, -Q and 2Q are placed on x-y plane at the coordinates (-a,+a), (+a,+a), (+a,-a) and (-a,-a) respectively. Calculate the E-field at the
origin.
Solution:
Q
y
5Q
r1'

r2'
x

r4'

r3'

- 2Q

-Q

Let r1 ~ r4 are the position vectors of the charges. And the position vector of the
observation point where E is to be determined is r=0.
Now E = E 1 + E 2 + E 3 + E 4 where
Q(r r3 )
Q(r r1 )
5Q(r r2 )
2Q(r r4 )
=
and
=
E1 =
,
E
,
E
=
E
2
3
4
3
3
3
3
4 0 r r1
4 0 r r2
4 0 r r3
4 0 r r4
r1 = aa x + aa y , r2 = aa x + aa y , r3 = aa x aa y

and r r1 = r r2 = r r3 = r r4 =
3

a2 + a2

and r4 = aa x aa y
3

= (2a 2 ) 3 / 2

Q (r r1 ) 5(r r2 ) (r r3 ) 2(r r4 )

3
4 0 r r1 3 r r2 3 r r3 3

r
4

Q
=
aa x aa y 5aa x 5aa y + aa x aa y 2aa x 2aa y
4 0 (2a 2 ) 3 / 2

E=

[
[

Q
5aa x 9aa y V/m
4 0 (2a 2 ) 3 / 2
Example 1.4 Two point charges -q and q/2 are placed at the origin and at (a, 0, 0)
respectively. Determine the position on x-axis where E vanishes.
Solution:
=

-q

q/2

(0, 0, 0)

(a, 0, 0)

P(x, 0, 0)

Let P be the position where the electric field E vanishes. Then the field due to q and q/2
are, respectively, given by:
q
q/2
E1 =
a x and E 2 =
ax .
2
4 0 x
4 0 ( x a) 2
E = E1 + E 2 = 0

1
1
+
= 0 x 2 = 2( x a) 2 x 2 4ax + 2a 2 = 0
2
2
x
2( x a )

or, x = a(2 2 ) = 3.41a or 0.59a

Among the two, the right answer is x=3.41a, because E will act only in opposite
directions for considered charge configuration if x>a.
Example 1.5 A line charge having a constant charge density l Cm-1 extends along the
y-axis from - to +. Calculate the E-field at a distance a along the positive x-axis. Also
calculate the potential V.
Solution:
y
y=+
dl=dy'
y'

r=y' ay

r- r= a ax - y' ay

r=a ax

x=a

r r = ( a 2 + y 2 )1 / 2

y= -
Now the charge in the small segment dl=dy' is ldy C. The field for this small charge is
+
+
l (aa x y a y )
l dy (r r )
l (r r )

E
dE =

=
d
y
=
3
4 r r 3
4 0 (a 2 + y 2 ) 3 / 2 dy
4 0 r r
0
+

l aa x

l y a y

dy
2

2 3/ 2
(a 2 + y 2 ) 3 / 2
4 0 ( a + y )
Assume y'=a tan dy'=a sec2 d and (a2+y'2)3/2=a3sec3. The limits of the
integration become = -/2 for
+
+
y'= - and =/2 for y'=.
y a y

aa x

E= l 2
d
y

d
y

(a 2 + y2 )3 / 2
4 0 (a + y 2 ) 3 / 2

= l
4 0

dy

+ / 2 2
+ / 2 a 2 sec 2 a x

a sec 2 tan a y
d
d

3
3
3
3
a sec
/ 2 a sec

/ 2

+ / 2
sin a y

l + / 2 cos a x
l
/2
/2
d

d =
sin / 2 a x + cos / 2 a y

4 0 / 2 a
a
/ 2
4 0 a

l
a x V/m
2 0 a

To calculate V using the above equation of E, we have to express the location of the point
where the field has been calculated as variable which lefts E =

l
a x V/m
2 0 x

Therefore the potential of a point, whose perpendicular distance from the line charge is x,
with respect to another point of distance x0 is
x
x
x
l
l

x
V = E dl = (
a x ) (dxa x + dya y + dza z ) =
dx = l ln
2 0 x
2 0 x
2 0 x0
x0
x0
x0

ln x + l ln x0
2 0
2 0

If we consider the potential of the point at x0 as zero, the potential about a line charge is
given by V =

l
ln x .
2 0

Example 1.6 A finite line charge of length L carrying uniform line charge l Cm-1 lies
along the z-axis. Determine E and V at a point P on the plane bisecting the line charge.
Solution:
z
z=+L/2
dl=dz'
r- r= a - z' az r r = ( 2 + z 2 )1 / 2
P
r=z' az
z= - L/2
r=a
z= -L/2
Now the charge in the small segment dl=dz' is ldz' C. The field for this small charge is
+
+L / 2
l ( a z a z )
l dz (r r )
l (r r )

dE =

E
=
d
z
=
dz
3
2
2 3/ 2
4 r r 3

z
4
(
+
)

4 0 r r
0
L / 2
0
+L / 2

l a

+L / 2

l z a z
dz
2

2 3/ 2
L / 2 4 0 ( + z )
L / 2 4 0 ( + z )
Assume z'= tan dz'= sec2 d and ( 2+z'2)3/2= 3sec3. The limits of the
integration become = -tan-1 (L/2 ) for z'= - L/2 and = tan-1 (L/2 ) for z'=L/2.
+L / 2
a

l + L / 2
z a z
E=
dz
dz

2
2 3/ 2
2
2 3/ 2

4 0 L / 2 ( + z )
L / 2 ( + z )

2 3/ 2

dz

tan ( L / 2 )
tan ( L / 2 ) 2 sec 2 a
2 sec 2 tan a z
d

3
3
3
3

sec

1
tan 1 ( L / 2 ) sec

tan ( L / 2 )
tan 1 ( L / 2 )
tan 1 ( L / 2 )

cos a
l
sin a z
=
d

4 0 tan 1( L / 2 )

tan 1 ( L / 2 )

= l
4 0

l
tan ( L / 2 )
sin tan ( L / 2 ) a + cos
4 0
l
L/2
=
a V/m
2 0 ( L / 2) 2 + 2

E=

tan 1 ( L / 2 )
tan 1 ( L / 2 )

l
az =
sin(tan 1 ( L / 2 ))a
2
0

The potential of the point P is

l
L/2
a (da + da + dza z )
2
2
2

(
L
/
2
)
+

V = E dl =

l
L/2
d
2 0 ( L / 2) 2 + 2

Assume =(L/2) tan, therefore d= (L/2) sec2 d and (L/2)2 + 2=(L/2)2 sec2.
l 1 1 cos
l
l
L/2
2 0 ( L / 2) 2 + 2 d = 2 0 cos ec d = 2 0 2 ln 1 + cos

1 ( L / 2) / ( L / 2) 2 + 2
= l
V = l ln
4 0 1 + ( L / 2) / ( L / 2) 2 + 2
4 0


l
=
4 0

( L / 2) 2 + 2 ( L / 2)

ln
2
2

( L / 2) + + ( L / 2)

( L / 2) 2 + 2 + ( L / 2)
V
ln
( L / 2) 2 + 2 ( L / 2)

Alternate Solution:
We can solve for V first and then determine E from E= -
V.
For the small segment dz', potential dV is
l dz
l L / 2
dz
V =
dV =

4 0 L / 2 z 2 + 2
4 0 z 2 + 2

Assume z'= tan, therefore dz'= sec2 d and z'2 + 2= 2 sec2.

dz
z 2 + 2

= sec d =

2
2
2
2
1 1 + sin 1 1 + z / z + 1 z + + z
ln
ln
ln
=

2 1 sin 2 1 z / z 2 + 2 2 z 2 + 2 z

L/2

( L / 2) 2 + 2 + ( L / 2)
l
l 1 z 2 + 2 + z

ln
=
ln
V =

4 0 2 z 2 + 2 z
4 0 ( L / 2) 2 + 2 ( L / 2)

L / 2

and

l ( L / 2) 2 + 2 + ( L / 2)

a
E = -V = -( a +
a + a z )V =
ln

z
4 0 ( L / 2) 2 + 2 ( L / 2)

l
4 0
l
=
2 0

2
2
2
2
ln ( L / 2) + + ( L / 2) ln ( L / 2) + ( L / 2) a

L/2
a
( L / 2) 2 + 2

Example 1.7 A semi-infinite line charge having a constant charge density l Cm-1
extends along the y-axis from 0 to +. Calculate the E-field and V at a distance a along
the positive x-axis.
Solution:

y
y=+
dl=dy'
r=y' ay

r- r= a ax - y' ay

y'=0

r r = ( a 2 + y 2 )1 / 2

x
r= a ax

The problem is similar to that given in example 4. In this case, however, the lower limit
of the integration is different, which is =0 as y'=0.
+ / 2 2

a sec 2 tan a y
l + / 2 a 2 sec 2 a x
E =
d
d

3
3
3
3
4 0 0 a sec
a sec

0
+ / 2
sin a y

l + / 2 cos a x
=
d

4 0 0
a
a
0

l
l
/2
/2
=
sin 0 a x + cos 0 a y =
(a x a y ) V/m
4 0 a
4 0 a
Calculation of V:
In calculating V using V = E dl , we have to replace a in equation of E by x to make

the position of P on the x-axis generalized. Assume the unit charge is moved from
infinity to P along the x-axis. Therefore dl=dx ax.
a

a
l
l

V = E dl =
(a r a y ) dxa r =
dx = l ln x = l ln a
4 0 x
4 0 x
4 0
4 0

Example 1.8 A semi-infinite line extending from - to 0 along the z-axis carries a
uniform charge distribution of 100 nC/m. Find the electric field intensity at point P(0, 0,
2). If a charge of 1C is placed at P, calculate the force acting on it.
Solution: Let us consider an elemental charge ldz at z= z. The position vector, r, of the
observation point P and that of the charge position, r, are given by
r=z az, r= -z az. Therefore r-r=(z+z) az.
Now using (8), E is given by
0 dz
l
E= l
az =
az
2
4 0 ( z + z )
4 0 z
9 10 9 100 10 9
a z = 450a z V/m
2
The force on 1C charge at P is F = qE =
=

1 10 6 450a z = 450a z N

Example 1.9 A line charge extends from x=0 to x=L and has a density l(x)=ax Cm-1.
Calculate the E-field at a point b (>L) on the x-axis.
Solution:

r=x' ax

r- r= (bx') ax

x
dx' x=L
x=b
r=b ax
|r- r|=(b-x'). Therefore the field dE at x'=b due to the charge l(x')dx'=ax'dx' is
L
L
ax (b x )a x
ax dx (r r )
ax (r r )

=
=
dE =

E
d
x
dx
3
3
3

x
4
(
)

4 0 r r
4
r

r
0
0
0
0
x=0

L
ax a x
x a x
a

d
x
=
dx
2
2

4
4
(
b

x
)
(
b

x
)

0 0
0
0
Assume y=b-x'. Therefore x'=b-y and dx'=-dy. The limits of the integration become y= b
for x'= 0 and y=b-L for x'=L. Then
b L
b L
b L
bL

(b y )a x
( y b)a x
a
a
a 1
b
E=

dy =
dy =
dy 2 dy a x
2
2

4 0 b
4 0 b
4 0 b y
y
y
b y

b L

a
1
a bL
L
b L
=
ln
+
ax V / m
ln y b + b
a x =

4 0
yb
4 0
b
(b L)

Example 1.10
A line charge has density and extends along the x-axis from a
to +a. Find the electric potential at a point r on the x-axis (r>a). Use your result to find
the E-field at r.
Solution:
r=x' ax
r- r= (r-x') ax

x
x= -a x=0 dx' x=a
x=r
r=r ax
|r- r|=(r-x'). Therefore the field V at x'=r due to the charge dx' is
a

a
a
dx

dV =
V =
dx =
dx =
ln(r x )

4 0 r r
4
r

r
4
(
r

x
)
4

0
0
0
a
a
a

r +a
ln
Volts
4 0 r a

Following the above equation of potential on the x-axis at any position r where (x'=r>a) is
given as function of x as:

x+a
V=
ln
Volts
4 0 x a

E on the x - axis where x = r > a, E = V = ( a x + a y + a z )V

x
y
z
=

1
1

a x
4 0 x a x + a

4 0

1
1

a x
r a r +a
Example 1.11
Two infinite line charges, one having charge distribution l C/m
and the other -l C/m, are separated by a distance a. Determine the electric field intensity
Therefore E at r on the x - axis is E =

E at a point P, which is located on the line joining the charges and which is a distance
from the positive line charge, for<a and >a.
Solution:
l
-l
l
-l

a
x

a
a

a
P

Following the result of example 5, E for an infinitely long line charge l at a

perpendicular distance r from the charge is (l/20) a, and that for -l is (-l/20 (a))(- a) and (-l/20 (-a)) a, respectively, for <a and >a.
l
l
l a
Therefore E at P for r<a is E = (
+
)a =
a and for r>a
2 0 2 0 (a )
2 0 r (a )
l
l
l a

E=(
)a =
a
2 0 2 0 ( a)
2 0 ( a )
Example 1.12
Determine electric field due to a ring charge of radius b at a point
on its axis. Assume the ring is uniformly charged with density l C/m.
Solution:
z
z
r=z az

r- r= z az - b a
dl=bd

r=b a
The field E at P due to charge of small segment dl is
2
2
l b( za z ba )
l dl (r r )
l b(r r )
dE =

=
d
=
E

3
0 4 r r 3
0 4 0 ( z 2 + b 2 ) 3 / 2 d
4 0 r r
0

l b 2
zd
bd
a

a
2
z
2 3/ 2

2
2 3/ 2
4 0 0 ( z + b )
0 (z + b )

The above equation shows that the field due to a segment has two components, one acts
along the z-axis and the other acts in the -direction. If we consider a similar segment
just opposite side of the present segment, the r-components of these two segments will
cancel each other. Therefore, when the entire ring is under consideration, the field only
acts in the z-direction. This leaves
b 2
2 l bz
zd
Qz
E= l 2
az =
az =
az
2 3/ 2
2
2 3/ 2
4 0 0 ( z + b )
4 0 ( z + b )
4 0 ( z 2 + b 2 ) 3 / 2
Here Q is the total charge on the ring.

Example 1.13
A charged ring of radius a carries a uniform charge distribution.
Determine the potential and E at any point on the axis of the ring.
Solution: Let us consider an elemental charge ldl=lad at P(a, , 0). The position
vector, r, of the observation point P and that of the charge position, r, are given by
r=z az, r= a a. Therefore r-r=(z az -a a) and |r-r|=(z2+a2)1/2.

Using (14), we obtain

2
l ad
l a
1
V ( z) =
=

2
2 1/ 2
2
4 0 0 ( z + a )
2 0 ( z + a 2 )1 / 2

l a
l az
d
a =
az

2
2 1/ 2 z
dz 2 0 ( z + a )
2 0 ( z 2 + a 2 ) 3 / 2
Example 1.14
A plane circular sheet of radius b has a uniform charge density s
Cm-2. Calculate the E-field and V at a point P which is a perpendicular distance a from
the centre of the sheet.
Now E = -V = -

Solution:
z
P
r- r= a az - a
r=a az

r= a
Let us consider a small surface ds=dd on the circular sheet placed on the z=0 xyplane. The electric field at P due to the charge on the surface ds is then given by
=b = 2
s ( a + aa z )

ds(r r )
s (r r )
dE = s

E
=
ds
=
d d

2
2 3/ 2
4 0 r r

= 0 = 0 4 0 ( + a )
s 4 0 r r
However, due to symmetry the radial component of the field is zero. Therefore,
s (aa z ) 2
=b
s r =b
s aa z
a a z
0

E=
d
d =
d =

2
2 3/ 2
2
2 3/ 2
2
2 0 r =0 ( + a 2 ) 3 / 2
4 0 (r + a )
2 0 ( + a )

r =0
=
0

Assume =a tan
d=a sec2 d and (2+a2)3/2=a3sec3. The limits of the
integration become = 0 for = 0 and =tan-1 b/a for =b. Then
r =b

 = tan 1 b / a

= tan b / a a 3 tan sec 2 a z

E= s
d = s cos
3
3

2 0 =0
2 0
a sec
=0

az

a
[
cos(tan 1 b / a) cos 0]a z = s
1a z = s
2 0
2 0 a 2 + b 2
2 0

a
1
a z V / m
2
2
a
+
b

This result indicates that the electric field is everywhere normal to the surface of the
sheet. Also one may observe that the field E at a perpendicular distance a on the ve zaxis, from the surface of the sheet, is same in magnitude to the previous one but it directs
in the ve z direction. Therefore, the field, of a sheet charge having +ve surface charge
density, directs outwards. However, it directs toward the sheet if the surface charge
density is ve.
Now using E with the position of the point on the z-axis as variable, we can find V from:

z
V = E dl = s 1
a z (da + da + dza z )
2 0
z 2 + b2
=

z
= s 1
dz = s z z 2 + b 2

2 0
2 0
z 2 + b2

s
2 0

[a

+ b2 a

Example 1.15
Show that the answer to example 6 tends to s/20 in the limit of
the sheet having an infinite radius. Also show that, the field due to two infinitely
extended sheet charge placed face to face only remains in the space between them.
Solution: In the final equation of E of question, if b is infinite, then the equation of the
field E reduces to E= (s/20) az.
Let us now consider two sheet charges, one having a surface charge density of + s C/m2
and the other having a surface charge density of - s C/m2, face to face as shown in the
following figure.

+ s
- s

The dashed arrows show the direction of the field due to the sheet with + s and the solid
arrows show the same for the sheet with - s. This particular system of sheet charge
results zero field outside the sheets; the field only remains in between the sheets which is
qual to (s/0) and will be directed from the sheet with + s to the sheet with - s.
Example 1.16 A thin annular disc of inner radius a and outer radius b carries a uniform
surface charge density s. Determine E at any point on the z-axis when z 0.

Solution: Let us consider an elemental charge sds=sdd at (, , 0). The position

vector, r, of the observation point P and that of the charge position, r, are given by
r=z az, r= a. Therefore r-r=(z az - a) and |r-r|=(z2+2)1/2.
Using (9), we obtain
b 2 d d
( za z a )
E= s 2
4 0 a 0 ( z + 2 ) 3 / 2
It is easy to prove that there would be no a component of E because of the symmetry
of the charge distribution with respect to P.
Therefore E =

s
4 0

For b , E =

b 2

(z
a 0

z
1
1
d d
za z = s 2
2
a
2
2 3/ 2
2 1/ 2
2 1/ 2 z
2 0 (a + z )
+ )
(b + z )

s z
1
a
2
2 1/ 2 z
2 0 (a + z )

Example 1.17
Charge is uniformly distributed over a strip which is infinite in
length and which has a width W with a surface distribution s C/m2. Determine the
electric field at a point P which is on a line perpendicular to the strip and which is a
distance D away from the center of the strip.
Solution:

z
r- r= D az - x' ax x P
r=D az
W
dx y

x
r= x' ax
sdx'
Sheet is assumed as the sum of some line charges of thickness dx' having a charge
distribution of sdx. The field for each line is determined first and then use the
superposition theorem to determine the required field.
Following the result of example 1.5, the field for an infinitely long line charge at P is give
by
l (r r )
s dx (r r )
dE =
=
2
2
2 0 r r
2 0 r r
dE =

s ( Da z x a x )dx
2 0 r r

s ( Da z x a x )dx
2

2
W / 2 2 0 {D + x )
W /2

E=

Let x'=D tan

dx'=D sec2 d and D2+x'2=D2 sec2.

( Da z x a x )dx

E = s
= s
2
2

2 0
W / 2 2 0 {D + x )
W /2

tan 1 (W 2 D )

[a

tan a x ]d =

tan 1 (W 2 D )

s
W
(tan 1
)a z
2D
0

Example 1.18
Charge is uniformly distributed over a WxW square sheet with a
surface distribution s C/m2. Determine V and E at a point P on the axis perpendicular to
the sheet and through its center.

Solution:
z
xP
R

r=z az

r - r= z az - x ax - y ay
r=x' ax+y' ay
dx'
y

dy'

sdx'dy'

x
Sheet is assumed as the sum of some elemental surface of area dxdy and having a charge
distribution of sdx'dy'. The field for each elemental surface is determined first and then
use the superposition theorem to determine the required field. The position of the
elemental surface is assumed at (x', y', 0) and the point P is taken at (0, 0, z).
The field dE at P due to the elemental surface charge is
dx dy (r r )
s W / 2 W / 2 ( x a x y a y + za z )dx dy
dE = s

E
=
3
4 0 W/ 2 W/ 2
( x 2 + y 2 + z 2 ) 3 / 2
4 0 r r
Let us integrate the above equation first w.r.t x' and then w.r.t y'.
Assume x'=(y'2+ z2)1/2 tan
dx'=(y'2+ z2)1/2 sec2 d and (x'2 + y'2 + z2)1/2=(y'2+ z2)1/2 sec .

( y a y + za z )dx
( x a x )dx
+ 2
2
2 3/ 2
(x + y + z )
(x + y + z )
(x + y 2 + z 2 ) 3 / 2
( y a y + za z ) cos d
( y a y + za z ) sin
sin d
cos
= 2
a x +
= 2
ax +
2 1/ 2
2
2
2 1/ 2
( y + z )
( y + z )
( y + z )
( y 2 + z 2 )
( y a y + za z ) x
( y 2 + z 2 )1 / 2
=
a
+
x
( y 2 + z 2 )1 / 2 ( x 2 + y 2 + z 2 )1 / 2
( y 2 + z 2 )( x 2 + y 2 + z 2 )1 / 2
( y a y + za z ) x
ax
=
+
( x 2 + y 2 + z 2 )1 / 2 ( y 2 + z 2 )( x 2 + y 2 + z 2 )1 / 2
Now

W /2

W / 2

( x a x y a y + za z )dx
2

2 3/ 2

( x a x y a y + za z )dx
( x 2 + y 2 + z 2 ) 3 / 2

2( y a y + za z )(W / 2)

( y + z 2 )[(W / 2) 2 + y 2 + z 2 ]1 / 2
2( ya y + za z )(W / 2) dy
W /2
Therefore, E = s
2
4 0 W / 2 ( y + z 2 )[( W / 2) 2 + y 2 + z 2 ]1 / 2
2

2(W / 2) s
=
4 0

W /2

W / 2

x =W / 2

( y a y + za z ) x

ax
= 2
+

2
2 1/ 2
( y 2 + z 2 )( x 2 + y 2 + z 2 )1 / 2 x =W / 2
(x + y + z )

( y a y + za z ) dy
( y 2 + z 2 )[(W / 2) 2 + y 2 + z 2 ]1 / 2

( y a y + za z ) dy

Now

( y 2 + z 2 )[(W / 2) 2 + y 2 + z 2 ]1 / 2
y a y dy
za z dy
= 2
+ 2
2
2
2
2 1/ 2
2
( y + z )[(W / 2) + y + z ]
( y + z )[(W / 2) 2 + y 2 + z 2 ]1 / 2
To solve the first integration, let us assume t = [(W / 2) 2 + y 2 + z 2 ]1 / 2 . Therefore

dt = y [( W / 2) 2 + y 2 + z 2 ] 1 / 2 dy and y 2 + z 2 = t 2 ( W / 2) 2
y a y dy
a y dt
2
= 2
2
2
2
2 1/ 2
( y + z )[(W / 2) + y + z ]
t ( W / 2) 2
If we put t = ( W/2) sin , dt = ( W/2) cos d then
y a y dy
a y dt
(W / 2) cos d
sec d
=
( y 2 + z 2 )[(W / 2) 2 + y 2 + z 2 ]1 / 2 t 2 ( W / 2) 2 = (W / 2) 2 cos 2 a y = (W / 2) a y
t
1
1 + sin
1
(W / 2)
=
ay =
ay
ln
ln
t
2(W / 2) 1 sin
2(W / 2)
1
(W / 2)
1+

1
(W / 2) + [( W / 2) 2 + y 2 + z 2 ]1 / 2
=
ln
ay
2(W / 2) (W / 2) [(W / 2) 2 + y 2 + z 2 ]1 / 2
To solve the second integration, let us assume y = [( W / 2) 2 + z 2 ]1 / 2 tan . Therefore
dy = [(W / 2) 2 + z 2 ]1 / 2 sec 2 d , y 2 + z 2 = [(W / 2) 2 tan 2 + z 2 sec 2 ] and

[( W / 2) 2 + y 2 + z 2 ]1 / 2 = [( W / 2) 2 + z 2 ]1 / 2 sec
za z dy
za z [(W / 2) 2 + z 2 ]1 / 2 sec 2 d
2
=
( y + z 2 )[(W / 2) 2 + y 2 + z 2 ]1 / 2 [( W / 2) 2 + z 2 ]1 / 2 sec [( W / 2) 2 tan 2 + z 2 sec 2 ]
za z d
za z cos d
=
=
2
sin
[(W / 2) 2 sin 2 + z 2 ]
[(W / 2) 2
+ z 2 sec ]
cos
Again if we put (W/2)sin=z tan, then (W/2)cos d=z sec2 d
za z dy
za z cos d
z 2 a z sec 2 d
2
=
=
( y + z 2 )[(W / 2) 2 + y 2 + z 2 ]1 / 2 [(W / 2) 2 sin 2 + z 2 ] ( W / 2) z 2 sec 2
=

az
a z d
1
W / 2

=
=
tan 1
sin a z
( W / 2) ( W / 2) ( W / 2)
z

1
(W / 2) y
tan 1
a
2
2
2 1/ 2 z
( W / 2)
z[(W / 2) + y + z ]

( ya y + za z ) dy

W /2

Now,

W / 2

( y + z )[(W / 2) 2 + y 2 + z 2 ]1 / 2
2

y =W / 2

1
(W / 2) + [(W / 2) 2 + y 2 + z 2 ]1 / 2
=
ln
ay +
2
2
2 1/ 2
2(W / 2) (W / 2) [(W / 2) + y + z ] y = W / 2

y =W / 2

(W / 2) y
tan 1
az

2
2
2 1/ 2
z[(W / 2) + y + z ] y =W / 2
( W / 2)

2
(W / 2) 2
tan 1
a
2
2
2 1/ 2 z
( W / 2)
z[(W / 2) + (W / 2) + z ]

Therefore, E =
=

2(W / 2) s
4 0

( ya y + za z ) dy

W /2

W / 2

( y + z )[(W / 2) 2 + y 2 + z 2 ]1 / 2
2

s
(W / 2) 2
tan 1
a
2
2
2 1/ 2 z
0
z[( W / 2) + (W / 2) + z ]

The solution of E, however, can be made simple by using the knowledge of symmetry
about the z-axis. In this case we could solve only
W /2 W /2
za z dxdy
E= s
2
4 0 W / 2 W / 2 ( x + y 2 + z 2 ) 3 / 2
Example 1.19
A spherical conducting shell of radius a carrying a charge Q
distributed uniformly over its surface. Determine the potential i) at an external and ii) at
an internal points.
Solution:
A
C
d
a

L
r
R

B D
Sphere is assumed as the sum of some rings of width ad. The potential for each ring is
determined first and then use the superposition theorem to determine the required
potential.
Potential at P due to the charged annular ring of radius CL=a sin and width AC= ad is
given by
s 2CL AC s 2a sin ad s 2a 2 sin d
=
=
dV =
4 0 r
4 0 r
4 0 r
2
2
2
2
2
2
But r =(CL) +(PL) =(asin) +(R-OL) =(asin) +(R-a cos)2=R2+a2-2Ra cos
2rdr=2Rasin d
asin d/r=dr/R
2adr
Qdr
dV = s
=
4 0 R
8 0 a R
Therefore when P is an external point V =
when it is an internal point V =

a+R

R+a

Qdr
Q
Q
R+a
=
r R a =
and
8 0 a R 8 0 a R
4 0 R
Ra

Qdr
Q
Q
a+R
=
r aR =
8 0 a R 8 0 a R
4 0 a
aR

Example 1.20
Show that the E-field produced by a spherical conducting shell of
radius r and carrying a charge Q is identical to that of a point charge Q placed at the
centre of the sphere for >r and is zero for <r.
Solution:
We can use E = V and the result of example 19 to determine E.

Q
a +
a )V =
a r which is
Therefore for >r E = V = ( a r +
r
rsin
r
4 0 r 2
identical to that of a point charge Q placed at the centre of the sphere.

And for <r E = V = ( a r +

a +
a )V = 0 .
r
r
rsin
Example 1.21
Determine E at a point due to a volume charge density. Assume
the volume charge density is v C/m3.
Solution:
ar
Say OP=r
OQ=s sin
P
QP=s cos
s dv=r2 sin d rdd
Q

O
r0

r
Volume charge distribution

O
If we consider the similar incremental volume dv at the similar position on the opposite
side of the line OP, it will show that the resultant field will act only in the OP direction or
in the radial direction ar. The magnitude of the field along the line OP due to the
incremental volume dv is only cos to that along OP. Now the magnitude of the field at
P due to vdv along ar is given by
dv
r 2 sin dr dd
dE = v 2 cos = v
cos .
4 0 s
4 0 s 2
r 2 = (r s cos ) 2 + ( s sin ) 2 = r 2 + s 2 2rs cos
r 2 + s 2 r 2
2rs
2
Similarly, s = (r r cos ) 2 + (r sin ) 2 = r 2 + r 2 2rr cos

cos =

r 2 + r2 s 2
sds
sin d =
2 rr
rr
2
2
r dsdr d r r
1 +

dE = v 2
8r 0
s 2

cos =

E=

v
8r 2 0

r + r

r 2 r2
=0 r=0 s=rr 1 + s 2
r0

r dsdr d

r +r

r
r 2 r 2
v 0 2
v (4 / 3)r 3

1
+
r
dsd
r
4
r
d
r
=
=

s 2
4r 2 0 r =0
4r 2 0
r =0 s = r r
Q
Q
= 2
E=
ar
4r 0
4 0 r 2
Example 1.22
The components of electric intensity in a certain region are given
Ky
Kx
and
Ey = 2
. Prove that the lines of force are the
by E x = 2
2 1/ 2
(x + y )
( x + y 2 )1 / 2
family of concentric circles.
Solution:
The ratio of Ex and Ey is given by
E x dx
y
x2
y2
=
=
xdx = -ydy
+ C1 = + C2 x 2 + y 2 = C
E y dy
x
2
2
2
2
where C, C1 and C2 are constants. The relation x + y =C is an equation of circle with
radius C and having center at the origin of the coordinate system. Thus we may have a
family of concentric circles with C as variable.
Example 1.23 An electron and a proton separated by a distance of 10-11 meter are
symmetrically arranged along the z-axis with z=0 as its bisecting plane. Determine the
potential and E at P(3, 4, 12).
Solution: The position vector, r, of the observation point P(3, 4, 12) is given by
r=3 ax +4 ay +12 az and |r|=13m. The dipole moment p=1.6x10-19x10-11 az=1.6x10-30az.
Now the potential at P can be obtained from (20) as
p ar
V=
= 7.865 10 23 V
2
4 0 r
And E is given by (19) as
p
E =
[2cosa r + sina ]
4 0 r 3
Now with the help of the coordinate system we have

v
4r 2 0

r0

x2 + y2
z 12
5
cos = =
and sin =
=
r 13
r
13
25
E = 5.04 10 [24a r + 5a ] V/m
New Example 1: A rod of length L has a uniform charge per unit length and a total
charge Q. Calculate the electric field and potential at a point P along the axis of the rod at
a distance d from one end.

y
x

P
d

Lets put the origin at P. The linear charge density and Q are related by

L = =

Q
L

y
dL=dx dQ = L dL= dx

r= x ax

d
Here r=0. Now from E =

E=

d +L

4 0

x a x
x3

Again from V =

V=

1
4 0

d +L

1
4 0

dx =

1
4 0
dx =

L (r r)

r r

L
dL we have

d +L

1
L
ax =
ax
4 0 x d
4 0 d (d + L)

L dL we have

r r
L

d+L
d +L
ln x d =
ln
4 0
4 0
d

New Example 2: A hollow cylinder of length 2L and radius a has its axis along the zdirection and is centered about the z=0 plane as shown in the following figure. Its outer
surface at =a has a uniform distribution of surface charge s. Determine electric field
and potential at z>L.

l= sdz

Solution: The cylinder is to broken up into several incremental hoops of height dz. The
line charge density for each hoop is given by l= s dz. Now, the electric field dE at z for

a single hoop at z can be obtained from the result of Example 1.12,

Qz
E=
a z , and is given in the following:
4 0 ( z 2 + b 2 ) 3 / 2
(z z ) s 2adz
dE =
az
2
2 3/ 2
4 0 ( z z ) + a

E =

(z z ) s 2adz

4 [(z z )

Let z-z=a tan

+ a2

3/ 2

az

dz= -a sec2 d and ( z z ) + a 2

(z z ) s 2adz
L 4 [(z z )2 + a 2 ]3 / 2 a z
0
L

E=

Hence

a
= s
2
2 0 ( z z ) + a 2

1/ 2

3/ 2

=a3 sec3.

s L
sin d a z = s cos

2 0 L
2 0

sa
1

2
2 0 ( z L ) + a 2

1/ 2

[(z + L)

1/ 2

+ a2

1
2

To determine voltage for a single hoop at z, let us use the result of example 1.13,
l a
V ( z) =
, which results
2
2 0 ( z + a 2 )1 / 2
s dz a
dV =
1/ 2
2 0 ( z z ) 2 + a 2

V =

sa L
dz

2 0 L [( z z ) 2 + a 2 ]1 / 2

Let z-z=a tan

Hence

dz= -a sec2 d and ( z z ) + a 2

a L
a
V = s sec d = s ln sec + tan
2 0 L
2 0
s a ( z L) +
=
ln
2 0

a ( z L) +
= s ln
2 0 ( z + L) +

1.8

[(z L)

+ a2

1/ 2

a ( z z ) +
= s ln
2 0

] ln ( z + L) + [(z + L)
a

[(z L)
[(z + L)

]
+a ]

+ a2

=a sec .

[(z z )
a

+ a2

+ a2

Electric Flux Density, D

The force acting on a test charge will move it along a certain path. This path is called the
electric flux line/field line. The flux lines have no real existence but they are a useful
concept to represent, visualize, and describe the electric field. The electric flux lines have
the following properties:
a) It must be independent of the medium,

b) Its magnitude solely depends upon the charge from which it originates,
c) If a point charge is enclosed in an imaginary sphere of radius r, the electric flux must
pass perpendicularly and uniformly through the surface of the sphere, and
d) The electric flux density D, the electric flux per unit area, is then inversely
proportional to r2.
The electric field E also satisfies the above properties except that its magnitude depends
upon the permittivity of the medium. Therefore, D and E are related as D=E. For air or
free space =0.

1.9

Gauss Law

Gauss law states that The electric flux coming out from any closed surface is equal to
the total charge in the volume enclosed by that surface.
ds

D
Assume an incremental surface dS over which flux density D is constant. If D makes an
angle with the normal to the surface, the number of flux lines d crossing the surface ds
is D ds cos = Dds. Integrating over the closed surface we have the mathematical
formulation of Gauss law:
(21)
D ds = Q or E ds = Q/ 0 .
s

To prove the Gauss law we may consider a point charge Q in free space at the origin of
the coordinate system and an imaginary sphere of radius r around it. Then,
1 Q
D= E=
a r and ds=r2 sin d d ar
4 r 2
Therefore,
2
Q
= D ds =
sin dd = Q .
4
s
0 0
If Q is expressed in terms of the volume charge density as Q = v dv , using the
v

divergence theorem we have from (21), D = v E = v / 0

In summary, Gausss law provides a convenient tool for evaluating electric field.
However, its application is limited only to systems that possess certain symmetry,
namely, systems with cylindrical, planar and spherical symmetry. The use of the above
equation (21) to determine E or D is easy if we are able to choose a suitable Gaussian
surface (any closed surface; generally cylindrical surface is considered for line charge,
rectangular or cylindrical surface is considered for surface charge and spherical surface is
considered for volume charge distribution) around the charge that satisfies the following
conditions:
1) D or E is everywhere either normal or tangential to the closed surface, so that
Dds becomes either Dds or zero, respectively.
2) D has the same value at all points on the surface where D is normal.

1.10 Gauss Law in Point form

z

y
x
Let P(x,y,z) be a point where the flux density is D(x,y,z) and it increases in the positive
direction of the coordinate axes.
D ds = D ds + D ds + D ds + D ds + D ds + D ds
s

front face

back face

left face

right face

top face

bottom face

Now only the x component of D contributes to the 1st and 2nd integrations. Similarly the
y component of D contributes to the 3rd and 4th integration, and the z component of D
contributes to the 5th and 6th integrations.
On the front face abcd, the x component of flux density is
x
x
Dx ( x +
Dx
, y , z ) D x ( x, y , z ) +
2
2 x
x, y,z
and on the back face efgh, the x component of flux density is
x
x
Dx ( x
, y , z ) D x ( x, y , z )
Dx
.
2
2 x
x, y, z
Thus the flux leaving the volume dv=xyz through the surface abcd is
xyz
D ds = D x ( x, y, z )yz +
D x x, y , z

2
x
front face
and the flux entering the volume dv through the surface efgh is
xyz
D ds = D x ( x, y, z )yz
Dx x, y, z

2
x
back face
On the left face abef, the y component of flux density is
y
y
D y ( x, y
, z ) D y ( x, y , z )
Dy
2
2 y
x, y, z
and on the right face cdgh, the x component of flux density is
y
y
, z ) D y ( x, y , z ) +
D y ( x, y +
Dy
.
x, y, z
2
2 y
Thus the flux entering the volume dv through the surface abef is
xyz
D ds = D y ( x, y, z )xz
Dy

2
y
left face
x, y, z
and the flux leaving the volume dv through the surface cdgh is
xyz
D ds = D y ( x, y, z )xz +
Dy

x, y, z
2
y
right face

On the top face bche, the z component of flux density is

z
z
D z ( x, y , z + ) D z ( x, y , z ) +
Dz x, y, z
2
2 z
and on the bottom face afgd, the z component of flux density is
z
z
D z ( x, y , z ) D z ( x, y , z )
Dz x, y, z .
2
2 z
Thus the flux entering the volume dv through the surface afgd is
xyz
D ds = D z ( x, y, z )xy
Dz x, y, z

2
z
bottom face
and the flux leaving the volume dv through the surface bche is
xyz
D ds = D z ( x, y, z )xy +
Dz x, y, z

2
z
top face
Therefore, the net flux leaving the volume dv in the x direction is
xyz
D ds + D ds = {D x ( x, y, z )yz +
Dx
}

2
x
x, y, z
front face
back face

xyz

D x x , y , z } = xyz D x
2
x
x
x, y, z
Similarly, the net flux leaving the volume dv in the y direction is
xyz
Dy
}
D ds + D ds = {D y ( x, y, z )xz +

x, y, z
2
y
right face
left face
{D x ( x, y, z )yz

xyz

Dy
} = xyz D y
x, y, z
x, y , z
2
y
y
and the net flux leaving the volume dv in the z direction is
xyz
D ds + D ds = {D z ( x, y, z )xy +
Dz x, y , z }

2
z

top face
bottom face
{D y ( x, y, z )xz

xyz

Dz
} = xyz D z x , y , z
2
z
z
x, y, z
Therefore net flux leaving the closed surface is

s D ds = ( x Dx + y D y + z Dz )xyz = ( x Dx + y D y + z Dz )dv
{D z ( x, y, z )xy

D ds
or,

dv

=(

Dx +
D y + Dz )
x
y
z

In case of volume dv tends to zero, following the divergence theorem, the right hand side
of the above equation becomes

D ds
lim

dv 0

dv

= D

Again from Gauss law

D ds = Q E ds = Q/ 0
s

(22)

D = v

E =

v
0

(23)

where v is the volume charge density. Equation (23) is the Gauss law in point form, one
of the Maxwells equations. In the charge free region E = 0
From equation (22)
s D ds = Q = v v dv = v ( D)dv
(24)
D

d
s
=
(

D
)
dv

Equation (24) is known as divergence theorem.

Equation (23) shows that the static electric field is not solenoidal. It is noted that a
divergenceless field is called a solenoidal field whose field lines always close upon
themselves. In the source free region, the electric field however becomes solenoidal.

1.11 Use of Gauss's Law to find the E-field

(a)

Spherical but non-point charge

A conducting sphere of radius a carries a charge Q on its surface. By symmetry the

resultant E-field must by spherically symmetric and hence can only have a radial
component. Take as a Gaussian surface (dotted) a sphere of radius r (>a) placed
concentric to the conductor. By symmetry the field E must be constant at all points on the
dotted surface and also normal to it. Flux through this surface
= field x surface area
= E.4r2
Total charge enclosed by surface is Q E.4r2=Q/0
hence E=Q/(40r2) V= Q/(40r) and zero at
For r<a (within conductor) charge contained within the Gaussian surface is zero (all
charge is on surface) E=0. Hence potential within sphere must be constant and equal
to surface potential = Q/(40a) Same result will be obtained if the sphere is hollow.

(b)

Infinite plane of charge (surface charge)

An infinite (or very large) sheet carries a uniform charge density s. By symmetry the
resultant E-field must have a direction normal to the plane and must have the same size at
all points a common distance from the plane. Take as a Gaussian surface a cylinder of
cross-sectional area A and height 2h.
Flux is only non-zero through ends of the cylinder. If field at cylinder ends is E then total
flux is 2EA.
As the charge enclosed is area x charge density = As, therefore, from Gauss's law
2EA=As /0 E=s /20

(c)

Infinitely long line charge

E

dS3

dS1

dS2

L
Assume an infinitely long line charge with uniformly distributed charge density of l
C/m. As seen earlier the resultant E-field must have a direction normal to the line and
have the same value at all points having a common distance from the line.
As the closed Gaussian surface, we choose a right cylinder of radius r and length L which
is coaxial with the line charge. Now from Gauss law
D ds = l L D dS 1 + D dS 2 + D dS 3 = l L
s

S1

S2

S3

2 L
LHS = 0 E dS 1 + E dS 2 + E dS 3 = 0 E dS 3 = 0 E r ddz = 2 0 LE
0
0
S

S2
S3
S3
1

Now, 2 0 LE r = l L E =

l
2 0

E=

l
a
2 0

1.12 Conductors in Electric Fields

In conductors, electrons are free to move. If the conductor is placed into an external
electric field E, a force F=-eE acts on each free electron.

The free electrons in the conductor will move in a direction opposite to the direction of
the electric field and soon pile up on the surface on one side of the conductor. The surface
on the other side, therefore, will be depleted of electrons and have a net positive charge.
These separated positive and negative charges on opposite sides of the conductor induce
their own electric field, which opposes the external field inside the conductor. The
resultant electric field inside the conductor can be found from the superposition of the
externally applied electric field and the induced electric field. When static equilibrium is
reached, the net electric field inside the conductor is exactly zero as shown in the
following figure. The zero fields inside a conductor imply that the charge density inside
the conductor is also zero as obtained by Gausss law.

Even in the absence of external electric field, an electric field will be set up in the
conductor due to the free electrons that exerts force on the electrons and making them
move away from one another. The movement will continue until all the charges reach the
conductor surface and redistribute themselves in such a way that both the charges and the
electric field inside the conductor vanish (v=0, E=0). Therefore, a solid conductor
carries all its excess charge on the surface. On the other hand, the potential difference
a

between any two points in a conductor Vab = Va Vb = E dl = 0 since E=0 inside the
b

conductor. This requires all points inside the conductor at the same potential.
As E=0 inside a conductor, the flux through Gaussian surface G is zero; and hence net
charge contained within G in Fig. (a) is also zero. A cavity in a conductor, completely
surrounded by conducting material as shown in Fig. (b), also is free of electric fields
unless the hollow region contains a net charge. A hollow conductor, therefore, shields its
interior from any outside electric fields; this is why we are safest inside our car during a
thunderstorm. It must also carry any excess charge on its outer surface to satisfy the
Gausss law for the Gaussian surface G. A hollow conductor has a charge on its inner
surface that is equal in magnitude and opposite in sign to any charge that may be
enclosed within the hollow region. This can be proved by considering a Gaussian surface
that lies within the conductor and that encircles the hollow region containing a charge as
shown in Fig. (c).

The presence of the charge on the conductor surface, however, sets an electric field
normally outward to the surface. This will be explained under the section Boundary
Conditions in details.

1.13 Dielectrics in Electric Fields

Dielectrics differ from conductors in that they have no free charges that can move
through the material under the influence of an electric field. Although the electrons are
unable to move far in a dielectric, a slight displacement of the positive and negative
charges in opposite directions (the positive charges move in the direction of the field and
the negative charges move in the opposite direction) happens when an electric field E is
applied. A dielectric, in which this charge displacement has taken place, is said to be
polarized. Each molecule of a polarized dielectric therefore acts as an electric dipole of
moment p which points in the same direction as the applied field. And the field produced
by these induced dipoles opposes the applied field. The total electric field within the
dielectric is thus less than the externally applied electric field.
The number of dipole moments per unit volume in a polarized dielectric is defined as the
polarization P which is proportional to the total electric field Em in the dielectric (Note
that the total electric field in the dielectric is different from the applied field.) and is given
by P=0e Em where e is a dimensionless constant known as the electric susceptibility of
the dielectric.
The effect of polarization is to produce accumulations of bound charge bv within the
dielectric and bs on the surface. Here, bv = P and bs = P a n (an is the unit
vector which is normally outward from the relevant surface) are, respectively, volume
and surface charge density. In the case of homogeneous and isotropic dielectrics where
the polarization of each molecule is to the same extent, the volume distribution of the
bound charges bv becomes zero as the positive and negative charges of dipoles
neutralize each other. The charges induced due to polarization thus left only on the
dielectric surface. The effect of polarized dielectric on the electric field may, therefore, be
realized equivalently by removing the dielectric material with the corresponding induced
dipoles in free space.
Suppose a dielectric slab is placed between the conducting plates of a parallel plate
capacitor. The upper and the lower plates of the capacitor are assumed to be charged with
+fs and -fs C/m2, respectively, as shown the following figure. The dielectric will be
polarized and bound charges will accumulate on the surfaces of the slabs, say bs C/m2.
We could now remove the dielectric conceptually, and modify the charge density of the
capacitor plates to +(fs-bs) and -(fs-bs) with free space in between them.
fs
bs

bs
fs

In absence of dielectric, there will an electric field Evac=fs/0 an in the space between the
charged plates as explained in example 1.15; in presence of dielectric material, the field

becomes Em=(fs-bs)/0 an. Now from the definition of dielectric constant/relative

permittivity, r=(field in free space)/(field in dielectric material)= Evac/Em, we have
fs bs
fs
fs
fs bs
fs
E
E m = vac

=
+
fs = 0
+ bs

r 0

r 0 0
r 0
= D m = 0 E m + P = 0 (1 + e )E m = 0 r E m = E m

or D vac = 0 E m + P D vac
because D is independent of media or material and r=(1+e).

If the field in which the dielectric is placed exceeds a certain limit, the atomic structure of
the dielectric breaks down and the electrons will come out of the atoms. These electrons
take part in avalanche breakdown mechanism of the material. The material will then
become conducting. The maximum electric field that a dielectric can withstand without
breakdown is the dielectric strength of the material.

1.14 Boundary Conditions for Electrostatic Field

Electromagnetic problem often involve media with different physical properties and
require the knowledge of the field quantities at the interface between two media (e.g.
vacuum and a dielectric, two different dielectrics etc). For instance, we may wish to
determine how E and D vectors change in crossing an interface between two media. Let
us consider an interface between two general media shown in the following figure.
t

s
h

For the E-field take the rectangular path as the closed loop. If its height h is made
infinitesimally small then the only contribution to the line integral is along the top and
bottom edges of length dL

E dl = E

(dL) + E 2 (dL) = E1t dL + E 2t dL = 0

D1t/1=D2t/2
(25)
leads to E1t=E2t
where Et is the tangential component of E.
Hence E1t=E2t; the tangential component of E is continuous across any interface.
In order to find a relation between the normal components of the fields at the interface,
let us construct a small cylinder with its top face in medium 1 and bottom face in medium
2. The faces have an area dS, and the height of the cylinder dt is vanishingly small.
Appling Gausss law to the cylinder, we have
D ds = ( D1 an1 + D2 an 2 )dS = (D n1 Dn2 )dS = lim v (tdS )
t 0

or D n1 Dn2 = lim v t = s 1 E n1 2 E n2 = s
h 0

(26)

where an1 and an2 are the unit vectors perpendicularly outwards to the top and bottom
surfaces of the cylinder, respectively.
In a special case where medium 2 is a perfect conductor and medium 1 is a dielectric,
electric field in the conductor E2t and hence E1t in eq.(25) become zero. From eq.(26) we
have Dn1=1En1=s as the electric field in a conductor is zero. It therefore shows that the
electric field at the interface between perfect conductor and insulator is normal to the
surface.
In the case where both media (1 and 2) are perfect dielectric, s in eq.(26) becomes zero
as no free charge can exist in the dielectric.

1.15 Calculation of Capacitance

A capacitor is formed by placing an dielectric material in between two conductors and its
capacitance is defined as the ratio of the charge on one of the conductors to the potential
difference between them i.e.
C=Q/V Farad.
(27)
If once the voltage difference between the conductors is calculated and the charge on one
of the conductor is known, capacitance for any arrangement of the conductors can easily
be obtained using eq. (27).
The capacitance can also calculate from the stored electrostatic energy of the system. It is
obtained from the work done in moving a small charge dQ against a potential difference
V and is given by dWe=VdQ=QdQ/C where V=Q/C.
Q
Q
Q2 1
= CV 2
(28)
Therefore We = dQ =
C
2C
2
0
1
1
1
1
CV 2 = QV = vVdv = ( D) Vdv
2
2
2
2
By vector identity we have VD=DV+V(D)
1
1
1
1
We = ( D) Vdv = [( VD) (D V )]dv = ( VD)dv (D V )]dv
2
2
2
2
Using divergence theorem we have
1
1
We = VDds (D V )]dv
2
2
As D varies inversely as square of distance (1/r2), potential varies as (1/r) and the surface
area increases as r2, it follows that the first term in We is a function of (1/r). In order to
consider all the fields in calculating We, we have to consider a very big Gaussian surface
with r, the first term in We thus becomes zero.
1
1
1
1 D2
2
We = (D V )]dv = (D E)dv = E dv =
dv
(29)
2
2
2
2
Therefore, C can also be obtained by using eq. (28) and (29).

Now We =

Example 1.24
An infinitely long cylinder with radius a has a uniform charge
2
distribution s C/m on its surface. Determine E and V at a perpendicular distance from
the center of the cylinder.
Solution: For >a

S3
S2
a
Gaussian surface of cylindrical shape with radius .

S1

Assume the Gaussian cylinder is of unit length. It is known that the electric field from a
cylindrical charge distribution acts in the radial direction . Therefore, using Gauss law,
we have
D ds = Q = s 2a D dS 1 + D dS 2 + D dS 3 = s 2a
S1

S
1

S2

S3

0 E dS 1 + E dS 2 + E dS 3 = l L 0 E dS 3 = l L
S2

S3

E ddz = L
l

S3

2 0 E = l L

sa
a
E = s a

Now for <a, the charge enclosed by the Gaussian cylinder is zero. It suggests that E is

E =

also zero for <a.

Example 1.25
Calculate E and V as a function of radial distance for a spherical
cloud of radius a and having a uniformly charge density v.
Solution: For r>a
Assume a Gaussian spherical surface of radius r around the charge cloud. For the
spherical charge distribution, the field E acts in the radial direction ar; now from Gauss
law
a 2
4 v a 3
2

d
s
=
Q

Q
=
dv
=
r
sin
drd
d
=
s
v v r=0 =0 =0 v
3
and

D ds = 0 E dS = 0
s

4 v a 3
3

=0 =0

4 0 r 2 E r =

Er =

E r r 2 sin dd = 4 0 r 2 E r

va3
3 0 r 2

E=

va3
ar
3 0 r 2

For ra

4 v r 3
s D ds = Q Q = v v dv = r=0 =0 =0 v r sin drdd = 3
r

and D ds = 0 E dS = 0
s

4 v r 3
4 0 r E r =
3
Calculation of V:
2

For ra

=0

Er =

=0

E r r 2 sin dd = 4 0 r 2 E r

vr
r
E = v ar
3 0
3 0

va3
a r (dra r + rda + r sin da )
2
3 0 r

V = E dl

V =

a3
a3
V = v 2 dr = v
3 0 r
3 0 r
r

at r=a, V =

va3
3 0 r

va
and for ra
3 0
2

V = E dl = E dl E dl =

va 2 r
E dl
3 0
a

va
v r
v a 2 r vr
V =

a r (dra r + rda + r sin da ) =

dr
3 0
3 0
3 0 a 3 0
a
2

va 2 vr 2
V=

3 0
6 0

va 2
vr 2 v a 2

=
(

) = v (3a 2 r 2 )
3 0
6 0
6 0
6 0

Example 1.26 Charge is uniformly distributed within a spherical region of radius a. An

isolated conducting spherical shell with inner radius b and outer radius c is placed
concentrically as shown in the figure. Determine E everywhere in the region.

Solution: For 0<r<a, the charge enclosed by a spherical Gaussian surface is given by
r 2
4 3
Q= v dv = v r 2 sin drdd =
r v
3
v
0 = 0 == 0
Now from Gausss law

D ds = Q
s

Dr

sin dd a r = Q

=0 =0

4 3
r
r
r
r v Dr = v E r =
v E =
var
3
3
3 0
3 0
For a r <b, therefore the charge enclosed by a spherical Gaussian surface is given by
a 2
4 3
Q= v dv = v r 2 sin drdd =
a v
3
v
r = 0 = 0 == 0

4 r 2 Dr =

Now from Gausss law

D ds = Q
s

D r

sin dd a r = Q

=0 = 0

4 3
a3
a3
a3
a v Dr = 2 v E r = 2 v E = 2 v a r
3
3r
3r 0
3r 0
For b r <c, electric field must be zero because of the presence of the conducting material
within this region. A charge (4a3v/3) must be induced on the surface at r=b to make
4 r 2 Dr =

the field within the conductor zero. The surface at r=c must also acquire a charge
(4a3v/3) to make the region charge neutral.
For r >c, the charge enclosed by a spherical Gaussian surface is given by
a 2
4 3
Q= v dv = v r 2 sin drdd =
a v
3
v
r = 0 = 0 == 0
Now from Gausss law

D ds = Q
s

Dr

sin dd a r = Q

=0 = 0

4 3
a3
a3
a3
a v Dr = 2 v E r = 2 v E = 2 v a r
3
3r
3r 0
3r 0
Example 1.27
A spherical volume charge distribution is given by v=0(1-r2/a2)
for ra and v=0 C/m3 for r>a. Find E as a function of radial distance r and also show that
the maximum value of E is at r=0.745 a.
Solution: For r>a
Assume a Gaussian spherical surface of radius r around the charge distribution. For the
spherical charge distribution, the field E acts in the radial direction ar; now from Gauss
law
a 2
8 0 a 3
r2 2

d
s
=
Q

Q
=
dv
=
(
1

)
r
sin
drd
d
=
s
v v r=0 =0 =0 0 a 2
15
4 r 2 Dr =

Now D ds == 0 E dS = 0
s

=0 = 0

8 0 a 3
2 0 a 3
4 0 r E r =
Er =
15
15r 2
For ra
2

E r r 2 sin dd = 4 0 r 2 E r

8 0 a 3
15

2 0 a 3
E=
ar
15r 2

r2 2
r3
r5
s D ds = Q Q = v v dv = r=0 =0 =0 0 (1 a 2 )r sin drdd = 4 0 ( 3 5a 2 )
and D ds = 0 E dS = 0
s

=0

=0

E r r 2 sin dd = 4 0 r 2 E r

0 r r 3
0 r r 3
Er =
(
) E=
(
)a r
0 3 5a 2
0 3 5a 2
To determine the position of maximum E, let us first differentiate |E| w.r.t r and then
equates the result with 0; we have
0 1 3r 2
( 2 ) = 0 9r 2 = 5a 2 r = 5a = 0.745a
0 3 5a
Example 1.28
A spherical volume of radius a has a volume charge density v=kr
where r is the radial distance and k=constant. Find E and V in the region 0r.
Solution:
The field for the given charge distribution is radial, i.e., in the direction of ar. We have to
consider a spherical Gaussian surface.
For r>a

 Q = v dv =

D ds = Q
s

krr

=0 = 0

4 0 r 2 E r = ka 4

ka 4
4 0 r 2

Er =

sin drdd = ka 4

r =0 =0 =0

and D ds = 0 E dS = 0
s

E R r 2 sin dd = 4 0 r 2 E R

E=

ka 4
ar
4 0 r 2

For ra

Q = v dv =

D ds = Q
s

and D ds = 0 E dS = 0
s

4 0 r 2 E r = kr 4

Er =

sin drdd = kr 4

r =0 =0 =0

kr 2
4 0

krr

=0 =0

E r r 2 sin dd = 4 0 r 2 E r

E=

kr 2
ar
4 0

Calculation of V:
For ra
r

ka 4
a r (dra r + rda + r sin da )
2
4 0 r

V = E dl

V =

ka 4
ka 4
V =
dr
=
2
4 0 r
4 0 r

ka 4
4 0 r

ka 3
at r=a, V =
and for ra
4 0
r

ka 3
V = E dl = E dl E dl =
E dl
4 0 a

a
r

V =

ka 3
kr 2
ka 3
kr 2
a r (dra r + rda + r sin da ) =

dr
4 0 a 4 0
4 0 a 4 0

ka 3 kr 3
V=

4 0 12 0

=
a

ka 3
kr 3
ka 3
ka 3 kr 3
)=
(

4 0 12 0 12 0
3 0 12 0

Example 1.29
A point charge Q is at the center of a spherical conducting shell of
an inner radius ri and an outer radius ro. Calculate E and V as functions of the radial
distance r.
Solution:
ri
ro

Calculation of E:

We know that the field E will act in the radial direction ar. Now, for r<ri, the Gaussian
surface shown by the dashed line encloses a total charge Q, therefore from Gauss law we
have
D ds = 0 E dS = 0

D ds = Q
s

=0

4 0 r 2 E r = Q E r =

4 0 r

E=

=0

E r r 2 sin dd = 4 0 r 2 E r

Q
4 0 r 2

ar

In the range ri r ro, the electric field E will be zero as a conducting media exists there.
To become E is zero, -ve Q charge will induce on the inner surface and +ve Q charge will
induce on the upper surface of the conducting shell. Hence no charge will be enclosed by
any Gaussian surface shown by the dashed line in this range of r.
For rro, the Gaussian surface will again enclose a charge Q, therefore the electric field in
this case is also given by
Q
E=
a r for rro.
4 0 r 2
Calculation of V:
For rro
r

V = E dl

V =

4 0 r 2

Q
4 0 r

V =

dr =

Q
4 0 r

a r (dra r + rda + r sin da )

Q

4 0 r

. As the electric field is zero in the range ri r r0, no work done

4 0 ro
will occur in bringing a charge from r= r0 to r= ri; the potential at ri is equal to that at ro.
Hence the potential at r= ri is given by
Q
V=
.
4 0 ro
And for rri

at r=ro , V =

ro

rio

ro

ri

V = E dl = E dl E dl E dl

V =
V=

Q
4 0 ro
Q

4 0 ro

ri

Q
4 0 r

Q
4 0 r

=
ri

a r (dra r + rda + r sin da ) =

Q
4 0 ro

+(

Q
4 0 r

Q
4 0 ri

)=

Q
4 0

Q
4 0 ro

ri

Q
4 0 r 2

dr

1 1 1
+ )
r0 r ri

Example 1.30 A metallic sphere of radius 10 cm has a surface charge density 10 nC/m2.
Calculate the electric energy stored in the system.
Solution: The field and potential for r10 cm is obtained from Gausss law as

D ds = Q

Dr r sin dd =
2

=0 = 0

sin dd

=0 = 0

0.01
0.01
0.01
0.01
0.01
D r = 2 s E =
sar V =
s a r dl =
s dr =
s
2
2
2
0r
r
0r
0 r
0 r
Now the potential for r10 cm is constant and is equal to that at r=10 cm: V=113 Volts.
Energy We is given by We=QV/2, where Q is given by

Q=

sin dd = 4r 2 s = 1.257 10 9 C . Therefore We=71.02 nJ

= 0 =0

Example 1.31
Consider two concentric spheres. The inner sphere has radius a,
and a charge Q is distributed uniformly over its outer surface. The outer sphere has radius
b, and a charge -Q is distributed uniformly over its inner surface. Determine the potential
difference between the two spheres and the capacitance of the structure.
Solution:
a

The fields E for this system in various regions are:

E=0
for r > b
=

V=

a r for a r b
4 0 r 2
E is also zero for r<a. Therefore the center and the outer surface of the inner sphere have
the same potential. Potential difference only exists in between the spheres which is
a
a
a
Q
Q
V = E dl =
a r (dra r + rda + r sin da ) =
dr
2
2
b
b 4 0 r
b 4 0 r

4 0 r

Now, C =

=
b

Q
4 0 a

Q
4 0 b

Q(b a )
4 0 ab

Q 4 0 ab
=
F
V (b a )

Example 1.32
A capacitor is constructed in the form shown in the following
figure. The capacitor consists of two large parallel plates 3 m apart, with a dielectric slab
1 m thick located midway between the plates. The electric field in region I is E=(s/0)ax
where s is the surface charge density on one plate. The surface area of each plate is A.
Determine the capacitance.

+
+
+
+
+

II

III
0

x
1

Solution:
The electric flux density in all regions is D=sax. The voltage difference between the
plates is
0

V03 = - E dl = E dl - E dl - E dl
3

0
s

a x dxa x - s a x dxa x - s a x dxa x

0
3 0
2
1
2 + 0

= s + s + s = s (
)
0
0 0
A s
A 0
Q
Now C =
=
=
V03
V03
(2 + 0 )
2

V03 =

Example 1.33
A system consists of two coaxial cylinders of radii a and b (b>a).
The inner cylinder is solid and has a charge + per unit length spread uniformly
throughout its volume. The outer cylinder has negligible thickness and carries a charge -
per unit length. Use Gausss law to find the E-field in the regions (i) <a, (ii) a<<b and
(iii) >b where is the radial distance from the axis of the cylinders. Sketch the form of
the dependence of the E-field on .
Solution:
The length of the coaxial cylinder is assumed to be large so that the E-field is radially
outwards, i.e. E acts in the direction a and constant over a cylindrical Gaussian surface.
Change per unit length is uniformly distributed within inner cylinder so there will be a
non-zero field for <a. Charge enclosed within a cylinder of radius (<a) is (2)/(a2).
This gives
1 2
2
s D ds = Q Q = a 2 and s D ds = 0 0 0 E ddz = 2E

2 0 E =

E =
2
2 0 a 2
a

For a<<b all of charge on inner cylinder is enclosed but no charge on outer cylinder. In
this case

D ds = Q Q =

1 2

and D ds = 0 E ddz = 2 0 E

2 0 E = E =

0 0

For >b charges on two cylinders cancel so E=0.

|E|

Example 1.34
A positive point charge Q is at the center of a spherical dielectric
shell of an inner radius ri and an outer radius ro. The dielectric constant of the shell is r.
Determine E, V, D and P as functions of radial position.
Solution:
ri

ro

=0r
=0

Let us consider a spherical Gaussian surface of radius rro as shown by the dashed line.
In this region E can be found using Gausss law as

D ds = Q now

D ds = 0 E dS = 0

Er =

=0 = 0

E r r 2 sin dd = 4 0 r 2 E r

E=

ar
4 0 r
4 0 r 2
and V at any point in this region is given by
r
r
Q
V = E dl
V =
a r (dra r + rda + r sin da )
2

4 0 r
2

V =

Q
4 0 r

dr =

Q
4 0 r
2

Q
4 0 r

Therefore, D=0E= Q/(4r ) ar and P=0 as it is free space.

In the region ri r ro(with in the dielectric) the electric field Em can be calculated from
the fact that D is equal in both region. Therefore, D=Dm=0r Em= Q/(4r2) ar and hence
Em= Q/(40rr2) ar. V at any point in this region is equal to V at r=ro (as obtained in the
previous section and is Q/(40ro)) plus the pd from ro to that point as

ro

ro

V = E dl = E dl E dl

V=

V =

4 0 ro
Q

4 0 ro

Q
Q
Q
dr
a r (dra r + rda + r sin da ) =

2
4 0 ro ro 4r 2
ro 4r

Q
4r

=
ro

Q
4 0 ro

+(

Q
4r

Q
Q
1 1
1
)=
[(1 ) +
]
r r0 r r
4ro
4 0

Now from Dm=0Em+P, we have P= Dm-0Em= Dm-0Dm/0r=Dm(1-1/r)= (1-1/r)

[Q/(4r2)] ar.
In the region rri, E can be found from the Gaussian surface of the sphere of radius r<ri as

D ds = Q

now

Er =

D ds = E dS =

E=

=0

=0

E r r 2 sin dd = 4 0 r 2 E r

ar
4 0 r
4 0 r 2
and V at any point in this region is equal to V at r=ri which is V= Q/(40)[(11/r)1/ro+1/(rri)] as obtained from the previous section plus the pd from ri to that point as
2

obtained from V = E dl . Therefore the potential at a point where rri,

r = ri

Q
1 1
1
1 1
1
V=
[(1 ) +
] E dl =
[(1 ) +
]
r r0 r ri ri
r r0 r ri
4 0
4 0
r

Q
1 1
1
r 4 0 r 2 a r (dra r + rda + r sin da ) = 4 0 [(1 r ) r0 + r ri ] r 4 0 r 2 dr
i
i
Q

Q
1 1
1
V=
[(1 ) +
]+
r r0 r ri
4 0
4 0 r

=
ri

Q
4 0

[(1

Q 1 1
1 1
1
) +
]+
( )
r r0 r ri
4 0 r ri

1 1
1
)+ ]
4 0
r r0 ri
r
2
In this region, D=0E= Q/(4r ) ar and P=0 as it is free space.
Example 1.35
The following figure gives a cross-sectional view of a co-axial
cable. The inner cylindrical conductor of radius a is charged of l C/m and the outer
cylinder of radius b is charged of -l C/m. Determine the capacitance between the
cylinders. Assume the dielectric between the cylinders has a permittivity .
Solution:
V =

[(1

)(

a
b

The field E outside the outer conductor is zero as the total charge enclosed by the
Gaussian surface considered around the outer conductor is zero. In the region arb, the
field for unit length of the inner conductor is given by Gauss law as

D ds = Q

Q = l

1 = 2

and

D ds = E dS = E ddz = 2 E
s

E =

l
2

E=

z =0

=0

l
a
2

The pd between the conductors is

r

V = E dl = E dl E dl

b
a

a
l
l

b
V = 0
a (da + da + dza z ) = 0
d = l ln = l ln
2
2
2
2 a
b
b
b
a

2
F/m
b
V
ln
a
Example 1.36
Solve the above problem considering that the space between two
cylinders is filled up with two different dielectrics (one is of permittivity 1 and attached
to the inner cylinder with a thickness t1 and the other is of permittivity 2).
C =

Solution:
The field E outside the outer conductor is zero as the total charge enclosed by the
Gaussian surface considered around the outer conductor is zero. In the region a(a+t1)
and (a+t1)b, the fields for unit length of the inner conductor are

E=

l
l
a and E =
a , respectively.
2 1
2 2

Now the pd between the conductors is

a
a+t
a+t
a

a l

l
Vab = E dl = E dl + E dl =
a dl +
a dl
2 2
b
b
b
a+t

a +t 2 1

a+t

a l

l
=
d +
d = l ln
2 2
2 1
b
a +t 2 1

a +t

a+t

l
+
ln
2 2
b

a
a+t
a+t
b
= l ln
+ l ln
ln
+ l ln
=
b 2 1
a
2 2 a + t
2 1 a + t 2 2
l
a+t
b
2 ln
=
+ 1 ln

2 1 2
a
a + t
2 1 2
a+t
b
Vab
2 ln a + 1 ln a + t
Example 1.37
Two wires having charges of l C/m and -l C/m are spaced D
meters apart. The radius of each wire is meter (<<D). Calculate the capacitance
between the wires.
Therefore the capacitance per unit length of the cable C =

Solution: Consider the wires are along y-axis on the x-y plane and assume the wire with
l charge is at x=0 and the other at x=D.
y

(0,0)

x=D

Assuming the wires are very long, the field at any generalized point P in between the
conductors, x meter away from the origin of the coordinate system using Gausss law can
be obtained as E = [

l
l
+
]a x .
2 0 x 2 0 ( D x)

Now determine pd between the conductors from

V = E dl =
D

[ 2

2 0 ( D x)

]a x (dxa x + dya y + dza z )

l
l

l
Dx
= [
+
]dx = l [ln x ln( D x)]
=
ln
2 0 x 2 0 ( D x)
2 0
2 0
x
D
D r

D
= l ln
0

0
l

C =

ln

. .

F/m

Example 1.38
A single core lead sheathed cable has a conductor diameter of 1 cm
and the two layers of different insulating materials each 1.1 cm thick. The relative
permittivities are 3.3 (inner) and 2.5 (outer). Calculate the potential gradient at the
conductor surface (maximum field in the inner insulator) when the potential difference
between the conductor and the lead-sheath is 65 kV.
Solution: Assume that the line charge on the inner conductor is l C/m. In the dielectrics,
maximum fields occur at the respective inner radii. Therefore, the maxm fields Emax1 and
Emax2 in inner and outer insulators are given by
E max 1 =

2 3.3 o

a
= 0.005

and E max 2 =

2 2.5 o

a .
= 0.016

However, if V1 and V2 are PDs across the insulators, then they are given as
0.005
0.005
la
l
V1 = E1 dl =
da =
ln(0.016 / 0.005)
2 3.3 o
= 0.016
= 0.016 2 3.3 o

= 0.00582E max1

0.016

V2 =

0.016

E 2 dl =

= 0.027

= 0.027

la
2 2.5 o

da =

l
2 2.5 o

ln(0.027 / 0.016)

= 0.0084E max2
Now at the interface of two insulators

E max 2 =

l
l 3.3 0.005
a =
a
2 2.5 o 0.016
2 2.5 o 0.016 3.3 0.005

.
E max 1 3.3 0.005
=
2.5 0.016
V1+V2=(0.00582+0.0084x0.4125)Emax1=65000 Emax1=7000538 V/m 70kV/cm
Example 1.39
Determine the maximum working voltage of a single core leadsheathed cable with conductor diameter 1 cm and sheath 5 cm inner diameter. Two
dielectrics are used:
Inner layer: maximum working stress 72 kV/cm, relative permittivity r=4.
Outer layer: maximum working stress 60 kV/cm, relative permittivity r=2.5.
Solution: Assume that the line charge on the inner conductor is l C/m and the thickness
of the inner dielectric is t m. In the dielectrics, maximum fields occur at the respective
inner radii. Therefore, the magnitudes of maxm fields Emax1 and Emax2 in inner and outer
dielectrics are given by
E max1 =

2 4 o

= 7200 kV/m and E max 2 =

= 0.005

2 2.5 o

= 0.005+ t

= 6000 kV/m
Putting the value of l=8 C/m as obtained from the equation of Emax1 into equation of
Emax2 we obtain t=0.0046 m.
If V1 and V2 are PDs across the dielectrics, then they are given as
0.005
0.005
la
l
V1 = E1 dl =
ln(0.0096 / 0.005)
da =
2 4 o
= 0.0096
= 0.0096 2 4 o

= 23.46kV.
0.0096

V2 =

E 2 dl =

= 0.025

0.0096

= 0.025

la
2 2.5 o

da =

2 2.5 o

ln(0.025 / 0.0096)

= 55kV.
Therefore the maximum working voltage of the cable is (23.46+55)=78.46 kV.
Example 1.40
The radius of the inner conductor of a co-axial cable is 0.40 cm.
Concentric layers of rubber (rr=3.2) with inner radius 0.40 cm and polystyrene (rp=2.6)
with inner radius p are used as insulating materials. Design the cable (i.e., determine p
and the inner radius of the outer cylindrical shell (outer radius of the polystyrene) o) that
is to work at a voltage rating of 20 kV. To avoid breakdown due to voltage surges, the
maximum electric fields in the insulating materials are not to exceed 25% of their
dielectric strengths. Assume that the dielectric strengths of rubber and polystyrene are 25
MV/m and 20 MV/m, respectively.

Solution: Assume that the line charge on the inner conductor is l C/m. In the dielectrics,
maximum fields occur at the respective inner radii. Therefore, the maxm fields Emax1 and
Emax2 in inner and outer insulators are given by

E max 1 =

2 3.2 o

and E max 2 =

= 0.004

2 2.6 o

a .
=p

Now using Emax1=0.25x25MV/m and Emax2=0.25x20MV/m in the above equations, we

have p=0.616 cm.
However, if V1 and V2 are PDs across the insulators, then they are given as
i
i
la
l
V1 = E1 dl =
d a =
ln( p / i )
2

3
.
2

3
.
2

o
o
=p
=p

l
2 o

1
1

3.2 ln( p / i ) + 2.6 ln( o / p ) = 20000V

Now finding the value of (l /20)=8x104 form the equation Emax1 we obtain o=0.832
cm from the above equation.

V1 + V2 =

Example 1.41
A cylindrical capacitor consists of an inner conductor of radius a
and an outer conductor whose inner radius is b. The space between the conductors is
filled with a dielectric of permittivity , and the length of the capacitor is L. Determine
the capacitance of the capacitor.
Solution: In the region ab, E can be found using Gausss law as E=Q/(2L) a. V
is then calculated from
a
a
a
Qa
Q
Q
b
d a =
V = E dl =
ln =
ln
2L
2L a
b
=b
=b 2L
Q 2L
=
b
V
ln
a
Example 1.42
A spherical capacitor consists of an inner conducting sphere of
radius Ri and an outer conductor with a spherical inner wall of radius Ro. The space in
between is filled with a dielectric of permittivity . Determine the capacitance. Assume a
charge Q is distributed on the surface of the inner sphere.
Solution: In the region Ri r Ro, E is found using Gausss law as E=Q/(2r2) ar. V is
then calculated from
R
ri
ri
Qa r
Q i
Q 1
1
V = E dl =
dra r =
=
( )
2
4r Ro 4 Ri R0
r = Ro
r = Ro 4r
C =

C =

Q
4
=
1
1
V
( )
Ri R0

1.16 Poissons and Laplaces Equations

Using E= - V and D=E in D=v, we have
V =
2V= - v/.

(30)

Equation (30) is known as Poissons which becomes Laplaces equation in absence of

free charge v. Now
2V 2V 2V
2V = 2 + 2 + 2 in Cartesian coordinate system
x
y
z
2V =

1
V
1 2V 2V
(
)+ 2
+

2 z 2

2V =

1 2 V
1

V
1
2V
(
r
)
+
(sin
)
+

r 2 r
r 2 sin
r 2 sin 2 2

in Cylindrical coordinate system and

in Spherical coordinate

system.
In most practical problems, the exact charge distribution is not known everywhere, so the
previous equations obtained from Coulombs law cannot be applied for finding the
potential and the field intensities. Sometimes the potentials of the conducting bodies may
be known and we wish to find the potential and field intensity in the surrounding space.
In such cases, the above differential equations must be solved subject to the appropriate
boundary conditions. These problems, therefore, are named by boundary value
problems.
Example 1.43
Two large parallel conducting plates are separated by a distance d
and maintained at potentials 0 and V0 as shown in the following figure. Assuming
negligible fringing effect at the edges, determine the potential at any point between the
plates and the surface charge densities on the plates when (a) the region between the
plates is empty and (b) the region between the plates is filled with continuous distribution
of charge having a volume charge density v= -0 y/d.
y

d
V0
0
Solution:
(a) We have to solve Laplace equation 2V=0.
2V 2V 2V
2
Now in rectangular coordinate system, V = 2 + 2 + 2
y
z
x
However, for this particular problem where V is only function of y, we have

2V
=0
y 2

V
= C1 V = C1 y + C 2
y
Now C1 and C2 are to be determined from the boundary conditions: V=0 ay y=0 and
V=V0 at y=d. Substituiting these conditions in the above equation of V, we have C2=0
and C1=V0/d. Therefore,
V
V = 0 y.
d
Now as the electric field and potential are related by E= -
V, threfore we have

V
V
V
V
V
ax +
ay +
a z ) = ( 0 y )a y = 0 a y .
x
y
z
y d
d
However, for parallel plate capacitor, |E| in between the plates is given by s/0.
V
V

Therefore s = 0 s = 0 0
0
d
d
As E is directed in the y direction, the upper plate is charged with positive charge i.e., s
and the lower plate is charged with negative charge i.e., -s.
(b) We have to solve Poissons equation 2V=-v/0.
2V 2V 2V
In rectangular coordinate system, 2 V = 2 + 2 + 2
y
z
x

E = -V = (

. However, for this particular problem where V is only function of y, we have

2V
= v
2
0
y

0 y 3
2V 0 y
V 0 y 2
=

=
+
C

V
=
+ C1 y + C 2
1
0d
y 2 0 d
6 0 d
y 2
C1 and C2 are to be determined from the boundary conditions: V=0 ay y=0 and V=V0 at
y=d. Substituting these conditions in the above equation of V, we have C2=0 and
V

C1 = ( 0 0 d ) . Therefore,
d 6 0
V

V = 0 y 3 + ( 0 0 d) y
6 0 d
d 6 0
Now as the electric field and potential are related by E= -
V, threfore we have
V

V
V
V
0 3
E = -V = (
ax +
ay +
az ) = (
y + ( 0 0 d ) y )a y
x
y
z
y 6 0 d
d 6 0

= (

0 2 V0 0
y +

d)a y .
2 0 d
d 6 0

At the lower plate

V

E = ( 0 0 d)a y . Using Gauss' s law at the lower plate 0 E Sa y = sl S

d 6 0
V0 0

d) = sl
d 6 0
At the upper plate

V
E = ( 0 d + 0 )a y .. Using Gauss' s law at the upper plate 0 E S (a y ) = su S
3 0
d

0 (

0(

V
0
d + 0 ) = su
3 0
d

Here sl and su are the surface charge densities on the lower and upper plates.
Example 1.44
A high voltage coaxial able consists of a single conductor of radius
Ri and a cylindrical metal sheath of radius 0 (0>i) with a homogeneous insulating
material between the two. Show that the potential distribution in the dielectric is

Vs ln( 0 / )
and hence determine the capacitance per unit length of the cable.
ln( 0 / i )
Assume Vs is the voltage applied on the inner conductor and the the outer metal sheath is
grounded.
Solution:
V=

i
V= Vs

0
V=0

We have to use Laplace equation in cylindrical coordinate system; therefore,

1

1 2
2
(
V)+ 2
V
V =0
2V = 0
+

2
z 2
1

As V is only function of , we have

(
V) = 0

C

V = C1
V = 1 Integrating this equation w.r.t we have

V = C1 ln + C 2
Now from the boundary conditions V = 0 at = 0 and V = Vs at = i we have
C1 =

Vs
Vs
and C 2 =
ln 0
ln( 0 / i )
ln( 0 / i )

V =

Vs
Vs
V ln( 0 / )
ln +
ln 0 = s
ln( 0 / i )
ln( 0 / i )
ln( 0 / i )

Now E = -V = -(

E=-

V ln( 0 / )

a +
a + a z ) s

z
ln( 0 / i )

Vs
Vs ln( 0 / )

a =
a
ln( 0 / i )
ln( 0 / i )

Vs
1
1
2

We = E dv =
2 v
2 ln( 0 / i )

2 2 L
0

Vs

dz
=

L
0 0 2
ln( / ) ln( 0 / i )
0
i

i
Here We is the energy stored in the cable of L meter length. If C is the capacitance of the
1

cable per unit length, then We is again given by We =

2

Vs
1
ln( 0 / i )
CLVs2 = L
2
ln( 0 / i )
C =

2
ln( 0 / i )

F/m

1
CLVs2 .
2

Example 1.45
A coaxial cable have a core diameter of 2 cm and an insulator
thickness of 0.5 cm for which r=3.5. Calculate the capacitance, voltage distribution in
the dielectric and the electric field at the core surface that is at a voltage of 1 kV with
respect to the outer shield.
Solution: Cross-sectional view of the cable is shown in the following:
0.5 cm

2 cm
V0=1 kV
r=3.5
In the insulator
1
1 2

2
(
V) + 2
V
+
V =0

2
z 2
Due to symmetry, V is only function of r, therefore
1

(
V) = 0
V = C1 V = C1 ln + C 2

For the boundary conditions, V = 0 at = 1.5 cm and V = 1 kV at = 1 cm we have
C 2 = 1 000 and C1 = 2466.3; therefore V = 2466.3 ln + 1000. Here in cm and V is in volts.

2466.3
E = V = -( a +
a + a z )V =
a

2V = 0

Therefore at the surface of the core E = 2466.3 Volts/cm

l
2466.3
2 = 2466.3 2
a l =

2
= 2466.3 2 0 3.5 = 0.48 C/m
Now for line charge E =

C=

l
V

0.48 10 6
= 0.48 10 9 F/m
1000

1.17 Theory of Images

In many practical problems, charge is placed near a conducting surface that induces
charges on it. The calculation of total electric field in such a region needs to consider
their effects also. For example, the effect of the earth on the fields from an open wire
transmission line cannot be ignored. Similarly, the field patterns of transmitting and
receiving antennas are greatly modified by the presence of the conducting bodies on
which they are mounted. To consider the influence of a nearby conductor on a field we
must know the charge distribution on the surface of the conductor that depends upon the
fields just above its surface. However, in the case of static fields, we know that (a) a
conductor forms an equipotential surface, (b) there are no fields inside an isolated

conductor, and (c) the fields are normal to the surface of a conductor. These observations
will help to quantify the charge distribution on the surface of a conductor and its
influence on the fields in the region.
In the calculation of potential and electric field of an electric dipole, we found that the
potential on the bisecting plane is zero, and the electric field is normal to the plane.
Therefore, the bisecting plane satisfies the requirements of a conducting plane. In other
words, if a conducting plane is inserted to coincide with the bisecting plane, the field
pattern of the dipole remains unchanged. If the negative charge below the conducting
plane is removed, the field distribution in the region above the plane still remains the
same, and the total charge induced on the surface of the conductor is Q as shown in the
following figure. Conversely, if we are given a point charge Q at a distance d above a
conducting plane of infinite extent, we can determine the potential and the electric field at
any point above the plane by ignoring the plane and imagining a charge Q at the same
distance away from the other side of the plane. However, these statements are only true
for a conducting plane of infinite extent and depth. For a curved surface, the imaginary
charge is neither equal in magnitude nor as far away on the other side of the conducting
surface. The imaginary charge Q is said to be the image of the real charge Q. Thus, in
the method of images, a conducting plane is temporarily ignored, and an imaginary
charge is placed behind the plane.

(a)

Point charge near an infinite grounded conducting plane

P

-Q

O(x, y)

P
Let PP be an earthed conducting infinite plane. Consider Q at a distance d on the right
side of the plan. Whatever may be the charge induced on the plane, it is always at zero
potential because it is earthed. It is obvious that if we remove the grounded conductor and
replace it by a point Q at a distance d to the left of the plane, then every point on the
plane will be equidistance from charges and will be at zero potential. Now these two
charges will give the proper solution of the problem. The charge Q is known as the
image of the charge Q. Now the potential at a point O(x,y) is given by

V =

Q
4 (d x) 2 + y 2

Q
4 (d + x) 2 + y 2

2
2

Q
1
1
(d x) + y

E = V =
a

a y
x
2
2
1
4 x (d x) 2 + y 2
y
(
d
+
x
)
+
y

2
2

(d + x) + y

Q
dx
d+x

=
+
2
2 3/ 2
4 {(d x) + y }
{( d + x) 2 + y 2 }3 / 2

2
2 3/ 2
{( d x) + y }

a
a x

y
y

2
2 3/ 2
{( d + x) + y }

Now for x=0 i. e., on the grounded conducting plane D=E=s

dQ
s =
2
2 (d + y 2 ) 3 / 2
(b)

Point charge near a grounded conducting sphere

P
r4
r3
r1

a
Q

r2 a

Q
P2

r
P3

P1

b
Let us try to find the position and magnitude of the image charge Q that will make the
potential on the spherical surface zero. The zero potential at P1 and P2 give us
Q
Q
+
= 0 and
d +a a+b

Q
Q
+
=0
d a ab

Solving these two equations gives Q= - aQ/d and b=a2/d.

Let us now calculate potential of P3 and hence first prove that the image charge Q and its
position as determined above are valid solution of this problem.
VP 3 =

Q
4r1

Q
Q
1
a/d

4r2 4 d 2 + a 2 + 2da cos

b 2 + a 2 + 2ba cos

Q
1
1

=0

2
2
4 d 2 + a 2 + 2da cos
a + d + 2ad cos

Since the original point charge and the image jointly satisfy the boundary condition, they
must give the correct field at every point in the space outside the sphere. The potential at
P is
=

Q
Q
1
a/d

2
2
2
2
4 r3 4 r4 4 d + r + 2dr cos
b + r + 2br cos

Q
1
a/d

2
2
4 d 2 + r 2 + 2dr cos

a
+
(
rd
/
a
)
+
2
rd
cos

VP =

Special Note: If the sphere is at a potential other than zero, we may still determine the
voltage or electric field by the method of images. We first replace the conducting sphere
with an image charge Q, as we did for a grounded sphere. This makes the surface
occupied by the sphere an equipotential. We next add a second image charge at the center
to raise the spherical surface to the required potential. If we are given a sphere of radius a
with a charge Q on it, we can replace the sphere by an image charge Q= - aQ/d at a
distance b=a2/d from the center, plus a charge (Q- Q) at the center.
(c)

Line charge near a parallel conducting cylinder

P3
r1

a
l

-l
P2

r2 a

P1

b
To make the surface equipotential, the image must be a parallel image line charge inside
the cylinder; and because of symmetry with respect to the line joining the charge l and
the center of the cylinder, the image line charge must lie somewhere on this line. If we
assume the image as -l, its position at a distance b from the center can be determined by
equating the potential of P1 with that of P1. We, therefore, have

l
r
l
r
l
r
l
r
ln 0 +
ln 0 =
ln 0 +
ln 0
2 0 d + a 2 0 a + b 2 0 d a 2 0 a b
a+b ab
a2
=
b=
d +a d a
d
Now the potential at P3 is

r
r
r
l
VP3 = l ln 0 +
ln 0 = l ln 0
2 0 r1 2 0 r2 2 0 r1

r1 = (d + a cos ) 2 + (a sin ) 2 = d 2 + a 2 + 2ad cos

and

r2 = (b + a cos ) 2 + (a sin ) 2 = b 2 + a 2 + 2ab cos =

a
d 2 + a 2 + 2ad cos
d

VP3 =

l
a
ln
2 0 d

(31)
Example 1.46
Determine the capacitance per unit length between two long,
parallel circular conducting wires of radius a. The axes of the wires are separated by a
distance D.

Solution:
The above problem can be realized by the following system:
2
-l
D

l
b

The equipotential surfaces of the wires can be realized if b=a2/d. Now if V1 and V2 are
the potentials at the surfaces of the wires 1 and 2, then using (31) we have

a
d
a
d
V1 = l ln = l ln
and V2 = l ln = l ln
2 0 d 2 0 a
2 0 d
2 0 a

V1 V2 =

l
d
ln
0 a

Again d = D - a 2 / d

C =

d=

l
V1 V2

0
d
ln
a

F/m

1
( D + D 2 4a 2
2

C is also given by C =

F/m
ln[( D / 2a) + ( D / 2a) 2 1] cosh ( D / 2a)
Example 1.47
A transmission line of radius a and having charge density l C/m is
situated at a height h above and parallel to the surface of the infinite ground plane.
Determine the capacitance per unit length of the line between the line and the ground
plane.
Solution:
1

a
l
h

l
h
Cg

Ground
h

Cg
-l

(a)

(b)

The problem is depicted in fig. (a). The conducting ground plane can be replaced by
considering an image line at a depth h below the ground plane. The problem then
becomes similar to that in Example 1.46 except D=2h. The capacitance between the
original line and the its image is therefore given by

C=

F/m
1
ln[(h / a) + (h / a) 1] cosh (h / a)
Now the capacitance between the line and the ground is, therefore,
2 0
C g = 2C =
F/m
cosh 1 (h / a)
Example 1.48 A point charge q is located at (0, 0, d) above the surface of a conducting
plane of infinite extent and depth. Determine the potential and electric field at P.Also
show that the total charge induced on the conducting surface is q.
Solution: The problem is depicted in the following figure. Following the theory of
images, to determine the fields at any position above the plane we have to place an image
charge q at (0, 0, -d) and ignore the existence of the plane.
Now the potential at P (x, y, z) is
q 1
1
2
2
2 1/2
2
2
2 1/2
V=
, where R1=[x +y +(z-d) ] and R2=[x +y +(z+d) ]
4 0 R1 R2
Now the electric field E is
z + d z d
1
q 1
1
1
E = V = x 3 3 a x + y 3 3 a y + 3 3 a z
4 0 R2 R1
R1
R2 R1
R2
2qd
On the surface of the conducting plane E= a z , where R=[x2+y2+d2]1/2.
4 0 R 3
2

Form the boundary condition at the interface of free space and conducting body we have
the surface charge density induced on the conducting surface must be equal to the normal
component of D. Therefore,
2qd
s = . Thus the total charge induced on the conducting surface of infinite extent
4R 3
2

dd
d
2qd
is Q = s ds = = qd
= q

2
2 3/ 2
2
2 3/ 2
4 =0 =0 ( + d )
s
=0 ( + d )

1.18 Equation of Continuity

Consider an arbitrary volume V bounded by a closed surface S contains a charge Q. If a
current I flow across the surface out of this region, the charge in the volume decrease at a
rate that equals the current. The current leaving the region is the total outward flux of the
current density vector through the surface. We, therefore, have
dQ
d
d
I=
= v dv Again I = J ds J ds = v dv
dt
dt v
dt v
s
s
Now with the help of divergence theorem we have

J ds = ( J)dv = dt
s

dv J =

d v
Continuity equation.
dt

(32)

As J = E , From eq. (32) we have

d
d v
d

v = v
+ v = 0
E = v

dt
dt
dt
The solution of the above equation is
v = v 0 e ( / )t C/m3. The initial charge v0 will decay to 1/e or 36.8% of its value in a
time equal to =/. The time constant is called the relaxation time. We are now in
position to explain why the charges introduced into the conducting material will move to
its surface and redistribute themselves in such a way as to make =0 and E=0 inside
under equilibrium conditions.
Now for dc current J = 0 that leads to J ds =0 . If the closed surface shrinks to a
s

point, we have the Kirchhoffs current law I = 0 . Putting J=E in the continuity
equation for dc current, we obtain
J = 0 E = 0 E + E = 0 . But for homogeneous medium =0,
therefore we obtain
E = 0 (V ) = 0 2V = 0
The equation 2V=0 states that the potential distribution within the conducting body
satisfies Laplaces equation as long as the medium is homogeneous and the current
distribution is time invariant.

1.19 Boundary Conditions for Current Density

Boundary condition for

tangential component

Boundary condition for

normal component

We study here how the current density J changes when passed through the interface
between two media of different conductivities. The boundary condition for tangential
components of static electric fields, Et1= Et2, relates Jt1 with Jt2, and the continuity
equation
J ds =0 determines the relation between the normal components of current densities
s

in two different media.

Using J=E in Et1= Et2, we have

J t1 1
=
(33)
1 2
Jt2 2
where 1 and 2 are respectively, the conductivities of media 1 and media 2.
Assuming the height of the pillbox shown in the above figure is small so that the
contribution from the current leaving or entering through the curved surface is negligible.
Now from the integration of the continuity equation over the closed surface s of the
pillbox with h0, we obtain
J n1a n sa n J n 2 a n sa n = 0 J n1 = J n 2 .
(34)
From (33) and (34) we obtain
2 J t1 J n1
J J
J /J
tan 1 1

=
t1 n 2 = 1 t1 n1 = 1
=
(35)
1 J t 2 J n2
J t 2 J n1 2
J t 2 / J n2 2
tan 2 2
Let us consider a special case where mediums 1 and 2 are respectively, poorly and highly
conducting. As is directly proportional to the conductivity, 1 will be very small
because 2 >>1. J in medium 1 is therefore almost normal to the interface. This leaves
the tangential component of J and hence that of E negligibly small in medium 1. Smaller
value of the tangential electric field in medium 1 leaves the tangential component of E in
medium 2 small as Et1= Et2. Now from (34), the normal component of E in medium 2
given by
J t1

En2 =

Jt2

1
E is also very small. This means that the electric field in a highly conducting
2 n1

medium is practically zero. As the medium 2 is a conducting body, there must exist a free
surface charge density at the interface we can be calculated from (26) as

s = Dn1 Dn 2 = 1 E n1 2 E n 2 = 1 E n1 2 1 E n1 = E n1 1 2 1 = J n1 1 2
2
2

1 2
Similar expression can be obtained in terms of Jn2 as

s = Dn1 Dn 2 = 1 E n1 2 E n 2 = 1 2 E n 2 2 E n 2 = E n 2 1 2 2 = J n 2 1 2 .
1
1

1 2

1.20 Leakage Resistance between Two Conductors

Whatever be the shapes of the conductors, the capacitance between the conductors is
given by
s D ds s E ds
Q
(36)
=
C= =
V - E dl - E dl
l

Now if the dielectric between the conductors is loosy, its resistance is given by
- E dl - E dl
V
R= = l
= l
I
J ds E ds
RC =

R=

(37)

(38)

Example 1.49
Two infinitely conducting parallel plates, each of cross-sectional
area A, are separated by a distance d. The potential difference between the plates is Vab.

If the medium between the plates is homogeneous and has a finite conductivity ,
determine the resistance of the region between the plates and the capacitance.
Solution: Let us consider the following problem configuration.

d 2V
The potential distribution may be expressed as
= 0.
dz 2
Integrating the above equation twice, we obtain V=az+b where a and b are constants and
can be obtained from the boundary conditions that V=0 at z=0 and V=Vab at z=d. Now
a=Vab/d and b=0. Therefore, the potential in the medium between the plates is given by
V=(Vab/d)z.
V
d V
E = V = ( ab z )a z = ab a z
d
dz d
Vab
J = E =
az
d
Now the current through a surface of area A normal to J is
V A
I = J ds = ab
d
V
d
Therefore, R= ab =
I
A

A
=
where is the permittivity of the medium.
Now following (38), we obtain C =
R
d
Example 1.50 The conductivity of the medium between two concentric metal spheres is
and the permittivity is . If the radius of the inner sphere is a and the inner radius of the
outer sphere is b, determine the resistance of the medium between the spheres.
Solution: Following the result of Example 1.42, Capacitance of the system is given by
4ab
C=
.
ba

(b a )
Now using (38), we have R =
=
.
C 4 ab
Example 1.51
A two-wire transmission line consists of wires of radius a. The
wires are separated by a distance D in a lossy media of permittivity and conductivity,
respectively, and . Determine the leakage resistance between the conductors for unit
length of the line.
Solution: The solution has two steps. In a first step, the capacitance per unit length of the
line between the wires is to be determined which is explained in Example 1.46. C is
given by

C=

ln[( D / 2a) + ( D / 2a) 2 1]

cosh ( D / 2a)

F/m.

Therefore the leakage resistance per unit length of the line is

R=

1
1
=
cosh 1 ( D / 2a ) =
ln[(D / 2a ) + ( D / 2a ) 2 1] m.
C

The total resistance of the line for a length of l meter is R/l.

Example 1.52
A coaxial cable consists of inner conductor of radius a and outer
conductor of inner radius b. The space between the conductors is a lossy media of
permittivity and conductivity, respectively, and . Determine the leakage resistance
between the conductors for unit length of the cable.
Solution: The solution has two steps. In a first step, the capacitance per unit length of the
cable is to be determined which is explained in Example 1.41. C is given by
2
C=
F/m.
b
ln
a
Therefore the leakage resistance per unit length of the cable is
R=

1
=
ln(b / a ) m.
C 2

The total resistance of the cable for a length of l meter is R/l.

1.21 Summary

The electric fields for any kind of charge distribution are found to obey the
following relations which are known as the fundamental postulates of
electrostatics:
E = v / E = 0.
In presence of source v, the electric field is not solenoidal but irrotational.
However, in a source free region the electric field is solenoidal as well as
irrotational.
No free charge stays within a conducting body; the electric field is thus zero
everywhere in it.
Electric field is always normal to the conducting surface which leaves the
conducting surface as equipotential.