Skip to content
Part IA Michaelmas, Lent, Easter Term

Divergence and Gauss's Theorem

The Divergence Operator

For a vector field F=(Fx,Fy,Fz)\mathbf{F} = (F_x, F_y, F_z):

F=Fxx+Fyy+Fzz\nabla \cdot \mathbf{F} = \frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z}

The result is a scalar field.


Physical Interpretation

F\nabla \cdot \mathbf{F} measures expansion: the rate at which flow “spreads out” from each point.

Vector field with positive divergence (source, flow outward) and negative divergence (sink, flow inward)

  • F>0\nabla \cdot \mathbf{F} > 0: Source (flow diverges)
  • F<0\nabla \cdot \mathbf{F} < 0: Sink (flow converges)
  • F=0\nabla \cdot \mathbf{F} = 0: Incompressible (volume preserved)

The Laplacian

Applying divergence to gradient of a scalar field:

(ϕ)=2ϕ=2ϕx2+2ϕy2+2ϕz2\nabla \cdot (\nabla \phi) = \nabla^2 \phi = \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2}

The Laplacian 2\nabla^2 is fundamental in PDEs:

PDEEquation
Laplace’s equation2ϕ=0\nabla^2\phi = 0
Heat equation2ϕ=1Dϕt\nabla^2\phi = \frac{1}{D}\frac{\partial\phi}{\partial t}
Wave equation2ϕ=1c22ϕt2\nabla^2\phi = \frac{1}{c^2}\frac{\partial^2\phi}{\partial t^2}

Example

For F=(x2,y2,z2)\mathbf{F} = (x^2, y^2, z^2):

F=2x+2y+2z=2(x+y+z)\nabla \cdot \mathbf{F} = 2x + 2y + 2z = 2(x + y + z)

At origin: F=0\nabla \cdot \mathbf{F} = 0.

On axis (x,0,0)(x, 0, 0): F=2x\nabla \cdot \mathbf{F} = 2x, increases with distance.


Gauss’s (Divergence) Theorem

For a closed surface SS enclosing volume VV: SFdS=V(F)dV\int_S \mathbf{F} \cdot d\mathbf{S} = \int_V (\nabla \cdot \mathbf{F})\,dV

The flux through a closed surface equals the volume integral of divergence inside.

Closed surface S enclosing volume V, with outward normal vectors and flux arrows through the surface


Proof Sketch (Cuboid)

Consider a cuboid VV: 0xa0 \leq x \leq a, 0yb0 \leq y \leq b, 0zc0 \leq z \leq c.

Faces perpendicular to xx:

  • x=0x = 0: flux =Fx(0,y,z)dydz= -\int F_x(0, y, z)\,dy\,dz (normal points inward)
  • x=ax = a: flux =Fx(a,y,z)dydz= \int F_x(a, y, z)\,dy\,dz

Combined: [Fx(a,y,z)Fx(0,y,z)]dydz\int [F_x(a, y, z) - F_x(0, y, z)]\,dy\,dz

By fundamental theorem of calculus: =0aFxxdxdydz= \int_0^a \frac{\partial F_x}{\partial x} dx\,dy\,dz

Sum over all six faces:

Φ=V(Fxx+Fyy+Fzz)dV\Phi = \int_V \left(\frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z}\right)dV


Proof Extension

To generalise from cuboid to any shape:

  1. Subdivide volume into small cuboids
  2. Apply theorem to each cuboid
  3. Sum: internal surfaces cancel (opposite normals)
  4. Result: theorem holds for whole volume

Requirement: Surface is smooth enough for subdivision.


Applications

Flux Calculation

Often easier to compute volume integral than surface integral.

Example: Flux of F=(x2,y2,z2)\mathbf{F} = (x^2, y^2, z^2) through unit sphere.

F=2(x+y+z)\nabla \cdot \mathbf{F} = 2(x + y + z)

Φ=V2(x+y+z)dV\Phi = \int_V 2(x + y + z)\,dV

By symmetry: xdV=ydV=zdV=0\int x\,dV = \int y\,dV = \int z\,dV = 0.

So Φ=0\Phi = 0.

Vector Area of Closed Surface

SdS=0\oint_S d\mathbf{S} = \mathbf{0}

(Apply divergence theorem to constant vector field.)


Solenoidal Fields

A vector field with F=0\nabla \cdot \mathbf{F} = 0 everywhere is solenoidal.

From divergence theorem: flux through any closed surface is zero.

Examples:

  • Magnetic field B\mathbf{B}: B=0\nabla \cdot \mathbf{B} = 0 (no magnetic monopoles)
  • Incompressible fluid velocity

Why These Theorems Matter

Gauss’s theorem converts between:

  • Surface integrals (flux)
  • Volume integrals (divergence)

This is powerful because:

  • Choose whichever is easier to evaluate
  • Derive conservation laws in fluid mechanics
  • Apply to Maxwell’s equations in electromagnetism