Vector Calculus: Gradient, Divergence, Curl & Theorems

Unit 5 — Vector Calculus (8 Hours)

Vector calculus deals with differentiation and integration of vector fields and provides the mathematical language for electromagnetism, fluid dynamics, heat transfer, and structural mechanics.

1. What Is Vector Calculus?

Vector calculus extends calculus to vector fields — functions that assign a vector to every point in space. It unifies differential and integral operations through powerful theorems that relate local properties (derivatives) to global properties (integrals over regions).

Why It Matters in Engineering

Electromagnetism: Maxwell's equations are written in vector calculus form
Fluid Dynamics: Velocity fields, vorticity, and flux calculations
Heat Transfer: Temperature gradients and heat flux
Structural Mechanics: Stress and strain fields
Gravitational Fields: Potential theory and field lines

2. Vector Differentiation

2.1 The Del Operator ( $\nabla$ )

The vector differential operator in Cartesian coordinates:

\nabla = \mathbf{i}\frac{\partial}{\partial x} + \mathbf{j}\frac{\partial}{\partial y} + \mathbf{k}\frac{\partial}{\partial z}

2.2 Gradient

Definition: For a scalar field $\phi(x, y, z)$ , the gradient is:

\nabla \phi = \mathbf{i}\frac{\partial \phi}{\partial x} + \mathbf{j}\frac{\partial \phi}{\partial y} + \mathbf{k}\frac{\partial \phi}{\partial z}

Properties:

$\nabla(\phi + \psi) = \nabla\phi + \nabla\psi$
$\nabla(\phi\psi) = \phi\nabla\psi + \psi\nabla\phi$

Geometrical Interpretation:

Points in the direction of maximum rate of increase of $\phi$
Magnitude $|\nabla\phi|$ equals the maximum rate of change
Perpendicular to level surfaces $\phi(x, y, z) = c$

Physical Interpretation:

If $\phi$ is temperature, $\nabla\phi$ is the temperature gradient (heat flows opposite to gradient)
If $\phi$ is electric potential, $\nabla\phi$ relates to electric field: $\mathbf{E} = -\nabla V$

2.3 Divergence

Definition: For a vector field $\mathbf{F} = F_1\mathbf{i} + F_2\mathbf{j} + F_3\mathbf{k}$ , the is:

\nabla \cdot \mathbf{F} = \frac{\partial F_1}{\partial x} + \frac{\partial F_2}{\partial y} + \frac{\partial F_3}{\partial z}

Physical Interpretation:

Measures the net outflow of flux from a point per unit volume
$\nabla \cdot \mathbf{F} > 0$ : Source (fluid emerging)
$\nabla \cdot \mathbf{F} < 0$ : Sink (fluid disappearing)
$\nabla \cdot \mathbf{F} = 0$ : Solenoidal (incompressible flow)

Properties:

$\nabla \cdot (\mathbf{F} + \mathbf{G}) = \nabla \cdot \mathbf{F} + \nabla \cdot \mathbf{G}$
$\nabla \cdot (ϕ F) = ϕ ($

2.4 Curl

Definition: The curl of a vector field $\mathbf{F}$ is:

\nabla \times \mathbf{F} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ F_1 & F_2 & F_3 \end{vmatrix}

Expanded:

\nabla \times \mathbf{F} = \mathbf{i}\left(\frac{\partial F_3}{\partial y} - \frac{\partial F_2}{\partial z}\right) - \mathbf{j}\left(\frac{\partial F_3}{\partial x} - \frac{\partial F_1}{\partial z}\right) + \mathbf{k}\left(\frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y}\right)

Physical Interpretation:

Measures the rotation or "swirl" of a vector field at a point
$\nabla \times \mathbf{F} = \mathbf{0}$ : Irrotational field
In fluid dynamics, curl of velocity = vorticity (twice the angular velocity)
In electromagnetism, $\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$ (Faraday's law)

Properties:

$\nabla \times (\mathbf{F} + \mathbf{G}) = \nabla \times \mathbf{F} + \nabla \times \mathbf{G}$
$\nabla \times (ϕ F)$

2.5 Directional Derivatives

Definition: The directional derivative of a scalar field $\phi$ in the direction of a unit vector $\hat{\mathbf{u}}$ is:

D_{\hat{\mathbf{u}}}\phi = \nabla\phi \cdot \hat{\mathbf{u}} = |\nabla\phi| \cos\theta

where $\theta$ is the angle between $\nabla\phi$ and $\hat{\mathbf{u}}$ .

Maximum Directional Derivative:

Occurs when $\hat{\mathbf{u}}$ is parallel to $\nabla\phi$ ( $\theta = 0$ )
Maximum value = $|\nabla\phi|$

Minimum Directional Derivative:

Occurs when $\hat{\mathbf{u}}$ is opposite to $\nabla\phi$ ( $\theta = \pi$ )
Minimum value = $-|\nabla\phi|$

3. Vector Integration

3.1 Line Integral

Definition: The line integral of a vector field $\mathbf{F}$ along a curve $C$ from point $A$ to point $B$ :

\int_C \mathbf{F} \cdot d\mathbf{r} = \int_C (F_1 \, dx + F_2 \, dy + F_3 \, dz)

where $d\mathbf{r} = dx\,\mathbf{i} + dy\,\mathbf{j} + dz\,\mathbf{k}$ .

Parametric Form: If $C$ is given by $\mathbf{r}(t)$ for $a \leq t \leq b$ :

\int_C \mathbf{F} \cdot d\mathbf{r} = \int_a^b \mathbf{F}(\mathbf{r}(t)) \cdot \frac{d\mathbf{r}}{dt} \, dt

Work Done: If $\mathbf{F}$ is a force field, the work done in moving a particle along $C$ is:

W = \int_C \mathbf{F} \cdot d\mathbf{r}

Circulation: For a closed curve $C$ :

\Gamma = \oint_C \mathbf{F} \cdot d\mathbf{r}

3.2 Surface Integral

Definition: The surface integral of a vector field $\mathbf{F}$ over a surface $S$ :

\iint_S \mathbf{F} \cdot d\mathbf{S} = \iint_S \mathbf{F} \cdot \hat{\mathbf{n}} \, dS

where $\hat{\mathbf{n}}$ is the unit normal vector to the surface.

Scalar Surface Integral:

\iint_S \phi \, dS

Flux: The surface integral $\iint_S \mathbf{F} \cdot \hat{\mathbf{n}} \, dS$ represents the flux of $\mathbf{F}$ through .

For a surface $z = f(x, y)$ projected onto region $R$ in the $xy$ -plane:

dS = \sqrt{1 + \left(\frac{\partial z}{\partial x}\right)^2 + \left(\frac{\partial z}{\partial y}\right)^2} \, dx \, dy

3.3 Volume Integral

Definition: The volume integral of a scalar field $\phi$ over a volume $V$ :

\iiint_V \phi \, dV

For a vector field $\mathbf{F}$ :

\iiint_V \mathbf{F} \, dV

4. Integral Theorems

4.1 Green's Theorem

Statement: Let $C$ be a positively oriented (counter-clockwise), piecewise smooth, simple closed curve in the plane enclosing a region $R$ . If $P(x, y)$ and $Q(x, y)$ have continuous partial derivatives on an open region containing , then:

\oint_C (P \, dx + Q \, dy) = \iint_R \left(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right) dA

Alternative Forms:

Area Formula:

\text{Area of } R = \frac{1}{2} \oint_C (x \, dy - y \, dx) = \oint_C x \, dy = -\oint_C y \, dx

Normal Form:

\oint_C \mathbf{F} \cdot \hat{\mathbf{n}} \, ds = \iint_R (\nabla \cdot \mathbf{F}) \, dA

4.2 Gauss's Divergence Theorem

Statement: Let $V$ be a solid region bounded by a closed surface $S$ with outward unit normal $\hat{\mathbf{n}}$ . If $\mathbf{F}$ has continuous partial derivatives in $V$ , then:

\iint_S \mathbf{F} \cdot \hat{\mathbf{n}} \, dS = \iiint_V (\nabla \cdot \mathbf{F}) \, dV

Physical Interpretation:

The total flux of $\mathbf{F}$ out of a closed surface equals the total divergence (source strength) inside the volume.
Foundation of conservation laws in physics.

Special Case: If $\mathbf{F} = \phi \nabla\psi$ (Green's first identity):

\iiint_V (\phi \nabla^2\psi + \nabla\phi \cdot \nabla\psi) \, dV = \iint_S \phi \frac{\partial\psi}{\partial n} \, dS

4.3 Stokes' Theorem

Statement: Let $S$ be an oriented piecewise smooth surface bounded by a simple, closed, piecewise smooth curve $C$ with positive orientation. If $\mathbf{F}$ has continuous partial derivatives on an open region containing $S$ , then:

\oint_C \mathbf{F} \cdot d\mathbf{r} = \iint_S (\nabla \times \mathbf{F}) \cdot \hat{\mathbf{n}} \, dS

Physical Interpretation:

The circulation of $\mathbf{F}$ around a closed loop equals the total curl (rotation) passing through any surface bounded by that loop.
If $\nabla \times \mathbf{F} = \mathbf{0}$ , the field is conservative and line integrals are path-independent.

Conservative Fields: $\mathbf{F}$ is conservative if and only if:

$\nabla \times \mathbf{F} = \mathbf{0}$ (irrotational)
$\oint_C \mathbf{F} \cdot d\mathbf{r} = 0$ for every closed curve

5. The Laplacian and Vector Identities

5.1 Laplacian Operator

For a scalar field $\phi$ :

\nabla^2 \phi = \nabla \cdot \nabla\phi = \frac{\partial^2\phi}{\partial x^2} + \frac{\partial^2\phi}{\partial y^2} + \frac{\partial^2\phi}{\partial z^2}

For a vector field $\mathbf{F}$ :

\nabla^2 \mathbf{F} = \nabla(\nabla \cdot \mathbf{F}) - \nabla \times (\nabla \times \mathbf{F})

5.2 Key Vector Identities

Identity	Formula
Gradient of sum	$\nabla(\phi + \psi) = \nabla\phi + \nabla\psi$
Divergence of curl	$\nabla \cdot (\nabla \times \mathbf{F}) = 0$

6. Coordinate System Formulas

6.1 Cylindrical Coordinates $(r, \theta, z)$

$x = r\cos\theta$ , $y = r\sin\theta$ , $z = z$

Gradient:

\nabla\phi = \frac{\partial\phi}{\partial r}\hat{\mathbf{r}} + \frac{1}{r}\frac{\partial\phi}{\partial\theta}\hat{\boldsymbol{\theta}} + \frac{\partial\phi}{\partial z}\hat{\mathbf{z}}

Divergence:

\nabla \cdot \mathbf{F} = \frac{1}{r}\frac{\partial(rF_r)}{\partial r} + \frac{1}{r}\frac{\partial F_\theta}{\partial\theta} + \frac{\partial F_z}{\partial z}

Curl:

\nabla \times \mathbf{F} = \frac{1}{r}\begin{vmatrix} \hat{\mathbf{r}} & r\hat{\boldsymbol{\theta}} & \hat{\mathbf{z}} \\ \frac{\partial}{\partial r} & \frac{\partial}{\partial\theta} & \frac{\partial}{\partial z} \\ F_r & rF_\theta & F_z \end{vmatrix}

Laplacian:

\nabla^2\phi = \frac{1}{r}\frac{\partial}{\partial r}\left(r\frac{\partial\phi}{\partial r}\right) + \frac{1}{r^2}\frac{\partial^2\phi}{\partial\theta^2} + \frac{\partial^2\phi}{\partial z^2}

6.2 Spherical Coordinates $(\rho, \phi, \theta)$

$x = \rho\sin\phi\cos\theta$ , $y = \rho\sin\phi\sin\theta$ , $z = \rho\cos\phi$

Gradient:

\nabla\psi = \frac{\partial\psi}{\partial\rho}\hat{\boldsymbol{\rho}} + \frac{1}{\rho}\frac{\partial\psi}{\partial\phi}\hat{\boldsymbol{\phi}} + \frac{1}{\rho\sin\phi}\frac{\partial\psi}{\partial\theta}\hat{\boldsymbol{\theta}}

Divergence:

\nabla \cdot \mathbf{F} = \frac{1}{\rho^2}\frac{\partial(\rho^2 F_\rho)}{\partial\rho} + \frac{1}{\rho\sin\phi}\frac{\partial(\sin\phi \, F_\phi)}{\partial\phi} + \frac{1}{\rho\sin\phi}\frac{\partial F_\theta}{\partial\theta}

Laplacian:

\nabla^2\psi = \frac{1}{\rho^2}\frac{\partial}{\partial\rho}\left(\rho^2\frac{\partial\psi}{\partial\rho}\right) + \frac{1}{\rho^2\sin\phi}\frac{\partial}{\partial\phi}\left(\sin\phi\frac{\partial\psi}{\partial\phi}\right) + \frac{1}{\rho^2\sin^2\phi}\frac{\partial^2\psi}{\partial\theta^2}

7. Engineering Applications

Domain	Application
Electromagnetism	Maxwell's equations: $\nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon_0}$ ,

8. Common Mistakes to Avoid

Mistake	Correction
Confusing $\nabla\phi$ (vector) with $\nabla \cdot \mathbf{F}$ (scalar)	Gradient of scalar = vector; Divergence of vector = scalar
Forgetting that curl of gradient is always zero	$\nabla \times (\nabla\phi) = \mathbf{0}$ — useful for checking conservative fields

9. Best Practices & Architecture

Always verify $\nabla \times \mathbf{F} = \mathbf{0}$ before claiming a field is conservative
Use potential functions $\phi$ where $\mathbf{F} = \nabla\phi$ to simplify line integral calculations
Choose the simplest surface for Stokes' theorem — any surface bounded by $C$ works

10. Practice Problems

Problem 1: Find the gradient of $\phi = x^2y + yz^2$ at the point $(1, 2, -1)$ .

Problem 2: Compute the divergence and curl of $\mathbf{F} = x^2y\,\mathbf{i} + yz^2\,\mathbf{j} + z^2x\,\mathbf{k}$ .

Problem 3: Find the directional derivative of $\phi = x^2 + y^2 + z^2$ at in the direction of the vector .

Problem 4: Evaluate $\displaystyle \int_C \mathbf{F} \cdot d\mathbf{r}$ where $\mathbf{F} = (x + y)\mathbf{i} + (y - x)\mathbf{j}$ and is the unit circle traversed counter-clockwise using Green's theorem.

Problem 5: Verify the divergence theorem for $\mathbf{F} = x\mathbf{i} + y\mathbf{j} + z\mathbf{k}$ over the unit sphere $x^2 + y^2 + z^2 = 1$ .

Problem 6: Verify Stokes' theorem for $\mathbf{F} = -y\mathbf{i} + x\mathbf{j} + z\mathbf{k}$ where $S$ is the upper hemisphere $x^2 + y^2 + z^2 = 1$ , .

Problem 7: Show that $\mathbf{F} = (2xy + z^3)\mathbf{i} + x^2\mathbf{j} + 3xz^2\mathbf{k}$ is conservative and find its scalar potential.

Problem 8: Evaluate $\displaystyle \iint_S \mathbf{F} \cdot \hat{\mathbf{n}} \, dS$ where $\mathbf{F} = z\mathbf{i} + x\mathbf{j} + y\mathbf{k}$ and is the surface of the cube bounded by ; ; .

11. Mini Project

Title: Electromagnetic Field Analysis Using Vector Calculus Theorems

Objective: Apply gradient, divergence, curl, and the integral theorems to analyze an electrostatic field configuration.

Given:

Electric potential field: $V(x, y, z) = \frac{k}{\sqrt{x^2 + y^2 + z^2}}$ where is constant

Tasks:

Compute $\nabla V$ and verify $\mathbf{E} = -\nabla V$ in spherical coordinates
Calculate $\nabla \cdot \mathbf{E}$ everywhere except the origin; interpret the result
Use the divergence theorem to find the total flux of $\mathbf{E}$ through a sphere of radius centered at the origin

Deliverable: A computational notebook (Python with SymPy/MATLAB) that symbolically computes all vector operations, visualizes the field using quiver plots and streamlines, validates each theorem numerically, and provides physical interpretations of divergence and curl in the electromagnetic context.

Quick Reference Formula Sheet

Concept	Formula
Gradient	$\nabla\phi = \mathbf{i}\frac{\partial\phi}{\partial x} + \mathbf{j}\frac{\partial\phi}{\partial y} + \mathbf{k}\frac{\partial\phi}{\partial z}$

End of Unit 5 — Vector Calculus

Unit 5 — Vector Calculus (8 Hours)

1. What Is Vector Calculus?

Why It Matters in Engineering

2. Vector Differentiation

2.1 The Del Operator (∇\nabla∇)

2.2 Gradient

2.3 Divergence

2.4 Curl

2.5 Directional Derivatives

3. Vector Integration

3.1 Line Integral

3.2 Surface Integral

3.3 Volume Integral

4. Integral Theorems

4.1 Green's Theorem

4.2 Gauss's Divergence Theorem

4.3 Stokes' Theorem

5. The Laplacian and Vector Identities

5.1 Laplacian Operator

5.2 Key Vector Identities

6. Coordinate System Formulas

6.1 Cylindrical Coordinates (r,θ,z)(r, \theta, z)(r,θ,z)

6.2 Spherical Coordinates (ρ,ϕ,θ)(\rho, \phi, \theta)(ρ,ϕ,θ)

7. Engineering Applications

8. Common Mistakes to Avoid

9. Best Practices & Architecture

10. Practice Problems

11. Mini Project

Quick Reference Formula Sheet

2.1 The Del Operator ( $\nabla$ )

6.1 Cylindrical Coordinates $(r, \theta, z)$

6.2 Spherical Coordinates $(\rho, \phi, \theta)$