Mpodule 2 Emw

Engineering Electromagnetics 17EC36 Module 2
ENGINEERING ELECTROMAGNETICS [17EC36]

B.E., III Semester, Electronics & Communication Engineering/
Telecommunication Engineering (VTU)
[As per Choice Based Credit System (CBCS) scheme]
Module-2
Gauss’s law and Divergence
Gauss’s law, Divergence, Maxwell’s First equation (Electrostatics), Vector operator
ai
) and
Divergence theorem
lek
Energy, Potential and Conductors
Energy expended in moving a point charge in an electric field, The line integral, Definition of
potential difference and potential, The potential field of point charge, Current and Current density,
Continuity of current.
Text Books:
Ba
1. W.H. Hayt and J.A. Buck, “Engineering Electromagnetics”, 7th Edition, Tata
ja
McGraw-Hill, 2009, ISBN-978-0-07-061223-5.
ite
Reference Books:
1. Ian John Krauss and Daniel A Fleisch, “ Electromagnetics with applications”, Mc Graw Hill.
av
2. N. Narayana Rao, “Fundamentals of Electromagnetics for Engineering”, Pearson.
PREPARED BY:
R
RAVITEJA BALEKAI
Asst Professor
ECE Dept, GMIT
Davangere 577006
Cell: +919739223504
Mail: [email protected] , [email protected]
Website: https://ravitejb.wixsite.com/ravitej
Raviteja B, Asst. Prof., Dept. of E&CE 1

Gauss’s Law
Statement: “The electric flux passing through any closed surface is equal to the total charge enclosed by
that surface”.
Mathematical representation of Gauss’s law
Consider any object of irregular shape as shown in figure. The total charge enclosed by the irregular
closed surface is Q coulombs. It may be in any form of distribution. Consider a small differential surface
dS at point P. The flux density at point P is D and its direction is such that it makes an angle θ with the
ai
normal direction at point P. The flux d passing through the surface dS is the product of the component of
D in the direction normal to the dS and d
lek
Mathematically this can be represented as
d = Dn dS
where Dn = Component of D in the direction of normal to the surface dS.

Ba
dS Closed irregular surface
ja
ɵ
ite
Q
av
Fig.1 Flux through irregular shaped closed surface

R
Dn = cos θ
d = cos θ dS
From the definition of dot product,
= cos θAB
d =

This is the flux passing through incremental surface area dS. Hence the total flux passing through the entire
closed surface is expressed as
= =
Now irrespective of the shape of the surface and charge distribution, total flux passing through the surface
is the total charge enclosed by the surface.
= =Q
Problems
ai
1. Given the electric flux density, D = 0.3r2 ar nC/m2 in free space
lek
a) Find E at point P(r =2 , θ=25° , ϕ=90°).
b) Find the total charge within the sphere r=3.
c) Find the total electric flux leaving the sphere r=4.
a)
D = 0.3r2 ar nC/m2
D = E εo
Ba
E= = ar =135.5 ar V/m
ja
b) = r=3
ite
Q=D = 305.40 C
av
c) =
= r2 sinθ dθ dφ
= 965.22 nC
R

Point Charge
ϕ r =a Ds
Q
ɵ
ai
lek
Fig.2 Application of Gauss’s law to the field of a point charge Q on spherical closed surface of radius ‘a’.
The electric flux density D is everywhere normal to the spherical surface and has a constant
magnitude at every point on it. Ba
The electric field intensity of the point charge has been found to be
E= - - - - - - - - - (1)
ja
Since D= E we have D = - - - - - - - - - (2)

ite
At the surface of the sphere, Ds=

av
The differential element of area on a spherical surface is in spherical coordinates

dsr = r dθ r sin θ dφ ar = a dθ a sin θ dφ ar = a2 sin θ dθ dφ ar
Ds . dS = . a2 sin θ dθ dφ ar (
R
Leading to the closed surface integral
sinθ dθ dφ = Q
= - - - - - - - - - (3)

This proves the Gauss’s law that Q coulombs of flux crosses the surface if Q coulombs of charge is
enclosed by that surface.
Infinite Line Charge
ai
L dS
lek
x -z
Ba
Fig. 3 Infinite line charge
Consider an infinite line charge of density C/m lying along z-axis from - ꚙ to ꚙ. Consider the Gaussian
surface as the right circular cylinder with z-axis as its axis and radius ‘r’ as shown in the figure. Length of
ja
the cylinder is ‘L’ mtr. The flux density at any point on the surface is directed radially outwards i.e in the
direction according to cylindrical coordinate system.
ite
As the line charge is along z-axis, there cannot be any component of in z direction . So has only
radial component.
av
Now, Q= - - - - - - - - - (1)
Q= + + = - - - - - - - - - (2)
R
= =0
Q= = . r dϕ dz = . r dϕ dz
Q = 2 r Dr L - - - - - - - - - (3)

Dr =
But =
Dr = (Due to infinite line charge)
E= V/m - - - - - - - - - (4)
ai
Coaxial Cable
lek
a b
L
Ba
ja
Fig. 4 Coaxial cable

ite
Consider a coaxial cable as shown in figure. The radius of inner conductor is ‘a’ while the radius of the
outer conductor is ‘b’. The length of the cable is ‘L’. The charge distribution on the outer surface of the
inner conductor is having density C/m2. The total surface area of the inner conductor is 2 aL.
av
Hence can be expressed interms of

R
= =
=2 a - - - - - - - - - (1)
Thus line charge density of inner conductor is C/m.
Consider the right circular cylinder of length L as the Gaussian surface. Due to the symmetry, has only
radial component.

Q = 2 r Dr L (From infinite Line charge) where a < r < b - - - - - - - - - (2)
Total charge on the inner conductor is to be obtained by evaluating the surface integral of the surface
charge distribution.
Q= dS = r dϕ dz
Q= . a dϕ dz = 2 a L - - - - - - - - - (3)
Equating (2) and (3)
ai
2 r Dr L = 2 a L
lek
Dr = - - - - - - - - - (4)
From eqn (1) =

Ba
D=
E= (a < r < b ) V/m - - - - - - - - - (5)

ja
This is same as obtained for infinite line charge. Every flux line starting from the positive charge on the
ite
inner cylinder must terminate on the negative charge on the inner surface of the outer cylinder. Hence the
total charge on the inner surface of the outer cylinder is
av
Qouter cylinder = - 2 a L
But Qouter cylinder = 2 b L

R
2 bL =-2 aL
=- - - - - - - - - - (6)

Gauss’s law applied to Differential volume element
D = Do = Dxo ax + Dyo ay +Dzo az
ai
Fig. 5 Differential volume element
lek
Consider a closed Gaussian differential surface in the form of rectangular box, which is a
differential volume element. The sides of this element are ∆x , and . Let us consider a point P
located in a rectangular coordinate system as shown in figure. The value of D at the point P may be
expressed in rectangular components,
Ba
Do = Dxo ax + Dyo ay +Dzo az - - - - - - - - - (1)
The components of Dxo ,Dyo and Dzo vary with distance in the respective directions.
ja
According to Gauss’s law,
ite
Q= - - - - - - - - - (2)
={ - - - - - - - - - (3)
av
Consider the front surface of the differential element. Though is varying with distance for small
surface like front surface it can be assumed constant.
R
And =
= Front .( ) - - - - - - - - - (4)
= Dx Front .( ) - - - - - - - - - (5)
Dx Front = Dx0 + [ Rate of change of Dx with x ] x [Distance of surface from P ]

Dx Front = Dx0 + - - - - - - - - - (6)
The point P is at the centre so distance of surface in x direction from P is ]
= Dx0 + ] - - - - - - - - - (7)
Consider the integral over back surface
= Back .( - - - - - - - - - (8)
ai
= Dx Back .( - - - - - - - - - (9)
lek
[ The flux is entering from back side and leaving from front in positive x direction hence is used
positive Back . While the surface is considered from point P is in negative x direction hence - is used
for expressing ]
Ba
Dx Back = Dx0 - [ Rate of change of Dx with x ] x [Distance of surface from P ]
Dx Back = Dx0 - - - - - - - - - - (10)

ja
= - Dx0 + ] - - - - - - - - - (11)
ite
Combining (7) and (11)
+ = - - - - - - - - - (12)
av
Similarly we can write,
+ = - - - - - - - - - (13)
R
+ = - - - - - - - - - (14)
From eqn (3)
= - - - - - - - - - (15)

= - - - - - - - - - (16)
Problems
1. In free space, let D = 8xyz4ax+4x2z4ay+16x2yz3az pC/m2. (a) Find the total electric flux passing
through the rectangular surface z = 2, 0 < x < 2, 1 < y < 3, in the az direction. (b) Find E at
P (2,−1, 3). (c) Find an approximate value for the total charge contained in an incremental sphere
located at P (2,−1, 3) and having a volume of 10−12 m3.
ai
D = 8xyz4ax+4x2z4ay+16x2yz3az pC/m2
lek
a) Q= Along z direction
Q=
Q=
Ba
Q = 1365 x C
b) E = =
ja
E at P (2, −1, 3)
ite
E= = (-146.4 ax+146.4ay-195.2az) v/m

av
c) Q= =
Q at P (2, −1, 3) = 10−12 m3

R
= = 8yz = - 648 x 10−12 = =0
= = = -1728 x 10−12
Q= = = - 2376 C

Divergence
Applying Gauss’s law to the differential volume element, we have obtained the relation,
Q= ΔV - - - - - - - - - (1)
This is the charge enclosed in the volume ΔV.
But by Gauss’s law - - - - - - - - - (2)
To apply Gauss’s law, we have assumed a differential volume element as the Gaussian surface, over which
ai
D is constant. Hence equations (1) and (2) can be equated in limiting case as ΔV →0.
lek
= ΔV
= = - - - - - - - - - (3)
This in general if is any vector then,

Ba
=
ja
This mathematical operation on is called a divergence and is given by
=
ite
div = - - - - - - - - - (4)
av
Divergence in different coordinate systems
. = div = (Cartesian)
R
. = div = (Cylindrical)
. = div = θ θ (Spherical)
θ θ θ

Problems
1. In each of the following parts, find a numerical value for div D at the point specified
a) D= (2xyz-y2)ax + (x2z-2xy)ay+x2yaz c/m2 at P(2, 3, -1).

b) D=2rz2sin2Φ ar + rz2 sin2ΦaΦ +2r2zsin2Φaz c/m2 at P(2, 1100, -1).
c) D=2r sinθ cosΦ ar + r cosθ cosΦ aθ - rsinΦaΦ c/m2 at P(1, 300, 500).
a) D= (2xyz-y2) ax + (x2z-2xy) ay+x2y az c/m2
ai
div =
lek
= = 2yz = = -2x = =0
div = = 2yz-2x At P(2,3,-1)
Ba
= 2 (3) (-1) – 2 (2) = -10 C/m3
b) D=2rz2sin2Φ ar + rz2 sin2Φ aΦ +2r2zsin2Φ az c/m2

ja
div =
ite
= = 4 z2 sin2Φ
= =
av
= =
R
div =
= 4 z2 sin2Φ + + At P (2,1100,-1)
= 3.53 - 1.532 + 7.064 = 9.062 C/m3

2
c) D=2r sinθ cosΦ ar + r cosθ cosΦ aθ - rsinΦ aΦ c/m
div = θ θ
θ θ θ
= =6
θ θ = θ =
θ θ θ θ θ
= =
θ θ θ
ai
div = θ θ
θ θ θ
lek
=6 + - At P (1,300,500)
θ θ
= 1.93 + 0.643 – 1.29 = 1.283 C/m3
Maxwell’s First Equation

Ba
The divergence of electric flux density is given by,
ja
div = - - - - - - - - - (1)
ite
According to Gauss’s law, it is known that
= - - - - - - - - -(2)
av
Expressing Gauss’s law per unit volume
=
R
Taking lim ΔV →0 i.e volume shrinks to zero,
= - - - - - - - - - (3)
But = at that point - - - - - - - - - (4)
The equation (4) gives the volume charge density at the point where divergence is obtained.

Equating (1) and (4),
div
.D= - - - - - - - (5)
This is volume charge density around a point. The equation (5) is called Maxwell’s first equation applied to
electrostatics. This is also called the Gauss’s law in point form or Gauss’s law in differential form.
Statement
“The divergence of electric flux density in a medium at a point (differential volume shrinking to zero), is
ai
equal to the volume charge density (charge per unit volume) at the same point ”.
lek
Divergence Theorem
From the Gauss’s law we can write,
Q=
Ba - - - - - - (1)
While the charge enclosed in a volume is given by,

ja
Q= - - - - - - (2)
ite
But according to Gauss’s law in the point form,
.D=
Q= - - - - - - (3)
av
Equating (1) and (3),
= - - - - - - (4)
R
The equation (4) is called divergence theorem.
Statement
“The integral of the normal component of any vector field over a closed surface is equal to the integral of
the divergence of this vector field throughout the volume enclosed by that closed surface ”.

Problems
1. Given that the field D = C/m2 . Find
a) Volume charge density b) The total electric flux leaving the surface of the spherical volume
of radius 2 m.
a) .D=
=0
.D = θ θ
ai
θ θ θ
θ Φ θ Φ
= = C/m3
lek
b) Ψ= =
θ Φ
= inθ dr dθ dϕ
Ψ=0C
Ba
2. Evaluate both sides of the divergence theorem for the field D = 2xyax + x2ay C/m2 and the
ja
rectangular parellelepiped formed by the planes x = 0 and 1, y = 0 and 2, and z = 0 and 3.
The divergence theorem states that
ite
Evaluating the surface integral first, we note that D is parallel to the surfaces at z = 0 and z = 3, so
av
D . dS = 0 there. For the remaining four surfaces,we have
= x=0 . (- dy dz ax ) + x=1 . (dy dz ax)

R
+ y=0 . (- dx dz ay ) + y=2 . ( dx dz ay )
= x=0 . (- dy dz ax ) + x=1 . (dy dz ax)
+ y=0 . (- dx dz ay ) + y=2 . ( dx dz ay )
= 12

= = + = 2y
= dx dy dz = 12
3. Given the field D = 6r sin ar + 1.5r cos aϕ C/m2, evaluate both sides of the divergence
theorem for the region bounded by r = 2, ϕ = 0, ϕ = π, z = 0, and z = 5.

ai
= r=0 . (-r dϕ dz ar ) + r=2 . (r dϕ dz ar )
lek
+ ϕ=0 . (- dr dz aϕ ) + ϕ= . ( dr dz aϕ )
= r=0 (-r dϕ dz ar )
+
Ba r=2 (r dϕ dz ar )
+ ϕ=0 (- dr dz aϕ )
ja
+ ϕ= ( dr dz aϕ )
= 225
ite
=
av
= + +0
= 12 sin - 0.75 sin

R
= dr dϕ dz
= 225

-r 2
4. Given D=30e ar - 2zaz c/m in cylindrical coordinates. Evaluate both sides of the divergence
theorem for volume enclosed by r=2, z=0, z=5.
+ +
Consider dS normal to ar direction which is for the side surface.
= r=2 . (r dϕ dz ar ) = 255.1
ai
= z=5 . (r dr d az ) = -40 π
lek
= z=0 . (r dr d az ) = 0
+ + = 129.43
=
Ba
= + + =- + –2
ja
= – dr dϕ dz
ite
= 129.43
av
R

Energy, Potential and Conductor
Work Done:
The electric field intensity is defined as the force on a unit test charge at that point at which we
want to find the value of E. Consider an electric field due to a positive charge Q. If a unit test positive
charge Qt is placed at any point in this field, it experiences a repulsive force and tends to move in the
direction of the force. But if a positive test charge Qt is to be moved towards the positive base charge Q
then it is required to be moved against the electric field of the charge Q. i.e against the repulsive force
exerted by charge Q on the test charge Qt. While doing so, an external source has to do work to move the
ai
test charge Qt against the electric field. This movement of charge requires to expand the energy. This work
done becomes the potential energy of the test charge Qt, at the point at which it is moved.
lek
Consider a positive charge Q and its electric field E. If a positive test charge Qt is placed in this
field, it will move due to the force of repulsion. Let the movement of the charge Qt is dl. The direction in
figure.
Ba
which the movement has taken place is denoted by unit vector , in the direction of dl. This is shown in
+ Qt
dL
ja
+Q
ite
Fig. 6
av
According to Coulomb’s law the force exerted by the field E is given by,
F = Qt E
R
But the component of this force exerted by the field in the direction of dl, is responsible to move the charge
Qt , through the distance dl. w.k.t the component of a vector in the direction of the unit vector is the dot
product of the vector with that unit vector. Thus the component of F in the direction of unit vector is
given by,
Fl = F. = Qt E . N

This is the force responsible to move the charge Qt through the distance dl, in the direction of the field.
To keep the charge in equilibrium, it is necessary to apply the force which is equal and opposite to the force
exerted by the field in the direction dl.
Fapplied = - Fl = - Qt E . N
In this case, the work is said to be done.
Thus there is expenditure of energy which is given by the product of force and the distance.
ai
Hence mathematically the differential work done by an external source in moving the charge Qt
through a distance dl, against the direction of field E is given by,
lek
dW = Fapplied x dl = - Qt E . dl
dW = - Qt . Joules
Ba
Thus if a charge Q is moved from initial position to the final position, against the direction of electric field
E then the total work done is obtained by integrating the differential work done over the distance from
initial position to the final position.
ja
W=
W= Joules
ite
The Line Integral

av
Consider that the charge is moved from initial position B to the final position A, against the electric field E
then the work done is given by,
R
W=
This is called the line integral, where gives the components of E along the direction .
Mathematical procedure involved in such a line integral is,
1. Choose any arbitrary path B to A.

2. Break up the path into number of very small segments, which are called differential lengths.
3. Find the component of E along each segments.

4. Adding all such components and multiplying by charge, the required work done can be obtained.
Thus line integral is basically a summation and accurate result is obtained when the number of segments
becomes infinite.
Consider an uniform electric field E. The charge is moved from B to A along the path shown in figure. The
path B to A is divided into number of small segments. The line integral from B to A can be expressed as
the summation of dot products.
W=-Q[ + + + ….._ ]
ai
But the electric field is uniform and is equal in all directions,
lek
= = = =
W=-QE.[ + + + …..+ ]
E2 E3
Ba
dL2 dL3 E4 A
ja
E1 dL4 E5 dL5
dL1
ite
Fig. 7
av
The sum of all such vectors is the vector joining initial point to final point.
R
W=-QE.
Key Point: The work done depends on Q, and does not depend on the path joining B to A.
This is true for nonuniform electric field as well.
The work done in moving a charge from one location B to another A, in a static, uniform or
nonuniform electric field E is independent of the path selected.

Differential length vector
Cartesian
= dx ax + dy ay + dz az
Cylindrical
= dr ar + r dφ a φ + dz az
Spherical
ai
= dr ar + r dθ aθ + r sinθ dφ aφ
lek
Problems
1. Find the work done in moving a charge of 2 C from (2, 0, 0) m to (0, 2, 0) m along the straight
line path joining two points, if the electric field is E= 12xax-4yay v/m.
Ba
W=
=
ja
= = 64 J
ite
2. Given the electric field E= (8xyzax + 4x2zay - 4x2yaz) v/m , find the differential amount of
work done in moving a 6nc charge a distance of 2μm, starting at P(2, -2, 3) and proceeding in
av
the direction aL a) - ax+ ay+ az b) ax - ay - az c) ax + ay
a) dW = - Qt E . dl
R
= - 6 x 10-9 (8xyzax + 4x2zay - 4x2yaz) . ( - ax+ ay + az ) (2 x 10-6)
Use Dot product operation
= - 12 x 10-15 ( - + - az )
At P (1, -2, 3)
dW = -149.34 x 10-15 J

b) dW = - Qt E . dl
= - 6 x 10-9 (8xyzax + 4x2zay - 4x2yaz) . ( ax - ay - az ) (2 x 10-6)

At P (1, -2, 3)
dW = 149.34 x 10-15 J
c) dW = - Qt E . dl
ai
= - 6 x 10-9 (8xyzax + 4x2zay - 4x2yaz) . ( ax+ ay) (2 x 10-6)
lek
= - 12 x 10-15 ( + )
At P (1, -2, 3)
Ba dW = 0
3. Calculate the work done in moving a 4 C charge from B (1, 0 ,0) to A (0, 2, 0) along the path
y = 2 – 2x, z = 0 in the field E a) 5 ax v/m b) 5 x ax v/m c) 5 x ax + 5 y ay v/m
ja
a)
W=
ite
W=
= -20 = 20 J
b)
av
W=
W=
R
= -20 = 10 J
c)
W=
W=
= -4 [ = - 30 J

Potential difference
The workdone in moving a point charge Q from point B to A in the electric field E is given by,
W=-Q - - - - - - (1)
If the charge Q is selected as unit test charge then from the above equation we get the workdone in moving
unit charge from B to A in the field E.
Thus “Workdone per unit charge in moving unit charge from B to A in the field E is called potential
difference between the points B to A”.
ai
Potential difference = V = - - - - - - - (2)
lek
NOTE: If B is the initial point and A is the final point then the potential difference is denoted as VAB
which indicates the potential difference between the points A and B and unit charge is moved from B to A.
Potential due to point charge

Ba
VAB = - - - - - - - (3)
Consider a point charge, located at the origin of a spherical coordinate system, producing E radially
ja
in all the directions as shown in the figure.
ite
Q rA A B
av
rB
R
Fig. 8 Potential due to point charge Q
Assuming free space, the field E due to a point charge Q at a point having radial distance r from
origin is given by
E= - - - - - - (1)

Consider a unit charge which is placed at a point B which is at a radial distance of rB from the
origin. It is moved against the direction of E from point B to point A. The point A is at a radial distance of
rA from the origin.
The differential length in spherical system is,
dL = dr ar + r dθ aθ + r sinθ dφ aφ - - - - - - (2)
Hence the potential difference VAB between points A and B is given by,
VAB = -
ai
VAB = - =- (dr ar + r dθ aθ + r sinθ dφ aφ )
lek
=
VAB =
Ba volts
Concept of Absolute Potential

ja
Instead of potential difference, it is more convenient to express absolute potentials of various points in
the field. Such absolute potentials are measured with respect to a specified reference position. Such a
ite
reference position is assumed to be at zero potential.
Most widely used reference which is used to develop the concept of absolute potential is infinity. The
potential at infinity is treated to be zero and all the potentials at various points in the field are defined with
av
reference to infinity.
Consider potential difference VAB due to movement of unit charge from B to A in a field of a point
R
charge Q. It is given by
VAB =
Now let the charge is moved from infinity to the point A i.e rB = , hence
VA = = - - - - - - (1)

Equation (1) is called potential of point A denoted as VA. This is also called Absolute potential of point A.
Similarly absolute potential of point B can be defined as
VB = - - - - - - (2)
This is the workdone in moving unit charge from infinity at point B.
Hence the potential difference can be expressed as the difference between the absolute potentials of the two
points.
ai
VAB = VA - VB - - - - - - (3)
lek
Thus absolute potential can be defined as,
“The workdone in moving a unit test charge from the infinity (or reference point at which potential is zero)
to the point, against the direction of the field”.
Ba
Thus, the absolute potential at any point which is at a distance r from the origin of a spherical system,
where point charge Q is located, is given by
V= - - - - - - (4)
ja
Potential due to several Point Charges

ite
Consider the various point charges Q1, Q2, Q3 …… Qn located at the distances r1, r2, r3 …… rn from the
origin as shown in the figure. The potential due to all these point charges, at point P is to be determined.
Use superposition principle.
av
R
Q1 r1
Q2 r2 P
Q3 r3
rn
Qn
Fig. 9 Potential due to several point charges

Consider the point charge Q1, the potential VP1 due to Q1 is given by,
VP1 =
The potential VP2 due to Q2 is given by,
VP2 =
ai
The potential VPn due to Qn is given by,
lek
VPn =
As the potential is scalar, the net potential at point P is the algebraic sum of the potentials at P due to
individual point charges, Ba
VP = VP1 + VP2 + …….. VPn
= + + …….
ja
ite
VP = volts
Problems
av
1. A point charge Q= 0.4nC is located at the origin. Obtain the absolute potential of A (2, 2, 3).
A (2, 2, 3)
R
V= Q= 0.4 nC(0, 0, 0)
= =
V= = 0.8719 volts

2. If Q = 0.4 nC is located at (2, 3, 3) then obtain the absolute potential of point A (2, 2, 3).
A (2, 2, 3)
VA = Q=0.4 nC (2, 3, 3)
= =1
VA = =3.595 volts
ai
3. If the point B is at (-2, 3, 3) in the above example. Obtain the potential difference between the
points A and B.
lek
VAB = VA - VB
Where VA and VB are the absolute potentials of A and B
Now VA =3.595 volts
A (2, 2, 3)
VB =
Ba Q=0.4 nC(2, 3, 3)
rB = Distance between point B and Q B (-2, 3, 3)

ja
= =4
VB = = = 0.8987 volts
ite
VAB = VA - VB = 2.6962 volts

av
4. If three charges 3μC, -4μC and 5μC are located at (0,0,0), (2,-1,3) and (0,4,-2) respectively. Find
the potential at (1,0,1) assuming zero potential at infinity.
R
Q1 (0,0,0) r1
Q2 (2,-1,3) r2 P(1,0,1)
Q3 (0,4,-2) r3
VP = V1 + V2 + V3

= + +
= = = =
= + + = 13.2 x 103 V
5. An electric field is expressed in rectangular coordinates by E = 6x2 ax + 6y ay + 4 az V/m.

Find a) VMN if points M and N are specified by M (2,6,-1) and N (-3,-3,2)
b) VM if V=0 at Q (4,-2,-35) c) VN if V=2 at P (1,2,-4)
ai
2
a) E = 6x ax + 6y ay + 4 az M (2,6,-1) and N (-3,-3,2)
VMN = - =-
lek
=- - - = -139 V
b) VM if V=0 at Q (4,-2,-35)
VMQ = - =-
Ba
M (2,6,-1)
=- - - = -120 V
ja
c) VN if V=2 at P (1,2,-4) (To calculate V at N , consider reference at point Q because V=0 at Q)
ite
VN = VNQ = - =-
=- - - = 19 V
av
Potential calculation when reference other than infinity

R
The expressions derived uptill now are under the assumption that the reference position of zero
potential is at infinity.
If any other point than infinity is selected as the reference than the potential at a point A due to point charge
Q at the origin becomes,
VA = +C
Where C = constant to be determined at chosen reference point where v=0.

Problem
1. A point charge of 6nC is located at the origin in free space, find potential of point P if P is located
at (0.2, -0.4, 0.4) and a) V= 0 at infinity b) V=0 at (1, 0, 0) c) V= 20V at (-0.5, 1, -1)
a) The reference is at infinity
VP = Q RP P
= = 0.6
ai
VP = = 89.87 V
lek
b) V=0 at (1, 0, 0) , Thus the reference is not at infinity. In such case potential at P is
Ba RP
P(0.2,-0.4,0.4)
VP = Q (0,0,0) RR R(1,0,0)
ja
Now VR at (1,0,0) is zero
VR =
ite
= =1
av
0= +C C = -53.92
VP = = 89.87 – 53.92 = 35.95 V

R

c) Now V = 20V at (-0.5, 1, -1) Let this is M (-0.5, 1,-1). The reference is not given as infinity.
VM = Rp P(0.2,-0.4,0.4)
= = 1.5 Q (0,0,0) RM
20 = +C C = -15.95 M(-0.5,1,-1)
VP = = 89.87 – 15.95 = 73.92V
ai
lek
Potential Gradient
Consider an electric field E due to a positive charge placed at the origin of a sphere
Ba
V=- =
The potential decrease as distance of point from the charge increases. It is known that the line integral of E
between the two points gives a potential difference between the two points. For an elementary length ΔL
we can write
ja
VAB = ΔV = -
ite
The rate of change of potential with respect to the distance is called the potential gradient.
= = Potential gradient
av
Relation between and V

R
Consider due to a particular charge distribution in space. The electric field and potential V is changing
from point to point in space. Consider a vector incremental making an angle θ with respect to the
direction of , as shown in figure.
To find incremental potential we use,
ΔV = -

=E = ΔL
Where = unit vector in the direction of ΔL
ΔL
P θ
ai
lek
Fig. 10 Incremental length at an angle θ
ΔV = - cos θ (using Dot product operation)
= Ba = Potential gradient
= - cos θ
ja
is maximum only when cos θ = -1 (θ = 180°) . This indicates that must be in the direction opposite
ite
to .
The maximum value of rate of change of potential with distance (dV / dL) is called gradient of V.
av
= - (grad V)
=- V
R
The gradient V in different coordinate systems
V= (Cartesian)
V= (Cylindrical)
V= θ (Spherical)
θ θ

Problem
1. Given the potential field in cylindrical coordinates V = r cos V, and point P at r=3m ,
=60° , z= 2m, find values at P for a) V b) E c) Magnitude of E d)
e) aN f) in free space
a) V = r cos at P ( r=3m , ϕ =60° , z= 2m)
V = 30 V
ai
b) =- V
lek
V= r cos r cos r cos
=-{ cos - - r z cos }
At P ( r=3m , ϕ =60° , z= 2m)

Ba
= - 10 + 17.3 + 24 V/m
ja
c) = = 31.23 V/m
ite
d) = | max = E = 31.23 V/m
e) E = -
av
= = = 0.32 - 0.55 -0.77

R
f) =?
D=E
= 8.854 x 10-12 { cos + + r z cos }

D= cos + + r z cos
. = div = =
r z cos
ai
=
lek
+ r cos
=
Ba
At P ( r=3m , ϕ =60° , z= 2m)
ja
= - 234 pC/m3
ite
av
R

Current and Current Density
The current is defined as the rate of flow of charge and it is measured in amperes.
I=
The current is considered to be the motion of the positive charges. The conventional current is due
to the flow of electrons, which are negatively charged. Hence the direction of conventional current is
assumed to be opposite to the direction of flow of the electrons.
ai
The current which exists in the conductors, due to the drifting of electrons, under the influence of
the applied voltage is called drift current. While in dielectrics, there can be flow of charges, under the
lek
influence of electric field intensity. Such a current is called the displacement current or convection current.
The analysis of such currents, in the field theory is based on defining a current density at a point in
the field. The current density is a vector quantity associated with the current and denoted as J.
Ba
The current density is defined as the current passing through the unit surface area, when the surface
is held normal to the direction of the current. The current density is measured in amperes per square meters
(A/ m2).
ja
Relation between I and J
ite
Consider a surface S and I is the current passing through the surface. The direction of current is
normal to the surface S and hence direction of J is also normal to the surface S.
av
J
R
dS
Fig. 11 J and dS are normal
I=

Continuity Equation
The continuity equation of the current is based on the principle of conservation of charge.
The principle states that, The charges can neither be created nor be destroyed.
Consider a closed surface S with a current density J, then the total current I crossing the surface S is
given by,
I= - - - - - - - - - (1)
The current flows outwards from the closed surface. It has been mentioned earlier that the current
ai
means the flow of positive charges. Hence the current I is constituted due to outward flow of positive
charges from closed surface S. According to principle of conservation of charge, there must be decrease of
lek
an equal amount of positive charge inside the closed surface. Hence the outward rate of flow of positive
charge gets balanced by the rate of decrease of charge inside the closed surface.
Let Qi = Charge within the closed surface

Ba
= Rate of decrease of charge inside the closed surface.
The negative sign indicates decrease in charge.

ja
Due to principle of conservation of charge, this rate of decrease is same as rate of outward flow of charge,
which is a current.
ite
I= = - - - - - - - - - (2)
av
This is the integral form of continuity equation of the current.
Using the divergence theorem, convert the surface integral in integral form to the volume integral.
= - - - - - - - - - (3)
R
= - - - - - - - - - (4)
But Qi =
= - - - - - - - - - (5)
For a constant surface, the derivative becomes the partial derivative.

= - - - - - - - - - (6)
If the relation is true for any volume, it must be true even for incremental volume ΔV.
=-
=- - - - - - - - - - (7)
This is the point form or differential form of the continuity equation of the current.
ai
Problem
– 4r ϕ mA/m2
lek
1. Given the vector current density J = 10 z
(a) Find the current density at P (r = 3, ϕ = 30◦, z = 2)
(b) Determine the total current flowing outward through the circular band r = 3, 0 < ϕ < 2π,
2 < z < 2.8.
a)
At P (r = 3, ϕ = 30◦, z = 2)
J = (10
Ba
z – 4r ϕ )x
J = (180 –9 )x A/ m2
ja
b) I= =
ite
.J=
av
–
=
R
. J = ( 30 r z + 8 sin ϕ cos ϕ ) x
I=
= ϕ
ϕ ϕ dr dϕ dz
I = 3.245 A

ASSIGNMENT
TOPIC: GAUSS'S LAW AND DIVERGENCE
1. State and prove and Gauss's law.

2. Show that s outer = (-a/b) s inner.
3. Calculate D in rectangular coordinates at point P(2,-3,6) produced by
a) a point charge Q= 55mc at Q(-2,3,-6). b) a uniform line charge L=20mc/m on the x-axis.
ai
c) a uniform surface charge density s=120μC /m2 on the plane z= -5m
4. Given the electric flux density, D= 0.3r2ar nc/m2 in free space. a) Find E at point P(2,250,900)
lek
b) Find the total charge within the sphere r=3.
c) Find the total electric flux leaving the sphere r=4.
Ba
5. In each of the following parts, find a numerical value for div D at the point specified
a) D= (2xyz-y2)ax + (x2z-2xy)ay+x2yaz c/m2 at P(2,3,-1).
b) D=2rz2sin2Φ ar + rz2 sin2ΦaΦ +2r2zsin2Φaz c/m2 at P(2,1100,-1).
c) D=2r sinθ cosΦ ar + r cosθ cosΦ aθ - rsinΦaΦ c/m2 at P(1,300,500).
ja
6. State Gauss's law in point form.
7. Determine an expression for the volume charge density associated with each D field following
ite
a) D= 5sinθ aθ + 5sinΦ aΦ c/m2 at P(0.5, /4, /4).

b) D= (4xy/z) ax + (2x2/z) ay +- (2x2y/z2) az.
av
c) D= z sinΦar + cosθ sinΦ aθ + rsinΦaz .

8. State and prove divergence theorem.
9. Given D=5r ar c/m2, prove divergence theorem for a shell region enclosed by spherical surface
R
at r=a and r=b (b>a) and centered at the origin.

10. Given D=30e-rar - 2zaz c/m2 in cylindrical coordinates. Evaluate both sides of the divergence
theorem for volume enclosed by r=2, z=0, z=5.

TOPIC: ENERGY, POTENTIAL AND CONDUCTOR
1. Find the work done in moving a charge of 2 C from (2,0,0)m to (0,2,0) m along the straight line
path joining two points, if the electric field is E= 12xax-4yay v/m.
2. Given the electric field E= (6xyzax + 2x2zay-4x2yaz)v/m , find the differential amount of work
done in moving a 3nc charge a distance of 5μm,starting at P(1,-2,3) and proceeding in the
direction aL a) -6/7 ax+2/7 ay+5/7az b) 6/7 ax- 4/7ay-5/7 az
3. Find V, E and D in free space a) V= Eoe-xsin( y/4) in cartesian at P(0,0,1) .
b) V= 100rcosΦ /(z2+1) at Q(4,500,3). c) V= 60sinθ/r2 at R(3,600,250)
ai
4. Derive the potential due to point charge.
5. Derive the potential due to several point charges.
lek
6. A point charge of 6nc is located at origin in free space, find potential of point P if P is located at
(0.2,-0.4,0.4) and a) V=0 at infinity b) V=0 at (1,0,0) c) V=20V at (-0.5,1,-1).
7. Derive an expression for J, v and velocity of the charge element.

Ba
8. With usual notations derive the expressions ᴧ.J= -∂ v /∂ t.
9. Find the current

a) J= 2x2ax 2xy3ay + 2xyaz A/m2 in outward direction of a cube of 1m, with one corner at the
ja
origin and edges parallel to the coordinate axes.
b) J= 10e-100 r A/m2 0.01≤ r ≤ 0.02 ,0 ≤ z ≤1m.
ite
c) J= 1 A/m2 r = 3m, 0 < θ < /6, 0 < Φ < 2 .

av
R

ai
lek
Ba
ja
ite
av
R

ai
lek
Ba
ja
ite
av
R

ai
lek
Ba
ja
ite
av
R

Mpodule 2 Emw

Uploaded by

Copyright:

Available Formats

Mpodule 2 Emw

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Mpodule 2 Emw

Uploaded by

Copyright:

Available Formats

Engineering Electromagnetics 17EC36 Module 2

ENGINEERING ELECTROMAGNETICS [17EC36]

Gauss’s law, Divergence, Maxwell’s First equation (Electrostatics), Vector operator

2. N. Narayana Rao, “Fundamentals of Electromagnetics for Engineering”, Pearson.

Raviteja B, Asst. Prof., Dept. of E&CE 1

Mathematical representation of Gauss’s law

where Dn = Component of D in the direction of normal to the surface dS.

Fig.1 Flux through irregular shaped closed surface

From the definition of dot product,

Raviteja B, Asst. Prof., Dept. of E&CE 2

Raviteja B, Asst. Prof., Dept. of E&CE 3

Since D= E we have D = - - - - - - - - - (2)

At the surface of the sphere, Ds=

The differential element of area on a spherical surface is in spherical coordinates

Leading to the closed surface integral

Raviteja B, Asst. Prof., Dept. of E&CE 4

Infinite Line Charge

Raviteja B, Asst. Prof., Dept. of E&CE 5

Dr = (Due to infinite line charge)

Fig. 4 Coaxial cable

Hence can be expressed interms of

Thus line charge density of inner conductor is C/m.

Raviteja B, Asst. Prof., Dept. of E&CE 6

Equating (2) and (3)

From eqn (1) =

E= (a < r < b ) V/m - - - - - - - - - (5)

But Qouter cylinder = 2 b L

Raviteja B, Asst. Prof., Dept. of E&CE 7

Gauss’s law applied to Differential volume element

D = Do = Dxo ax + Dyo ay +Dzo az

Dx Front = Dx0 + [ Rate of change of Dx with x ] x [Distance of surface from P ]

Raviteja B, Asst. Prof., Dept. of E&CE 8

Dx Front = Dx0 + - - - - - - - - - (6)

The point P is at the centre so distance of surface in x direction from P is ]

Consider the integral over back surface

Dx Back = Dx0 - - - - - - - - - - (10)

Combining (7) and (11)

Similarly we can write,

From eqn (3)

Raviteja B, Asst. Prof., Dept. of E&CE 9

E= = (-146.4 ax+146.4ay-195.2az) v/m

Q at P (2, −1, 3) = 10−12 m3

= = 8yz = - 648 x 10−12 = =0

Raviteja B, Asst. Prof., Dept. of E&CE 10

This is the charge enclosed in the volume ΔV.

But by Gauss’s law - - - - - - - - - (2)

This in general if is any vector then,

Divergence in different coordinate systems

Raviteja B, Asst. Prof., Dept. of E&CE 11

a) D= (2xyz-y2)ax + (x2z-2xy)ay+x2yaz c/m2 at P(2, 3, -1).

a) D= (2xyz-y2) ax + (x2z-2xy) ay+x2y az c/m2

div = = 2yz-2x At P(2,3,-1)

b) D=2rz2sin2Φ ar + rz2 sin2Φ aΦ +2r2zsin2Φ az c/m2

= 3.53 - 1.532 + 7.064 = 9.062 C/m3

Raviteja B, Asst. Prof., Dept. of E&CE 12

= 1.93 + 0.643 – 1.29 = 1.283 C/m3

Maxwell’s First Equation

According to Gauss’s law, it is known that

Expressing Gauss’s law per unit volume