Search Continuum Mechanics Website

# Vector Calculus

## Differentiation With Respect To Time

Differentiation with respect to time can be written in several forms.

$\qquad \; \; \text{velocity} \qquad = \qquad {d{\bf x} \over dt} \qquad = \qquad \left( {dx_1 \over dt} , {dx_2 \over dt} , {dx_3 \over dt} \right) \qquad \ = \qquad \dot{\bf x} \qquad = \qquad \dot{x}_i \qquad = \qquad x_{i,t}$
$\text{acceleration} \qquad = \qquad {d{\bf v} \over dt} \qquad = \qquad \left( {dv_1 \over dt} , {dv_2 \over dt} , {dv_3 \over dt} \right) \qquad \ = \qquad \dot{\bf v} \qquad = \qquad \dot{v}_i \qquad = \qquad v_{i,t}$
$\text{acceleration} \qquad = \qquad {d^2{\bf x} \over dt^2} \qquad = \qquad \left( {d^2x_1 \over dt^2} , {d^2x_2 \over dt^2} , {d^2x_3 \over dt^2} \right) \qquad \ = \qquad \ddot{\bf x} \qquad = \qquad \ddot{x}_i \qquad = \qquad x_{i,tt}$
One can use the derivative with respect to $$\;t$$, or the dot, which is probably the most popular, or the comma notation, which is a popular subset of tensor notation. Note that the notation $$x_{i,tt}$$ somewhat violates the tensor notation rule of double-indices automatically summing from 1 to 3. This is because time does not have 3 dimensions as space does, so it is understood that no summation is performed.

### Differentiation of a Vector

Suppose $${\bf v} = (5t^2, \sin t, e^{3t})$$. Then $$\dot{\bf v}$$ equals

$\dot{\bf v} = (10t, \cos t, 3e^{3t})$

### Helix Example

The position of an ant crawling around and up a pipe is given by $${\bf x} = (2 \cos t, 2 \sin t, 5t)$$. The velocity, $${\bf v}$$, equals

$\dot{\bf x} = {\bf v} = (-2 \sin t, 2 \cos t, 5)$
and the acceleration is

$\ddot{\bf x} = \dot{\bf v} = {\bf a} = (-2 \cos t, -2 \sin t, 0)$
which always points toward the center of the pipe.

## Differentiation With Respect To A Coordinate

Suppose you want to differentiate a function, $$f(x,y,z)$$, with respect to $$y$$. This is written as

${\partial f \over \partial y} \quad \qquad \text{or} \qquad \quad {\partial f \over \partial x_2} \quad \qquad \text{or} \qquad \quad f_{,2}$
where the comma is common tensor notation for a derivative.

In the more general case, differentiation with respect to $$x_j$$ is (yes, this is a gradient)

${\partial f \over \partial x_j} \qquad \qquad \text{or} \qquad \qquad f,_j$
Differentiation of a vector, $${\bf v}$$, is

${\partial {\bf v} \over \partial x_j} \qquad \qquad \text{or} \qquad \qquad \left( {\partial \, v_x \over \partial x_j} , {\partial \, v_y \over \partial x_j} , {\partial \, v_z \over \partial x_j} \right) \qquad \qquad \text{or} \qquad \qquad v_{i,j}$
Differentiation of a tensor, $$\boldsymbol{\sigma}$$, is

${\partial \boldsymbol{\sigma} \over \partial x_k} \qquad \qquad \text{or} \qquad \qquad \sigma_{ij,k}$
As with vectors, every component of a tensor is differentiated.

### Differentiation of a Vector

Suppose $${\bf v} = (3x^2 - 2y, z^2+x, y^3 - z)$$. Then $$\partial {\bf v} \over \partial y$$ equals

${\partial {\bf v} \over \partial y} = (-2, 0, 3y^2)$

## Gradient

The gradient of a function, $$f({\bf x})$$, is written, $$\nabla f({\bf x})$$, and is a vector. It is formed by differentiating the function with respect to each coordinate.

$\nabla f({\bf x}) = \left( {\partial f \over \partial x_1} , {\partial f \over \partial x_2} , {\partial f \over \partial x_3} \right)$

### Tensor Notation

The gradient can also be written as $${ \partial f \over \partial x_i}$$, or simply as $$f_{,i}$$.

### Gradient Example

Suppose $$f({\bf x}) = 3x^2 - 2 y z^2$$. Then the gradient is

$\nabla f({\bf x}) = ( 6x, -2 z^2, -4 y z )$
The gradient of a scalar function tells how much the function increases along each global coordinate. In the above example, $$f$$ increases at the rate of $$6x$$ along the $$x$$ axis, $$-2 z^2$$ along the $$y$$ axis, and $$-4 y z$$ along the $$z$$ axis.

Coincidentally, the gradient also gives the direction, or orientation, in space that corresponds to the greatest rate of increase. The following example, in 2D space, demonstrates this.

### 2nd Gradient Example

Take for example the paraboloid, $$f(x,y) = 2 x^2 + y^2$$. The gradient at $${\bf x} = (5,3)$$ is

$\nabla f(x,y) \quad = \quad ( 4x, 2y ) \quad = \quad (20, 6)$
Therefore, at $${\bf x} = (5,3)$$, $$f$$ is increasing at the rate of 20 along the $$x$$ axis, and at the rate of 6 along the $$y$$ axis. $$20{\bf i} + 6 {\bf j}$$ also corresponds to the direction in the $$x,y$$ plane along which $$f$$ will increase the most quickly.

Gradients of vectors can also be computed. The result will be a 2nd order tensor. For example, the gradient of a velocity field is written as $$\nabla {\bf v}$$. Writing this in tensor notation $$v_{i,j}$$ shows more clearly that the result is a 2nd order tensor because of the presence of the $$i$$ and $$j$$ subscripts. Gradients arise in mechanical deformation and heat conduction applications. Mechanical strains are related to gradients of displacements and heat conduction is related to the gradient of the temperature distribution.

## Divergence

The divergence of a vector is a scalar result, and the divergence of a 2nd order tensor is a vector. The divergence of a vector is written as $$\nabla \cdot {\bf v}$$, or $$v_{i,i}$$ in tensor notation. It is computed as

$\begin{eqnarray} \nabla \cdot {\bf v} & = & ( {\partial \over \partial x} {\bf i} + {\partial \over \partial y} {\bf j} + {\partial \over \partial z} {\bf k} ) \cdot( v_x {\bf i} + v_y {\bf j} + v_z {\bf k} ) \\ \\ & = & {\partial v_x \over \partial x} + {\partial v_y \over \partial y} + {\partial v_z \over \partial z} \end{eqnarray}$

### Tensor Notation

As stated above, the divergence is written in tensor notation as $$v_{i,i}$$. It is very important that both subscripts are the same because this dictates that they are automatically summed from 1 to 3. They can in fact be any letter one desires, so long as they are both the same letter.

### Divergence Example

If $${\bf v} = (3x^2 - 2y, z^2+x, y^3 - z)$$, then the dot product is

$\nabla \cdot {\bf v} \quad = \quad {\partial \over \partial x}(3x^2 - 2y) + {\partial \over \partial y}(z^2 + x) + {\partial \over \partial z}(y^3 - z) \quad = \quad 6x - 1$
The divergence of velocity vectors often arises in the discussion of incompressibility and conservation of mass.

## Curl

The curl of a vector is the cross product of partial derivatives with the vector. Curls arise when rotations are important, just as cross products of vectors tend to do. Rotations of solids automatically imply large displacements, which in turn automatically imply nonlinear analyses. And this is why one seldom comes across curls... because most analyses are linear.

Curls are calculated as follows.

$\begin{eqnarray} \nabla \times {\bf v} & = & ( {\partial \over \partial x} {\bf i} + {\partial \over \partial y} {\bf j} + {\partial \over \partial z} {\bf k} ) \times ( v_x {\bf i} + v_y {\bf j} + v_z {\bf k} ) \\ \\ & = & \left| \matrix { {\bf i\;} & {\bf j\;} & {\bf k\;} \\ {\partial \over \partial x} & {\partial \over \partial y} & {\partial \over \partial z} \\ v_x & v_y & v_z } \right| \\ \\ & = & ({\partial \, v_z \over \partial y} - {\partial \, v_y \over \partial z}) {\bf i} + ({\partial \, v_x \over \partial z} - {\partial \, v_z \over \partial x}) {\bf j} + ({\partial \, v_y \over \partial x} - {\partial \, v_x \over \partial y}) {\bf k} \end{eqnarray}$

### Tensor Notation of Curls

The curl of a vector is written in tensor notation as $$\epsilon_{ijk} v_{k,j}$$. It is critical to recognize that the vector is written as $$v_{k,j}$$ here, not $$v_{j,k}$$. This is because the curl is $$\nabla \times {\bf v}$$, not $${\bf v} \times \nabla$$.

An easy way to get the tensor notation right is to think of $$\nabla \times {\bf v}$$ as $$\epsilon_{ijk} \nabla_j v_k$$ and note the order of the subscripts. Of course, this reduces to the correct result: $$\epsilon_{ijk} v_{k,j}$$.

As with cross products, the fact that $$j$$ and $$k$$ both occur twice in $$\epsilon_{ijk} v_{k,j}$$ dictates that both are automatically summed from 1 to 3. The term expands to

$\matrix { \epsilon_{ijk} v_{k,j} & = & \epsilon_{i11} v_{1,1} & + & \epsilon_{i12} v_{2,1} & + & \epsilon_{i13} v_{3,1} & + & \\ & & \epsilon_{i21} v_{1,2} & + & \epsilon_{i22} v_{2,2} & + & \epsilon_{i23} v_{3,2} & + & \\ & & \epsilon_{i31} v_{1,3} & + & \epsilon_{i32} v_{2,3} & + & \epsilon_{i33} v_{3,3} }$

### Curls Using Tensor Notation

To obtain the yth component of a curl, set $$i$$ equal to 2 in the above equation.

$\matrix { \epsilon_{2jk} v_{k,j} & = & \epsilon_{211} v_{1,1} & + & \epsilon_{212} v_{2,1} & + & \epsilon_{213} v_{3,1} & + & \\ & & \epsilon_{221} v_{1,2} & + & \epsilon_{222} v_{2,2} & + & \epsilon_{223} v_{3,2} & + & \\ & & \epsilon_{231} v_{1,3} & + & \epsilon_{232} v_{2,3} & + & \epsilon_{233} v_{3,3} }$
All subscripts are now specified, and this permits evaluation of all alternating tensors. All of them will equal zero except two, leaving

$\epsilon_{2jk} v_{k,j} \; = \; v_{1,3} - v_{3,1} \; = \; \frac{ \partial \, v_x}{\partial z} - \frac{\partial \, v_z}{\partial x}$
which is again consistent with the determinant result (as it must be). Results for the xth and zth components are obtained by setting $$i$$ equal to 1 and 3, respectively.

### Curl Example - Rotating Disk

Consider a disk rotating about the $$z$$ axis such that

$\begin{eqnarray} x & = & X \cos (\omega \, t) - Y \sin (\omega \, t) \\ y & = & X \sin (\omega \, t) + Y \cos (\omega \, t) \\ z & = & Z \end{eqnarray}$
where $${\bf X}$$ is the vector of original coordinates of each point at $$t = 0$$, and $${\bf x}$$ is the vector of that point's coordinates at any other time, $$t$$.

Note that this is common in Continuum Mechanics to use $${\bf X}$$ as the position vector at $$t = 0$$, the so-called reference configuration, and $${\bf x}$$ for the position vector following any translations, rotations, and deformations, the so-called current configuration.

The velocity vector is

${\bf v} = {\partial {\bf x} \over \partial t} = \left( \begin{eqnarray} & - \omega X \sin (\omega \, t) - \omega Y \cos (\omega \, t) \\ & + \omega X \cos (\omega \, t) - \omega Y \sin (\omega \, t) \\ & 0 \end{eqnarray} \right)$
which simplifies to

${\bf v} = ( - \omega \, y, \omega \, x, 0 )$
making the curl of the velocity vector relatively simple to compute.

$\nabla \times {\bf v} = ( 0, 0, 2 \omega )$
As stated above, the curl is related to rotations. It turns out that $$\nabla \times {\bf v}$$ gives the axis of rotation, and $$\frac{1}{2} | \nabla \times {\bf v} |$$ is the rotational rate. So $$\frac{1}{2} ( \nabla \times {\bf v} )$$ gives

$\frac{1}{2} ( \nabla \times {\bf v} ) = ( 0, 0, \omega )$

## Laplacian

The Laplacian is the divergence of the gradient of a function. It often arises in 2nd order partial differential equations and is usually written as $$\nabla^2 \! f({\bf x})$$. It can also be written in the less popular, but more descriptive form of $$\nabla \cdot \nabla f({\bf x})$$. Its definition is

$\nabla^2 \! f({\bf x}) \equiv {\partial^{\,2} \! f({\bf x}) \over \partial \, x^2} + {\partial^{\,2} \! f({\bf x}) \over \partial \, y^2} + {\partial^{\,2} \! f({\bf x}) \over \partial \, z^2}$

### Tensor Notation

The Laplacian is written in tensor notation simply as $$f,_{ii}$$ where the two $$i$$ indices means that they are automatically summed from 1 to 3.

### Laplacian Example

Determine the Laplacian of $$f({\bf x}) = 2 x^3 y - z \sin(y)$$.

Start by calculating the gradient of $$f({\bf x})$$.

$\nabla f({\bf x}) = \left( 6 x^2 y, 2 x^3 - z \cos(y), - \sin(y) \right)$
And the divergence of the gradient (which is the Laplacian after all) is

$\nabla^2 \! f({\bf x}) = \nabla \cdot \nabla f({\bf x}) = 12 x y + z \sin(y)$

## Derivatives of Vector Products

Differentiation of vector products (dot, cross, and diadic) follow the same rules as differentiation of scalar products. For example, the derivative of a dot product is

${d \over dt} ( {\bf a} \cdot {\bf b} ) = {d \, {\bf a} \over dt} \cdot {\bf b} + {\bf a} \cdot {d \, {\bf b} \over dt}$
while the derivative of a cross product is

${d \over dt} ( {\bf a} \times {\bf b} ) = {d \, {\bf a} \over dt} \times {\bf b} + {\bf a} \times {d \, {\bf b} \over dt}$
and the derivative of a diadic product is

${d \over dt} ( {\bf a} \otimes {\bf b} ) = {d \, {\bf a} \over dt} \otimes {\bf b} + {\bf a} \otimes {d \, {\bf b} \over dt}$

### Dot Product Derivative Example

Suppose $${\bf a} = (5t, \sin t, e^t)$$ and $${\bf b} = (t^2, \sin t, 6t)$$, then

${\bf a} \cdot {\bf b} = 5 t^3 + \sin^2 t + 6t e^t$
and the derivative is

${d \over dt} ( {\bf a} \cdot {\bf b} ) = 15 t^2 + 2 \sin t \cos t + 6 (t + 1) e^t$
Applying the differentiation product rule gives the same result.

$\begin{eqnarray} {d \, {\bf a} \over dt} \cdot {\bf b} + {\bf a} \cdot {d \, {\bf b} \over dt} & = & (5, \cos t, e^t) \cdot (t^2, \sin t, 6t) + (5t, \sin t, e^t) \cdot (2t, \cos t, 6)\\ \\ & = & 15 t^2 + 2 \sin t \cos t + 6 (t + 1) e^t \end{eqnarray}$

## Textbooks

### Thank You

Thank you for visiting this webpage. Feel free to email me if you have questions.

Also, please consider visiting an advertiser above. Doing so helps to cover website hosting fees.

Bob McGinty
bmcginty@gmail.com

### Entire Website - $12.95 For$12.95, you receive two formatted PDFs (the first for 8.5" x 11" pages, the second for tablets) of the entire website.

Click here to see a sample page in each of the two formats.

 Email address to receive PDFs

### Math Chapter - $4.99 For$4.99, you receive two formatted PDFs (the first for 8.5" x 11" pages, the second for tablets) of the complete math chapter.

Click here to see a sample page in each of the two formats.

 Email address to receive PDFs

Copyright © 2012 by Bob McGinty