Vector Calculus Solution Manual
Vector Calculus Solution Manual
o 0
RddzdR 27r loa 2RJb
2
- R2 dR
47r[ - (b
2
- R
2
)3/2
47r(b
3
- (b
2
- a
2
)3/2) /3.
The proportion of the apple removed is this volume divided by 47rb
3
/3,
which is 1 - (1 - a
2
/b
2
)3/2.
6.6 The limits on the coordinates are 7r /3 () 27r /3, < 27r, r a,
where a is the radius of the Earth. The volume V is
V = r
27r
/3 r
27r
r r2 sin () dr d d() = 27r [r3 /3] [- cos = 27ra
3
/3,
17r/3 10 10
so exactly half of the Earth's volume is less than 30 away from the Equator.
6.7 In spherical polar coordinates the divergence of e", is zero, using (6.23).
The curl, from (6.25), is
1 a sin () 1 ar cot () 1
V x e", = -.---er - --eB = --er - -eB.
r sm () a() r ar r r
6.8 The formula u . Vu = V(luI
2
/2) - u x (V x u) is used to find u . Vu.
Since u is a unit vector, its magnitude is constant so V(luI
2
/2) = O. Using
(6.17), V xu = ez/R, so u x (V x u) = (0,1,0) x (0,0, I/R) = (I/R,O,O).
Therefore u . Vu = -eR/ R.
174 Vector Calculus
6.9 Recall that the definition of the Laplacian of a vector field is '\7
2
v = VV .
v - V x V x v. Using the formulae (6.14) and (6.15), the R component of
VV v is
which can be expanded to give
1 8VR VR 2 8
2
vR 1 8vI/> 1 8
2
vI/> 8
2
vz
R 8R - Ii + 8R2 - R2 8 + R 88R + 8R8z'
Computing the R component of V x V x v using (6.17) gives
1 8 ( 8VI/ 1 8
2
vR 8
2
vR 8
2
vz
R28 vI/> + R 8R - R2 82 - 8z2 + 8R8z'
Subtracting these two quantities gives the R component of '\7
2
v:
2 1 8VR VR 8
2
vR 2 8vI/> 1 8
2
vR 8
2
vR
('\7 V)R = R 8R - R2 + 8R2 - R2 8 + R2 82 + 8z2 .
Note that this cannot be the Laplacian of the R component of v, since it
involves vI/>' In fact ('\7
2
V)R and '\7
2
(VR) are related by
2 2 VR 2 8vI/>
('\7 V)R = '\7 (VR) - R2 - R2 8 .
Solutions to Exercises for Chapter 7
7.1 The definition Lij = ~ ej states that Lij is the cosine of the angle between
the ~ and ej unit vectors. Referring to Figure 7.1, the cosine of the angle
between ~ and el is cos 0, and the cosine of the angle between e; and e2 is
the same. The cosine of the angle between ~ and e2 is cos(1r/2-0) = sinO,
and the cosine of the angle between e; and el is cos(1r /2 + 0) = - sin O.
These results agree with the matrix in (7.3).
7.2 u is a vector, so ui = Lijuj. To show V u is a scalar we need to compute
its value in the dashed frame:
( )
' 8u: 8 8uj 8Xk
V'u =-8' =-8 ,(LijUj) = L ij -
8
-8"
Xi Xi Xk Xi
using the chain rule. Now using (7.8),
( )
' 8Uj 8Uj 8Uj
Vu =Lij -
8
Lik = Ujk-8 = -8 = V u.
Xk Xk Xj
Solutions 175
Thus V . u has the same value in the dashed and undashed frame, so V . u
is a scalar.
7.3 Since a and b are vectors, a ~ = L;kak and bj = Ljmbm. Note that the
suffices are carefully chosen to avoid repetition. Then
This agrees with the transformation rule (7.13), so a;b
j
is a second-rank
tensor.
7.4 Note that Tij can be written as Tij = aibj where a is the position vector
(XI,X2) and b = (x2,-xd which was shown to be a vector in Example
7.3. Therefore Tij is a tensor by the result of the previous question. Al-
ternatively, the fact that Tij is a tensor can be confirmed by following the
method of Example 7.3.
7.5 Since is a scalar, / = . The transformation rule for Tjk is
, EP'
Tjk = a x j a x ~
so Tjk is a second-rank tensor.
7.6 If Tij is a tensor, Tij obeys the rule
Setting j = i in this formula,
So the value of Tii is the same in the dashed and undashed frames, i.e. Tii
is a scalar.
7.7 In suffix notation, the divergence theorem (5.1) becomes
iii ; dV = ffs Ujnj dS.
This result holds if Uj is replaced by T Ij , T2j or T
3j
, giving the required
result (7.16).
7.8 If Qijkl is a tensor of rank four,
7.9 Q;jkl obeys the above rule, so by setting k = j,
so Qijjl obeys the transformation law for a tensor of rank two.
176
Vector Calculus
7.10 Given that Uiai is a scalar,
so - UjLij = O. If this holds for any ai, then = Lijuj, so Ui is a
vector.
7.11 B
rs
is an anti-symmetric tensor, so
= . L. .. B = -L. .. B =
Hence B is also anti-syr
7.12 If Brs is anti-symmetric, Ors = -Osr so Or,' = -Orr and hence Bn = O.
7.13 A ijk has the properties A ijk = .1# and = -Aikj . Repeatedly apply-
ing these rules alternately gives
So any element of A ijk is equal to minus itself, hence all elements are zero.
7.14 (a) Aij = fijkBk, so
so Aij obeys the transformation rule for a second-rank tensor. Since fijk is
anti-symmetric with respect to any two indices, Aij is antisymmetric.
(b) Multiply through by fijm:
making use of Exercise 4.4(c). Hence Bm = fijmAij/2.
7.15 The most general isotropic fourth-rank tensor is aijkl = Mijl5kl + J.U5ik6j1 +
v6i/6j k, from Theorem 7.4. From (4.12), the difference between products of
6ij terms can be written in terms of a products of fijk terms. For example, if
A = 1, J.L = -1, v = 0, we have the isotropic tensor aijkl = bijbkl - 8ik6jl =
filmfmjk
7.16 Since bij and fijk are isotropic tensors, the combination aijklm = bijfklm
is an isotropic fifth-rank tensor. Since any two of the five suffices can be
chosen for the 6, for example aijklm = 6ik fjlm or aijklm = 6i/fjkm, the total
number of different tensors of this type is the number of ways of choosing
two objects from five, which is 5!/2!(5 - 2)! = 10. So there are at least ten
different components in the most general isotropic fifth-rank tensor.
Solutions
177
7.17 Following the approach of Section 7.4.2, the kinetic energy of a volume
element dV is P Iv1
2
dV/2 and v = n x r, so
E = 1/2 IIi P Ivl2 dV
= 1/2 IIi PijkfljrkilnJ1Ir m dV
= 1/2 IIi p(fljfljrkrA' - fljrjflkrk) dV
= 1/2 III P (6jkr2 - rjrk) fljflk dV
1jk flj flk /2.
Solutions to Exercises for Chapter 8
8.1 The size of the body L is a length and the diffusivity k has units of
length
2
/time. The only combination of these which has the units of time is
L
2
/k. Therefore the time for heat to diffuse through a body is proportional
to the square of the size of the body.
(a) Assume that the mammoth is the same shape as the chicken, made of
the same material (so the wooliness is ignored) and 20 times the length of
the chicken. The defrosting time is therefore 20
2
x 6 hours which is 100
days.
(b) Assume that to cook properly, a certain temperature must be reached
in the interior. The time required for this is proportional to L2. The mass
M is proportional to L
3
, so L ex: M
1
/
3
Therefore the cooking time should
be proportional to the two-thirds power of the mass: t ex: M2/3.
8.2 Using (8.4),
8p
8t
8E
= oV-
8t
oV (V x B -/1-0j)/EO/1-0 from (8.7)
-Vj
since the combination div curl is always zero.
8.3 Gauss's law says that the total flux of electric field through the surface is
the total charge within the surface divided by EO. Taking the surface to
be the surface of a sphere of radius r, with area 41Tr2, Gauss's law gives
41Tr2 Er = Q / EO, which agrees with the result of Example 8.4.
178
Vector Calculus
8.4 Take the curl of (8.12):
{)V x E
V x (V x B) == /-Lofo {)t .
Expanding the r.h.s. and using (8.10) and (8.11) gives
2 {)2 B
- V' B == -/-Lofo {)t
2
'
so B obeys exactly the same wave equation as E.
8.5 Given the electric field E == Eof(kx -wt) == Eof(u) where u == kx -wt,
the magnetic field B can be found using (8.11):
VxE V x (Eof(u))
Integrating with respect to t and changing the sign to find B gives
1 1
B == -k x Eof(k . x - wt) == -k x E.
w w
Hence the magnetic field B is perpendicular to the electric field E.
8.6 The energy can be written w == B . B /2 + E . E /2c
2
, so the rate of change
of energy is
{)w
at =
B . {)B + -.!:..E . {)E
{)t c
2
{)t
= -B V x E + E . V x B using (8.11) and (8.12)
== -V (E x B).
This gives the conservation law
{)w
at+V.P==O
where P = Ex B is known as the Poynting vector, representing the energy
flux of the electromagnetic wave.
8.7 (a) If the material is in equilibrium then the stress tensor must be symmet-
ric (8.16), so a == 0 and c == b. Also, {)Pij/{)Xj == O. For the first row (i = 1)
this is satisfied. The second row (i = 2) gives m x ~ l l + bX2 = 0, so m == 2
and b = -2. From the third row, CX3 + x ~ l = 0, so n = 2 and c == -2.
(b) For the surface Xl == 1 the normal is (1,0,0) (the normal points in
the direction of the side which is exerting the force). Thus the force is
Fi == Pijnj dB == Pi1 dB. Since P ll = P 21 = 0 and P 31 = X2, the force is
only in the X3 direction. Its magnitude is
F3 == 11 11 X2 dX2 dX3 == 1/2.
Solutions
8.8 The Lh.s. of (8.22) is
-!f"k [V x v]
2 tJ k
=
1 aVm
-"2fijkfklm aXI
1 aVm
-"2 (6il
6
jm - oimojd aXI
! (aVi _ aVj) = Sij.
2 aXj aXi
8.9 Applying the condition aPijjaXj = 0 and using (8.21),
179
8.10 (a) From the definition of the strain tensor, Ell = aX2, E12 = E2l =
axI/2 + bXl, E22 = -2bx2 and the other components of Eij are zero.
(b) Using (8.21), Pu = A(a - 2b)X2 + 2/LaX2, P12 = P2l = 2/L(axI/2 + bXl),
P22 = A(a - 2b)X2 - 4/Lbx2, P33 = A(a - 2b)X2 and other components of Pij
are zero.
(c) Applying the equilibrium condition aPijjaXj = 0, this is identically
satisfied for i = 1 and i = 3, but for i = 2, /La + 2/Lb + A(a - 2b) - 4/Lb = O.
This is satisfied if a = 2b.
8.11 (a) The term V (pu) in (5.9) can be written u V(po + PI) + (Po + pdV u.
Now since Po is a constant and terms involving products of u and PI can be
neglected, this simplifies to Po V u, so (5.9) becomes apI/at + Po V u = O.
(b) In the Navier-Stokes equation (8.28) the product term u . Vu can be
ignored and there is no body force b or viscosity /L, so the remaining terms
are (Po + pd (au j at) = - V (Po + pd The term involving PI and u can be
ignored, Po is constant and PI = apI, so poau j at = -a V Pl
(c) The velocity u can be eliminated from (8.35) and (8.36) by taking the
time derivative of (8.35):
2
at
2
= - Po V . at = a V . V PI = a V' PI'
Hence the density perturbation obeys the wave equation (8.13). Physically,
the waves are sound waves travelling through the fluid.
Alternating tensor 70
Ampere's law 93, 135
Angular momentum 13
Bernoulli's equation 149
Circulation 27
Closed curve 26
Combinations of grad, div and curl 76
Conductivity tensor 126
Conservation of mass 85
Conservative 28, 51, 61, 93
Coordinate curves 99
Coordinate surfaces 99
Coriolis force 13
Cross product 9, 70
Curl 58, 74, 106
Curvilinear coordinates 99
Cylindrical polar coordinates 107
Del squared 56, 76
- of a vector, 77
dij 68
Determinant 71
Diffusion equation 133, 148
Directional derivative 50
Divergence 53, 74, 105
- theorem, 83
Dot product 4, 66
Dummy suffix 66
Electric field 134
Electromagnetic waves 137
Electromagnetism 134
Electrostatics 135
ijk 70
Faraday's law 135
Field 17
Fluid 85
Fluid mechanics 145
Flux 31
Free suffix 65
Gauss's law 135
Gauss's theorem 83
Gradient 48, 74, 104
Green's identities 89
Green's theorem 96
Green, George 89
Heat equation
Heat transfer
Hooke's law
133
132
143
Incompressible 86
Inertia tensor 128
Integral
- line, 25
- surface, 31
- volume, 39
Integration 21
- by parts, 23
- by substitution, 22
Irrotational 60
Isosurface 48
Isotropic 123, 127
Jacobian 102
Kronecker delta 68
Lame's constants 143
Index
Laplace's equation 56, 87, 133, 135
Laplacian 56, 76, 105
- of a vector, 77
Level surface 48
Line element 25
Line integral 25
Line, equation of 12
Lorentz force 13
182
Magnetic field 134
Magnitude 4
Mass conservation
Maxwell's equations
85
134
Newtonian fluid 146
Normal vector 33
Operator 50
Order 119
Orthogonal 101
- matrix, 116
Partiai derivative
Plane, equation of
Poisson's equation
Potential 52
Potential theory
Pressure 145
45
7
135
135
Quotient rule 120
Rank 119
Right-handed 3
Rotation 11, 61
- of coordinates, 115
Scalar 1, 117
- field, 17
Scalar product 4
Scalar triple product 14, 71
Scale factor 100
Solenoidal 56
Vector Calculus
Spherical polar coordinates 110
Stokes's theorem 91
Strain 143
Stress tensor 140
Suffix notation 65
Summation convention 66
Surface element 32
Surface integral 31
Symmetric tensor 122
Taylor series 47
Tensor 68, 119
- isotropic, 123
- symmetric, 122
Transpose 73
Vector 1,117
- addition, 2
- components, 3
- field, 17
- magnitude, 4
- normal, 33, 50
Vector product 9
Vector triple product 16
Viscosity 146
Volume element 39
Volume integral 39
Vorticity 61
Vorticity equation 148
Wave equation 137
Zero vector 3