Divergence ($\nabla \cdot$)

For $F=(F_1,F_2,F_3)$ in Cartesian coordinates,

For a continuously differentiable vector field, the divergence is the limiting outward flux per unit volume through a small closed surface. Opposite faces of a small rectangular box almost cancel; the remaining differences in the three components give the sum of partial derivatives above. The divergence theorem extends this local balance to a whole region.

For $F=(x,y,0)$, the divergence is $1+1+0=2$. In a planar unit square, the right edge contributes outward flux $1$ and the top contributes $1$, while the left and bottom contribute zero. The total $2$ equals the integral of the planar divergence over area $1$.

If $F$ is a velocity field, positive divergence means local volume expansion, not automatically creation of mass. The continuity equation gives $D\rho /Dt=-\rho \nabla \cdot F$, where $D/Dt$ is the material derivative. For this field, a uniform initial density evolves as $\rho (t)=\rho (0)e^{-2t}$ while material volumes grow as $e^{2t}$. Their product stays constant.

Horizontal flow convergence is the negative horizontal divergence of a velocity field. Positive convergence reduces the horizontal area of a material patch; for constant-depth linear shallow water equations, the surface equation is ${\eta }_t=HC$, so convergence raises the free surface and divergence lowers it. This is a fluid-flow notion, distinct from convergence of a numerical sequence.

The divergence defined through weak derivatives: pair a locally integrable vector field $J$ with a smooth compactly supported test function $\phi$ using the displayed integration-by-parts identity. A piecewise smooth field can have a surface contribution where its normal component jumps. Thus zero classical divergence on each side does not imply zero weak divergence across the interface.

The divergence with respect to velocity coordinates, with position held fixed. In $d$ dimensions, $\mathbf{A}(\mathbf{x},\mathbf{v})=\mathbf{B}(\mathbf{x})-a(\mathbf{x})\mathbf{v}$ has ${\nabla }_{\mathbf{v}}\cdot \mathbf{A}=-da$. A kinetic conservation law must retain this term for velocity-dependent forces.

For the Levi-Civita connection of a Riemannian metric, the divergence of $X$ is ${\mathrm{div}}_{g}X=\mathrm{tr}(Y\mapsto {\nabla }_{Y}X)$. In coordinates it equals $(\det g{)}^{-1/2}{\partial }_i((\det g{)}^{1/2}{X}^i)$. On an oriented manifold it satisfies ${\mathcal{L}}_X{\omega }_g=({\mathrm{div}}_{g}X){\omega }_g$. The codifferential on a one-form satisfies $\delta \theta =-{\mathrm{div}}_g({\theta }^{♯})$, where $♯$ is the musical isomorphism; this follows by integration by parts or the Hodge star formula.

For a vector field directed radially from a fixed axis and depending only on distance $\rho$ from that axis, cylindrical coordinates give $\nabla \cdot [f(\rho ){\mathbf{e}}_{\rho }]=f\prime +f/\rho$. Zero divergence gives $f=A/\rho$ on any positive radial interval. Unless $A=0$, this field is singular on the axis, so applications of the divergence theorem must use a domain avoiding it or account separately for that singularity.

For a locally integrable vector field, this formula defines its divergence as a distribution. For a piecewise smooth field, a jump of its normal component across an interface adds a surface measure to the smooth divergences. A continuous normal component eliminates that measure, as follows by integration by parts on the two sides.

The inverse-square radial vector field has zero ordinary divergence at every nonzero radius but flux $4\pi$ through every centred sphere. It is locally integrable in three dimensions. For a smooth compactly supported test function $\eta$, integrate by parts outside a ball of radius $ϵ$. The inward-oriented inner boundary gives $-\int \mathbf{F}\cdot \nabla \eta ={\int }_{S^2}\eta (ϵ\hat{\mathbf{r}})d\Omega$, which tends to $4\pi \eta (0)$. This proves the displayed distributional divergence, with the three-dimensional Dirac delta function representing the point source.

For $\mathbf{u}=u_r{\mathbf{e}}_r+u_{\theta }{\mathbf{e}}_{\theta }$,

For a scalar field $\varphi$ and vector field $F$,

A total divergence has the form $\nabla \cdot \mathbf{F}$. The divergence theorem converts its volume integral into a boundary flux, which vanishes when the field has suitable decay or boundary conditions.

Adding $dF(q,t)/dt$ to a mechanical Lagrangian changes its action by $F(q_f,t_f)-F(q_i,t_i)$. Fixed endpoint variations therefore leave its bulk Euler-Lagrange equations unchanged. For fields, adding ${\partial }_{\mu }{K}^{\mu }$ changes the action by a boundary flux. The same conclusion holds for compactly supported variations or appropriate fixed boundary data. If $K$ depends on field derivatives, those boundary data must also control their variations. Such terms can still change canonical momenta, boundary observables, and representatives of Noether currents; equivalence of bulk equations does not mean that every boundary-value problem is identical.

For a twice continuously differentiable vector field, commuting partial derivatives and antisymmetry of the Levi-Civita symbol give