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=∂Fx∂x+∂Fy∂y+∂Fz∂z\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 equation∇2ϕ=0\nabla^2\phi = 0
Heat equation∇2ϕ=1D∂ϕ∂t\nabla^2\phi = \frac{1}{D}\frac{\partial\phi}{\partial t}
Wave equation∇2ϕ=1c2∂2ϕ∂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: ∫SF⋅dS=∫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: 0≤x≤a0 \leq x \leq a, 0≤y≤b0 \leq y \leq b, 0≤z≤c0 \leq z \leq c.

Faces perpendicular to xx:

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

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

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

Sum over all six faces:

Φ=∫V(∂Fx∂x+∂Fy∂y+∂Fz∂z)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: ∫x dV=∫y dV=∫z dV=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