05.04.2013 Views

UNIT 11 VECTOR DIFFERENTIAL CALCULUS - IGNOU

UNIT 11 VECTOR DIFFERENTIAL CALCULUS - IGNOU

UNIT 11 VECTOR DIFFERENTIAL CALCULUS - IGNOU

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Vector Calculu~<br />

<strong>UNIT</strong> <strong>11</strong> <strong>VECTOR</strong> <strong>DIFFERENTIAL</strong><br />

<strong>CALCULUS</strong><br />

Structure<br />

<strong>11</strong>.1 Introductio~i<br />

Objectives<br />

<strong>11</strong>.2 Scalar and Vector Fields<br />

<strong>11</strong>.3 Vector Calculus<br />

<strong>11</strong>.3.1 Limit and Continuity<br />

<strong>11</strong>.3.2 Differentiability<br />

<strong>11</strong>.3.3 Applications of Derivatives<br />

<strong>11</strong>.4 Directional Derivative and Gradient Operator<br />

<strong>11</strong>.5 Divergence of a Vector Field<br />

<strong>11</strong>.5.1 Physical Interpretation<br />

<strong>11</strong>.5.2 Formulae on Divergence of Vector Functions<br />

<strong>11</strong>.6 Curl of a Vector Field<br />

<strong>11</strong>.6.1 Physical Interpretation<br />

<strong>11</strong>.6.2 Rotation of a Rigid Body<br />

<strong>11</strong>.6.3 Formulae on Curl, Divergence and Gradient<br />

<strong>11</strong>.7 Summary<br />

<strong>11</strong>.1 INTRODUCTION<br />

In physical problems we often come across quantities such as temperature of a liquid,<br />

distance between two points, density of a gas, velocity a~id acceleratio~i of a particle or<br />

a body, tangent to a curve and normal to a surface. In Unit 10, you have learnt that<br />

physical quantities can be categorised either as a scalar or as a vector. You might have<br />

noticed that some physical quantities, whether scalars or vectors, are variable. That is,<br />

their values are not constant or static but change with the change in variable. For<br />

example, the density of a gas, which is a scalar quantity, changes from place to place<br />

and is different at different times. Similarly, tangent to a curve may have different<br />

directions at different points of the curve. This variable character of scalars and vectors<br />

give rise to scalar functions and vector functions. Further, you may notice that at<br />

different positions the temperature or velocity of a body does not remain same. The<br />

distribution of temperature or velocity is therefore defined at each point of a given<br />

domain in space which leads to the idea of scalar fields aiid vector fields. We shall<br />

discuss about scalar functions and scalar fields, vector functions and vector fields in<br />

Section <strong>11</strong>.2.<br />

In Section <strong>11</strong>.3, we shall extend, in a very simple and natural way, the basic concepts<br />

of differential calculus to vector-valued functions. We shall also discuss about<br />

physically and geometrically important concepts related to scalar and vector fields<br />

namely, directional derivatives, gradient, divergence and curl in Section <strong>11</strong>.4 and give<br />

their applications.<br />

We shall introduce the vector operator V in Sectioii <strong>11</strong>.5 and give the physical<br />

interpretation of the divergence of a vector field and some basic formulas involving it.<br />

Finally, the concept of curl of a vector field and its invariance is discussed in Section<br />

<strong>11</strong>.6. Formulas involving curl, divergence, curl and gradient, divergence and curl and<br />

Laplacian operator, v2, are also developed here.<br />

Objectives<br />

After studying this unit you should be able to<br />

* define a scalar function, a scalar field, a vector function and a vector field,<br />

* state coilditions for the coilvergence of a sequence of vectors,<br />

* define limit, continuity and differentiability of vector functions,<br />

* differentiate sum, difference and products iilvolvi~lg vectors,


* describe the notion of directional derivative and compute directional Vector DLlltmatial Calculus<br />

den va tives,<br />

* define and compute gradient of a scalar field and divergence and curl of vector<br />

fields,<br />

* interpret physically the gradient, divergence and curl of a vector,<br />

* define conditions for solenoidal and irrotational vector fields, and<br />

* solve problems on application of del operator and product rules involving the<br />

del operator.<br />

<strong>11</strong>.2 SCALAR AND <strong>VECTOR</strong> FIELDS<br />

Let us first talk about the scalar field.<br />

Scalar Fields : Consider the distance of a point P from a fixed point Po which will be a<br />

real number. As we vary the point P its distance from Po also changes. It depends only<br />

on the location of point P in space and may be regarded as function f (P). If we<br />

consider cartesian coordinate system in space and take coordinates of fixed point Po as<br />

(x,, yo, zo) and the variable point P as (x, y, z), then f is given by the well-known<br />

f(P) = f(x.y,z) = \/(x - x")~ + (y -<br />

+(z -1")'.<br />

Next, consider a room fitted with an air-conditioner (A.C.). Once A.C. is switched on<br />

for its cooling effect; the temperature of room falls down. Now, if we put off the A.C.,<br />

the temperature starts rising up till it reaches the room temperature. The temperature<br />

further rises up if we now switch on the A.C. for its heating effect.<br />

Thus, the temperature of the room depends on the switching system of the air<br />

conditioner.<br />

In this case, the temperature of the room can be considered as a function of switching<br />

system of the A.C.<br />

In both the examples taken above, you may note that distance, as well as, temperature<br />

give us only the magnitude and not the direction. Hence both are scalar quantities.<br />

These quantities depend either on the position of P or on the switching system of the<br />

A.C. In both the situations, we get a function. Functions of this type are called Scalar<br />

Functions. It may also be noted here that distance or temperature functions do not<br />

depend on the choice of coordinate system or brand of A.C., but only on the physical<br />

situations such as actual distance or actual duration of switching on the A.C.<br />

Formally, we give the following definition of a scalar.<br />

Definition : A Scalar function is a function which is defined at eacli polnt of certai~z<br />

region (domain) in space and whose values are real numbers depending only on the<br />

points in space but not on particular choice of the coordinate system.<br />

In most of the applications, the domaill of a scalar function is a curve, a surface, or a<br />

three-dimensional region in space.<br />

In the case of temperature of a room fitted with A.C., the do~nain of temperature<br />

function is the set of points on the regulator of A.C., which controls the temperature of<br />

the room.<br />

The function f associates with each poiut in domain D a scalar (a real number) and we<br />

say that a scalar field is obtained in D. More formally, we hgve the following<br />

definition :<br />

Definition : If i) be a function which associates a unique scalar with eacli point in a<br />

given region, then $ is called a scalar field function, or simply a scalar field.<br />

In a plane, for instance, the equation of a curve is given by f (x, y) = constant, i.e.,<br />

these are curves along which f has a constant value for all points in the x y - plane.<br />

Similarly, a surface may be given by i) (x, y, z) = constant, i.e., these are surfaces for<br />

which i) has a constant value for all points in space. In these examples scalar functions .<br />

f and i), having plane and space as their respective domains, are scalar'fields.<br />

Some more examples of scalar fields are the density of the air of the earth's atmosphere<br />

and the pressure within a region through which a compressible fluid is flowing.


Vector Calculus<br />

We can represent a scalar field by a formula as well as pictorially. For a pictorial<br />

representation, we may use the curves and surfaces. The pictorial representation of<br />

scalar fields are maps showing physical geography of a region (indicating hills, lakes,<br />

land above or below sea level, etc.).<br />

In the same manner when we assign a vector to each point of a certaill region we obtain<br />

a vector field. Let us now talk about vector fields.<br />

Vector Fields<br />

Consider a curve in a plane or in space. At each point of the curve we can draw a<br />

Figure <strong>11</strong>.1 : Tangent Vectors of a curve Figure <strong>11</strong>.2 : Normal Vectorsofa surface<br />

A<br />

We can assign to each point P of the curve, a tangent vector t (P) (see Figure. <strong>11</strong> .I).<br />

Similarly at each point of a surface, we can draw a nornlal (see Figure <strong>11</strong>. 2). These<br />

llormals inay have different directions at different points of the surface. Thus to each<br />

A<br />

point Q of the surface, a nonnal vector n ( Q ) lnay be assigned.<br />

Moreover, a vector may depend on one or more than one indepeildent scalar variables.<br />

The velocity of a particle, for instance, depelids on the position of the particle as well<br />

as time. The same is true for position vector and acceleration of a particle.<br />

We now give the following definition :<br />

Definition : If to each point P of n certain region G in spnce a vector V (P) is<br />

assigrted, then V (P) is called n vector function.<br />

There are many examples of vector functions ill physics. For example velocity,<br />

acceleration, and forcc are all vector functions. The gravitational force exerted by the<br />

sun on a unit mass is also a variable vecto~, depending on the position of the mass, and<br />

thus represents a vector function.<br />

The collectio~l of all such vector fuilctio~ls V (P) is called a vector field on G. More<br />

precisely, we give the following definition :<br />

Definition : IfF be a function which assigns to each point x in its domain, a vector<br />

F (x) is called a vector field function or a vector field.<br />

A physical exaillple of vector field is given by the particles of a fluid under flow.<br />

At any instant the velocity vector V (P) of a rotatilig body collstitutes a vector field,<br />

called the velocity field of rotation. If we introduce a cartesian coordinate system<br />

having the origin on the axis of rotation, then<br />

A L A<br />

V(x,y,z) = LI) x (xi + y j + zk),<br />

where r, y, z ate the coordinates of any point P of the body in a plane perpendicular to<br />

the axis of rotaltion and o is the rotation vector or constant angular vclocity of the body<br />

(see Figure <strong>11</strong>.3).<br />

Next consider a particle A of inass M, which is placed at a fixed point Po and let a<br />

particle B of mass m be free to take up various positions P in space (see Figure <strong>11</strong>.4).


' axis<br />

- 1<br />

of rotation<br />

Figure <strong>11</strong>.3. Velocity Field of a body rotating with consmat angular wlocity O in<br />

the positive (counter clockwise) di~ction<br />

I b<br />

Then particle A attracts particle 8. According to Newton's law of attractionlgravitation,<br />

the corresponding gravitational force p is directed from P to Po and its magnitude is<br />

1<br />

proportional to -, where r is the distance between P and Po. Then<br />

r2<br />

GMm<br />

IpI=--p<br />

where G is the gravitatio~ial constant. Hence p defi~ies a vector field in space.<br />

0<br />

Vector Differential Cak<br />

Figure <strong>11</strong>.4 Some of the vecton of gravitational field<br />

By now you must have clearly understood what we mean by scalar functions, scalar<br />

fields, vectors functions and vector fields. You can test your knowledge by attempting<br />

the following exercise :<br />

E 1.<br />

Which of the followillg are scalar functioils, scalar fields, vector fu~~ctiorls and vector<br />

fields ?<br />

(a)<br />

(b)<br />

(c)<br />

(d)<br />

(e)<br />

The gravitational force on a particle of mass Mat distance r due to<br />

another particle of mass m.<br />

A collstant force applied to a particle.<br />

The temperature at every point of a mass of heated liquid.<br />

Potential of an electric charge placed at the origin.<br />

The force on a unit charge placed at a point P due to an electric charge<br />

e placed at the origin.<br />

You have already learnt about the basic concepts of limit, continuity, differentiability<br />

and partial differentiation of scalar fuiictio~ls in Units 2 and 5 of Block-I. We shall<br />

now, in the next section, illtroduce these basic concepls of calculus for vector functiolls<br />

in a simple and natural way.<br />

<strong>11</strong>.3 <strong>VECTOR</strong> <strong>CALCULUS</strong><br />

We begin by introducing the idea of limit and continuity of a vector function. We<br />

know that the limit concept is a basic and useful tool for the analysis of functions. It<br />

enables us to study derivatives, improper integrals and other important features of<br />

W ball inlrodu[c Ihc mlcepls of limit and co~~tinuity<br />

~unctions, e s<br />

manner without going into the intricaciof delta, epsilon method.<br />

* A 4 1 !-.-!A -...I<br />

c#.-4:...-:4 --<br />

~II<br />

a most informal


Vector Calcu<strong>11</strong>~<br />

the length and direction of I. more precisely, we define \;mite of v-c*o* raluca<br />

functions in te- of the familiar limits of real-valued functions in the following way.<br />

A A A<br />

Definition : Let f ( t ) - fi ( t ) i + f2 ( t ) j + f3 ( t ) k be a vector valued function oft<br />

defined in some reighbourhood of a (possibly except at a). The limit off ( t ) as t<br />

approaches the number a is the vector 1 #the limit of I f ( t ) - 1 1 as t approaches a is<br />

zero. In symbols<br />

l==lim f(t)elim If(t)-<strong>11</strong> -0<br />

I-a r-a<br />

Observe that f ( t ) having las a limit means that the components off have the<br />

corresponding components of 1 as limits. In other words, if<br />

then<br />

f )<br />

= f )<br />

- A b<br />

A<br />

i + f2(t) j + f3(t).2 and 1 ll i + 123 + l3 i<br />

lim f(t) = 2- lirn fi = <strong>11</strong>, lim f2 - I2 and lirn J; = I3<br />

This equivalence says that we may calculate limits of vector- valued functions<br />

componentwise, i.e., one component at a time.<br />

The limit of a vector valued functions can also be defined in the usual 6, E lnani~er as<br />

we do for the real-valued functions in the following way :<br />

A vector function f (t) of a real variable t is said to tend vector 1 as t approaches a, if<br />

to any pre-assignedpositive number E, however small, there corresponds a povitve<br />

number 6 such that<br />

I f (t)-I(c.rwhen)t-a1


In view of the equivalence ia equation (1 1.2) above, we may say that f ( t ) is<br />

continuous at t = a if each component off is continuous at t = a. Thus we may test a<br />

vector function for continuity by applying our knowledge of real-valued functions to<br />

each component of$<br />

Also f (t) is continuous function if it is continuous at every point of its domain.<br />

Just as in the case of real valued functions the sum, the difference, scalar product and<br />

vector product of two continuous vector functions are also continuous. We shall not be<br />

proving these results here. You can check them yourself,<br />

Let us consider the following example :<br />

Example 2 :<br />

Discuss the continuity of the function<br />

Solution :<br />

The function f (t) is continuous at every value oft > 0 because each colnponellt is<br />

continuous fort > 0. However, f is discontinuous fort s 0 because -, 1 the first<br />

di<br />

component off, is not defined fort s 0.<br />

You may now try the following exercises :<br />

E2<br />

Find lim f ( t ) if<br />

1-0<br />

b)<br />

A A A<br />

f(t) = (efsint)i + (efcost)j - e'k<br />

t A 1-cost 1<br />

C) At) = - i + J + i<br />

sin t t<br />

At what values oft are the following vector functions f (t) continuous.<br />

A A A<br />

a) f(t) - (cost)i+ (sint)j + k<br />

In mechanics, if the position vector of a particle is given by r (t) and if we wish to find<br />

its velocity, we will have to differentiate r(t) with respect to time. Similarly, if the<br />

potential of an electric charge is given, then to determine force due to a,n electric<br />

charge, we shall have to take a recourse to differentiation. We now discuss the<br />

~'tC-,,t:,+:~~ nf 9 vector function.<br />

Vector Diiltcrentid Calculus


0<br />

* -<br />

Figure <strong>11</strong>.5 : Derivative of a vector<br />

hnction<br />

<strong>11</strong>.3.2 Differentiability<br />

We define the derivative of a vector-valued function f (t) at a point t = a by the same<br />

type of limit equation as we use for scalar functions. Thus<br />

j'(a) = h--0 ~ i m f ( ~ h + ~ ) - f ( ~ ) ,<br />

provided the limit on the right exists. We then expect that f is differentiable at t = a iff<br />

each of its components is differentiable at t = a. In this connection we prove the<br />

following fqsult :<br />

Theorenl 1 :<br />

A A A<br />

A vector function f ( t ) = f, ( t ) i +fi ( t ) j +f3 ( t ) k is differentiable at t = a iff<br />

each of its component function is differentiable at t = a. If this condition is met, then<br />

Proof :<br />

Consider the difference quotient<br />

f(a+h)-f(a)<br />

h<br />

..<br />

!(a) =i,(a):+ha)j +h(a)P<br />

The left hand side of equation (<strong>11</strong>.6) has a limit as h + 0 iff each component on the<br />

right hand side has a limit as h -, 0. The ith component on the right has a limit iff fi is<br />

differentiable at a.<br />

Finally, if each component is differentiable at a, then taking the limit h + 0 in equation<br />

(<strong>11</strong>.6) of each of the quotient, we get equation (<strong>11</strong>.5), thus proving the result.<br />

I<br />

Note that the differential coefficient f ( a ) is itself a vector and is called the derivative<br />

of f (t) at t = a. From the definition it is clear that every derivable vector fu~lctio~~ is<br />

continuous. Consider the followi~lg example :<br />

Example :<br />

Obtain the derivative off( r ) = ( sid t ) ; + (In t ); + tan- ' (3t ) ;<br />

Solution :<br />

The function and its derivative are defined at every positive value of r.<br />

\<br />

We now give the geometrical interpretatipn of the derivative at a vector valued<br />

function.<br />

.Geometrical Representation of Derivative<br />

Draw the vector f (t) forvalues of the indepellJellt variable t in some interval<br />

co~ltai~ling t and t + A t from the same initial point 0. The<strong>11</strong> the locus of head of<br />

nti~ig f for different values oft traces out a space curve (See<br />

and OQ = f (t + A t), then<br />

f(t+At)-f(t) = OQ-OP = PQ = AS,(say).<br />

a lim<br />

dt A,-oAt


The direction of uis the limiting directioq of of A/: But as Q tends top, PQ<br />

d t At<br />

tends to the tangent line at P. Hence the dimtion of gis along the tangelit to the<br />

d t<br />

space curve traced out by P. Lets denote the length of the arc of this curve from a<br />

fixed point on it upto P. Then the magnitude of gis given by<br />

d t<br />

since the ratio IL?fl, chOrdpQ -1 a, At-0,<br />

As arc PQ<br />

Thus derivative of a vectorfunction represetzts a vector whose direction is tangent to<br />

ds<br />

tlre space curve traced by tlre vectorfurtction arzd the magrzitude is -, whm s is the<br />

d t<br />

arc length from a Fred point on the curve to the variable point representing tlre<br />

vector function.<br />

You may note here that a vector will change if either its magnitude changes or<br />

direction changes or both direction and magnitude changes. Jn this regard. the<br />

following results may be remembered :<br />

df<br />

(a) The necessary and sufficient condition for f (t) to be constant is - = 0.<br />

d t<br />

(b) The necessary and sufficient condition for f (t) to have constant magnitude is<br />

(c)<br />

1 he necessary and sufficie~it co~iditioll forf (t) to have constailt/uni€orm<br />

directioii is<br />

The familiar rules of differentiation of real functions yield corresponding rules for<br />

differentiating vector functions; for example,<br />

(i) (cf)' = cf (Caconstant).<br />

I I I<br />

(ii) (u v) = u * v<br />

du<br />

(iii) ( u f ) = - f + u (u is a scalar function of t).<br />

dl dl<br />

I I I<br />

(iv) (u.v) = u .v + u.v<br />

I I I<br />

(v) (uxv) = u XV+UXV<br />

In (v) above the order of the vectors must be carefully observed, as cross multiplication<br />

of vectors is not commutative.<br />

The chain rule of differentiation is also valid for vector valued functioils. That is, if<br />

f (t) is a differentiable functioii of t, aild t = g (s) is a differentiable function of s, then<br />

the composite function f (g(s)) is a differentiable function of s and<br />

We call write equation (1 1.7) in the form<br />

Vector DUlcreathl Calculus


Vector Calculus The chain rule given for vector functions by equation (1 1.8) is an immediate<br />

consequence of the chain rule for scalar functions that applies to the components f,, f,<br />

and f3.<br />

Consider the following example :<br />

Example 4 :<br />

~x~ress~intermsofsif~ (t)=; + sin(t+l)J? + et+ 'i and<br />

ds<br />

t - g (s) - 2-1.<br />

Solution :<br />

From the chain rule we have<br />

We would obtain the same result if we first substitute t = g ( s ) - 2 - 1 in the<br />

formula for f (t) and then differentiate w.r.ts.<br />

And now a few exercises for you.<br />

E4<br />

Find the derivative of the vector function f (t) in each case and give the<br />

domain of the derivative :<br />

(a) f (t) = eZAt + A t e- A j<br />

(b) f(t) = i + 3j-k<br />

(c) f(t) -;sin- '2r +J? taa13t + k -<br />

As we have already mentioned, not all vector functions are functions of one variable.<br />

The velocity of a fluid particle in motion is a function of time and position. The<br />

position vector of a fluid particle at any time and at any position is a vector function of<br />

four variables x, y, z and t. Thus, in any physical problem involving a vector function<br />

of two or more scalar variables, we may be required to find the partial derivatives of


this vector function. In other words, we may be required to find the derivative of the Vtar ~ i ~ ~ ~ C.I~I,,~ t i . 1<br />

vector function w.r.t. one scalar variable treating the other scalar variables as constant.<br />

Partial derivatives can be calculated for vector functions by applying the rules we<br />

already know for differentiating vector fuactions of a siagkscalar variable.<br />

If a vector function f (u, v) be a differentiable function of iwo scalar variables U, v<br />

given' in the component form as<br />

then partial derivatives off w.r.t. u and v are denoted by<br />

respectively and are defined as<br />

?<br />

%,<br />

af 1 af2 1 ah A<br />

-=- r+-J+-k<br />

au au a u a u<br />

and a f afl t af2 1 ah A - - --r+-J+-k<br />

av av a v a v<br />

Similarly, aZ J a2 fl A a2 f2 1 aZh A<br />

- = -- i + ~ j + ~ k<br />

a u2 a u2 a u<br />

etc. are the second order partial derivatives.<br />

a2f a2fl: a2j+ aZf3 A<br />

-=- r+- j+auav<br />

auav auav auav k<br />

The derivatives a2f a2j :<br />

auav avau<br />

are called mixedpartial derivatives and they are equal if each of them is a continuous<br />

function.<br />

Physically, a gives the rate of change 0ffw.r.t. to u at a given point (u,v) in space.<br />

a u<br />

Thus, partial derivatives<br />

of a function f (x, y, z, t) give us the rate of change off in the directions of x, y, z axes at a<br />

given instant and ugives the rate of change off with respect to time at a given point in<br />

at<br />

space.<br />

Consider the following example :<br />

Example 5 :<br />

A A<br />

4<br />

Find the first order partial derivatives of r ( t, , t2 ) = a cos tl i + a sin tl J + t2 k<br />

Solution :<br />

A A<br />

a r<br />

We have - = - a sin tl i + a cos'tl j<br />

a tl<br />

Note that r (tl, t2) is a position vector. It represents a cylinder of revolution of radius a,<br />

having the z- axis as axis of rotation.


Vector Calculus<br />

F i <strong>11</strong>.6 Pammebic<br />

Repmentation<br />

of a curve<br />

X Y<br />

Figore <strong>11</strong>.7 : Parametric<br />

Representation of Straight lie<br />

Figure <strong>11</strong>.8 : Representatiota<br />

of the tangent to a m e<br />

You may now try the following exercise :<br />

E L For each of the vector functionf, find the first partial derivatives w.r.t. x,y,z.<br />

A A<br />

(a) f=xyi+yzj<br />

i<br />

A A<br />

(b) f = #i- e-"j<br />

(c) f - ~ ~ ; + ~ 2 3 + ~ i<br />

You know that curves occur in many considerations iq calculus as well as in physics,<br />

for example, as paths of moving particles. Let us consider some basic facts about<br />

kurves in space as an important application of vector calcuhs about which we are<br />

going to talk in our next section.<br />

1 1.3.3 Applications of Derivatives<br />

The simplest application of vector calculus is some basic facts about curves in space.<br />

Given a cartesian coordinate system, we map represent a curve C by a vector futiction~<br />

A A A<br />

r(t) = x(t)i + y(t)j+ z(t)k.<br />

Here to each value of the real variable t, there corresponds a point of C having<br />

position vector r (to). (see Figure <strong>11</strong>.6).<br />

For example, any straight line L can be represented in the fonn<br />

r(t) = a + tb,<br />

whete a and b are contant vectors and line L passes through the point A with position<br />

vector r = a and has the direction of b (see Figure <strong>11</strong>.7).<br />

The vector function<br />

A A<br />

r(t) = acosti + bsintj<br />

represents an ellipse in the xy-plane with centre at origin and axes in the directions of<br />

x and y axes.<br />

Further, if a curve C is represented by a continuously differentiable vector functioi~<br />

r (t), where t is any parameter, then the vector<br />

drr - lim r(t+At)-r(t)<br />

dl At--0 At<br />

has the direction of the tangent to the curve C at P (t) (see Figure <strong>11</strong> 3).<br />

Thus, the position vector of a point on the tangent is the sum of the gosition vector r<br />

of a point. P on the curve and a vector in the direction of the tangent. Hence the<br />

parametric representation of the tangent is<br />

.


dr<br />

where both r and -depend on P and the parameter o is a real variable.<br />

dt<br />

Let us now consider the following example :<br />

Example 6 :<br />

If a (t) be a variable unit vector, show that<br />

do<br />

(i) - is a vector normal to a.<br />

d t<br />

Solution :<br />

(i)<br />

da<br />

(ii) - is a unit normal vector to a, 0 being the angle through which a turns.<br />

d0<br />

Since a (t) is a unit vector,<br />

Differentiating Eqn. (1 1.9) w.r.t tot, we get<br />

da da<br />

So a and - are at right angles, i.e., - is a vector normal to a.<br />

dt dt<br />

(ii) Let OP = a and OQ = a + A a be two neighbourirrg values of the given<br />

vector making an angle A0 with each other (Figure <strong>11</strong>.9).<br />

Then<br />

and<br />

da A a<br />

- = lim -<br />

d0 ~e-oA8<br />

Since Aa is nonnal to a in the limiting position when A0 -. 0, therefore<br />

do<br />

- is normal to a<br />

d0<br />

Also<br />

da<br />

Hence - is a unit vector normal to a<br />

d0<br />

Let us now look into some of the applications of derivatives to dynamics.<br />

Let the position vector of a point moving on a curve be given by 4t). Its<br />

displacement in time At is<br />

Since the velocity V of the moving point is the rate of change of its<br />

displacement w.r.t to time, therefore<br />

Again, the acceleration of the point, being the rate of change of velocity, is<br />

given by<br />

0 a P<br />

Figure lL9


Vector Calculus<br />

f' *<br />

72<br />

1~<br />

Consider the following examples :<br />

Example 7:<br />

Show that if r = a sin ot + b cos ot, where a and b are constants, then<br />

Solution :<br />

We haver=asino t + bcosot<br />

Differentiating w.r.t. tot, we get<br />

d2 r 2 dr<br />

- = - o r and r x - = - a) (a x b).<br />

d t2 dt<br />

and dZr 2 2<br />

-<br />

dt2<br />

= - a w sin o r - bw cos ot<br />

Also dr<br />

r x - = (asinot + bcoswt) x (aocosot-bosinot)<br />

dt<br />

2<br />

= -(wsin2wt + ocos wt)a x b<br />

= - w (a x b).<br />

(-: a x a = 0 and b x b = 0).<br />

Let us take up another example.<br />

Example 8 :<br />

A particle P moves on a disk towards the edge, the position vector being<br />

r(t) = tb,<br />

where b is a unit vector, rotating together with the disk with constant angular velocity<br />

o in the counter clockwise direction. Find the acceleration a of P.<br />

Solution :<br />

Since the particle is rotating with constant angular velocity a), therefore b is of the<br />

form<br />

The position vector of particle P is<br />

A<br />

J .r(t) = tb<br />

-<br />

Differentiating equatiou (1 1.1 1) w.r.to t, we get<br />

X<br />

V=i=b+tb<br />

Obviously b is the velocity of P relative to the disk and t b is the additional velociry<br />

Fiun <strong>11</strong>.10 : Motion in due fo the rotation (see Figure <strong>11</strong>.70)).<br />

Exnmple 8<br />

Differentiatingequation (<strong>11</strong>.12) once more w.r.t to t, we obtain


..<br />

a X v = 2 b + t b (<strong>11</strong>.13) Vector Diiemntial Calculus<br />

. .<br />

In the last term of equation (1 1.13), using (<strong>11</strong>.10), we have b = - w2 b. Hence the<br />

accelerationt b is directed towards the centre of the disk and is called the Centripetal<br />

Acceleration due to the rotation.<br />

The most interesting term in equation (1 1.13) is 2 b, which results from the<br />

interaction of the rotation of the disk and the motion of P on disk. It has the direction<br />

of b, i.e., it is tangential to the edge of the disk and it points in the direction of<br />

You may now try the following exercises :<br />

A particle moves along the curve x - e ' , y = 2 cos 3t, z = 2 sin 3t, where t<br />

is the time variable. Determine its velocity and acceleration at t = 0.<br />

A particle moves so that its position vector is given by<br />

A A<br />

r = cos ot i + sin otj.<br />

Show that the velocity V of the particle is perpendicular to r and r x V is a<br />

constant vector.<br />

Find the Coriolis acceleration when the particle moves on a disk towards the<br />

edge with position vector<br />

where b is a unit vector, rotating together with the disk with the constant<br />

angular speed o in the anti- clockwise sense.


fl -<br />

If we cobider a scalar field f (x,y,z) in space, then we know that<br />

?f?.f?f<br />

ax9ay'az<br />

are the rates of change off in the directions of x, y and z coordinate axes. It seem<br />

unnatural to restrict our attention to these three directions and you may ask the ~iatural<br />

question- How to find the rate of change off in any direction ? The answer to this<br />

question leads to the notion of directional derivative which we shall try to aliswer i<strong>11</strong> thk<br />

next sectian.<br />

<strong>11</strong>.4 DIRECTIONAL DERIVATIVE AND GRADIENT<br />

OPERATOR<br />

Let us consider a scalar field in space given by the scalar function f (P) = f (x, y, z),<br />

where we have chosen the point P in space. Let us choose any direction at P, say given<br />

by vector b. Let C be a ray from P in the direction of b and let Q be a point on C,<br />

whose distance from P is s (see Figure <strong>11</strong>.<strong>11</strong>).<br />

The limit<br />

f(Q) -f(P) 7<br />

S<br />

as .+o<br />

af =lim (Q+P)<br />

if it exists'is called the directional derivative of ttie scalar function f at P in ttie<br />

J b<br />

direction of b.<br />

P In this way, there can be infinitely many directio~ial derivatives off at P, each<br />

corresponding to a certain direction. An iilterestiilg question arises - Can we represent<br />

r p <strong>11</strong>-<strong>11</strong> : Ilircctionrl any such directio~ial derivative in terms of some derivative or derivatives off at P ?<br />

tkrivalive<br />

The answer is 'yes' and it is achieved as follows :<br />

Let a cartesian coordiiiate system be given. Let a be the positio~i vector of P relative to<br />

the origin of this system. Then any point on ray C, or ray C itself, can be represented in<br />

form.<br />

A A A<br />

r(s) = x (s)i +y(s)j + z(s)k = a + s b (s r 0)<br />

Now a is the derivative off [x (s), y (s), z (s)] with respect to arc-length s of my C.<br />

as<br />

Hence assuming that f has continuous first partial derivatives and applying the chain<br />

rule, we obtain<br />

~ = af dx + af dy + af dz,<br />

as axds ayds azds<br />

dx d~ dz<br />

where - and - are evaluated at s = 0.<br />

ds ' d s ds<br />

Also from equation (1 1.15), we have<br />

dr dxA dyA dzA<br />

- = -i+-j+-k=b.<br />

ds ds ds ds<br />

This suggests that we introduce the vector<br />

aft afA afA<br />

gradfn-r+-j+-k<br />

ax ay az<br />

and write equation (<strong>11</strong>.16) in the form of a scalar product<br />

- - iagradf<br />

as


The vector grad f is called the gradient of the scalar function f. More precisely, we<br />

give the following definition :.<br />

Definition :<br />

The vector finction<br />

is called the gradient of the scalar finction f and is written as grad$ viz,<br />

A2f+;2f+ii!f<br />

grad f = i ax ay az'<br />

Here we have assumed that scalar function f is a continuously differentiable function.<br />

Thus the directional derivative of the scalar point function f along the direction of<br />

vector b can be written as<br />

--<br />

In other words, the directional derivative i!f is the resolved pan of grad f in the<br />

as<br />

direction of b.<br />

We see that the gradient of a scalar field f is obtained by operating onfby the vector<br />

operator<br />

This operator is denoted by the symbol V (read as "del" or nabla) and it operates<br />

distributively.<br />

In terms of V, we write<br />

The operator V is also known as differential operator or gradient operator and ( V f (<br />

gives the greatest rate of change sf f.<br />

The interest and usefulness of introducing this symbol V lies in the fact that it can be<br />

formally assumed to have the character of a vector and as such it facilitates the<br />

manipulations with expressions involving differential operators. Thus formally, V f<br />

being product of a vector V by a scalar f is a vector.<br />

Before we take up the properties of gradient of a scalar field and operator V, we take<br />

up a few examples to illustrate how directional derivatives are calculated.<br />

Example 9':<br />

Solution :<br />

Find the directional derivative 2fof (x, y, z) = 2$ + 3y2 + 2 at 'the point P<br />

as<br />

A ?<br />

(2, 1,3) in the direction of the vector a = i - 2 k.<br />

Vector DiereaIl.l Cdculos


Vector Calculus<br />

:. Unit vector in the direction of a is<br />

Therefore,<br />

The minus sign indicates that f decreases in the direction under consideration.<br />

Let us take another example.<br />

Example 10 :<br />

Find the directional derivative off (x, y, z) = 2 y2 2 at the point (1, 1, -1) in the<br />

di~ction of the tangent to the cuwe<br />

~=k,~=2sint + 1,z-t-cost,-1 r t r 1.<br />

Solution :<br />

A A A<br />

:.AtP(1,1,-l),(gmdj) = 2i + 2j - 2 k<br />

Now any point on the given curve has position vector<br />

:. Direction of tangent t to the given curve is<br />

The point (1, 1, -1) on the curve corresponds to the value t = 0.<br />

:. Required direction derivative


You may now try the following exercises :<br />

E 10<br />

E <strong>11</strong><br />

E 12<br />

Find the directional derivative o f2<br />

@ @ A<br />

direction of 2 1 - 25 + k.<br />

Find the direction in which the directional derivative of<br />

f(x,y) = (2 - Y2)/(~y)at(1,1)iszero.<br />

Find the directional derivative of 4 x z3 - 3 2 y2 2 at (2, -1, 2) alpng the<br />

z-axis.<br />

We shall now discuss some of the important properties of gradient of a scalar field<br />

functiov.<br />

Consider a differentiable scalar function f (x, y, z) in space. For each constant c the<br />

equation<br />

f(x,y,z) = c - constant<br />

represents a surface S in space. Then, by letting c assume all values, we obtain a family<br />

of surfaces; which are called level surfaces of the function f. Since, by the definition of<br />

a function, 6ur function f has a unique value at each point in space, it follows that<br />

through each.point in space there passes one, and only one, level surface off.<br />

If $ (x, y, z) denotes the potential, the surface<br />

is called an equipotential surface. The potential of all points on this surface is equal to<br />

the constant C.<br />

Important geometrical characterisation of the gradient of a scalar funktion f is in tenris<br />

of a vector normal to a level surface or an equipotential surface. This property can also<br />

be used in obtaining the normal to a given surface at a given point. We shall now take<br />

up this property.<br />

Vector Dillercotid Cdcnlus


Gradient as Normal Vector to Surfaces<br />

Let P be a point on the level surface<br />

{see Figure <strong>11</strong>.12).<br />

where 6f is the difference in values off at Q and P.<br />

Hence if Q lies on the same level surface as P. Vf6r = 0.<br />

(by differential calculus) (<strong>11</strong>.17)<br />

n ?<br />

Moreover, let Vf = ( V f ( n, where n is a unit vector normal to the surface.<br />

I<br />

A<br />

Let Q be a point on a neighbouring level surface f + 6f and let 6 n be the<br />

perpendicular distance along PN between the two surfaces. Then the rate of change<br />

1<br />

off nonnal to the surface 3<br />

= lim V f. =, using (1 1.17)<br />

6n-0 6n<br />

af<br />

Hence the magnitude of Vf is equal to -. Thus the gradient of a scalar field f is a<br />

an<br />

vector normal to the surface f = constant and having a magnitude equal to the rate of<br />

change off along this normal.<br />

Let us take up an example for the beV.er understanding of what we have discussed<br />

above.<br />

1


~xam~le <strong>11</strong> :<br />

Find a unit vector normal to surface<br />

2 y + 2 x z - 4 at the point (2, -2,3).<br />

SsWtion :<br />

Letf - xZY + 2xz 5 4.<br />

Then a vector normal to the surface is<br />

Hence a unit vector normal to the surface<br />

1 A A A l? 2A 2A<br />

= - ( -2i+4j+4k) = --i+-j+-k.<br />

6 3 3 3<br />

somi of the vector fields occuring in physics and engineering are given by vector<br />

functions which can be obtained as the gradients of suitable scalar funcltions. Such a<br />

scalar function is then called a potential function or potential of the corresponding<br />

vector field. The use of potentials simplifies the investigation of those vector fields<br />

considerably. To understand this let us consider the gradient of a potential due to an<br />

electric charge.<br />

P I1 Gradient of a potential due to an electric charge<br />

Consider an electric charge e placed at the origin.<br />

Figure <strong>11</strong>.13<br />

The potential of this charge at a point P (x, y, z) is 4 where OP = ; (Figure <strong>11</strong>.13).<br />

r<br />

The force on a unit charge &aced at P will be<br />

in the direction OP.<br />

The component of the force in the direction of x will be<br />

Similarly, the component of the force inj and z directions can be obtained as<br />

ez<br />

--<br />

ey and - - ? respectively.<br />

?


Vector Calculus The resultant force may then be written as<br />

This result also follows from equation (1 1.1 8).<br />

e<br />

If we deilote the potential - by I$, then the force on the particle is<br />

r<br />

A a$. ; * + 1 * + k -, i.e., grad I$.<br />

ax ay. a z<br />

The physical definition of potential shows that this will he true for all cases.<br />

Similarly, the gradient of the potential due to a iluillber of charges placed at various<br />

points will give force due to these charges.<br />

In the fluid flow, the gradient of the velocity potential will give the velocity at a point.<br />

Also the gradient of gravitational potential will give the force due to the gravity.<br />

In Section <strong>11</strong>.3, we discussed the differelltiability of a vector function. Now we take<br />

the derivative along a curve, which is expressed in tenlls of gradient of the function<br />

representing the curve at a point and tangeilt to curve at that point.<br />

P <strong>11</strong>1 Derivative along a curve<br />

A<br />

For a curve f (x, y, z) = coilstant in a space, the derivative 1 . V f along a curve is zero<br />

A<br />

if and only if the function f is constalit along the curve, t beiug the unit tangent to the<br />

curve.<br />

P IV The length and direction of grad f are independent of the particular choice of<br />

cartesian coordinates<br />

We shall give the proofs of P-<strong>11</strong>1 and P-IV in Appendices I and I1 respectively.<br />

You may now try the following exercises :<br />

E 13<br />

Find grad J", where r is the distance of any point from the origin.<br />

E 14<br />

E 15<br />

Find the directiol~al derivative of the hhction x y2 + y 2 + z x' along the<br />

3<br />

tallgel<strong>11</strong> to thecurvex = t,y = t2 ,z = t at (1, 1, 1)and at (-1, 1, -1).<br />

Find a unit vector normal to the surface x3 +y3 + 3 x y z = 3 at the point<br />

(132, -1).


E 16 vector DIffercntial calculus<br />

Find the gravitational potential when the force of attraction between two<br />

particles, distance r apart, is proportional to -(r/r3)<br />

In Section <strong>11</strong>.4, we have developed a tool for determining the rate of change of a scalar<br />

field in space and time in the form of directional derivative and gradient. A natural<br />

question arises - How fast a vector field varies in a given region ? Another question,<br />

which may be asked is - Can we extend the analysis of gradient to a vector field ?<br />

The answer to the second question is that it is not possible to extend the analysis of<br />

gradient to a vector field. For a vector field, we can find two types of derivatives -<br />

(i) derivatives involving rate of change of a vector component in its own<br />

direction, called the divergence.<br />

(ii) rate of change of components in directions other than their own, called the<br />

curl. Let us first study the concept of divergence of a vector field.<br />

<strong>11</strong>.5 DIVERGENCE OF A <strong>VECTOR</strong> FIELD<br />

We have seen that starting with a scalar point fu~lctionf, we can construct a vector<br />

point function grad f. But, what happens if we sta with a vector point function. If the<br />

constructed fu~lction is a scalar point function, the TI we call it the divergence of a vector<br />

point function. More precisely, we give the follow\ng definition :<br />

Definition :<br />

Iff (x, y, z) be rfny given continuously differentiable vectorfunction, then the function<br />

is called the divergence off or divergence of the vectorfield defined by f and is<br />

denoted by divf.<br />

hi terms of operator V, we write<br />

The symbol "V." is pronounced as "del dot".<br />

Let us consider a cartesian coordinate system Oxyz. Let f have the scalar field<br />

components fi,fZ, f3 along the directions of x, y, z axes respectively, so that<br />

:. divf = V.f =<br />

- -


Vector Calculus<br />

v -<br />

%-,-<br />

P<br />

While writing in the above form, we have the understanding that<br />

a fl<br />

in the dot product i - - (L (:: fi) means the partial derivative -, etc. It is a<br />

ax<br />

convenient notation that is being used.<br />

From the definition and the notation used, i.e., V .f, it is clear that the divergence of a<br />

vector field function is itself a scalar field. Thus we can construct a scalar field from a<br />

vector field by taking its divergence.<br />

The meaning of divergence of a vector field is indicated in the name itself. div f is a<br />

measure of how much the vector field f diverges (or spreads out) from a point.<br />

We shall now give the physical interpretation of divergence. But before, that we would<br />

like to mention that the function div f is a pointfunction. By a point function we mean<br />

that the value of div f is independent of the particular choice of coordinates, i.e., its<br />

value is invariant w.r.to coordinate transformation. Learner interested in knowing the<br />

details about the invariance of the divergence may see Appendix-<strong>11</strong>1.<br />

1 1 .S.l'Physical Interpretation<br />

We consider the motion of a compressible fluid in a region R having no sources or<br />

sinks in R, i.e., no points in R at which fluid is produced or disappears.<br />

Let p (x, y, z, t) be the density of fluid and v = v (x, y, z, t ) be the velocity of fluid<br />

particle at a point (x) y) z) at time t.<br />

Let V = p v, then V is a vector having the same direction as v and a magnitude<br />

1 V 1 = p 1 vl. It is known as 'flux'. Its direction gives the direction of the fluid flow and<br />

its magnitude gives the mass of the fluid crossing per unit time a unit area placed<br />

perpendicular to the flow.<br />

If the unit area is placed with its normal at an angle 8 to the flow (see Figure <strong>11</strong>.14),<br />

then the nlass of the fluid crossing it per unit time is<br />

4 1. cos 8 v p(1 .cos8)v = Vcos8<br />

Unit normal<br />

Figurc <strong>11</strong>.14 : Flux<br />

= Resolved part of V in the direction<br />

of the normal to the area.<br />

,<br />

Consider a small fixed rectangular parallelepiped of sides bx, 6y, 62 parallel to the<br />

coordinate axes as shown in Figure <strong>11</strong>.15 enclosing a point P (x, jr, z).<br />

--<br />

A<br />

A A<br />

Let V = V,i+V,j+ V,k<br />

Figure <strong>11</strong>.15 :<br />

The velocity component parallel to y-axis at any point of the face ABCD.


omitting powers of 6y higher than one.<br />

1 - Vy(x,y,z)+ j 6 ay, ~<br />

Thus the mass of fluid that lnoves out of the face ABCD in time 6t<br />

slimilarly, the mass of the fluid that enters through the face A' B' C'D' in time 6t.<br />

Thus the net mass of the fluid that moves out through the faces ABCD and A' B' C' D'<br />

perpendicular to y-axis.<br />

Similarly, considering the other two pairs of faces, we see that the total mass of fluid<br />

flowing out of the parallelopiped iq time 6t<br />

The volume of the parallelopiped is 6x. 6y .6 z. On taking the limit, when<br />

6x, 6y, 6 z, 6 tall tend to zero, the amount of fluid per unit time that passes through a<br />

point P (x, y, z).<br />

a v, a vy a v,<br />

+ + - -<br />

ax ay az<br />

- div V<br />

Thus div V = div ( p v) gives the net rate of fluid outflow per unit volume per unit time<br />

at a point of the fluid.<br />

The outflow will cause a decrease in the density of the fluid inside the parallelopiped,<br />

say- 6p in time 6t. Thus loss of mass per unit time per unit voluine at a point<br />

E -<br />

at'<br />

Equating the loss of mass to thq outflow, we get<br />

I<br />

-<br />

at<br />

=, div(pv) + * = 0<br />

at<br />

This import'ant relation is called the condition for conservation of mass or the<br />

continuity equation of a compressible fluid flow.<br />

(<strong>11</strong>.19)<br />

Similarly, we can discuss the flow of electricity or the flow of heat or flow of particles<br />

from a radioactive source or water flowing into a drain. Thus, in gener~zl, ifF is any<br />

vector field defined at all points in a given region, then the divergence of F at any<br />

point represents the flu. per unit volume out of the volume dV enclosing the point, as<br />

dV is made smaller and smaller, i.e., as dV + 0.<br />

You may note that the net rate of fluid out flow ht a point P is positive i.e., div V > 0<br />

when the fluid has the tendency to diverge away from P, but if the fluid flows towards<br />

the point, then div V < 0. Thus a point ofpositive divergence means that there is a net<br />

outflow )'?om that point. Similarly a point of negative divergence implies a net inward<br />

I<br />

Vector ,Merentid Calculus


Vector Calculus<br />

If we consider the steady fluid motion of a<br />

p is constant, then equation (<strong>11</strong>.19) becomes<br />

div v = 0<br />

i.e., the rates of outflow and inflow are equal for any given volume at any time, i.e., the<br />

amount of the material in a volume remains constant.<br />

We know that for a magnetic field, the lines of force are closed -they neither flow out<br />

of a point nor into a point. Thus for a magnetic field B, we have<br />

div B = 0<br />

Thus there exists vector fields where divergence is zero. Such vector fields are called<br />

divergencefree or solenoidal. We give below the formal definition of solenoidal vector<br />

field.<br />

Definition :<br />

A vector field F is called divergencefree or solenoidal in a given region iffor all<br />

points in that region<br />

V.F = 0<br />

Thus magneticfield or velocily of a steady flow of compressiblefiid are examples of<br />

solenoidal vector fields.<br />

We now apply the concepts discussed in this section to some examples.<br />

Example 12 :<br />

Solution :<br />

Find the divergence of the vector<br />

A = 2y; + ~~Zi-2xzj.<br />

From definition,<br />

a a a<br />

divA = -(2y) +'- (-2x2) + -(2yz)<br />

ax ay a z<br />

Example 13 :<br />

5 2xy + 0 + 2y<br />

= 2y(x+l)<br />

Showthatdiv(grad$) =n(n+~)$-~.<br />

Solution : , ,<br />

Let r be the distance of a point P (x, y, z) from a fixed poi<strong>11</strong>tA (x,, y, 2,).<br />

:. r = \lfx-xo~2 + (-Y-yo)2 + (~-10)~<br />

:. P - [(x-x0)2 + ( y - ~ + ~ (z-s)2]y ) ~<br />

Now<br />

-<br />

% i-<br />

- n<br />

grad(P) 2 ; &[(X-X~)~ + (y-Yo)2 + (Z-Z~)~]~<br />

n


Hence the result.<br />

Example 14 :<br />

Solution :<br />

~howthatthevector~=(x+3~);+ (y-3z)~+(x-2z)~issolenoidal.<br />

We know thatA is solenoidal if div A = 0.<br />

a a a<br />

Now div A = -(x+3y) + -Cy-3z) + -(x-22)<br />

ax ay a z<br />

= 1 + 1 - 2<br />

= 0.<br />

HenceA is solenoidal vector.<br />

Example 15 :<br />

A rigid body is rotating about a fixed axis with a constant angular speed o. The<br />

velo~ity vector field V of the rigid body at any point r is given by V = o x r. Show<br />

that V is a divergence free vector.<br />

Solution :<br />

If V is a divergence free vector, then<br />

div V = 0.<br />

Let z-axis be the axis of rotation for the rigid body.<br />

If r is the position vector of any particle P of the rigid body, then<br />

Vector Difltrtntial Calculus<br />

2(x-d2


Vector Calculus<br />

:. V = Velocity of P = o x r<br />

A A A A<br />

= ok x (xi + y j + z k)<br />

A A<br />

= oxj - o yi<br />

Now, by definition,<br />

a a a<br />

div V = - (-my) + - ( ox) + -(O)<br />

ax ay a z<br />

Hence velocity vector V is a divergence free vector.<br />

How about trying a few exercises now.<br />

E 17<br />

A A A<br />

Ifr = xi+yj+zk,showthat<br />

E 19<br />

(a) divr=3<br />

(b)<br />

(c)<br />

(d)<br />

div ( r/r3 ) = 0<br />

div ( r 4 ) = 3 4 + r. grad 4, where 4 is a scalar function ofx, J; z.<br />

A 2<br />

div ( r ) =<br />

lKjz7<br />

The gravitational forcep of aFraction of two particleh is the gradient of the<br />

scalar function f ( x, y, z ) = -. Show that for r > 0,p is a solenoidal vector.<br />

r<br />

Determine the electric field and charge distribution corresponding to potential<br />

2 2 a3<br />

4 = a? add 4 a(-a + -). r<br />

(Hint : Electric field, E = - V 4 and charge distribution p = E V . E).


Vector Calculus<br />

1 1.6.1 physical Interpretation<br />

Let F be a continuosly differentiable vector field.<br />

Then<br />

A aF,-aF,<br />

CurlF=VxF=i - dF, dF, Y<br />

( a a-.)+~(X-xF)+i(dd".$$)<br />

Let us consider the z-component of V x f, i.e.,<br />

Now ( V x F ), will be always positive if<br />

-<br />

dF a F x<br />

(i) increases and - decreases or<br />

ax ay<br />

dF dFX a F,, aFx<br />

(ii) > -, when both - and - are positive.<br />

ax ay ax<br />

-<br />

a))<br />

A A<br />

The projection of F on xy-plane will be OA = F, i + F,, j. At point A (x, y, O), the<br />

dF<br />

y-componetlt of F increases by the factor dx and the x-component of F decreases<br />

ax<br />

a Fx<br />

by the factor - dy due to displacement (h, dy, 0) in A. Hence at the point B (x + &,<br />

ay<br />

Y + dy, O),<br />

and the resultant displacement is AB, which has turned left. If we give a further<br />

displacement, F,, would increase and F, would decrease, giving a further resultant BC.<br />

Thus in going from point A to C, the field vector has rotated anti-clockwise.<br />

We can relate this rotation of F to the z-component of Curl F. From the right-hand<br />

rule, the direction of ( V x F ), will be along the ziaxis. The magnitude of ( V x F ),<br />

tells us how the magnitude of the field vector F changes as it rotates.<br />

We can extend this argument to the x and y compol~ents of V x F and can say that the<br />

x and y components of V x F represent the rotation about x and y-axes, respectively.<br />

Thus curl of a vector field gives us an idea of its rotation about an axis, It is termed ns<br />

VORTEXfield. The directiorl of V x F (i.e., curl F) is along the axis about which the<br />

vector field F rotates (or curls) most rapidly and I V x F I is a measure of speed of<br />

this rotation.<br />

The sense of rotation (clockwise or anti-clockwise) is determined by the right-hand rule.<br />

The curl of a vector functions plays an important role in many applications. Its .<br />

significance will be explained in more detail in Unit 12. At present we confine<br />

ourselves to some simple examples.<br />

<strong>11</strong>.6.2 Rotation of a Rigid Body<br />

A rotation of a rigid body B in space can be simply addpniquely described by a vector<br />

o. The direction of o is that of the axis of rotation and is such that the rotation appears<br />

clockwise if we look from the initial point of o to its terminal point. The magnitude of<br />

o is equal to the angular speed o ( 0) of the rotation, i.e., the linear (or tangential)<br />

speed of a point of B divided by its distance from the axis of rotation<br />

(see Figure <strong>11</strong>.16).<br />

Let P be any point of the body B and let d be the distance of P from the axis of<br />

rotation. Let r be the positi6n vector of P referred to some origin 0 on the axis of<br />

rotation. Then<br />

d = I r I sin y, where y is the angle between o and r.


0V<br />

Figure <strong>11</strong>.16 : Rotation of a rigid body<br />

From above and the definition of vector product, the velocity vector V of P is given by<br />

Let the axis of rotation be along the z-axis and we choose right-handed cartesian<br />

A<br />

coordinates such that w = o k<br />

Then<br />

A A A A<br />

V = (<strong>11</strong> x r = (wk) x (xi + y j + zk)<br />

(Here (I) is a co~~$tant)<br />

Hence in the case of rotation of a rigid body, tile curl of the velocity field Itas the<br />

direction of the axis of rotation and its magnitude equals twice the angular speed o of<br />

the rotation.<br />

We can also say that the curl of a velocity vector field at a point of a rotating body<br />

represents iwice the rate at which a material element occupying the point is rotating. -<br />

This result justifies the tenn rotation which is sometimes used for curl. We so~neti~nes<br />

write rot f instead of curl$<br />

It may be remarked that a similar velocity field is obtained by stirring coffee in a cup.<br />

We now take up another property of curl.<br />

For any twice continuously differentiable scalar field f,<br />

Vtctor Differential Calculus


vector Calculue and<br />

= 0 + 0 + 0 ( '.'f is twice continuously differentiable)<br />

= 0<br />

Hence ifa vector function is the gradient of a scalar function, its curl is the zero vector.<br />

Since curl characterises the rotation in a field, we may say that gradient fields<br />

describing a motion are irrotational.<br />

If a vector field occurs which is the gradient of a scalar field and if such a field occurs<br />

not as a velocity field, then it is usually called conservative.<br />

We can now give the formal definition of irrotational vector field.<br />

Definition :<br />

A field that has a vanishing curl everywhere is called an irrotational (or<br />

conservative) field i.e., ifcurlA = 0, the vector fieldA is called irrotational.<br />

Let us now take up a few examples.<br />

Example 16 :<br />

In electrostatics, the force of attraction (or repulsion) between two particles of<br />

oppasite (or like) charges Ql and Q, is<br />

where K = Q, Q, /(4 x E) and E is the dielectric constant. Is p an irrotational field ?<br />

Solution :<br />

p is an irrotation field, if curl p = 0.<br />

. .<br />

Now<br />

Then, by definition,<br />

curl p =<br />

A A A<br />

i j k<br />

a/ax a/a y a/a z<br />

Kx K_y Kz<br />

(x2+y2+2)3R (y2+y2+t)3/2 (2+y2+t)3/2


Let us take up another example.<br />

Example 17 :<br />

Solution :<br />

= 0 + 0 + Q= 0 . Hencep is an irratatio~lal field.<br />

Show that the vector field defined by<br />

F = 2xy2; + 22j+ 3xzY2i<br />

is irrotational. Find a scalar ~otential u such that F = grad u.<br />

By def.,<br />

=o+o+o=o<br />

Hence F is irrotational and hence F can be expressed as grad u.<br />

Also<br />

F = 2xyil + $2; + 3xzy2i<br />

Comparing the two expressions for F, we get<br />

a u a u a u<br />

- = 2xy2, -<br />

ax ay = x 22, = 32y.2<br />

NOW au a u a u<br />

du = - d x + -dy + - d z<br />

ax ay a z<br />

You may now try the fallowing questions and see whether you have underst .ood the<br />

concepts given in this sectiod.<br />

E 20<br />

FindCurlF,whereF = ~(1?+y~+3-3xyz).


I Vecbr Calculus<br />

E 21<br />

A fluid notion is given by<br />

A A A<br />

q = (y+z)i + (z+x)j + (x+y)k.<br />

Is this motion imtational ? If so, find the velocity pote~ltial.<br />

You know that the operator V is a vector operator. You can use this operator to prove<br />

some formulas on gradient, divergence and curl. We shall now state these formulas in<br />

the next section.<br />

<strong>11</strong>.6.3 Formulae on Gradient, Divergence and Curl<br />

You already know from Sections <strong>11</strong>.4, <strong>11</strong>.5 and 1 1.6 that ,<br />

(i) The V operator, when operated on a scalar field f gives rise to a vector field<br />

V f (gradient of j).<br />

(ii) Scalar product of V with a vector field F gives a scalar field, divergence of F<br />

or (V . F).<br />

(iii) Vector product of V with a vector field F gives a vector field, curl For<br />

(V x F).<br />

We now give the following formulas :<br />

(1) grad(+9)=4gradV+9grad+.<br />

(<strong>11</strong>) grad(A.B)=A x curlB +B x curlA +(A.V)B + (B.V)A.<br />

(In) grad (div A) = Curl Curl A + v2 A.<br />

(IV) div(A x B) = B.curlA-A.curlB<br />

(V) div (curl A) = 0<br />

(VI) div CfgraG )= f v2 + + v f. V g<br />

(W) curl($A)=(grad+) x A t+arlA.<br />

(WI)curl(A x B) =A div B-BdivA +(B.V)A -(A .V)B.<br />

For the proofs of formulas I-WI see Appendix-V.


Using the above formulas, you may now try the following exercises : v,tor MCPC.~~ Calculos<br />

E 24<br />

If a is a constant vector and r denotes the position vector of any point in space<br />

and iff = ( a x r ) S, show that<br />

divf - 0 and curl f = (n+2)r"a-nP"-2(a.r)r.<br />

Show that<br />

<strong>11</strong>.7 SUMMARY<br />

(a) div[(rxa)xb]=-2(a.b)<br />

@) grad[r,a,b] =a x b<br />

(c) curl(rxa)=-2a<br />

(d) div(rxa)=O<br />

(e) grad(a.r)=a,<br />

where a and b are a constant vector and r is the radius vector.<br />

We will now summarise the result of this unit.<br />

* A function which is defined at each point of a certain region in space and<br />

whose values are real numbers, depending on the points in space but not on<br />

particular choice of coordinate system, is called a scalar fectwn.<br />

* If I$ be a function which associates a unique scalar with each point in a given<br />

region, then 4 is called a scalarfield.<br />

* If to each point P of a certain region in space, or if to each set P of variables of<br />

a certain sets of variables, a vector V (P) is assigned, then V (P) is called a<br />

vector function.<br />

* If F be a function which associates a unique vector with every point in a given<br />

region, then F is called a vector fwM.<br />

* A vector function f (t) of a real variable t is said to tend to its limit I as t<br />

approaches a iff (t) moves into the position occuped by I as t -. a.<br />

* A vector'function f (t) is said to be continuous at t = a if it is defined in some<br />

neighbourhood of a and<br />

* The vector function f (t), of scalar variable t, is said to be differentiable at a<br />

. point t if the limit ,<br />

lim f(t+At)-f(t)<br />

At-0 At


Vector Calculu<br />

exists and it is then denoted by gf dt or f ' ( t ).<br />

The derivative of a vector function represents a vector in the direction of the<br />

tangent to the curve traced by the vector function.<br />

If a vector fu~ictio~i is expressed iii terms of its comnponeiits, then the limit,<br />

continuity and differentiability of the function exists provided limit, continuity<br />

and differentiability respectively of each component exists.<br />

The partial derivative gives the rate of chahge off (u , v) w.r.to u at a<br />

a u<br />

given point (u, v) in space.<br />

Thevector<br />

is called the gradient of scalarf; where f is a contin~rousl~ differentiable<br />

function. Here the symbol V is read as del and represents vector operator<br />

The directional derivative of a scalar function f along the vector b can be<br />

written as<br />

* The gradient of a scalar field f is a vector normal to the level surface<br />

f = constaiit and has a magnitude equal to the rate of change off alo~ig this<br />

normal.<br />

* If F (x, y, z) be any given continuously differentiable vector function, -<br />

then the function<br />

is called the divergence of F and is written as div F or V . F.<br />

* The divergence of a vector field represents the net outward flux per unit time<br />

at any point of the vector field.<br />

A vector field F is called divergence free or solenoidal field in a given region<br />

if for all points in that kgion div F = 0.<br />

If F (x, y, z) be any continuously differentiable vector function, then the<br />

function<br />

is cilled the curl of F and is denoted by curl F or V X F.<br />

The &rl of a vector field represents a rotation about an axis.<br />

A field that has a vanishing curl everywhere is called imtational (or<br />

conservative) field, i.e., if curlA = 0, thenA is called irratational.<br />

<strong>11</strong>.8 SOLUTIONS1 ANSWERS<br />

I<br />

El (a) Vector function . , -<br />

@) Vector Function<br />

-<br />

" (c) Scalar field<br />

/-<br />

-


(d) Scalar function<br />

(e) Vector function<br />

E2 (a) The limits of the components of f (t) as t -, 0 are<br />

(b)<br />

Lim ef = 1 , Lim t ef = 0<br />

r+O I--0<br />

n<br />

:. Lim f (t) = i<br />

r + 0<br />

The limits of the components off (t) as t -4 0 are<br />

n n<br />

:. Lirn f (t) = j - k<br />

1-0<br />

(c) The limits of the components off (t) as t 0 are<br />

I--0<br />

t 1 - cos t 0 + sin t<br />

Lim - = 1, Lim = Lim - = 0, Lim 1 = 1<br />

, , s1n t t-o t r-o 1 1-0<br />

E3 (a) The function<br />

A n<br />

:. Lim f (t) = i + k<br />

1-0<br />

n n<br />

f( t ) - (COS t) l + (sin t) j + k<br />

is cot~tinuous at every value oft, because each colnpollel~t is co~~ti~~uous<br />

, for all<br />

t .<br />

(b) The function<br />

A<br />

f (t) = d i + cos - j+<br />

1 +t<br />

log <strong>11</strong> +rl i<br />

1<br />

is continuous at all t except t = -1, because the function cos -is not defined<br />

l+t<br />

at t = -1 and at all other values of r each component is continuous.<br />

E4 (a) me derivative of f (I) = e2' 7 + t e-'<br />

is f' (t) = 2 2' I+ - re-');<br />

The function f (t) and its derivative f' (t) are defined at every value oft.<br />

(b)<br />

(c)<br />

He& f' (t) = 0 ,became f (t) is a constant vector. Domain is all t<br />

2 3<br />

Here f ' ^ l n<br />

(t) =<br />

r n j + ~ + g ~ . l - ; i . k<br />

Domain : For (21 < 1 and t 0.<br />

ES From the chain-rule, we have<br />

. -<br />

- d d du<br />

-f dx 0 - p 4 dx<br />

vector Diffemtirl C~rolos


Vector Calculus<br />

This result is valid for (x + 2) * 0<br />

E6<br />

Alsa I r x V = cos cllt sin o)t 0<br />

- w sin ot w cos ot 0<br />

E9<br />

- - l sin (42TZX)<br />

. 2 (X + 1) A 3 /- I<br />

2 d z ,<br />

n L\X + 1)<br />

= - i sin (X + 1) + j<br />

-,- . 1;<br />

+j l+(l?+~r+ly<br />

" 1<br />

, + k- 2<br />

1+(x+ 1) (x + 2)<br />

The position vector of the particle at any time t is<br />

A<br />

r (t) = e-' i + 2 cos 3t j^ + 2 sin 3t i !<br />

d, . A<br />

Vel. - V (t) = -r (l) - -I' 1 - 6 sin 3tj + 6 ms 3t i<br />

dt<br />

A A<br />

Hence V (0) = - i + 6 k<br />

d<br />

Accel =a (t) = - V(t) = e-' 1 - 18 ms 13 - 18 sin 3t i<br />

dt<br />

Here<br />

A A<br />

r = cos wt i + sin wt j<br />

:. V- -<br />

* dt<br />

-<br />

:. r . V = - w sin wt cos wt + w sin wt cos wt = 0<br />

:. V is perpendicular to r<br />

dr = -w sin wt 1 + w cos wt j<br />

A A A<br />

i j k<br />

= i (u, cos2 wt + o sin2 or = o i, which is independent of t.<br />

)<br />

Hence r x V is a constant vector.<br />

Here r (t) = ? b<br />

:. Vel=V=r=ab+?b<br />

Here 2 2 b is the velocity of any point P relative to the disk and b is the additional<br />

velocity due to rotation.<br />

-2b+4tb+?b<br />

Also be#cause of rotation, b is of the form


:. 6 = - o sin o t I + o cos wtj<br />

. . 2 A 2 2<br />

and b=-w coswti-o sinwtj=-o b<br />

Hence the Coriolis acceleration = 4t 6, which is due to the interaction of rotation of<br />

the disk and the motion of P on the disk. It is in the direction of 6, i.e., tangential to<br />

the edge of the disk and it points in the direction of rotation.<br />

El0 Here f = 2 + J + 4Vz<br />

:. gradf - ~r i + 4yzi + (2y + 4x2) 3 + + i<br />

A A A<br />

Herea = 2i - 2 j + k<br />

:. Required directional derivative<br />

The minus sign indicates that & decreases in the direction under consideration.<br />

.,<br />

xL - yL<br />

Here f (x , y) = -<br />

V<br />

:. grad f = -<br />

A A<br />

:. (gradf) l,l = 2i - 2j A A<br />

Let a = a, i + a2 j be the direction in which directional derivative off is zero.<br />

Hence the required direction<br />

a1 i + a1 j<br />

-\~;lf+af<br />

:. grad f = xi+xj+i!li<br />

ax ay az<br />

a l + 1<br />

= -(i + j )<br />

= (6al) fl


Vector Calculus<br />

Nnw<br />

:. Directional detivative alone z-axis<br />

El3 Here r - q2 + y2 + 2 , where(x, y, z) is any point in space<br />

a p A ap<br />

Hence grad f" - - i +-<br />

ax ay<br />

- m #-7 [x i + y j + zj)<br />

El4 Here parametric represertation of the curve is<br />

If i' is the unit tangent to the curve, then<br />

At (1,1,1) for the given curvet - 1<br />

At (-1,1,-1) for the given curve t - -1<br />

A A A<br />

Now f-xy2+y2+&<br />

:. (gnuij) - (J + a) i + (2 + ~ry)<br />

.'. (grad fl(l,l, 1) -<br />

:. Directional derivative at (1, 1, 1) along the tangent<br />

3 + (2 + 2yz) i<br />

I<br />

A , , - 3(i + j + k)<br />

1 : A A<br />

(1 + 2j + 3k) 1<br />

F<br />

i


Also Directional derivative at (-1 1, -1) along the tangent<br />

A A A<br />

= (3i-j- k) .- m<br />

=-<br />

:IS Letf = x3 + y3 + 3xyz - 3<br />

A vector normal to the surface is<br />

1<br />

2<br />

(3+2-3) = -<br />

m<br />

(7 - $+ 3;)<br />

= (3.12 + 3.Z2.(-1)) 7 + (3.2 + 3.1 . (-1)) j + (3.1.2)<br />

A A A<br />

--3i+9j+6k .<br />

- -3 i + 9j + 6k 1<br />

A<br />

A A A<br />

-- (-i + 3j +<br />

3 m -m<br />

A A -<br />

a z<br />

(a) divr-@+*+- - 1 +1+1= 3<br />

ax ay az<br />

U;)<br />

Vector Dilferenti.1 Calculus


(d) Hen lrl zd2+y2+2<br />

A A A


-<br />

'.' div p - 0, :. p is a solenoidal vector.<br />

El9 When4 - a?<br />

When - a (-2 + $)<br />

Charge distribution - E<br />

V . (=iX)<br />

E = E + - - 5<br />

E2O HereF-V(r'+y'+P-3xyt)<br />

:., cud F - curl [V (2 + y' + 2 - 3xyt)]<br />

(9) A(?)] +


Vector Calculus<br />

Hence the fluid motion is imtational.<br />

Let $ be the velocity potential<br />

A<br />

Comparing the two expression for q , we get<br />

--lydr+afx+zdy+xdy+xdz+y&]<br />

- [(ydx + xdy) + (zdr + xdz) + (zdy + ydz)]<br />

= -d (xy + yz + zx)<br />

:. $ = - (xy + yz + uc) , which is required<br />

E22 Here V-= r = - =<br />

Curl v+=<br />

-<br />

A r x z<br />

7'<br />

42+y2+t<br />

A A<br />

I j k<br />

a/ ax aiay .aiaz<br />

5 1 z<br />

A


1<br />

1<br />

!<br />

I<br />

:. div f - div a x r S ( )<br />

- r S curl a - a - curl (r S)<br />

. a is a constant vector, :. curl a = 0 *<br />

i<br />

curl (rrn\ = I ajar<br />

a<br />

NOW cud ((a x I] I"] - Z I x - ((a XI) $1<br />

ax<br />

a<br />

Now -((axr)S]-nS-l%(axr)+S ax-<br />

i ax ax % ( ::)<br />

I<br />

I + A a<br />

I :. ix-[(axr)S)- x ~ $ - ~ ~ x ( u x ~ ) + s ~ x ( u x ~<br />

ax<br />

i<br />

A A A<br />

= n . #'-2<br />

I x [(i. r) a - ( i. a) r]<br />

A A<br />

I<br />

+S[(i. i) a - (;.a) E ]<br />

I<br />

I - n - P-~. A? a - n . f-2 (xi- (I) r<br />

I<br />

I<br />

a<br />

+~a-$(1.a)C<br />

. . ix-[(a~r)f]-n$-~a(2+~*+2)<br />

dx<br />

Hence the result<br />

E24 (a) div [(r x a) x b]<br />

- div [(r b) a - (a . b) r]<br />

-nf'%(r-a)+3Pa-f<br />

- r" a (n+2) - n S-* (a r) r<br />

Vector DiRerential Calculus


vector Calculus<br />

= div ~ b ~ + ~ b ~ a-(a-b) + z &<br />

- [( )<br />

and<br />

- A A A A A A A A A<br />

= xi + yj + zk, a = ai i +a2j + a3k, b bl i + 6j + b3 k.<br />

grad [r 9 a , bl<br />

- V [r . (a x b)]<br />

A - i (a bk-am<br />

6); - component + Y (a x b)j-wrnp,nt + (a x b) - amp,mn]<br />

A A<br />

A A A<br />

= 2 ((?* a) i - (i i) a<br />

- a - 34<br />

= -2a<br />

div(r x a) - a . curl r - r curl a<br />

= I(0-0) +j (0-0) + C^ (0-0)<br />

= 0<br />

and curl a - 0 (-: a is a constant vector )<br />

div (r x a) - 0<br />

grad(a r) - V (alx + a2y + ajz),<br />

A A A<br />

ifa- al i + azj+a3k<br />

.-. gwd (a . r) - 2<br />

A A A<br />

r - xi + yj + zk<br />

A a<br />

ic (OLx + a2y + a3z 1


I<br />

[<br />

I I<br />

I<br />

APPENDIX - I : Derivative along a curve<br />

Let a curve C in space be given.by<br />

f( x, y, z ) - constant<br />

If we take a fued point A on the given curve C, from which we measure the arc length s of<br />

the curve, then along the curve, f (x , y, z) is a function of s.<br />

I The position vector r (s) of any point P on the curve at a distance s from A is given by<br />

I<br />

Let Q be a neighbouring point of P, such thatA Q - s + 6 s. Then from calculas, by the<br />

[ definition of differential of a scalar hndion, we have<br />

I<br />

jdx ! ^ d z<br />

ds +k-<br />

ds<br />

ds +'<br />

A<br />

where t denotes unit tangent to the curve.<br />

Thus the derivative along a curve is tip dot product of the unit tangent to the curve and<br />

gradient of the function, defining the curve.<br />

Now if the function f is constant along the curve, then its gradient is zero. This shows that<br />

A<br />

'The derivative (t. grad f) along a curve is zero if and only if, the function f is constant<br />

along the curve'.


Vector a Calculus<br />

APPENDIX - I1 : Invariance of length and direction of =ad f<br />

- -- -<br />

The grad f is given by<br />

A A A<br />

where i , j , k are unit vectors in the direction ofx, y & z - axes of a cartesian coordinate<br />

system, say Oqz and (r, y, z) are the coordinates of any point P referred to Ow (see<br />

Figure <strong>11</strong>.12).<br />

Since relation (1) involves partial derivatives, which depend on the choice of the<br />

coordinates, hence the result is not obvious.<br />

Now, by the definition of a scalar function, the value off at a point P depends on the<br />

location of P but is independent of the coondinates. Also the arc length s of ray C is<br />

independent of the choice of coordinates. . b<br />

Hence<br />

Q-p<br />

contains only the quantities which are independent of the choice of coordinates.<br />

Also<br />

where y is the angle between grad f and k. From above, 2L is maximum when cos y = 1,<br />

as 1<br />

i.e., y = 0 and then - I grad f I. This shows that thelength and the direction of y df is<br />

as<br />

independent of the choice of coordinates.


vector Caleulue<br />

((be;) ax + (;.$I +<br />

a Y<br />

. a .<br />

Thus the expression for div f in terms of ,? ,I , k<br />

-2; x<br />

'ax1, .<br />

-2 ;.I;.;) Lf ax + (;.<br />

(since is dimctional derivative alongxr axis, i.e., along i<br />

ax<br />

direction)<br />

- p [(i.?)<br />

- i. ax<br />

"X<br />

i ax . + ( J. ~<br />

A A A<br />

A A A<br />

= exp~ssion of div f in terms of i , j , k.<br />

Thus div f is independent of the choice of coordinates, i.e., div f is essentially a point<br />

- .<br />

function.


APPENDIX - IV : Proof of Formulae on Divernence<br />

I)<br />

- kV.f<br />

- k div f<br />

-<br />

-f.grad$ + $ V .f<br />

$ div f + f. grad $<br />

If $ is twice differentiable scalar function, then<br />

a * + a f i a*<br />

By definition, div ( grad 41 ) - -<br />

ax (ax) a, (a,)+ a~(a,)<br />

The expression on the right hand side is called the Laplacian of scalar fundion $ and<br />

is denoted as v2 $.<br />

Thw.div ( grad a$ ) - V' $.<br />

vcdsr Dur~cwlid Cddw


vector Calculus<br />

APPENDIX -V : Proofs of Formulae on Gradient,<br />

Divergence and Curl<br />

(* ddt))=( ax)<br />

NOWAX ix- A.Z<br />

InterchangingA and B in (2), we get<br />

From (I), (2) and (3), we get<br />

A<br />

:-(~.i)dB<br />

ax<br />

=A x curl B +(A.V)B<br />

grad (A.B) =A x curlB + B x curlA + (A.V) B + (R.V)A<br />

<strong>11</strong>1) curl ( curl A ) = V x ( V x A )<br />

= V (V.A) - (V.V)A,<br />

(applying the triple vector product formula and taking care that the variable follows<br />

the del operator).<br />

* grad(div A) = curl (curlA) + v2 A


div ( curl A ) =<br />

x ~ x ax , ~ . ~ - x a ax ; xA - (applying solar . triple property product) of<br />

= (curl A ).B - (curl B).A.<br />

- k.- + i.- + j.-<br />

A aLA A a L ~ aLA aLA<br />

because ): k . ~<br />

x ay axay ayaz axaz<br />

" a 2 ~<br />

& similar expression for j . - axaz*<br />

[r*+'j*+i")]<br />

= div f i ax ay az


curl(A x B) = V x (A x B)<br />

Assume V = VA + VB , where we consider only A as variable while operating VA<br />

and only B is to be considered as variable when operating by VB. Since this takes<br />

care of the product as far as differentiation is concerned and the VA and V5 can be<br />

treated as vectors and rules of multiplication of vectors can be applied.<br />

:. cud (A x B)= V x (A x B)<br />

=VAx(AxB)+VBx(AxB)<br />

(VA.B)A - (VA.A) B + (VB.B)A - (VB.A) B<br />

(B.VA)A - B (VA.A) +A (VB.B) - (A.VB) B<br />

= (B.V)A - B (V.A) + A (V.B) - (A.V) B<br />

=A div B - B divA + (B.V)A - (A.V)B<br />

Note that you should exercise special care to avoEd inadvertent mistakes. For this<br />

reason, in the first two steps, the variable vectors A, B are both written to the right of<br />

del operator. In the next step, we have brought all those vectors which can be treated<br />

as constant to the left of the operator. Since<br />

V .A - VA .A + VB .A - VA .A + 0<br />

. B = V5 . B, therefom in the final step VA & VB have been<br />

andsi&larly~<br />

replaced by V as th<br />

is makes no difference.<br />

1<br />

i<br />

i<br />

1

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!