Imprtant Problems EE339
Imprtant Problems EE339
1017
2.56 8.854 x10 12
9
t o e 441 x10 t
MLP0402
Plot rhov vs time
for a charge placed in silver
4-2
grid on
subplot(2,1,2)
loglog(t,rho)
xlabel('t (s)')
ylabel('rhov (C/m^3)')
grid on
Fig. P4.02
P4.3: A current density is given by J = e-.01ta A/m2. Find the charge density after 10
seconds if it has an initial value of zero.
v 1 2 0.01t
e
2e 0.01t
t
t
0.01t
d v 2e
dt; v 200e 0.01t C
gJ
At t = 0, v = 0 = 200 + C, C = -200.
So we have
v 200 e 0.01(10) 1 19
C
m3
4-3
12
11.8 8.854 x10
237.4 x109 sec
4
4.4 x10
Q
60C
s
; area=4 r 2 1.257 x10 3 m 2 ; s
47.8 mC 2
3 2
m
area
1.257
x
10
m
(c)
2. Wave Fundamentals
P4.5: A propagating electric field is given by
(a) Determine the attenuation constant, the wave frequency, the wavelength, the
propagation velocity and the phase shift. (b) How far must the wave travel before its
amplitude is reduced to 1.0 V/m?
(a)
;
45o
s
4
(b)
1
; z 460m
100
E ( z ) 100e0.01z 1; 0.01z ln
P4.6: A 10.0 MHz magnetic field travels in a fluid for which the propagation velocity is
1.0x108 m/sec. Initially, we have H(0,0)=2.0 ax A/m. The amplitude drops to 1.0 A/m
after the wave travels 5.0 meters in the y direction. Find the general expression for this
wave.
H( y , t ) H o e y cos t y a x
m
A
H(0,0) 2.0a x ; H o 2.0, 0
m
H ( y 5) 1 H o e y 2.0e (5) ; solving we get 0.14
Finally,
A
H ( y, t ) 2.0e0.14 y cos 20 x106 t 0.2 y a x
m
A
m
4-4
P4.7: MATLAB: Modify the simple wave program in MATLAB 4.1 to include
attenuation. Generate a plot for the case where the amplitude is 4 V/m, the attenuation
constant is 0.001 Np/m, and the frequency is 1 MHz. Take your snapshot in time at 0
seconds, and let your phase shift be 0.
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
M-File: MLP0407
This is a modification of ML0401 that
includes attenuation. It plots a wave (in vac)
versus position (in z direction) for a fixed time.
Wentworth, 1/19/03
Variables:
atten
Eo
f
omega
t
phi
phir
c
lambda
B
E
z
clc
clear
attenuation constant
wave amplitude (V/m)
frequency (Hz)
angular frequency (rad/sec)
time snapshot (sec)
phase constant (degrees)
phase constant (radians)
speed of light in vacuum (m/s)
wavelength (m)
phase constant (1/m)
electric field intensity
position
% Initialize Variables
atten=0.001;
Eo=4;
f=1e6;
t=0;
phi=0;
phir=phi*pi/180;
c=2.998e8;
lambda=c/f;
B=2*pi/lambda;
omega=2*pi*f;
4-5
% Perform Calculation
z=0:4*lambda/100:4*lambda;
E=Eo*exp(-atten.*z).*cos(omega*t-B*z+phir);
% Generate the Plot
plot(z,E)
axis('tight')
%sets axes min & max data values
grid
xlabel('z(m)')
ylabel('E(V/m)')
Fig. P4.7
P4.8: MATLAB: Modify the traveling wave program in MATLAB 4.2 to include
attenuation. Use the parameters from problem P4.7, except for the fixed time of course.
%
%
%
%
%
%
%
%
%
%
%
%
%
M-File: MLP0408
This program illustrates a traveling wave including
attenuation. It modifies ML0402.
Wentworth, 1/19/03
Variables:
atten
Eo
f
omega
t
attenuation (Np/m)
wave amplitude (V/m)
frequency (Hz)
angular frequency (rad/sec)
time snapshot (sec)
4-6
%
%
%
%
%
%
%
phi
phir
c
lambda
B
E
z
phase constant
phase constant
speed of light
wavelength (m)
phase constant
electric field
position
(degrees)
(radians)
in vacuum (m/s)
clc
clear
(1/m)
intensity
% Initialize Variables
atten=0.001
Eo=4;
f=1e6;
t=1;
phi=0;
phir=phi*pi/180;
c=2.998e8;
lambda=c/f;
B=2*pi/lambda;
omega=2*pi*f;
% Perform Calculation
z=0:4*lambda/100:4*lambda;
E=Eo.*exp(-atten.*z).*cos(omega*t-B*z+phir);
% Generate a Reference Frame
plot(z,E)
axis([0 4*lambda -2*Eo 2*Eo])
grid
xlabel('z(m)')
ylabel('E(V/m)')
pause
% Make the Movie
t=0:1/(40*f):1/f;
for n=1:40;
E=Eo.*exp(-atten.*z).*cos(omega*t(n)-B*z+phir);
plot(z,E)
axis([0 4*lambda -2*Eo 2*Eo])
grid
4-7
title('General Wave Equation');
xlabel('z(m)');
ylabel('E(V/m)');
M(:,1)=getframe;
end
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.
4-8
Fig. P4.10
P4.11: Considering Figure 4.7, suppose the area of a single loop of the pair is 100 cm 2,
and the magnetic flux density is constant over the area of the loops but changes with time
t
as B Bo e a z , where B = 4.0 mWb/m2 and = 0.30 Np/sec. Determine V at 1, 10,
o
B
dB
gdS;
Bo e t a z ;
t
dt
Vemf 2 Bo e t S
1m
VR 2 Bo Se 2 0.30 4 x10 100cm
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
t
We have
i2
N1
N
i1
v2 2 v1
N 2 and
N1
4-9
N2
v1
N
v1
v2 N1
Z1 , Z 2
2
i1
i2 N1
N1
N i1
2
N
Z1 1
N 2
Z1
Z2
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.
4 0.04m)
1 l
1m
R
3.5m
7
A 5.8 x10 0.0005m 2
0.04
0.04
dB
Wb
dB
Wb
1 2 a z ; Vemf
gdS 1 2 dx dy 1.6mV
dt
ms
dt
ms 0
0
1.6mV
I ind
0.46 A
3.5m
(note that this answer has no time dependence)
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):
I
H
4
a z
z
2
2
a
I
2
0.46
1
-a z
2
2a
1
A
a z 2.59a z ;
m
2 0.02
4-10
So
B 4o H 13
Wb
az .
m2
P4.14: The mean length around a nickel core of a transformer like the one shown in
Figure 4.11 is 16 cm, and its cross sectional area is 1 cm 2. There are 30 turns on the
primary side and 45 on the secondary side. If the current on the primary side is 1.0
sin20x106t mA, (a) calculate the amplitude of the magnetic flux in the core in the
absence of the output windings. (b) With the output windings in place, calculate i2.
(a)
Vm
l
; Vm N1i1 ; R
R
r o A
(b)
i2
1cm
2
m2
100cm
14 x10 9 Wb
N1
30
2
i1
1sin 20 x106 t mA sin 20 x106 t mA
N2
45
3
P4.15: A triangular wire loop has its vertices at the points (2, 0, 0), (0, 3, 0) and (0, 0, 4),
with dimensions in meters. A time-varying magnetic field is given by B = 4t ay Wb/m2 (t
in seconds). If the wire has a total distributed resistance of 2 , calculate the induced
current and indicate its direction in a carefully drawn sketch.
B
gdS,
t
B
Wb
4a y 2
t
ms
1
MN
S MN
2
MN
Vemf
M 2a x 3a y , N 2a x 4a z
S 6a x 4a y 3a z m 2
Vemf = -(4)(4)=-16V
Vemf 16V
I
8A
R
2
Fig. P4.15
4-11
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,
h
o I
m o I
Vemf
u
B
g
d
L
1
a
a
g
dz
a
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 Bgu y dxa z 3 2 3 0.2m 1.8V
m
s
1.8V
I
0.18 A
10
(clockwise when viewed from the +z axis)
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 weve assumed dS = -dSaz to get iind and Vemf as shown. Our approach will be to
find , then Vemf = -d/dt.
4-12
BgdS 1a z g d d a z
r
2 d r 2
0
d
dr
1
2 r
2 r 2V
dt
dt
r
Vemf 2V
I
2V
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.
u B gdL Bo d
a
Bo a 2
2
Bo a 2
u B gdL Bo 0 d 2
a
Fig. P4.19a
Fig. P4.19b
4-13
P4.20: In Figure 4.16, replace the rectangular loop with a circular one of radius a and
rederive Vemf.
Bo a 2 sin ;
Vemf
d flux
dt
Bo a 2 sin
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 gdL
u
d
a a
dt
Bo l
a Bo a z gd a
2
0
l
Vemf
1
rad
Wb
2 Vs
Vemf 120
0.1 2 0.06m
2
s
m
Wb
68mV
Fig. P4.21
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 B o = 100. mWb/m2,
a = 6.0 cm, and R = 50. , determine I if the bar rotates at 1.0 revolution per second.
d
a a , dL d a
dt
Figure P4.22 indicates one of the paths for the circulation integral.
Vemf
u B gdL; u
4-14
a
Vemf
a Bo a z gd a Bo d
0
Bo a 2
2
Bo a
1 rev
rad
Wb
1 2 100 x10 3 2
R
2R
2 s
rev
m
I 22.6 A
I
Vemf
m
0.06
1 Vs A
50 Wb V
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 V emf 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:
BgdS Boa z gdxdya z
Now, notice that x and y are not independent and are in fact related: x=y+w
So we have
y yw
y
1
Bo dxdy Bo y w dy Bo y 2 wy
2
y 0 x 0
0
4-15
d
dy
Wb
m
Bo y w
Bo y w u y 0.100 2 0.1m 0.1m 2
dt
dt
m
s
40mV
Vemf
Vemf
Alternate Method:
Vemf u B gdL u y a y Boa z gdxa x
1
y
2
Vemf u y Bo
dx u y Bo w y
1
w y
2
Fig. P4.24
5. Displacement Current
P4.25: Suppose a vector field is given as
A 3x 2 yz 3a x .
Verify that the divergence of the curl of this vector field is equal to zero.
verify that g A 0
ax
A
x
2
3 x yz 3
ay
az
3 x 2 yz 3 a 3 x 2 yz 3 a
y
z
y
z
z
y
0
0
9 x 2 yz 2a y 3x 2 z 3a z
g A
9 x 2 yz 2 3x 2 z 3 9 x 2 z 2 9 x 2 z 2 0
y
z
4-16
P4.26: Suppose a vector field is given by
A 2 cos a z .
Verify that the divergence of the curl of this vector field is equal to zero.
1
1
1
g A
sin
2 cos
1
2 sin 2sin 0
P4.27: A pair of 60 cm2 area plates are separated by a 2.0 mm thick layer of ideal
dielectric characterized by r = 9.0. If a voltage v(t) = 1.0 sin (2x103t) V is placed across
the plates, determine the displacement current.
12
2
dv
S 9 8.854 x10 F m 60cm 10mm 1m
id C ; C
239 pF
dt
d
cm 100cm
2mm
dv
V
2 x103 cos 2 x103 t
dt
s
3
id 1.5cos 2 x10 t A
P4.28: Plot the loss tangent of seawater ( = 4 S/m and r = 81) versus log of frequency
from 1 Hz to 1 GHz. At what frequency is the magnitude of the displacement current
density equal to the magnitude of the conduction current density?
tan
4 S m
889 x106
2 f ( Hz ) 81 8.854 x1012 F m
f
4-17
end
loglog(f,tand)
xlabel('frequency (Hz)')
ylabel('loss tangent')
grid on
Fig. P4.28
P4.29: A 1.0 m long coaxial cable of inner conductor diameter 2.0 mm and outer
conductor diameter 6.0 mm is filled with an ideal dielectric with r = 10.2. A voltage v(t)
= 10.cos(6x106 t) mV is placed on the inner conductor and the outer conductor is
grounded. Neglecting fringing fields at the ends of the coax, find the displacement
current between the inner and outer conductors.
C
Q
, so Q Cv (t )
V
DgdS D a
2 l
ln b a
so Q
2 l Vo cos t
ln b a
g d dza 2 l D Q
2 l Vo cos t Vo cos t
2 l ln b a ln b a
so
D Vo sin t
a J d
t
ln b a
and D
Vo cos t
a .
ln b a
Vo sin t 1 2 l
2 l Vo sin t
id J d gdS
d dz
ln b a 0
ln b a
0
4-18
Now evaluate id with the given parameters:
2 6 x106 rad s 10.2 8.854 x10 12 F m 1m 0.010V sin t C As
id
ln 0.006 0.002
FV C
id 97 sin 6 x106 t A
5e 2 zt a 10te 2 zt a
E o
y
y
x
y
z
z
t
5e 2 zt
0
0
10t 2 zt
10
dH
e dta y , H= te 2 zt dta y
o
o
udv uv vdu
2 zt
This integral is solved by parts
where we let u t and dv e dt.
We arrive at:
10t 2 zt
10 2 zt
A
H
e
e ay
2
4o z
m
2o z
V
2
m
1
, f 2 MHz ,
, 50m, u p 100 x106
; r 9
m
H
E Eo sin t y a x o
t
Eo
Eo
H
sin t y dt
cos t y a x
o
o
Eo
H( y, t )
cos t y a x
o
A
H ( y , t ) 0.796 cos 4 x106 t 0.126 y a x
m
Inserting the given values we find:
Eo 100
4-19
H( z, t ) 20.sin 108 t z ax
A
m
1
x108 50 MHz
2
8
c
3x10
2 rad
6m,
8
f 0.5 x10
3 m
D
H
(no J c )
t
2 f x108 , f
ax
H
x
H o sin t z
ay
az
H sin t z a
y
y
z
z o
0
0
H o cos t z a y o
E
t
Ho
H
cos t z dta y o sin t z a y
o
o
Ho
V
E( z , t )
sin t z a y
o
m
so
E( z , t ) 7.5sin x108 t
3
dE E =
kV
z a y
m
P4.33: Find the instantaneous expression for E for the magnetic field of problem P4.6.
From P4.6 we have H 2.0e 0.14 y cos 20 x106 t 0.2 y
as
H H o e y cos t y
A
ax
m . We can write this wave
A
ax
m , and can reinsert the values at the end.
c
E
up
,
.
r
t Assuming a nonmagnetic medium, we have
First we have
or r = 9. Now we evaluate the curl of H:
ax
a y a z
H
H o e y cos t y a z .
x
y z
y
y
0
H o e cos t y 0
H
4-20
H o y
H o y
E
e sin t y a z
e cos t y a z .
56,
6
Ho
2 0.2
251.
So
E 56e 0.14 y sin 20 x106 t 0.2 y a z 251e 0.14 y cos 20 x10 6 t 0.2 y a z
V
.
m
E e y a z A b cos t y c sin t y ,
we can let
where b=cosy, c=siny, and tany=c/b. Now, since b/c = 251/56, then y=12.6 and A=257.
Thus
V
E( y, t ) 257e0.14 y cos 20 x106 t 0.2 y 12.6o a z .
m
This result is easier to obtain using the phasor approach of the next section.
P4.34: Given, at some point distant from a point source at the origin in free space,
V
E(r , t ) 8.0 cos 9 106 t r a ,
m
find frequency, phase constant and H(r,t).
2 f , f 4.5MHz ,
E
B
H
o
t
t
c
2
rad
66.7m,
0.094
f
4-21
E
1
rE a , E Eo cos t r ,
r r
2 f 3 x108 ; f 48MHz
m
c
u p 0.3 x108
; r 100
s
r
r
rad
10
c
m
D
E
H
t
t
1
1
H
H o sin t 10 a z H o sin t H o cos t a z
In the far-field, we will ignore the first term (has 1/) leaving us with
H H o cos t a z
Ho
dE E
cos t dta z
Now,
Ho
E
sin t a z .
4-22
E( , t ) 226sin(3 x108 t 10 )a z
V
.
m
u
c
m
2
rad
1.73x108 , p 34.6m,
0.181
s
f
m
r
f 5MHz , u p
E 9 sin t y a x 4 cos t y a z o
9
4
cos t y a x
sin t y a z
o
o
Inserting the given values we arrive at:
H
t
H ( y, t ) 41cos t y a x 18sin t y a z
mA
m
j B s e j t .
Re B s e j t Re
t
Re B s e j t Re
B s e j t
t
t
Re B s e j t Re B s cos(t ) jB s sin(t )
t
t
B s cos(t ) B s sin(t )
t
Now expanding the right side we have:
Re
B s e jt Re jB s e j t Re j B s cos t j sin t
t
So since we have :
j B s e jt ,
Re B s e j t Re
B s e j t Re
t
then
V
m.
4-23
Re B s e j t
t
Re B
j t
jB s e jt
e Re
t
P4.38: Derive the differential phasor form of (a) Gausss Law, and (b) Amperes Circuit
Law.
gD v , gD x, y, z , t v x, y, z, t
vs e jt ; Re gD s e jt Re vs e jt ;
gRe Ds e j t Re
gDs vs
H
D
E
t
t
Re H s e jt
H s j E s
Re E s e j t ;
t
Re H s e jt Re j E s e jt
Eo j y
Eo
e a x and H ( y, t )
cos t y a x .
o
o
Plugging in the values from P4.31, we find
A
H ( y , t ) 0.796 cos 4 x106 t 0.126 y a x .
m
Hs
ay
az
jH e j z a
o
y
y
z z
0
0
jH o e j z a y j 2 H o e j z a y H o e j z a y j o E s
z
H o e j z a y
H o j z
Ho
V
Es
j
e a y , E( z , t )
sin t z a y
j o
o
o
m
Plugging in the values from P4.32, we find
Hs
4-24
kV
E( z , t ) 7.5sin 108 t z a y
.
3
m
20e j z a 10e j z a
Es
x
y
x
y
z
z
z
10e j z 20e j z
0
E s j 20e j z a x j 10e j z a y j o H s
Hs
j 20 j z
j 10 j z
20 j z
10 j z
e ax
e ay
e ax
e ay
jo
j o
o
o
x106
rad
;
10.47 x10 3
8
u p 3 x10
m
Inserting the appropriate values we then find:
mA
H s 53e j z a x 27e j z a y
,
m
up
mA
m
ax
Es
x
j 4e j y
ay
az
j 9 e j y a x j 2 4 e j y a z
y
z
0
9e j y
j 9 e j y a x 4 e j y a z jo H s
j9 j y
4 j y
9 j y
4 j y
Hs
e ax
e az
e ax j
e az
jo
jo
o
o
Inserting the appropriate values we then find:
V
.
m
4-25
up
c
m
rad
1.73 x108 ,
0.181
s
up
m
r
H ( y, t ) 41cos t y a x 18sin t y a z
mA
m
P4.43: MATLAB: In MATLAB 4.4, a polar plot of the phasor corresponded to a location
on a sine wave for a particular time, at a fixed position in space. You can also make a
polar phasor plot for a snapshot in time, where you change position. Modify MATLAB
4.4 to provide an animation of the phasor versus sine wave as you change the position.
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
M-File: MLP0443
This program modifies ML0404, generating a
position plot animation synchronized with
a polar plot.
Wentworth, 1/19/03
Variables:
Eo
E
E
f
theta
theta1
T
t
t1
clc
clear
clf
4-26
% Perform Calculation
z=0:lambda/100:2*lambda;
theta=-beta*z;
E=Eo*cos(theta);
% Generate the Plot
subplot(211),plot(z,E,0,Eo,'ro');
subplot(212),polar(0,Eo,'ro');
pause
%Make the Movie
z1=0:lambda/50:2*lambda;
for n=1:100
theta1(n)=-beta*z1(n);
E1(n)=Eo*cos(theta1(n));
subplot(211),plot(z,E,z1(n),E1(n),'ro');
subplot(212),polar(theta1(n),Eo,'ro');
M(:,1)=getframe;
end
M-File: MLP0444
Same as MLP0443 but add attenuation.
This program modifies ML0404, generating a
Run the
4-27
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%
clc
clear
clf
4-28
z1=0:lambda/50:2*lambda;
for n=1:100
theta1(n)=-beta*z1(n);
E2(n)=Eo.*exp(-atten.*z1(n));
E1(n)=E2(n).*cos(theta1(n));
E2(n)=Eo.*exp(-atten.*z1(n));
subplot(211),plot(z,E,z1(n),E1(n),'ro');
subplot(212),polar(theta1(n),E2(n),'ro');
M(:,1)=getframe;
end