Nayfeh Solutions
Nayfeh Solutions
b. Thus charge density p = 0 for p> be
a « But dq = oda = ordr® so
le- 3]
(pe, + W)/2 4
° an
al rodrd@ _ oW '
eae f J 7 = Ter + It doesn't depend on py
° o/2 ° °
‘0
a bax 2
2.17 a) o= tre, ee « Note that x = rytan@, dx = rysec*6d6
and r = r,sec8, thus ¢ = 7
1 ome,
15terete aera: ones eta ast remnant tn =n rhe earn eam ener
don r+ constant
b) In the £ >> r, limit we have @ = -
1 2ne, 1
c) For ry, ry << £, we have Vjy = r,) - (ry) =
i . This is, as we would expect, the same as the
°
potential difference of two points in the presence of an
infinite line of charge.
eG
30° rsind a6"
+
2.18 a) We use E = -V where V =
spherical coordinates. Thus E = t(2ar~>cos® + br@2) +
8(ar~3sine).
b) We use a = VeE. Consider first ? # 0.
°°
z_ 1 38 (2 ec)
WE = “5 oe (EL) * reine 96 (2198 Eg
Pea ie -1 toe oe
€, "(28 (lar “cos® + b) + 755 gq (ar sin’)
= -2ar “cos® + 2ar “cosé = 0
Thus there is no charge at f # 0. We now inspect the equation
as @ +0. Also we note that Ve#/r? = -V*(1/r) = 4n6(#) where
from example 1.4, the Dirac 6-function 6(f) = 0 for r # 0 and
[s(2)43F = 1 if volume includes %= 0. This would have been
obvious if we noticed % = b/r is a point charge potential. We
should also recognize (a/r2)cos® as a dipole potential of dipole
moment 4nega, in the 2 direction and placed at r = 0. See
16example 2.9.
©) We use the following in one step: 2= df = §.(-¥6) = -v7s.
(See Ex. 3-20). z
é
€
oe
a °o >
2.19 a) Using Gauss' law § Bead = nt ives LA “get.
Q
b) Using Gauss' law gives i, = om r
0
“2
€) Take pf} = -p, gives the cavity, thus #) = 7-27". Adding
+e + MY
gives B= Ey + Ey = 52 zoe
0
ee ES ea = ly
a) @=-ffedt =—22 f az= =a 2 r' cose’. One
3e, 70 36. 70
° oo °
can also find the potential by finding the potential due to the
>
sphere and the cavity. (Same procedure used for E.)
2.20 We notice that this problem has symmetry about the y axis,
so it is natural to set up the problem as in example 2.9, with
the axis redefined: : P
17-L @-
The potential at point P is clearly ¢, Tre, & da}. Now,
we note that
Vr? + (2/2)? = recosa, x" = fx? + (2/2)? + recosa
We are to assume £ b. We also have the boundary conditions:
* :
Pease x gives o3(*) Seba ~E pcos
2) At p = bwe have (hb) = &(b) and D,(b)-B = 4(h) +8
3) At p = a we have (a) = (a) oe constant
Thus we write the following expansions for
A c
+ $6, 0 = VV + (S- Sb.
eV ¢ bo)cosd, %, = V ¢ & p)cosé. We now apply
these conditions. Conditions two and three give
gb) »
e(h - a/b”) = e,(Ey + A/b?), (Bd + A/D) = (Clb = BG)
V, = (A/a + ba)cosé + V, or A/a = -ba. Solving these equations
simultaneously sive
49B= 28,b2/(a2(1 - k) - b%(k + 1))
A= -28,b2a?/(a(1 - k) - b%(k + 1)
wEjb'(a2 - b? + k(a? + b2))/(a2(1 - k) + b*(K + 1)
b) 946) = Dea. From (a) of = ~e(B ~ A/a”)cosd = -2eBcos¢ with
B given in (a).
4.15 a) First consider the case €, > €)+ The electric field
will exert a torque on the cylinder due to the induced
polarization: the cylinder will orientate itself to create zero
net torque. The two such positions are with the axis of the
cylinder aligned (1) and perpendicular (2) to the field. In (1)
the cylinder will be in a stable equilibrium. In (2), the
cylinder will be in unstable equilibrium. For a thin disk, its
axis will be perpendicular to the field in a stable equilibrium.
b) For the case €, < €2, we have stable equilibrium when the
axis of the cylinder is perpendicular to the field, and the axis
of the thin disk parallel to the field.
4.16 a) Since the polarization of the material is permanent,
then the permeability of the unpolarized material is
irrelevant. The field inside the cavity is dictated by the
polarization only. We can use the polarization charge
ea te
technique. Now pp, = -vep = 0, Oye Ee Per = Pcos@ on the
surface of the cavity. ‘There are no other charges because the
material is infinite. ‘The electric field corresponding to this
50charge is uniform inside the cavity and it is equal to
E= P/3e,. [See Exs. 4.2, 2.17, and 3.6-]
b) When the polarization is not permanent, then we have to
solve the problem in an external field with the permitivity e of
the material becoming relevant. We have already solved a
similar problem in Ex. 4.8. In this example we have a
dielectric sphere of permitivity ¢, place in an infinite medium
+ s
of permitivity €, and with an electric field E = Ej2. Now we
take ey €, = €, and E, = P/(e- €,). Thus from tq. 4.80
we have
: : _ 3 7
“sphere — + 2€, ‘o”
changing this result to the present case gives
>
eee eeere
2
+ 2e Ce - €,)
z =
hole — Ce,
4.17 a) See Ex. 3.5 to get # = Asse 47/44?
b) The potential is (x,y,z) = (1/4ne,) {a/(x - d)* +
Vey. vhus
yrs at) ee q(x + ay? + y+
= -(36/ax)| =
oe for x=0
“n
(alan) (x = d)/(Cx = a)? + y2 + 24) 91? +
Ge + a CCx ay? + yh & 22)! HL
-(2qd/4m) A/a? + y* + 22)3/2) ow ~Cqd/2my/(o? + d2)3/2, Now
=
om da where da = 2npdp, thus
0
5122 2,2
y- ad yp _2nede ge tt We do the
Bre, (0 a2)? M0 Cp? + a3
integrat?on by letting u = p* + d*, hence du = 2pdp and the
ae 242 ‘i ee
force becomes 7—— f Stee ao | -—15%
0 0 Qu (p24 de) 0 lene d
which is the same as in (a).
4.18 Because the potential has to be finite at r = 0 and r =
and because the normal component of P on the surface of the
sphere P, = P,cos@ involves cos@ only, we keep only the cos®
terms in the expansion of the potential inside and outside %)
and $5:
+9) = Ajrcosé,
Taking ) = ) at r= R gives A)R = B/R2, Matching Dj,
at r= R, we get
Sse ooo 3
(ek, + P)en = e,Eyen or -(Ay + Po/eg) = 2B,/R°-
Solving for Ay and By we get Ay = P,/3e,, By = R’P,/3e,, hence
2
= R°P,cos@/3e,r° +
= Pyreos8/3€,5
4.19 We use expansions for the potentials of the form of
Eq. 3-65. The potential inside the cylinder @, has to be finite
52at p = 0, and the potential outside the cylinder $) has to
vanish at p = *. Moreover because at p= a, the polarization is
Pycoss, then we will keep the cos$ terms only. ‘Thus
@) = Apcosng, 4. Atcos¢/p
Taking ¢) = 6) at p = a gives Aja = A'/a. Matching Din = Yon at
2n
8 a 5
p= a gives (c.f) + Py)en = e,£yen or -(e,A + Py) =
egAt/a*. Solving for A and At gives: A -P,/2e, and
2
At = -a*p,/2e,, and hence
PoP a’P
4, = - go cose, ~~ Tap con
Px a’p . .
[pcos + gsing]
2e,e
4.20 We use the results of section 4.7 to introduce the charge
q' and q” as images of q, and the charges -q' and -q" as images
of the charge
» as shown in the figure.
&5€,Thus the images are: an image dipole $' in material €, at
distance d, and an image dipole $" in material €, at the
location of the real dipole. ‘The moments of the image dipoles
are
& -€ 2e.
2a) ‘f 2
ee re Dl De econ Pla
te it %
4.21 a) , or
* X+Ax
ae 1 - a
Co. ee Je xn ‘dx
x
a n xt+Ax =
Xp’ x nn ecko
=F] = Ge + wo" - x] = SP. The operations of
averaging and differentiation are indeed interchangeable.
+ * > is > >
4.22 a) Since f = x°x then E = x*x and D = ek since
1 Coulomb. At one face of the differential cube, we have,
+ :
cE = ex2x, and at the opposite face we have
oro
u
* >
e(E + (8E/dx]Ax. Thus the net flux of
54Bee = AxsAydz = 2xe AxdyAz. This net flux must be equated to
the total charge inside the cube, that is p¢ AxAyAz = 2xe AxAyAz
or p = 2xe+C/m.
:
b) pe = VoD. Thus p¢ = €8x2/ax = 2ex C/m>.
4.23 a) Because of the symmetry, the fields and the potentials
> +
will depend on x only. ‘Thus VeD = p and VeD = 0 give
> a > «
D) = (ox + b)x for x ¢d, Dy = fx for d = + pk/kT = 1.3 x 107%.
b) B= N 2 = 7 x 1077 c/n?
©) = p, thus p of the
block is N pV = 4.95 x 1079 C/m
5-3 a) Baipore = (p/4mR?) (2e0s er + sind 8}. ‘The induced
dipole is 8, = af = (pa/4ne,k>) (2cosor + sindd}
b) u = Bef = - P22 cose + sine8}).—P-{2coser + sindd}
a) 3
me SR 4neR
= -a( 25)? (4eos?9 + sin®e) = -a(p/4negk?)7(1 + 3eos70)
4ne.R
562 3,2 7 2
: poa(2n)/(4ne R°) J (1 + 3cos“@)sinede
U(8)da
c) = =-
Jaa
7
(4m)
= -0p2/ar2e2n6
5-4 a) Let 6; and 0) be the angles py and py make with a line
joining the molecules. Then we know the interaction energy
= (1/413) (308) *5)(By*2) ~ By eB). Now
Pitt = cos®), pyr = c088,
Pit by = sin®,cos$,sin8,cosd, + sin®)sing,sind,sing, +
€080,cos8, = sind,sin8,(cos(4, - $)) + cos®,cos9,.
U= (-P1Py/4ne,R>)(2cos8 ,cos0.y ~ sind,sin8,(cos($, - $)))«
i exp(-U/kT)
Bb) P(81,8946),6)) = —————— 2 om
7
f sind 140! fsindsdey, J 4g) f ds,exp(-u/kt)
° ° ° °
c) Let $ = ¢) - ¢)
1 Qu
Jind 49, Jatne, 485 ! dé U exp(-U/kT)
p =
w Ta
fsino 40, Jsind,d0, [ do exp(-U/kT)
° ° °
Taking U/kT << 1 then exp (-U/kT) * 1 - U/kT. Let us call
Denominator = D and Numerator = N of the above integral
. 7 20 3
D= tear fsind de, J do(1 + pypy/(4ne RKT) x
(2cos8,cos8y - eanestneacet))
Qn
Note that [ cos¢do = 0, and Joind cos0d0 = 0, then D = 81
" oon 2n 0
N= fsind 49, if sind,d8, ! ds[U-U7/kT]. Note that we have
° °
already calculated the fire term in this integral, so we need
to calculate the second term. Remember that U* is proportional
57to 4cos20,cos*8, - 4cos8,cos#,sin8,sin@,cos¢ +
sin?@,sin26,cos”. The middle term goes to zero since
fTcos sae = 0, but [™ cos*@sined® = 2/3, and fsin2ed0 = 4/3,
therefore N = 4+2/3+2n + 4/3+4/3+n = 485/9 = 10n/3 so
pe p? pe
_ 16n/3 Pe fie
ee cam, 226. oy nee ep
lone RRT = 24 eR KT
0 °
5.5 a) Using Eq. 5-17, we write py = 2a Eo/[1-(2a/4negR3)] «
Solving for 2a = Po! Laney? Py + EQ).
b) The molecular polarizability for all R is given by Eq. 5-17
at = 2a/[1-(2a/4negR?) ]-
5.6 Consider the figure where it shows the atoms and the
oe >
external field. From Eq. 5.14 § = a(E, + E") where E' is an
A B
——
additional electric field produced by each atom at the site of
4
for 9 = 90° we get E'(8 = 90°) =
the other. Using Eq. 2+
-p2/AnegR?. Thus &
= -aE + BY) /4ne, 3 Solving for £', and
hence for p and & gives
58* +
fr = - —_,_—_"______ , p=
4ne R[L + a/dne R>] 1 + a/4neR
° 0 3.
at = la/(l + a/ane,k?)
5.7 Using Eqs. 5.9 and 5.21 we get a +a’ ta =
L :
ane + Po(coth n=) where n= pyE/kT. At high temperatures
ue have a = 4neyx? + po2/3KT. ‘The average dipole moment is
= Chnegr? + ng2/3kT) EG
5.8 using Hq. 5.36 we write K - 1 = 3/6(i-T,). Since K is very
large, then K - 1 = K, and we have I/k = 4(1 «+ The points
ra)
(i/", T): (0035, 133) and (.0105, 153) define the Line and we
met 6 = 1.05 x 1079 and Ty = 123°K.
5.9 the ferroelectric condition is given by Eq. 5.34:
Na/3e, = 1. sut from Hq. 5.5 we have
ha/3e, “= 1V/(K + 2) = 290090 which is different from 1;
thus it is not ferroelectric. ror tne liquid we have
ha/3e, = .367 which is not ferroelectric. In this case
5.10 In Eq. 4.4 we nave o, = - (K-1) of/K. For ferroelectric
naterials x >> 1, hence o, = o,-
59CHAPTER 6
6.1 See example 6.1.
6.2 (4) = Aexp(-ar
2 2
a) From Poisson equation we have V#= -p/eor p= ~e Vo.
Note @ has only radial dependence, thus
Peed ee. tor Gb ertoun san
shaceep(-at?> (net). Thus the density is
0 = 3Ae,ar(4-3ar)exp(-ar’).
b) v= (1/2) fotdv = (3/2) Ae a frexp(-2ar?)(4-3ar>) dv
= 4m (3/2) A2ega fexp(-2ar9)(4-3or yr Jar.
We know fx"lexp(-x) dx = T(n) with T(n) the gamma function
0)
T(nt1) = n(n), and T(n+l) = n!, if n is a positive integer.
3 ang put the integral in the form of the f integral
Ane) A’
we get V * ras 1(4/3).
oe ae caer ae
6.3 We use Eq. 6.9: 0) = us P,Q, where the coefficients of
ik
potential are: Pj; = eae Pjq = Waneyt for j # k since the
spheres are identical and they are at the same distance from
each other. Grounding sphere | changes the charge on it but
does not affect the charge on the other two. Thus 6) = 0 =
60(qy/a) + (q/2) + (g/t). Solving for q) gives q) = -2qa/2. Now,
ungrounding sphere 1 and grounding 2 gives %) = 0 = (q)/&) +
2a
e
(ag/a) + (q/2) which gives qy = - = (1 - 24). Now ungrounding
2, and grounding 3 gives % = 0 = (a\/2) + a2/%) + (q3/a)- This
pol ves| for qo 2 se i =) te
i e ie
to allowing the charge on each sphere to adjust to the external
“until equilibrium” refers
potential. You do not have to keep grounding spheres.
6-4 Since the force between the spheres is given by F = -8U/ar
where U is the electrostatic energy of the system, and r is the
distance between then, then we need to calculate U as a function
of r. We use the coefficients of potential to solve the
problem. We use Eqs. 6.8 and 6 olde
PLO P 4,
; HOt Pyy!ys P50, + Poot, where
My
te oe Dil 202
j
Py, = P92» and Py) = Pyy because of the symmetry. Now when
© = qy, % = 0, and ) = V. These equations give V = Py,9q and
y = Po14, = Pyyqy vinich can be solved for Py, = V/qy+ ‘the
2
energy of the systen is 1 = 1/0] Py,, and hence from the force
ap,
rs laze
Yo aq (=) we get dP),/dr = 2F/ay-+ In the second
F
operation, we have *) = Pyyqy + Pyydy, and V = Pyyqy + Pyyay
which gives Py) = Y(qy - a2)/aq, and hence Uy =p Py (ay + 49) +
Pyya,4y+ From the force between the spneres we can calculate
aP,o/dr as follows:
ap
1
+ 03) a
614Py/dr = -F(2qf + 43)/ 439+
Now that we know Pj), Pyy and their derivatives with r we can
calculate the properties of the last operation. In the last
operation sphere 1 is grounded, thus %, = U, while the charge on
the second does not change, therefore 0) = Pi Q + Pid, oF
Qy = -a9€ay - 47)/q,+ The force between them can be calculated
from F = ~U3/r where Uz = My Py (QT + 43) + Py2Qazy OF
F = qp(2ay - 43) Flaps
6.5 a) See Ex. 6.8. By symmetry Py, = P97 = P33 = Pays Now
=i
4ne a
°°
% = Pij4+ thus Py] = I/énegae By Symmetry: Py2 = Pa] =
Pyy = Pay = Pay = Pay = W4negt = Pog = P32+ Also by symmetry:
Pig = P31 = Pan = Pou = MW4ne,(v28).
c) (1) When we connect 1 and 2 then the charge merely splits
because of symmetry, that is q) = 4/2, q) = @/2+ (2) When we
connect 1 and 3, then because spheres 1 and 3 are symmetric with
respect to sphere 2, the charge splits equally q, = 4/4, and
a3 = 4/4. (3) Similarly there is symmetry in the last
operation and we get q, = 4/8 and q, = q/8. Thus the final
charges are q/8, 4/2, q/2, q/8-
6.6 a) See example 6.6.
b) From above we have C = 4ne,ab/(b - a). Now if b - a=
d << a,b then C* 4ne,a’/d = Ae,/d with & = 4nd the surface
area of the sphere.
626.7 Let us call the fields in the regions r < a, a< 4 + + .
b p9, we have a field as if all the current is at
uot
p20: B= 5
+
b) By symmetry we have A= A,2 and A, = A(p) only. Thus we
>
oJ
solve V7A,
: V2 oo :
i) For p< oy: Noa, = pap % ap “ce 7 Oe ches
A, C £n p+ D. Now A, must be finite at p = 0; therefore
C= 0. Moreover we take the potential to be zero so we set
>
D= 0. Thus A= 0 for p < p)-
95ii) For 9; p 9: A, = Cy inp +d B-- 5; but from
Ampere's law we have
ulz
o
Qn
2n p+ De
.
The constant D, can be determined by matching A at p = p9-
2
ior tp ut
ot | Mo 1 o
ee
eo eee
96.
uniformly magnetized sphere we have m = f Mdv =
+ + os
= uo(M +H). But f =m x B, thus
Sut + ty =f aaty th x + x
> +
ray Ma sina where a is the angle between M and H.
9-9 Let us label the three regions r < Ry, Ry = (Rez)12, where I, = Vj2 is the
amplitude of the current through the source. The source voltage
ts taken of the form V= Vet. Thus = RIZ/(1 + u?R?C?).
e) We need to calculate the current in the resistor. First let
143us calculate the impedance of the circuit and the current in the
inductor and hence in the source for wLC = 1.
S 2p2 = vetut a
Zo = (R+ t0l)/(1 + w?R2C?), T= voet/z. Vp =v - vy
V - iwLT = v1 - iwt/z,)
~iVRwC. Thus Ip = Vp/R = -iVuC.
Thus the average power =1/)R IR Ip = "ye ve wrc?.
| i . 1
13.25 2) = tl - 2 2) = goede
b) From example 13.7, cosh(ty) = 1 + 2/22 =
1 - (wLc - 1)%/2u2Lc. Then we have the following frequency
range: 0 < (wLC - 1)7/wLe < 4.
144CHAPTER 14
14.1 a) Gauss' law between the plates where 2 is normal to
q. * qu
* - 2 s + aD 0. -
he 7 pd a .
them gives D = oz = 7~ sin(ut)z, thus J) = 3¢ A c08(wt)Z
u
The displacement current 1) = J,-A = qyucos(ut)-
The conduction current Ig = dq/dt = qq weos( wt) « Observe that
the two are equal.
b) Since Ip is in the 2 direction then B= Bb
e)
=
ap
1
vx(B, 9) = (uq,w/A) cos (wt)2 > (oB,) (aq, 0/A)eos (ut)2
qu
pe, cy) 0
= 24+ tt). Buc B#@ = hi = 0;
By wy > x cos(wt) we 8 at p= 0, then C = 0;
9,8
Ta Pcos(ut) d+
therefore B =
14.2 We are given V = 250 sin(377t) volts. (See Ex. 6-5)
tee dQ/dt where Q= CV. We know the capacitance of the
coaxial cylinders from Eq. 6-36 is
gw me. 2 (25m) (6-7_x 8.85 107!28/m)
Tn(b/ay ~ 1n(67 +5) :
145LLL LLL LLL LLL LL LO NTC CC CCC TTT Rte:
Thus C= 1.0 * 107 Therefore I, = dQ/dt = 9.6 x 1075
:
cos(¥77t). In = f Seeda =e f Baad. But inside the capacitor
>
= O/(2nepL). Since Q/L = (2meV/2n(b/a), then =
# V/pln(b/a) )6- Now dE/dt = 250/(pln(b/a))377 cos(377t). Note:
250(377) cos 377t 2
eln(b/a)
cos(377t). We see that I) = Tg-
fda = 2nLp then t) = P= (9.6 x 107)
14.3, We are given V = 150 sin(500t) volts. (See Ex. 6.6.)
I, = dg/dt where Q = 4neV ba/(b ~ a), therefore I, = 7.1 x 10>
z)eos(500t). Now the displacement current is
if r
ec (fb us. we b= 2 then SE ee
at 2 at 2
: fuer) 4ner
so Ip = I,/(4nr') fda = 1,/(4mr*)4mr* = 1, so Ip = 1
Tp
z. z
14.4 From Eq. 14.17 we have v°(') - ue % (5) = 0
B ae B
a) E = (Ac) exp(x-ct), B = A exp(x-ct). We use a single
notation.
E 3
ee 2 :
x (8) = (1) A expGr-ety, 5 (B) = (£2) A exp(x-ct)
ax” ae
Using ve = 1/c* we get
2 E 32 E c c
v°(B) - pe =s (8) = (1) 4 exp(x-ct)-(1) A exp(x-ct) = 0.
ot
146b) E = (Ac)&n(x + ct), B= -A gn(x + ct)
E Ac 1 a2 E e
(8) = (-4) sarge and Sy (8) = aC)
ax
ze
|
en
1
eo) x Oy
ot (x +
(x + ct)
14.5 a) We are given E = Ey sin(wt - kz). For free space we
have €= €g- Sob = ch = cf = cy sin(ut - kz). We
Ore = oro
integrate Maxwell's Equation
ab *
vxi=- 2 where vx 8 = = = KE eos(ue - kz):
t
B-fouxdae = - KES sim(ut - ke) = - X sin(wt - kz)
2 E
t- Be. - 2 sin(wt - kz) = -e E cx sin(ut - kz).
we oo
Mo
b) We are given ii = H fexp(-i(ut + ke))+ Integrating Maxwell's
*
Equation V x ie Je + a where V x i = iky Hyexp (-iCut + kz))
and 7 O we get:
t
Be fixate = - « ¥ Hexp(-iCut + kz))
B ks
sof "eee? Hyexp(-i(ut + kz).
147