#09 (09/18/2024)

Vector Analysis

Scalar product (inner product, dot product)

The scalar product between the two vectors, a and b, is defined as

a·b = a b cos θ,
(1)
where θ is the angle between the two vectors.
innerpro.jpg
The following properties can be derived based on the definition of Eq.(1).
  1. a ·b = b ·a
  2. c a ·b = a ·c b = c (a ·b)
  3. a ·(b + c ) = a ·b + a ·c
It is possible to derive the scalar product in terms of the components of a vector by introducing the base vectors, ex and ey. Note that

ex·ex
=
ey·ey
=
1,
ex·ey
=
ey·ex
=
0.
(2)
Using Eq.(2), the scalar product of two vectors in components is
a·b
=
(ax ex +ay ey )·(bx ex +by ey )
=
ax bx ex ·ex + ax by ex ·ey + ay bx ey ·ex + ay by ey ·ey
=
ax bx + ay by.
(3)
The definition of Eq.(1) is a preferred form as it is independent of the coordinate system.

Vector product (cross product)

The vector product is defined between two vectors in 3-D as

a×bS n,
(4)
where S is the area of the parallelogram spanned by the two vectors and n is the unit vector whose direction is in the right hand system as shown in the figure.
crosspro.jpg
The following properties can be derived based on the definition of Eq.(4).
  1. a ×b = −b ×a
  2. c a ×b = a ×c b = c (a ×b)
  3. a ×(b + c ) = a ×b + a ×c
Note that a×a = 0 (no area by two identical vectors).
The relationship among 3-D base vectors is summarized as

ex×ey
=
ey×ex
=
ez
ey×ez
=
ez×ey
=
ex
ez×ex
=
ex×ez
=
ey,
(5)
and

ex ×ex = ey ×ey = ez ×ez = 0.
(6)
Using the relationship above, the vector product is expressed in components as
a×b
=
(ax ex +ay ey +az ez )×(bx ex +by ey +bz ez )
=
ax bx ex×ex+ax by ex×ey+ax bz ex×ez+ …
=
=




ex
ey
ez
ax
ay
az
bx
by
bz




.
(7)

Equation of plane (application of dot product)

eqplane.jpg
The relationship among x0 (fixed point), x (moving point) and n (normal, fixed vector) is

n ·(xx0) = 0.
(8)
So the equation of a plane is expressed as
n ·x
=
n ·x0
=
c.
(9)
Note that c represents the shortest distance between the plane and the origin.

Triple products

There are two types of vector products, scalar triple product and vector triple product, whose main usage is found in vector analysis.

Scalar triple products, a·(b×c)

a·(b×c) = a·S n = S a cosθ = S h represents the volume of the parallelepiped spanned by the three vectors, a, b and c. Note the cyclic relationship below.
a·(b×c) = b·(c×a) = c·(a×b)
triplepro.jpg

Vector triple products, a ×(b×c)


a ×(b×c) = (a·c)b − (a·b)c.
Note that the right hand side is a combination of b and c both of which are inside the parentheses in the left hand side. The one that is in the middle (b) comes first.
Exercise
  1. a×( b×c) +b×( c×a) +c×( a×b) = ?
  2. (a×b)·(c×d) = ?

Curves, Surfaces and Volumes

If the position vector, R, is expressed by one parameter, u, a set of all the points of R varying u forms a curved line. If R is expressed by two parameters, u and v, as R(u,v), a set of all the points of R by varying both u and v forms a curves surface. (obvious statement for R(u,v,w)).

Arc length

linelength.jpg
The length of a curved line from u = u to u = u + du can be computed by carrying out integration of ds (line segment)

dR
=
R(u + d u ) − R(u)
=
R(u) + d R

du
du + …− R(u)
 ∼ 
dR

d u
du
=
eu d u,
(10)
where
eu dR

d u
.
Therefore,

ds
=


 

dR ·dR
 
=


 

eu d u ·eu d u
 
=


 

eu ·eu
 
d u.
(11)
and

s =
u=u1

u=u0 


 

eu ·eu
 
d u.
(12)
Example Compute the length of y = x2 from x = 0 to x = 1.
lengthexample.jpg
Solution As
R =



x
y
z




=



u
u2
0




,
it follows

eu =



1
2u
0




,
and

ds =

 

eu ·eu
 
d u =

 

1 + 4 u2
 
du,
so

s =
1

0 


 

1 + 4 u2
 
du
Let the computer carry out this integral.
Enter the following
Sqrt[1 + 4 u^2] u 0 1

exactly (case sensitive, square bracket) in the boxes then press the query button.
Integrate with respect to
from to .

Surface area element

surfacearea.jpg
To obtain the surface area expressed by two parameters, u and v, one can first compute the surface element spanned by the two vectors shown above:

R(u + d u, v)− R(u, v)
=
R(u, v) + R

u
duR(u, v)
=
eu du,
(13)
R(u, v + dv) − R(u, v)
=
R(u, v) + R

v
dvR(u, v)
=
ev dv,
(14)
where
eu
R

u
,
(15)
ev
R

v
.
(16)
Remembering that the area of a parallelogram spanned by two vectors, a and b, is

ds
=
  ⎛
 ⎥


det


eu ·eu ,
eu ·ev
ev ·eu ,
ev ·ev



 
du dv
=


 

E FG2
 
du dv.
(17)
where
E
=
eu ·eu
F
=
ev ·ev
G
=
eu ·ev.
Alternatively, recall that the magnitude of a ×b is the area of the parallegram spanned by a and b, Equation (17) can be also written as

d s = | eu ×ev | du dv.

Example of surface integral

As an example of surface integrals, we attempt to compute the surface area of a sphere x2 + y2 + z2 = a2. By setting x = u and y = v, a position vector on the sphere is expressed as

R(u, v) =




u
v


 

a2u2v2
 





,
(18)
so
eu
=






1
0
u



 
[ˉ(a2u2v2)]






(19)
ev
=






0
1
v



 
[ˉ(a2u2v2)]






(20)
and

eu ·eu
=
1+ u2

a2u2v2
ev ·ev
=
1+ v2

a2u2v2
eu ·ev
=
uv

a2u2v2
.
(21)
Hence,

dS
=


 

(eu ·eu) (ev ·ev) − (eu ·ev)2
 
du dv
=
a




a2u2v2
du dv.
(22)
The surface area of a semi-sphere (radius = a) is

S
=


a du dv




a2u2v2
=
a

r dr dθ




a2r2
1
=
a
a

0 
r dr




a2r2



0 
1 dθ
=
a
a

0 
r dr




a2r2
2
=
a2.
(23)
Alternative approach
spherical.jpg
A bit faster approach than the one above is to use the spherical coordinate system. Noting that the equation, x2 + y2 + z2 = a2 is automatically satisfied if

x
=
a sin ϕcos θ,
y
=
a sin ϕsin θ,
z
=
a cos ϕ.
(24)
Hence we can set
R =



a sin ϕcos θ
a sin ϕsin θ
a cos ϕ,




.
(25)
It follows
eϕ
=




a cos ϕcos θ
a cos ϕsin θ
a sin ϕ




(26)
eθ
=




a sin ϕsin θ
a sin ϕcos θ
0




,
(27)
and
E
=
a2
(28)
F
=
a2 sin 2 ϕ
(29)
G
=
0
(30)
dS
=


 

EFG2
 
dθdϕ
(31)
=
a2 sin ϕdθdϕ,
(32)
thus
S
=

π/2

0 



0 
a2 sin ϕdθdϕ
=
a2
π/2

0 
sin ϕdϕ
=
a2.
(33)

z = f(x, y)

If the equation of a curved surface is given as
z = f(x, y),
(34)
the area element can be expressed as
d S =

 

1 + fx2 + fy2
 
dx dy,
(35)
where
fx
f

x
,
fy
f

y
.
(36)
Proof

ex
=
R

x
=




1
0
fx




ey
=
R

y
=




0
1
fy




,
E
=
ex ·ex
=
1 + fx2,
F
=
ey ·ey
=
1 + fy2,
G
=
ex ·ey
=
fx fy,
EFG2
=
(1 + fx2)(1 + fy2) − (fx fy)2
=
1 + fx2 + fy2.
(37)

Volume

volumeelement.jpg
The volume element of the parallelepiped spanned by the three vectors above can be expressed as3

dv
=
  ⎛
 ⎥
 ⎥






eu ·eu, eu ·ev, eu ·ew
ev ·eu, ev ·ev, ev ·ew
ew ·eu, ew ·ev, ew ·ew




 
du dv dw
=
|eu·(ev×ew)| du dv dw.
(38)
Example: Volume of sphere
spherical2.jpg

R =



x
y
z




=



r sin ϕcos θ
r sin ϕsin θ
r cos ϕ




,
(39)
so

er =



sinϕ cos θ
sinϕ sin θ
cos ϕ




,   eϕ =



r cosϕ cos θ
r cosϕ sin θ
r sin ϕ




,   eθ =



r sin ϕ sin θ
r sin ϕ cos θ
0




(40)
Noting that
eϕ ×eθ =



r2  sin2 ϕ cos θ
r2  sin2 ϕ sin θ
r2  sin ϕ cos ϕ




(41)
it follows
er ·eϕ ×eθ = r2 sin ϕ.
(42)
so
V
=



dV
=



0 

π

0 

a

0 
r2 sin ϕdr dϕdθ
=



0 
dθ
π

0 
sin ϕdϕ
a

0 
r2 dr
=
a3

3
.
(43)

Summary

  1. Length of Arc

    dl
    =


     

    eu ·eu
     
    du
    =
    |eu | du.
  2. Area

    d s
    =
      ⎛
     ⎥


    |
    eu ·eu,
    eu ·ev
    eu ·ev,
    ev ·ev
    |
     
    du dv
    =
    |eu ×ev| du dv.
  3. Volume

    d v
    =
      ⎛
     ⎥
     ⎥






    eu ·eu,
    eu ·ev,
    eu ·ew
    ev ·eu,
    ev ·ev,
    ev ·ew
    ew ·eu,
    ew ·ev,
    ew ·ew




     
    du dv dw
    =
    |eu ·(ev ×ew) | du dv dw.


Footnotes:

1
u = r cos θ,     v = r sin θ,     du dv = r dr dθ.
(44)
2 Let
a2r2t,
(45)
then − 2 r dr = dt and as r → 0, ta2 and as ra, t → 0 so that

a

0 
r dr




a2r2
=
0

a2 
r

t
dt

(− 2 r)
= 1

2

a2

0 
dt

t
= a.
(46)
3 Let a=(a1, a2, a3), b=(b1, b2, b3), c=(c1, c2, c3).
a·(b×c) = (−a3 b2 + a2 b3) c1 + (a3 b1a1 b3) c2 + (−a2 b1 + a1 b2) c3 =



a1
a2
a3
b1
b2
b3
c 1
c2
c3




.
On the other hand,




(a, a)
(a, b)
(a, c)
(b, a)
(b, b)
(b, b)
(b, a)
(b, b)
(b, b)




=(−a3 b2 c1 + a2 b3 c1 + a3 b1 c2a1 b3 c2a2 b1 c3 + a1 b2 c3)2 =



a1
a2
a3
b1
b2
b3
c 1
c2
c3




2





 
.



File translated from TEX by TTH, version 4.03.
On 17 Sep 2024, 13:41.