Skip to content

Mathematical Foundations

For a general vector field,

F=Fxax+Fyay+Fzaz.\mathbf{F}=F_x\mathbf{a}_x+F_y\mathbf{a}_y+F_z\mathbf{a}_z.

Use the coordinate system that matches the symmetry of the source.

SystemCoordinatesBest forDifferential length
Cartesian(x,y,z)(x,y,z)planes, rectangular boundariesdl=dxax+dyay+dzazd\mathbf{l}=dx\mathbf{a}_x+dy\mathbf{a}_y+dz\mathbf{a}_z
Cylindrical(ρ,ϕ,z)(\rho,\phi,z)lines, cylinders, coaxial structuresdl=dρaρ+ρdϕaϕ+dzazd\mathbf{l}=d\rho\mathbf{a}_\rho+\rho d\phi\mathbf{a}_\phi+dz\mathbf{a}_z
Spherical(r,θ,ϕ)(r,\theta,\phi)points, spheres, radial fieldsdl=drar+rdθaθ+rsin⁡θdϕaϕd\mathbf{l}=dr\mathbf{a}_r+r d\theta\mathbf{a}_\theta+r\sin\theta d\phi\mathbf{a}_\phi

Coordinate systems chosen to match source symmetry.

Coordinate systems chosen to match source symmetry.

Useful surface elements are

dSρ=ρ dϕ dz,dSr=r2sin⁡θ dθ dϕ.dS_\rho=\rho\,d\phi\,dz,\qquad dS_r=r^2\sin\theta\,d\theta\,d\phi.

Useful volume elements are

dV=dx dy dz,dV=ρ dρ dϕ dz,dV=r2sin⁡θ dr dθ dϕ.dV=dx\,dy\,dz, \qquad dV=\rho\,d\rho\,d\phi\,dz, \qquad dV=r^2\sin\theta\,dr\,d\theta\,d\phi.

In Cartesian coordinates,

∇=ax∂∂x+ay∂∂y+az∂∂z.\nabla=\mathbf{a}_x\frac{\partial}{\partial x} +\mathbf{a}_y\frac{\partial}{\partial y} +\mathbf{a}_z\frac{\partial}{\partial z}. ∇V=∂V∂xax+∂V∂yay+∂V∂zaz.\nabla V= \frac{\partial V}{\partial x}\mathbf{a}_x+ \frac{\partial V}{\partial y}\mathbf{a}_y+ \frac{\partial V}{\partial z}\mathbf{a}_z.

The gradient is normal to a constant- VV surface. In electrostatics,

E=−∇V,\mathbf E=-\nabla V,

so the electric field points in the direction of the steepest decrease in potential.

∇⋅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}.

Equivalently, divergence is the limiting outward flux per unit volume:

∇⋅F=lim⁡ΔV→01ΔV∮ΔSF⋅dS.\nabla\cdot\mathbf F =\lim_{\Delta V\to0} \frac{1}{\Delta V}\oint_{\Delta S}\mathbf F\cdot d\mathbf S.

Positive divergence indicates a local source and negative divergence a local sink. In electromagnetics, ∇⋅D=ρv\nabla\cdot\mathbf D=\rho_v identifies volume charge density as the source density of electric flux.

∇×F=∣axayaz∂∂x∂∂y∂∂zFxFyFz∣.\nabla\times\mathbf{F}= \begin{vmatrix} \mathbf{a}_x & \mathbf{a}_y & \mathbf{a}_z \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ F_x & F_y & F_z \end{vmatrix}.

Equivalently,

∇×F=(∂Fz∂y−∂Fy∂z)ax+(∂Fx∂z−∂Fz∂x)ay+(∂Fy∂x−∂Fx∂y)az.\nabla\times\mathbf{F}= \left(\frac{\partial F_z}{\partial y}-\frac{\partial F_y}{\partial z}\right)\mathbf{a}_x+ \left(\frac{\partial F_x}{\partial z}-\frac{\partial F_z}{\partial x}\right)\mathbf{a}_y+ \left(\frac{\partial F_y}{\partial x}-\frac{\partial F_x}{\partial y}\right)\mathbf{a}_z.

The component normal to a small surface is the limiting circulation per unit area:

(∇×F)⋅an=lim⁡ΔS→01ΔS∮ΔCF⋅dl.(\nabla\times\mathbf F)\cdot\mathbf a_n =\lim_{\Delta S\to0} \frac{1}{\Delta S}\oint_{\Delta C}\mathbf F\cdot d\mathbf l.

A small paddle wheel placed in the field therefore rotates only where the curl is nonzero. An electrostatic field is irrotational, ∇×E=0\nabla\times\mathbf E=0, whereas time-varying fields obey

∇×E=−∂B∂t,∇×H=J+∂D∂t.\nabla\times\mathbf E=-\frac{\partial\mathbf B}{\partial t}, \qquad \nabla\times\mathbf H=\mathbf J+\frac{\partial\mathbf D}{\partial t}.

Curl as circulation around a small cell.

Curl as circulation around a small cell.

IntegralFormMeaning
Line integral∫CF⋅dl\int_C \mathbf{F}\cdot d\mathbf{l}field component along a path
Flux integral∫SF⋅dS\int_S \mathbf{F}\cdot d\mathbf{S}field crossing a surface
Volume integral∫Vf dV\int_V f\,dVtotal amount inside a volume

For a closed curve or surface, use ∮\oint.

∮SF⋅dS=∫V∇⋅F dV.\oint_S \mathbf{F}\cdot d\mathbf{S} =\int_V \nabla\cdot\mathbf{F}\,dV.

To see why, consider a differential rectangular volume dV=dx dy dzdV=dx\,dy\,dz. The net flux through its two faces normal to the xx-axis is

dΨx=[Fx(x+dx)−Fx(x)]dy dz≃∂Fx∂x dV.\begin{aligned} d\Psi_x & =\left[F_x(x+dx)-F_x(x)\right]dy\,dz \\ & \simeq \frac{\partial F_x}{\partial x}\,dV. \end{aligned}

Adding the corresponding yy- and zz-face contributions gives dΨ=(∇⋅F)dVd\Psi=(\nabla\cdot\mathbf F)dV. When a finite volume is divided into small cells, fluxes through common internal faces cancel in opposite pairs. Only the outer boundary remains, and the limiting sum gives the divergence theorem.

Applied to Gauss’s law, it converts a finite-region statement into a local field equation:

∮SD⋅dS=∫Vρv dV=∫V∇⋅D dV.\oint_S\mathbf D\cdot d\mathbf S =\int_V\rho_v\,dV =\int_V\nabla\cdot\mathbf D\,dV.

Because this equality holds for every volume,

∇⋅D=ρv.\boxed{\nabla\cdot\mathbf D=\rho_v}. ∮CF⋅dl=∫S(∇×F)⋅dS.\oint_C \mathbf{F}\cdot d\mathbf{l} =\int_S (\nabla\times\mathbf{F})\cdot d\mathbf{S}.

The contour direction and surface normal follow the right-hand rule. For a fixed contour, applying the theorem to Faraday’s law gives

∮CE⋅dl=−ddt∫SB⋅dS,∫S(∇×E)⋅dS=−∫S∂B∂t⋅dS.\begin{aligned} \oint_C\mathbf E\cdot d\mathbf l & =-\frac{d}{dt}\int_S\mathbf B\cdot d\mathbf S, \\ \int_S(\nabla\times\mathbf E)\cdot d\mathbf S & =-\int_S\frac{\partial\mathbf B}{\partial t}\cdot d\mathbf S. \end{aligned}

Since the surface is arbitrary,

∇×E=−∂B∂t.\boxed{\nabla\times\mathbf E=-\frac{\partial\mathbf B}{\partial t}}.

The same conversion applied to the Ampere-Maxwell law yields

∇×H=J+∂D∂t.\boxed{\nabla\times\mathbf H=\mathbf J+\frac{\partial\mathbf D}{\partial t}}.

Stokes theorem orientation for a boundary curve and surface normal.

Stokes theorem orientation for a boundary curve and surface normal.

Integral quantityDifferential operator
outward flux through a closed surfacedivergence
circulation around a closed curvecurl

The usual conversion path is

integral law⟷divergence theorem or Stokes’ theoremdifferential law.\mathrm{integral\ law}\quad \overset{\text{divergence theorem or Stokes' theorem}}{\longleftrightarrow} \quad \text{differential law}.