Ampere’s Circuit Law

P3.14: A pair of infinite extent current sheets exists at z = -2.0 m and at z = +2.0 m.
The top sheet has a uniform current density K = 3.0 ay A/m and the bottom one has K
= -3.0 ay A/m. Find H at (a) (0,0,4m), (b) (0,0,0) and (c) (0,0,-4m).

1
We apply H  K  a N ,
2
(a) H   3a y   a z   -3a y   a z  0
1 1
2 2
(b)
H   3a y   a z   -3a y   a z
1 1
2 2
A
 3a x
m
(c) H = 0 Fig. P3.14

P3.15: An infinite extent current sheet with K = 6.0 ay A/m exists at z = 0. A

conductive loop of radius 1.0 m, in the y-z plane centered at z = 2.0 m, has zero
magnetic field intensity measured at its center. Determine the magnitude of the
current in the loop and show its direction with a sketch.

Htot = HS + HL
H S  K  a N   6a y   a z  3a x
1 1 A
2 2 m
For the loop, we use Eqn. (3.10):
I
H az
2a
where here
I I
HL   a x   a x
2a 2
(sign is chosen opposite HS).
So, I/2 = 3 and I = 6A.
Fig. P3.15

P3.16: Given the field H = 3y2 ax, find the current passing through a square in the x-y
plane that has one corner at the origin and the opposite corner at (2, 2, 0).

Referring to Figure P3.6, we evaluate the circulation of H around the square path.
b c d a

 H dL  I enc     
a b c d
b

 3 0 a x dxa x  0
2

a
c

 3y a dya y  0
2
x
b
d 0

 3  2  a x dxa x  12 dx  24

2

c 2
a


d
0

So we have Ienc = 24 A. The negative

Sign indicates current is going in the
-az direction.

P3.17: Given a 3.0 mm radius solid wire centered on the z-axis with an evenly
distributed 2.0 amps of current in the +az direction, plot the magnetic field intensity H
versus radial distance from the z-axis over the range 0 ≤  ≤ 9 mm.

Figure P3.17 shows the situation along with the Amperian Paths. We have:
 H dL  I enc , where H  H a and dL   d a ; 2 H   I enc
This will be true for each Amperian path.

 2
I I I 2
AP1: I enc   J dS, J = 2 a z , I enc  2   d   d  2
a a 0 0
a
I
So: H   a for   a
2 a 2
I
AP2: Ienc = I, H  a for   a
2

Fig. P3.17a Fig. P3.17b

% MLP0317
% generate plot for ACL problem
 a=3e-3; %radius of solid wire (m)
I=2; %current (A)
N=30; %number of data points to plot
rmax=9e-3; %max radius for plot (m)
dr=rmax/N;

for i=1:round(a/dr)
r(i)=i*dr;
H(i)=(I/(2*pi*a^2))*r(i);
end
for i=round(a/dr)+1:N
r(i)=i*dr;
H(i)=I/(2*pi*r(i));
end
plot(r,H)
xlabel('rho(m)')
ylabel('H (A/m)')
grid on

P3.18: Given a 2.0 cm radius solid wire centered on the z-axis with a current density J
= 3 A/cm2 az (for  in cm) plot the magnetic field intensity H versus radial distance
from the z-axis over the range 0 ≤  ≤ 8 cm.

We’ll let a = 2 cm.

 H dL  I enc , where H  H a and dL   d a ; 2 H   I enc
AP1 ( < a): I enc   J dS   3 a   d  d a z  2 3
and H   2a for   a
a3
AP2 ( > a): Ienc = 2a3, so H  a for   a


The MATLAB plotting routine is as follows:

% MLP0318
% generate plot for ACL problem

a=2; %radius of solid wire (cm)

N=40; %number of data points to plot
rmax=8; %max radius for plot
(cm)
dr=rmax/N;

for i=1:round(a/dr)
r(i)=i*dr;
H(i)=r(i)^2;

Fig. P3.18
end
for i=round(a/dr)+1:N
r(i)=i*dr;
H(i)=a^3/r(i);
end
plot(r,H)
xlabel('rho(cm)')
ylabel('H (A/cm)')
grid on

P3.19: An infinitesimally thin metallic cylindrical shell of radius 4.0 cm is centered

on the z-axis and carries an evenly distributed current of 10.0 mA in the +az direction.
(a) Determine the value of the surface current density on the conductive shell and (b)
plot H as a function of radial distance from the z-axis over the range 0 ≤ ≤ 12 cm.

I 10mA mA mA
(a) K s    39.8 ; so K  40a z
2 a 2  0.04m  m m
(b) for  < a, H = 0. For  > a we have:
I
H a
2
The MATLAB routine to generate the plot is as follows:
% MLP0319
% generate plot for ACL problem
a=4; %radius of solid wire (cm)
N=120; %number of data points to plot
I=10e-3; %current (A)
rmax=12; %max plot radius(cm)
dr=rmax/N;

for i=1:round(a/dr)
r(i)=i*dr;
H(i)=0;
end
for i=round(a/dr)+1:N
r(i)=i*dr;
H(i)=100*I/(2*pi*r(i));
end
plot(r,H)
xlabel('rho(cm)')
ylabel('H (A/m)')
grid on

Fig. P3.19b
 P3.20: A cylindrical pipe with a 1.0 cm
wall thickness and an inner radius of 4.0
cm is centered on the z-axis and has an
evenly distributed 3.0 amps of current in
the +az direction. Plot the magnetic field
intensity H versus radial distance from
the z-axis over the range 0 ≤  ≤ 10 cm.

For each Amperian Path:

Fig. P3.19a

 H dL  I enc , where H  H a and dL   d a ; 2 H   I enc

Now, for  < a, Ienc = 0 so H = 0.
I
For a <  < b, I enc   J dS, where J  a z and dS = d  d a z
 b  a2 
2

I
 2
 2  a2 I   2  a2 
  b2  a 2  a
I enc   d   d  I 2 2 , H  a
0
b a 2  b 2  a 2  

Fig. P3.20a Fig. P3.20b

% MLP0320
% generate plot for ACL problem

a=4; %inner radius of pipe (cm)

b=5; %outer radius of pipe(cm)
N=120; %number of data points to plot
I=3; %current (A)
rmax=10; %max radius for plot (cm)
dr=rmax/N;
aoverdr=a/dr
 boverdr=b/dr

for i=1:round(a/dr)
r(i)=i*dr;
H(i)=0;
end

for i=round(a/dr)+1:round(b/dr)
r(i)=i*dr;
num(i)=I*(r(i)^2-a^2);
den(i)=2*pi*(b^2-a^2)*r(i);
H(i)=100*num(i)/den(i);
end

for i=round(b/dr)+1:N
r(i)=i*dr;
H(i)=100*I/(2*pi*r(i));
end

plot(r,H)
xlabel('rho(cm)')
ylabel('H (A/m)')
grid on

P3.21: An infinite length line carries current I in the +az direction on the z-axis, and
this is surrounded by an infinite length cylindrical shell (centered about the z-axis) of
radius a carrying the return current I in the –az direction as a surface current. Find
expressions for the magnetic field intensity everywhere. If the current is 1.0 A and
the radius a is 2.0 cm, plot the magnitude of H versus radial distance from the z-axis
from 0.1 cm to 4 cm.

I
 H dL  I enc ; for 0<  a, H 
2
a and for   a, H  0.

Fig. P3.21a
 The MATLAB routine used to
generate Figure P3.21b is as follows:
% MLP0321
% generate plot for ACL
problem
clc
clear

a=2; %inner radius

of cylinder(cm)
Fig. P3.21b N=80; %number of data
points to plot
I=1; %current (A)
rmax=4; %max radius for plot (cm)
dr=rmax/N;

for i=1:40
r(i)=.1+(i-1)*dr;
H(i)=100*I/(2*pi*r(i));
end

for i=40:N
r(i)=i*dr;
H(i)=0;
end

plot(r,H)
xlabel('rho(cm)')
ylabel('H (A/m)')
grid on

P3.22: Consider a pair of collinear cylindrical shells centered on the z-axis. The inner
shell has radius a and carries a sheet current totaling I amps in the +az direction while
the outer shell of radius b carries the return current I in the –az direction. Find
expressions for the magnetic field intensity everywhere. If a = 2cm, b = 4cm and I =
4A, plot the magnitude of H versus radial distance from the z-axis from 0 to 8 cm.

 H dL  I enc ; for 0<  a, H = 0;

I
for a<  b, H  a ;
2
and for   b, H  0.
The MATLAB routine used to generate Figure P3.22b is as follows:
% MLP0322
% generate plot for ACL problem
a=2; %inner radius of coax (cm)
b=4; %outer radius of coax(cm)
N=160; %number of data points to plot

Fig. P3.22a Fig. P3.22b

I=4; %current (A)

rmax=8; %max radius for plot (cm)
dr=rmax/N;
aoverdr=a/dr
boverdr=b/dr

for i=1:round(a/dr)
r(i)=i*dr;
H(i)=0;
end

for i=round(a/dr)+1:round(b/dr)
r(i)=i*dr;
H(i)=100*I/(2*pi*r(i));
end

for i=round(b/dr)+1:N
r(i)=i*dr;
H(i)=0;
end

plot(r,H)
xlabel('rho(cm)')
ylabel('H (A/m)')
grid on

P3.23: Consider the toroid in Figure 3.55 that is tightly wrapped with N turns of
conductive wire. For an Amperian path with radius less than a, no current is enclosed
and therefore the field is zero. Likewise, for radius greater than c, the net current
enclosed is zero and again the field is zero. Use Ampere’s Circuital Law to find an
expression for the magnetic field at radius b, the center of the toroid.

 H dL  I enc ;
Within the toroid, H = Ha, so
 H dL   H a    b
da  2 bH .
Then, Ienc by the Amperian path is: Ienc = NI.
NI
H  a .
2 b

4. Curl and the Point Form of Ampere’s Circuital Law

P3.24: Find  A for the following fields:
a. A = 3xy2/z ax
b. A = sin2 a – 2 z cos a
c. A = r2sin ar + r/cos a

Ax A 3xy 2
(a)   3xy 2 z  a x 
6 xy
a y  x az  ay  az
z y z 2
z
A A 1 A 
(b)    sin 2 a    2 z cos a    a  a     A     a z
z z     
3 z cos 
2
1
  2 cos a   0a   2sin  cos a z
az 
 
  2 cos a  3 z cos   2sin  cos   a z
(c)
 r 
   r 2 sin  a r  a 
 cos  
1  A  1  1 Ar  1    rA  Ar 
      a
r sin     r  sin   
a a
r  r  
r

r   cos  
1
1   r2  1  2

r sin  
ar    a 
r r  cos   r 
 r sin   a
 sin   2 
 ar    r cos   a
sin  cos 
2
 cos  

P3.25: Find J at (3m, 60, 4m) for H = (z/sin) a – (2/cos) az A/m.

 1 Az  A A  1 A
 H  a      z  a  a
   z     z
  sin   1 2  z cos 
 a     a  a
cos 
2
 sin  cos    sin 2  z
Now find J by evaluating H at the given point:
A
J  10a   13a  0.89a z 2
m

P3.26: Suppose H = y2ax + x2ay A/m.

a. Calculate  H dL around the path A  B  C  D  A , where A(2m,0,0),
B(2m,4m,0), C(0,4m,0) and D(0,0,0).
b. Divide this  H dL by the area S (2m*4m = 8m2).
c. Evaluate H at the center point.
d. Comment on your results for (b) & (c).

B C D A
(a) Referring to the figure, we evaluate  H dL  I enc     
A B C D
B 4

   x2 dy  4  4   16C
A 0 x2
C 0

   y2 dx  32C
B 2 y 4
D 0

   x2 dy  0
C 4 x 0
A 2

   y2 dy  0
D 0 y 0

So we have  H dL  16C Fig. P3.26

(b) dividing by S = 8m2, we have -2 C/m2

(c) Evaluating the curl of H:
 A A 
 H   y  x  a z   2 x  2 y  a z , and at the center point (x = 1 and y = 2) we
 x y 
have
C
  H center  2 2 a z
m
(d) In this particular case,   H   H dL S , even though S is of appreciable size.
P3.27: For the coaxial cable example 3.8, we found:
 I
for   a, H a ,
2 a 2
I
for a<  b, H a ,
2
I  c2   2 
for b<  c, H=   a ,

2 c 2  b 2    
and for c< , H =0.
a. Evaluate the curl in all 4 regions.
b. Calculate the current density in the conductive regions by dividing the current
by the area. Are these results the same as what you found in (a)?

1   I2  I
(a)  H   a  2 a z for   a
2  z
   2 a  a
1   I 
H   a z  0 for a    b
   2 
1   I c    
2 2
I
 H    az  a z for b<  c
   2  c  b  
2 2
  c 2
 b 2

 
  H = 0 for   c

 
(b) for   a, S   a 2 , J  I S a z  2 a z
I
a
I
for b<  c, S    c 2  b 2  , J  az
  c2  b2 
Comment: H = J is confirmed.

P3.28: Suppose you have the field H = r cos  a A/m. Now consider the cone
specified by  = /4, with a height a as shown in Figure 3.56. The circular top of the
cone has a radius a.
a. Evaluate the right side of Stoke’s theorem through the dS = dSa surface.
b. Evaluate the left side of Stoke’s theorem by integrating around the loop.

1   1   rH 
(a)  H  
r sin   
 sin  H  ar 
 r r
a

ar derivative:
 1    cos 2   sin 2  
 r sin  cos   r  cos2   sin 2   ; sin  H  ar  ar
 r sin    sin 
 Fig. P3.28
a derivative:
1   r cos  
2

a  2 cos  a
r r

So,   H 
 cos 2   sin 2  
a r  2cos  a
sin 
Now we must integrate this over the a surface:
  cos 2   sin 2   
   H dS  
 a  2 cos  a  r sin  drd  a 
sin 
r
 
2a 2
  2r sin  cos  drd 2  sin  cos       rdr  d  2 a
2

4 0 0
2

 H dL   r r a cos  a ad a  a 2 cos      d  2 a

2
(b) 2
4 0

Clearly in this case the circulation of H is the easiest approach.

5. Magnetic Flux Density

P3.29: An infinite length line of 3.0 A current in the +ay direction lies on the y-axis.
Find the magnetic flux density at P(7.0m,0,0) in (a) Teslas, (b) Wb/m2, and (c) Gauss.

I 3A mA
H a   a z   68 a z
2 2  7m  m
 H  A   Wb  9 Wb
B  o H   4 x107  68  a z    86 x10 a z  86 x109 Ta z
 m  m   HA  2
m

B   86 x109 Ta z 
10, 000G
 860 x106 Ga z
T

P3.30: Suppose an infinite extent sheet of current with K = 12ax A/m lies on the x-y
plane at z = 0. Find B for any point above the sheet. Find the magnetic flux passing
through a 2m2 area in the x-z plane for z > 0.
 1
H  K  aN ;
2

B
o
K  aN 
 4 x10
H m 
7
A  6 Wb
12a x  a z   7.54 x10 ay
2 2  m  m2
This is valid at any point above the sheet.

 Wb 

Now,    B dS  B S =  -7.54x10-6 2 a y  2m2  a y   15Wb
 m 


P3.31: An infinite length coaxial cable exists along the z-axis, with an inner shell of
radius a carrying current I in the +az direction and outer shell of radius b carrying the
return current. Find the magnetic flux passing through an area of length h along the z-
axis bounded by radius between a and b.

I o I
H a , B  a ,
2 2
o I a I
For a <  < b,    B dS   d  dza  o ln    a h
b

2  2
o Ih  b 
 ln   Wb
2 a

6. Magnetic Forces
P3.32: A 1.0 nC charge with velocity 100. m/sec in the y direction enters a region
where the electric field intensity is 100. V/m az and the magnetic flux density is 5.0
Wb/m2 ax. Determine the force vector acting on the charge.
m Wb Wb
F  q  E  u  B  ; u  B  100 a y  5 2  500 az
s m ax sm
 V Wb Vs  mN
F  109 C 100 a z  500 az   400nNa z
 m sm Wb  VC

P3.33: A 10. nC charge with velocity 100. m/sec in the z direction enters a region
where the electric field intensity is 800. V/m ax and the magnetic flux density 12.0
Wb/m2 ay. Determine the force vector acting on the charge.

 V m Wb 
F  q  E  u  B   10 x109 C  800 a x  100 a z  12 2   4  Na x
 m s m ay
 
P3.34: A 10. nC charged particle has a velocity v = 3.0ax + 4.0ay + 5.0az m/sec as it
enters a magnetic field B = 1000. T ay (recall that a tesla T = Wb/m2). Calculate the
force vector on the charge.

 Wb 
F  q  u  B   10 x109 C  3a x  4a y  5a z  1000 2 a y 
m
 s m 
The cross-product:
a x a y az 
3 4 5   5000a x  3000a z

 0 1000 0 
Evaluating we find: F = -50ax + 30az N

P3.35: What electric field is required so that the velocity of the charged particle in the
previous problem remains constant?

dv
F  ma  m  0 (constant velocity)
dt
F  q  E  u  B   0;

E  u  B = -  3a x  4a y  5a z 
m Wb V V
1000 2 a y  3000 a z  5000 a x
s m m m
kV
E  5a x  3a z
m

P3.36: An electron (with rest mass Me= 9.11x10-31kg and charge q = -1.6 x 10-19 C)
has a velocity of 1.0 km/sec as it enters a 1.0 nT magnetic field. The field is oriented
normal to the velocity of the electron. Determine the magnitude of the acceleration
on the electron caused by its encounter with the magnetic field.

F  ma  q  u  B  ;
q
a u  B  ;
m

a
quB  1.6 x10

19

C  1000 m 
s
109 Wb 2
m  175 x103 m
m  9.11x10 kg 
31 s2

P3.37: Suppose you have a surface current K = 20. ax A/m along the z = 0 plane.
About a meter or so above this plane, a 5.0 nC charged particle is moving along with
velocity v = -10.ax m/sec. Determine the force vector on this particle.
 1 1 A
H K  a N  20a x  a z  10a y
2 2 m
B  o H  10  4 x10  a y Wb 2  40 x107 a y Wb 2
7
m m
 Wb 
F  qu  b   5 x109 C   10a x  40 x107 a y 2   0.63 pNa z
m
 s m 

P3.38: A meter or so above the surface current of the previous problem there is an
infinite length line conducting 1.0 A of current in the –ax direction. Determine the
force per unit length acting on this line of current.

0
Wb
F12   I 2 dL2  B1  I 2  dxa x  40 x107 a y
L
m2
F12  I 2   L   40 x107 a z  ;
F12 N
 40 x107 a z  12.6 az
L m

P3.39: Recall that the gravitational force on a mass m is F  mg, where, at the earth’s
surface, g = 9.8 m/s2 (-az). A line of 2.0 A current with 100. g mass per meter length
is horizontal with the earth’s surface and is directed from west to east. What
magnitude and direction of uniform magnetic flux density would be required to
levitate this line?

 m  Ns 2 kg 
Fg  mg; for 1 m Fg  100 g   9.8 2   a z     0.98 Na z
 s   kg m 1000 g 
1m
F  IdL  B   F
0
g  0.98 Na z

By inspection, B = Bo(-ax)
L

 Idya y  Bo  a x   ILBoa z
o

 Wb 
  2 A 1m   Bo 2   2 Bo  N 
 m 
The unit conversion to arrive at
Newtons is as follows:
Am Wb Vs W J Nm
N
m 2 Wb VA Ws J

So we have Bo = 0.490 Wb/m2, and

B = 0.490 Wb/m2 (-ax) (directed north) Fig. P3.39
P3.40: Suppose you have a pair of parallel lines each with a mass per unit length of
0.10 kg/m. One line sits on the ground and conducts 200. A in the +ax direction, and
the other one, 1.0 cm above the first (and parallel), has sufficient current to levitate.
Determine the current and its direction for line 2.

F12 o I1I 2 4 x107 H m 200 A

Here we will use  ay  I 2a z  4 x103 I 2a z
L 2 y 2 0.01m
2
  0.10 kg m   9.8 m s 
F mg 2 Ns N
  0.98
L L kgm m
So solving for I2:
0.98
I2   245 A in the -a x direction.
4 x103

P3.41: In Figure 3.57, a 2.0 A line of current is shown on the z-axis with the current in
the +az direction. A current loop exists on the x-y plane (z = 0) that has 4 wires
(labeled 1 through 4) and carries 1.0 mA as shown. Find the force on each arm and
the total force acting on the loop from the field of the 2.0 A line.

o I1
F12   I 2 dL2  B1 ; B1 = a
2
A  B : dL2  5da , a  a  0
C  D : a  a  0
3
o I1  I I 3
B  C : F12  I 2  d  a   a  o 1 2 ln   a z
5
2 2 5

Fig. P3.41

4 x107 H m
 2 A 103 A ln   a z  204 pNa z
3
F12 
2 5
So for B to C: F12 = -0.20 nN az
Likewise, from D to A: F12 = +0.20 nN az
P3.42: MATLAB: Modify MATLAB 3.4 to find the differential force acting from
each individual differential segment on the loop. Plot this force against the phi
location of the segment.

%MLP0342
%modify ML0304 to find dF acting from the field
% of each segment of current; plot vs phi

clear
clc
I=1; %current in A
a=1; %loop radius, in m
mu=pi*(4e-7); %free space permeability
az=[0 0 1]; %unit vector in z direction
DL1=a*2*(pi/180)*[0 1 0]; %Assume 2 degree increments
%DL1 is the test element vector
%F is the angle phi in radians
%xi,yi is location of ith element on the loop
%Ai & ai = vector and unit vector from origin
% to xi,yi
%DLi is the ith element vector
%Ri1 & ri1 = vector and unit vector from ith
% point to test point

for i=1:179
phi(i)=i*2;
F=2*i*pi/180;
xi=a*cos(F);
yi=a*sin(F);
Ai=[xi yi 0];
ai=unitvector(Ai);
DLi=(pi*a/90)*cross(az,ai);
Ri1=[a-xi -yi 0];
ri1=unitvector(Ri1);
num=mu*I*cross(DLi,ri1);
den=4*pi*(magvector(Ri1)^2);
B=num/den;
dFvect=I*cross(DL1,B);
dF(i)=dFvect(1);
end
plot(phi,dF)
xlabel('angle in degrees')
ylabel('the differential force, N')

Fig. P3.42

P3.43: MATLAB: Consider a circular conducting loop of radius 4.0 cm in the y-z
plane centered at (0,6cm,0). The loop conducts 1.0 mA current clockwise as viewed
from the +x-axis. An infinite length line on the z-axis conducts 10. A current in the
+az direction. Find the net force on the loop.

The following MATLAB routine shows the force as a function of radial position
around the loop. Notice that while there is a net force in the -y direction, the forces in
the z-direction cancel.

% MLP0343
% find total force and torque on a loop of
% current next to a line of current

% variables
% I1,I2 current in the line and loop (A)
% yo center of loop on y axis (m)
% uo free space permeability (H/m)
% N number of segements on loop
% a loop radius (m)
% dalpha differential loop element
% dL length of differential section
% DL diff section vector
% B1 I1's mag flux vector (Wb/m^2)
% Rv vector from center of loop
% to the diff segment
% ar unit vector for Rv
% y,z the location of the diff segment

clc
clear

% initialize variables
I1=10;
I2=1e-3;
yo=.06;
uo=pi*4e-7;
N=180;
a=0.04;
dalpha=360/N;
dL=a*dalpha*pi/180;
ax=[1 0 0];

% perform calculations
for i=1:N
dalpha=360/N; Fig. P3.43
alpha=(i-1)*dalpha;
phi(i)=alpha;
z=a*sin(alpha*pi/180);
y=yo+a*cos(alpha*pi/180);
B1=-(uo*I1/(2*pi*y))*[1 0 0];
Rv=[0 y-.06 z];
ar=unitvector(Rv);
aL=cross(ar,ax);
DL=dL*aL;
dF=cross(I2*DL,B1);
dFx(i)=dF(1);
dFy(i)=dF(2);
dFz(i)=dF(3);
end
plot(phi,dFy,phi,dFz,'--k')
legend('dFy','dFz')
Fnet=sum(dFy)

Running the program we get:

Fnet = -4.2932e-009
>>
So Fnet = -4.3 nN ay

P3.44: MATLAB: A square loop of 1.0 A current of side 4.0 cm is centered on the x-y
plane. Assume 1 mm diameter wire, and estimate the force vector on one arm
resulting from the field of the other 3 arms.
% MLP0344 V2
%
% Square loop of current is centered on x-y plane.
Viewed
% from the +z axis, let current go clockwise. We want
to
% find the force on the arem at x = +2 cm resulting
from the
% current in arms at y = -2 cm, x = -2 cm and y = +2 cm.

% Wentworth, 12/3/03

% Variables
% a side length (m)
% b wire radius (m)
% I current in loop (A)
% uo free space permeability (H/m)
% N number of segments for each arm
% xi,yi location of test arm segment (at x = +2 cm)
% xj,yj location of source arm segment (at y = -2 cm)
% xk,yk location of source arm segment (at x = -2 cm)
% xL,yL location of source arm segment (at y = +2 cm)
% dLi differential test segment vector
% dLj, dLk, dLL diff vectors on sources
% Rji vector from source point j to test point i
% aji unit vector of Rji

clc;clear;

a=0.04;
b=.0005;
I=1;
uo=pi*4e-7;
N=80;

for i=1:N
xi=(a/2)+b;
yi=-(a/2)+(i-0.5)*a/N;
ypos(i)=yi;
dLi=(a/N)*[0 -1 0];
for j=1:N
xj=-(a/2)+(j-0.5)*a/N;
yj=-a/2-b;
dLj=(a/N)*[-1 0 0];
Rji=[xi-xj yi-yj 0];
aji=unitvector(Rji);
 num=I*CROSS(dLj,aji);
den=4*pi*(magvector(Rji))^2;
H=num/den;
dHj(j)=H(3);
end
for k=1:N
yk=(-a/2)+(k-0.5)*a/N;
xk=-(a/2)-b;
dLk=(a/N)*[0 1 0];
Rki=[xi-xk yi-yk 0];
aki=unitvector(Rki);
num=I*CROSS(dLk,aki);
den=4*pi*(magvector(Rki))^2;
H=num/den;
dHk(k)=H(3);
end
for L=1:N
xL=(-a/2)+(L-0.5)*a/N;
yL=(a/2)+b;
dLL=(a/N)*[1 0 0];
RLi=[xi-xL yi-yL 0];
aLi=unitvector(RLi);
num=I*CROSS(dLL,aLi);
den=4*pi*(magvector(RLi))^2;
H=num/den;
dHL(L)=H(3);
end
H=sum(dHj)+sum(dHk)+sum(dHL);
B=uo*H*[0 0 1];
F=I*CROSS(dLi,B);
dF(i)=F(1);
end
Ftot=sum(dF)
plot(ypos,dF)

Running the program:

Ftot =

7.4448e-007
Fig. P3.44
So Ftot = 740 nN
P3.45: A current sheet K = 100ax A/m exists at z = 2.0 cm. A 2.0 cm diameter loop
centered in the x-y plane at z = 0 conducts 1.0 mA current in the +a direction. Find
the torque on this loop.

  m  B;
 
m  ISa N  103 A    0.01m  a z  314 x109 Am 2a z
2

1 1 A A
B   o H; H  K  a N  100a x    a z   50a y ;
2 2 m m
Wb
B  50oa y ;
m2
  m  B = -20 pNma x

P3.46: 10 turns of insulated wire in a 4.0 cm diameter coil are centered in the x-y
plane. Each strand of the coil conducts 2.0 A of current in the a direction. (a) What
is the magnetic dipole moment of this coil? Now suppose this coil is in a uniform
magnetic field B = 6.0ax + 3.0ay + 6.0az Wb/m2, (b) what is the torque on the coil?


(a) m  NISa x  10  2 A   0.02m 
2
  25.1mAm a
2
z

(b)   m  B  25.1mAa z   6a x  3a y  6a z 
Wb
 0.151a y  0.075a x Nm
m2
  75a x  151a y mNm

P3.47: A square conducting loop of side 2.0 cm is free to rotate about one side that is
fixed on the z-axis. There is 1.0 A current in the loop, flowing in the –az direction on
the fixed side. A uniform B-field exists such that when the loop is positioned at  =
90, no torque acts on the loop, and when the loop is positioned at  = 180 a
maximum torque of 8.0 N-m az occurs. Determine the magnetic flux density.

At  = 90°, m  ISa N  0.0004 Am 2a x . Also, since   m  B  0, B is in direction of

m, and therefore B = ±Boax.
At  = 180°, m  ISa N  0.0004 Am2a y , and   m  B  8 x106 a z Nm.
Therefore, B = -Boax and mBo = 8x10-6, so
8 x106 mWb mWb
Bo   20 2 , and B  20 2 a x .
0.0004 m m

7. Magnetic Materials
P3.48: A solid nickel wire of diameter 2.0 mm evenly conducts 1.0 amp of current.
Determine the magnitude of the magnetic flux density B as a function of radial
distance from the center of the wire. Plot to a radius of 2 mm.
 I 1A kA
J a  a z  31.8 2 a z
a  1x103 m 
2 z 2
m

 H dL  I enc   J dS
I I
2 
2 H    d  d  2  2
a a
I   I
for   a H  a ; B  r o 2 a
2 
2 a 2 a
I I
for   a H  a ; B  o a
2 2

% MLP0348
% generate plot for ACL problem

a=2e-3; %radius of solid wire (m)

I=1; %current (A)
N=30; %number of data points to plot
rmax=4e-3; %max radius for plot (m)
dr=rmax/N;
uo=pi*4e-7;
ur=600;

for i=1:round(a/dr)
r(i)=i*dr;
B(i)=(ur*uo*I/(2*pi*a^2))*r(i);
end
for i=round(a/dr)+1:N
r(i)=i*dr;
B(i)=uo*I/(2*pi*r(i));
End
rmm=r*1000;
plot(rmm,B)
xlabel('rho(cm)')
ylabel('B (Wb/m^2)')
grid on

Fig. P3.48

8. Boundary Conditions
P3.49: A planar interface separates two magnetic media. The magnetic field in media
1 (with r1) makes an angle 1 with a normal to the interface. (a) Find an equation for
2, the angle the field in media 2 (that has r2) makes with a normal to the interface,
in terms of 1 and the relative permeabilities in the two media. (b) Suppose media 1
is nickel and media 2 is air, and that the magnetic field in the nickel makes an 89
angle with a normal to the surface. Find 2.

H1  H1N a N  H1T aT ; H 2T aT  H1T aT ; B1N a N  r1o H1N a N

 r1
B2 N a N  B1N a N  r1o H1N a N  r 2 o H 2 N a N ;  H 2 N  H
 r 2 1N
H 2T H1T
tan  2   ;
H2N  r1
H1N
r 2
r 2 H1T r 2
tan  2   tan 1
r1 H1N r1
 r 2 
 2  tan 1  tan 1 
  r1 
 1 
 2  tan 1  tan 89   5.5
 600  Fig. P3.49

P3.50: MATLAB: Suppose the z = 0 plane separates two magnetic media, and that no
surface current exists at the interface. Construct a program that prompts the user for
r1 (for z < 0), r2 (for z > 0), and one of the fields, either H1 or H2. The program is to
calculate the unknown H. Verify the program using Example 3.11.

% M-File: MLP0350
%
% Given H1 at boundary between a pair of
% materials with no surface current at boundary,
% calculate H2.
%
clc
clear
% enter variables
disp('enter vectors quantities in brackets,')
disp('for example: [1 2 3]')
ur1=input('relative permeability in material 1: ');
ur2=input('relative permeability in material 2: ');
a12=input('unit vector from mtrl 1 to mtrl 2: ');
F=input('material where field is known (1 or 2): ');
Ha=input('known magnetic field intensity vector: ');

if F==1
ura=ur1;
urb=ur2;
a=a12;
else
ura=ur2;
urb=ur1;
a=-a12;
end

% perform calculations
Hna=dot(Ha,a)*a;
Hta=Ha-Hna;
Htb=Hta;
Bna=ura*Hna; %ignores uo since it will factor out
Bnb=Bna;
Hnb=Bnb/urb;
display('The magnetic field in the other medium is: ')
Hb=Htb+Hnb

Now run the program (for Example 3.11):

enter vectors quantities in brackets,
for example: [1 2 3]
relative permeability in material 1: 6000
relative permeability in material 2: 3000
unit vector from mtrl 1 to mtrl 2: [0 0 1]
material where field is known (1 or 2): 1
known magnetic field intensity vector: [6 2 3]

ans =
The magnetic field in the other medium is:
Hb = 6 2 6

For a second test, run the program for problem P3.52(a).

enter vectors quantities in brackets,

for example: [1 2 3]
relative permeability in material 1: 4
relative permeability in material 2: 1
unit vector from mtrl 1 to mtrl 2: [0 0 -1]
material where field is known (1 or 2): 1
known magnetic field intensity vector: [3 0 4]

ans =
The magnetic field in the other medium is:
Hb = 3 0 16
P3.51: The plane y = 0 separates two magnetic media. Media 1 (y < 0) has r1 = 3.0
and media 2 (y > 0) has r2 = 9.0. A sheet current K = (1/o) ax A/m exists at the
interface, and B1 = 4.0ay + 6.0az Wb/m2. (a) Find B2. (b) What angles do B1 and B2
make with a normal to the surface?

(a) B N 1  4a y
(b) B N 2  B N 1  4a y
(c) BT 1  6a z
BT 1 2
(d) HT 1   az
 r1o o
3
(e) HT 2  a z (see below)
o
(f) BT 2  r o HT 2  27a z
Wb Fig. P3.51
(g) B 2  4a y  27a z
m2
Now for step (e):
1
a 21   H1  H 2   K ;  a y   H T 1a z  H T 2a z   ax
o
1 1 1 2 3
  HT 1  HT 2  a x  a x ; HT 2  HT 1  ; HT 2   
o o o o o

 B   BT 2 
Angles: 1  tan 1  T 1   56 ;  2  tan 
1
  82
 BN 1   BN 2 

P3.52: Above the x-y plane (z > 0), there exists a magnetic material with r1 = 4.0 and
a field H1 = 3.0ax + 4.0az A/m. Below the plane (z < 0) is free space. (a) Find H2,
assuming the boundary is free of surface current. What angle does H2 make with a
normal to the surface? (b) Find H2, assuming the boundary has a surface current K =
5.0 ax A/m.

(a)
BN 2
(1) H N 1  4a z , (2) B N 1  16oa z , (3) B N 2  B N 1  16oa z ,(4) H N 2   16a z
o
A
(5) HT 1  3a x , (6) HT 2  HT 1  3a x , (7) H 2  3a x  16a z
m
 HT 2 
 2  tan 1    10.6
 HN2 
(b)
Now step (6) becomes a21   H1  H2   K , shere a21 = az.
Let’s let H 2  Aa x  Ba y  16a z , then a 21   3a x  4a z  Aa x  Ba y  16a z   5a x
Solve for A and B:
 ax a y az 
 0 0 1   Ba x   A  3 a y  5a x ; so A = 3 and B = 5

3  A  B 12

A
Finally, H 2  3a x  5a y  16a z .
m

P3.53: The x-z plane separates magnetic material with r1 = 2.0 (for y < 0) from
magnetic material with r2 = 4.0 (for y > 0). In medium 1, there is a field H1 = 2.0ax
+ 4.0ay + 6.0az A/m. Find H2 assuming the boundary has a surface current K = 2.0ax
– 2.0az A/m.

8 o
(1)H1N  4a y ,(2)B1N  8oa y ,(3)B 2 N  B1N  8oa y ,(4)H 2 N  a y  2a y
4 o
(5)H1T  2a x  6az ,(6)a21   H1  H2   K, let HT 2  H xa x  H zaz ,
so  a y    2  H x  a x   6  H z  a z   2a x  2a z ,
ax ay az
0 1 0   H z  6  a x   2  H x  a z  2a x  2a z , so H z  8, H x  4
2  Hx 0 6  Hz
A
(7)H 2T  4a x  8a z ,  H 2  4a x  2a y  8a z
m

P3.54: An infinite length line of 2 A current in the +az direction exists on the z-axis.
This is surrounded by air for  ≤ 50 cm, at which point the magnetic media has r2 =
9.0 for  > 50 cm. If the field in media 2 at  = 1.0 m is H = 5.0a A/m, find the sheet
current density vector at  = 50. cm, if any.

Method 1:
I1
From just the line of current we would have H1  a  1a .
2 1
Now, since HTOT  5a  H1  H 2 , then H 2  4a is the contribution from the sheet
current.
I2 I2 8 A
H2  a  4a , so I 2  8 A, then K  az   8 az
2 1 2 a 2  0.5 m
A
K  8 az
m
Method 2:
From I1 at boundary we have
2 a
H1   2a , but 5a at  = 1.0m corresponds to 10a at  =0.5m since H
2  0.5 
varies as 1/. So
a    2a  10a   K ,  a    8a   8a z  K
A
K  8 az
m
Faraday’s law and EMF
P4.9: The magnetic flux density increases at the rate of 10 (Wb/m 2)/sec in the z
direction. A 10 cm x 10 cm square conducting loop, centered at the origin in the x-y
plane, has 10 ohms of distributed resistance. Determine the direction (with a sketch)
and magnitude of the induced current in the conducting loop.

dB Wb
 10 2 a z
dt ms
dB Wb
Vemf    dS    10 2 a z dxdya z
dt ms
Wb Vs
Vemf  0.1  0.1V
s Wb
0.1V
I= I   10mA
10
I=10mA clockwise (when viewed
From +z) Fig. P4.9

P4.10: A bar magnet is dropped through a conductive ring. Indicate in a sketch the
direction of the induced current when the falling magnet is just above the plane of the
ring and when it is just below the plane of the ring, as shown in Figure 4.22.

Refer to Figure P4.10.

When the north pole first goes through the loop, flux is increasing and the current
induced to oppose this change in flux is as shown.

When the south pole is exiting the loop, flux is decreasing and the current induced
acts to oppose this change in flux.

Fig. P4.10

P4.11: Considering Figure 4.7, suppose the area of a single loop of the pair is 100
cm2, and the magnetic flux density is constant over the area of the loops but changes
with time as B  Bo e   t a z , where Bo = 4.0 mWb/m2 and  = 0.30 Np/sec. Determine
VR at 1, 10, and 100 seconds.

B dB
Vemf   N  dS;    Bo e   t a z ; Vemf  2  Bo e   t S
t dt
2
 1m  0.30t
VR  2  Bo Se   t  2  0.30   4 x103 100cm 2    e  24 x106 e 0.30t
 100cm 
at t = 1 sec, VR = 17.8 V
at t = 10 sec, VR = 1.20 V
at t = 100 sec, VR = 2.25x10-18 V

P4.12: Sometimes a transformer is used as an impedance converter, where impedance

is given by v/i. Find an expression for the impedance Z1 seen by the primary side of
the transformer in Figure 4.11 that has a load impedance Z2 terminating the
secondary.

N1 N
We have i2  i1 and v2  2 v1
N2 N1
 N2 
 v1  2
v1 v2  N 1   N 2 
Z1  , Z 2     Z1
i1 i2  N1   N1 
 N i1 
 2 
2
N 
Z1   1  Z 2
 N2 

P4.13: A 1.0 mm diameter copper wire is shaped into a square loop of side 4.0 cm. It
is placed in a plane normal to a magnetic field increasing with time as B = 1.0 t
Wb/m2 az, where t is in seconds. (a) Find the magnitude of the induced current and
indicate its direction in a sketch. (b) Calculate the magnetic flux density at the center
of the loop resulting from the induced current, and compare this with the original
magnetic flux density that generated the induced current at t = 1.0 sec.

We find the distributed resistance of the loop and work the problem assuming this
resistance is lumped in one spot as shown in the figure.
(a) The induced current is Vemf divided by the distributed resistance of the wire loop.
1 l 1m 4  0.04m) 
R   3.5m
 A 5.8 x10   0.0005m 2
7

0.04 0.04
dB Wb dB Wb
 1 2 a z ; Vemf    dS  1 2  dx  dy  1.6mV
dt ms dt ms 0 0
 1.6mV
I ind   0.46 A (note that this answer has no time dependence)
3.5m

Fig. P4.13

(b) The field at the center of the loop from a single arm of the loop is found from Eqn.
(3.7):
a
I  z  I 1  0.46  1 A
H  2 2   a z    -a z    a z   2.59a z ;
4   z   2   a 2 2a 2 2  0.02  m

Wb
So B  4 o H  13 az .
m2

P4.16: Referring to Figure 4.23, suppose a conductive bar of length h = 2.0 cm

moves with velocity u = -1.0 m/s a towards an infinite length line of current I = 4.0
A. Find an expression for the voltage from one end of the bar to the other when 
reaches 10 cm and indicate which end is positive.
In Figure P4.16, an imaginary circuit has been chosen. For the chosen circulation
direction, we have the sign for Vemf as shown. Then,
 o I   m  o I
h
Vemf    u  B  dL    -1a   a  dza z   1   dz ,
 2   s  2 0
Vemf  160nV . Therefore, the bottom of the bar is positive.

Fig. P4.16
P4.17: Suppose we have a conductive bar moving along a pair of conductive rails as
in Figure 4.12, only now the magnetic flux density is B = 4.0ax + 3.0az Wb/m2. If R =
10. , w = 20. cm, and uy = 3.0 m/s, calculate the current induced and indicate its
direction.

 Wb  m 
Vemf    B u y dxa z    3 2  3   0.2m   1.8V
 m  s 
1.8V
I  0.18 A (clockwise when viewed from the +z axis)
10

P4.18: The radius r of a perfectly conducting metal loop in free space, situated in the
x-y plane, increases at the rate of (r)-1 m/sec. A break in the loop has a small 2.0
ohm resistor across it. Meanwhile, there exists a magnetic field B = 1.0 az T.
Determine the current induced in the loop, and show in a sketch the direction of flow.

Here we’ve assumed dS = -dSaz to get iind and Vemf as shown. Our approach will be to
find , then Vemf = -d/dt.

   B dS   1a z    d  d a z 
r
 2   d    r 2
0

d dr  1 
 2 r  2 r    2V
dt dt r 
Vemf  2V
2V
I  1A, clockwise as shown
2

Fig. P4.18

P4.19: Rederive Vemf for the rectangular loop of Figure 4.16 if the magnetic field is
now B = Boaz.

We see in Figure P4.19a that  u  B  dL  0 for the 1  2 and 3  4 line sections.

For the 2  3 section we have:
d l  d a
3

    u  B  dL; B  Boa z , dL  d  a  , u    a

2
dt dt

Fig. P4.19a Fig. P4.19b

  Bo a 2
0

  u  B  dL   Bo a  d   2
 Bo a 2
a

  u  B  dL   Bo   d  
0
2
, so for the 2  3 section, the contributions to

Vemf cancel. This will also be the case for the 4  1 section, and therefore Vemf = 0;
no current is induced.

P4.20: In Figure 4.16, replace the rectangular loop with a circular one of radius a and
rederive Vemf.

   B dS   Boa y dS = Boa y  a 2   sin  a x  cos  a y   Bo a 2 cos 

d flux
  Bo a 2 sin  ;
dt
d flux
Vemf     Bo a 2 sin 
dt

P4.21: A conductive rod, of length 6.0 cm, has one end fixed on a grounded origin
and is free to rotate in the x-y plane. It rotates at 60 revolutions per second in a
magnetic field B = 100. mT az. Find the voltage at the end of the bar.

 rev  2 rad  rad

   60    120
 s  rev  s
Vemf    u  B  dL
d
u a   a
dt
l
 Bo l 2
Vemf    a  Boa z  d  a  
0
2
1 rad  Wb  2 Vs
Vemf  120  0.1 2   0.06m 
2 s  m  Wb Fig. P4.21
 68mV

We can confirm the sign by observing that a positive charge placed in the middle of
the bar would move to the ungrounded end by the Lorentz force equation.

P4.22: Consider the rotating conductor shown in Figure 4.24. The center of the 2a
diameter bar is fixed at the origin, and can rotate in x-y plane with B = Boaz. The
outer ends of the bar make conductive contact with a ring to make one end of the
electrical contact to R; the other contact is made to the center of the bar. Given Bo =
100. mWb/m2, a = 6.0 cm, and R = 50. , determine I if the bar rotates at 1.0
revolution per second.
 d
Vemf    u  B  dL; u
dt
a   a , dL  d  a 
Figure P4.22 indicates one of the paths for the circulation integral.
a a
 Bo a 2
Vemf    a  Boa z  d  a     Bo   d  
0 0
2
Vemf  Bo a 2
1  rev  rad  3 Wb  2  1  Vs A
2 
I   1  2 100 x10 0.06m   
R 2R 2  s  rev  m   50  Wb V
I  22.6  A

Fig. P4.22

P4.23: A Faraday Disk Generator is similar to the rotating conductor of P4.22, only
now the rotating element is a disk instead of a bar. Derive an expression of the Vemf
produced by a Faraday Disk Generator, and using the parameters given in problem
4.22, find I.

Worked exactly as P4.22.

P4.24: Consider a sliding rail problem where the conductive rails expand as they
progress in the y direction as shown in Figure 4.25. If w = 10. cm and the distance
between the rails increases at the rate of 1.0 cm in the x direction per 1.0 cm in the y
direction, and uy = 2.0 m/sec, find the Vemf across a 100.  resistor at the instant when
y = 10. cm if the field is Bo = 100. mT.

First we modify the figure so that the top rail is horizontal and all the spreading occurs
via the bottom rail. As before, our approach will be to find  and then d /dt. We
have:
   B dS   Boa z dxdya z
Now, notice that x and y are not independent and are in fact related: x=y+w
So we have
 y yw y
1 
  Bo   dxdy  Bo   y  w dy  Bo  y 2  wy 
y 0 x 0 0 2 
d dy  Wb   m
Vemf     Bo  y  w    Bo  y  w  u y   0.100 2   0.1m  0.1m   2 
dt dt  m   s
Vemf  40mV

Alternate Method:
Vemf    u  B  dL    u y a y  Boa z  dxa x
1
 y
2
Vemf  u y Bo  dx  u y Bo  w  y 
1
w y
2

Fig. P4.24
Uniform Plane Wave and Power Transmission

1. A wave with  = 6.0 cm in air is incident on a nonmagnetic, lossless liquid media.

In the liquid, the wavelength is measured as 1.0 cm. What is the wave’s frequency (a)
in air? (b) in the liquid? (c) What is the liquid’s relative permittivity?

up c 3x108 m s
(a) f     5GHz
  0.06m
(b) the frequency doesn’t change with the media (the wavelength does) so f = 5 GHz
(c)
 1 m c
u p   f   5 x109   0.01m   5 x107 
 s s r
2
 3x108 
 r   8 
 36
 0.5 x10 
_____________________________________________________________________
2. Given  = 1.0x10-5 S/m , r = 2.0, r = 50., and f = 10. MHz, find , , , and .

  jr o   j r o     j
jr o

  j r  o
jr o  j 2 10 x106   50   4 x107   j 3948
  j r  o  1x105  j 2 10 x106   2  8.854 x1012   1x105  j1.11x103
Inserting these into the expressions for  and ,
  9.4 x103  j 2.1 1 m ,   9.4 x103 Np m ,   2.1 rad m ,   1880e j 257 
_____________________________________________________________________
3. Suppose in free space, H(x,t) = 100.cos(2x107t – x + /4) az mA/m. Find E(x,t).

H s  0.100e j  x e j a z , a P  a x ,    4 
E s   a P  H s  120 a x  0.100e j  x e j a z  12 e  j  x e j a y
E  12 cos t   x    a y
Since free space is stated,
2 2
   2 30 rad m
 c f
and then
 2  V
E  12 cos  2 x107 t  x  ay
 30 4 m
_____________________________________________________________________
4. A 100 MHz wave in free space propagates in the y direction with an amplitude of 1
V/m. If the electric field vector for this wave has only an az component, find the
instantaneous expression for the electric and magnetic fields.
 rad
From the given information we have   2 f  200 x106 and
s
 2 rad
  ,
up 3 m
 2  V
or E( y, t )  1cos  200 x106 t  y  az .
 3  m
Now to find H.
1 1 1  j y
Es  1e j y a z , H s  a P  Es  a y 1e j y a z  e ax
 120 120
So
1  2  A
H  y, t   cos  200 x106 t  y  ax
120  3  m
or
 2  mA
H  y, t   2.7 cos  200 x106 t  y  ax .
 3  m
_____________________________________________________________________
5. In a lossless, nonmagnetic material with r = 16, H = 100 cos(t – 10y) az mA/m.
Determine the propagation velocity, the angular frequency, and the instantaneous
expression for the electric field intensity.

 c 3x108 m
up     0.75 x108
 r 16 s

  u p    0.75 x108  10   7.5 x108

rad
s
H ( y, t )  100 cos  7.5 x108 t  10 y  a z
mA
m
 j y
H s  0.100e a z ,
120
Es  a P  H s  a y  0.100e j y a z  3 e j y a x
r
E( y, t )  9.4 cos  7.5 x108 t  10 y  a x
V
m
_____________________________________________________________________
6. In a media with properties  = 0.00964 S/m , r = 1.0, r = 100., and f = 100. MHz,
a 1.0 mA/m amplitude magnetic field travels in the +x direction with its field vector in
the z direction. Find the instantaneous form of the related electric field intensity.

 mA   x
 e cos t   x  a z ; H s  H o e e a z
 x  j  x
H  1
 m 
Es  aP  Hs  a x  Hoe x e j xa z   Hoe x e j xa y

j j 2 100 x106  100   4 x107 

   2664e j 30 
  j 0.00964  j 2 100 x10  8.854 x10 
6 12
 1
 j   j   14.8  j 25.7
m
Finally,
E( x, t )  2.66e 15 x cos  200 x106 t  26 x  30  a y
V
m
_____________________________________________________________________
7. In seawater, a propagating electric field is given by E(z,t) = 20.e-z cos(xt –
z + 0.5) ay V/m. Assuming ’’=0, find (a)  and , and (b) the instantaneous form of
H.

For seawater we have r = 72,  = 5, and r = 1.

So: jo  j 7.896, j r o  j 0.004
j
  1.257e j 44.98 
  j
  j   j   4.441  j 4.445 1 m
1
    4.4
m
V V
Es  20e  z e   z e j 0.5radians a y
 20e  z e   z e j 28.6 a y
m m
1 1  20 A
H s  a P  Es  a z  20e  z e   z e j 28.6 a y  e  z e   z e j 28.6 a x
   m
H( z, t )  15.9e 4.4 z cos  2 x106 t  4.4 z  28.6  45  a x
A
m
or with appropriate significant digits:
H( z, t )  16e 4.4 z cos  2 x106 t  4.4 z  16  a x
A
m
_____________________________________________________________________
8. For Nickel ( = 1.45 x 107, r = 600), make a table of , , , up, and  for 1Hz,
1kHz, 1MHz, and 1 GHz.

For Ni we have = 1.45x107S/m, r = 600

     f    f  Hz  600  4 x107 1.45x107   34.35x103 f ( Hz )
= 1/
 j 45
 2 e   18.08 x106 f ( Hz )e j 45 

c m
up   12 x106
r r s
Table
f(Hz)= 1 103 106 109
(Np/m) 185 5860 185x103 5.9x106
(rad/m) 185 5860 185x103 5.9x106
 18e 
j45º
570ej45º 18ej45ºm 0.57ej45º
 5.4mm 170m 5.3m 170nm
up(m/s) 12x106 12x106 12x106 12x106
9. A semi-infinite slab exists for z > 0 with  = 300 S/m, r = 10.2, and r = 1.0. At
the surface (z = 0),
E(0,t) = 1.0 cos( x 106t) ax V/m.
Find the instantaneous expressions for E and H anywhere in the slab.

The general expression for E is: E( z, t )  1.0e  z cos  x106 t   z  a x

V
m
j  j  x106  4 x107   j 3.948
j  j  x106  10.2   8.854 x1012   j 284 x10 6
Here,  (i.e. it is a good conductor), so
1
   f   24.3  
m
 j 45
  2 e  0.115e j 45 

So now we have
E( z, t )  1.0e 24 z cos  x106 t  24 z  a x
V
m
To find B we’ll work in phasors.

1 1 1
Es  1e z e j z a x , H s  a P  Es  a z 1e z e j z a x  e z e j z a y
  
e 24 z cos  x106 t  24 z  45  a y
1 A
H( z, t ) 
0.115 m
H( z, t )  8.7e 24 z cos  x106 t  24 z  45  a y
A
m
_____________________________________________________________________
10. A 600 MHz uniform plane wave incident in the z direction on a thick slab of
Teflon (r = 2.1, r = 1.0) imparts a 1.0 V/m amplitude y-polarized electric field
intensity at the surface. Assuming  = 0 for Teflon, find in the Teflon (a) E(z,t), (b)
H(z,t) and (c) Pav.


E(0, t )  1cos 2  600 x106  t   z a y  V
m
V
E( z , t )  1e  z cos t   z  a y
m
Teflon:  = 0 so  = 0,
 2  600 x106  rad
and      r  8
2.1  18.2
c 3x10 m
(a) E( z , t )  1cos 1.2 x109 t  18.2 z  a y
V
m
1 2.1 V
(b) H s  a P  Es  a z 1e j z a y ,
 120 m
H ( z , t )  3.8cos 1.2 x109 t  18.2 z  a x
mA
m
 1 1 2.1
2
mW
(c) Pavg  a z  1.9 2 a z
2 120 m
_____________________________________________________________________
11. A 200 MHz uniform plane wave incident on a thick copper slab imparts a 1.0
mV/m amplitude at the surface. How much power passes through a square meter at
the surface? How much power passes through a square meter area 10. m beneath the
surface?

mV 1 Eo 2
f  200MHz , Eo  1 , Pavg 
m 2 
 j 45 Np
Cu:   2 e ,    f   214 x103 , so   5.22e j 45 m
 m
1 10 
3 2
W
Pavg  3
 96 2 ; P  Pavg S  96W
2 5.22 x10 m
Now at 10 m beneath the surface, we have
V
E ( z  10 m)  Eo e  (10  m )  103 e  (214 x10 )(10  m )  118 x106
3

m
1 118 x10 
6 2
W
Pavg  3
 1.3 2 ; P  1.3W
2 5.22 x10 m
_____________________________________________________________________
12. Given E(z,t) = 10.cos(t-z)ax - 20.cos(t-z-45)ay V/m, find the polarization
and handedness.

The field can be rewritten as E(z,t) = 10.cos(t-z)ax + 20.cos(t-z-45-180°)ay

or E(z,t) = 10.cos(t-z)ax + 20.cos(t-z+135°)ay
_____________________________________________________________________
13. Given

E( z, t )  Exo cos t   z  ax  E yo cos t   z    ay ,

we say that Ey leads Ex for 0 <  < 180, and that Ey lags Ex when –180 <  < 0.
Determine the handedness for each of these two cases.

For 0 <  < 180°, we have LHP

For 180° <  < 360°, we have RHP

_____________________________________________________________________
14. Suppose a UPW in air carrying an average power density of 100 mW/m2 is
normally incident on a nonmagnetic material with r = 11. What is the time-averaged
power density of the reflected and transmitted waves?

 1  1
120  
;    0.537
11
1  o  120; 2 
11  
 1  1
 11 
  1   0.463
mW
  Pavg  28.8
r i 2
Pavg
m2
2
1 Exo mW
P 
t
  2 Pavg
i
11  71.2 2
2 2
avg
m
_____________________________________________________________________
15. A UPW in a lossless nonmagnetic r = 16 media (for z < 0) is given by
E(z,t) = 10.cos(t-1z)ax + 20.cos(t-1z+/3)ay V/m.
This is incident on a lossless media characterized by r = 12, r = 6.0 (for z > 0). Find
the instantaneous expressions for the reflected and transmitted electric field
intensities.

Eis  10e j1z a x  20e j1z e j 3a y

Ers  10e j1za x  20e j1z e j 3a y
120 12
1   30; 2  120  120 2
16 6
 
  2 1  0.700;   1    1.70
2  1
Ers  7e j1z a x  14e j1z e j 3a y
  V
E(rz ,t )  7 cos t  1 z  a x  14 cos  t  1 z   a y
 3 m
 j 2 z  j 2 z j 3  j 2 z
Es  10 e
t
a x  20 e e a y , or Es  17e
t
a x  34e j2 z e j 3a y , so
  V
Et( z ,t )  17 cos t   2 z  a x  34 cos  t   2 z   a y .
 3 m
_____________________________________________________________________
16. The wave Ei = 10.cos(2x 108t - 1z) ax V/m is incident from air onto a copper
conductor. Find Er, Et and the time-averaged power density transmitted at the
surface.

For copper we have

 2 j 45
2  2 e 
2

where  2   f 2 2   108  4 x107  5.8 x107   151x103

Np
 2
m
so 2  3.7e j 45 m
22 2
We find   1, and  =  2  19.6 x106 e j 45
1   2 1
So E = -10.cos(2x 10 t + 1z) ax V/m
r 8

V
, and Et  196e  2 z cos t   2 z  45  a x
V
Ets  196  e  2 z e  j2 z e j 45 a x
m m
1 196 x10 V m 
6 2
W
P t
 cos  45  a z  3.7 2 a z .
2  3.7 x10  
3
avg
m
17. A wave specified by Ei = 100.cos(x107t-1z)ax V/m is incident from air (at z < 0)
to a nonmagnetic media (z > 0,  = 0.050 S/m, r = 9.0). Find Er, Et and SWR. Also
find the average power densities for the incident, reflected and transmitted waves.

rad 2  rad
1  120,    x107 so f  5 x106 Hz, 1    0.105
s 1 c m
In this problem we find in medium 2 (z > 0) that  = 0.0025 and  = 0.05. These
values are too close to allow for simplifying assumptions. Using (5.13) and (5.31),
we calculate:
Np rad
 2  0.969 ,  2  1.019 , 2  28.1e j 43.6  .
m m
Then,
  1 
  2 1  0.898e j174 , SWR   18.6,   1    0.141e j 40.8
2  1 1 
V
Eis  100e  j 1z a x
m
V V
Ers  100e  j 1z a x  89.8e  j 1z e j174 a x ,
m m
so Er ( z , t )  89.8cos  x107 t  0.105 z  174  a x
V
.
m
V V
Ets  100 e j2 z a x  14.1e j 2 z e j 40.8 a x ,
m m
so Et ( z , t )  14.1cos  x107 t  1.02 z  40.8  a x
V
.
m
14.1  j 43.6  j 2 z j 40.8 A A
H ts  e e e a y  0.502e  j 2 z e  j 2.8 a y
28.1 m m
P t  14.1 0.502  cos  40.8  2.8  a z  2.6 2 a z
1 W
2 m
100   13.3 W a
2

P 
i

2 120 
z
m2
  89.8 
2
W
Pr   10.7 2  -a z 
2 120  m
(check: 13.3 W/m2 = 10.7 W/m2 + 2.6 W/m2)
_____________________________________________________________________
18. A 100 MHz TM polarized wave with amplitude 1.0 V/m is obliquely incident
from air (z < 0) onto a slab of lossless, nonmagnetic material with r = 25 (z > 0). The
angle of incidence is 40. Calculate (a) the angle of transmission, (b) the reflection
and transmission coefficients, and (c) the incident, reflected and transmitted fields.

(a) The material parameters in this problem are the same as for P5.48. So, once again
we have t = 7.4°. Also, 1 = 2.09 rad/m and 2 = 10.45 rad/m.
(b)
 cos t  1 cos i
TM  2  0.589
2 cos t  1 cos i
 22 cos i
 TM   0.318
2 cos t  1 cos i
(c)
Incident:
Eis  1e j1.34 x e j1.60 z  cos 40 a x  sin 40 a z 
V
Ei ( z , t )   0.766a x  0.643a z  cos t  1.34 x  1.60 z 
m
1  j1.34 x  j1.60 z A
H is  e e ay
120 m
mA
H i ( z, t )  2.65cos t  1.34 x  1.60 z  a y
m
Reflected:
Ers  0.589e j1.34 x e j1.60 z  cos 40 a x  sin 40 a z 
V
Er ( z , t )   0.452a x  0.379a z  cos t  1.34 x  1.60 z 
m
0.589  j1.34 x  j1.60 z A
H rs  e e ay
120 m
mA
H r ( z , t )  1.56 cos t  1.34 x  1.60 z  a y
m
transmitted:
Ets  0.318e j1.35 x e j10.4 z  cos 7.4 a x  sin 7.4 a z 
V
Et ( z , t )   0.315a x  0.041a z  cos t  1.35 x  10.4 z 
m
mA
H t ( z, t )  4.22 cos t  1.35 x  10.4 z  a y
m
_____________________________________________________________________