Solenoidal vector field

A solenoidal vector field has zero divergence. Magnetic fields are solenoidal because $\nabla \cdot \mathbf{B}=0$.

For a solenoidal vector field $J$ with compact support and suitable integrability, apply the divergence theorem to $x_i{x}_{j}J$. It gives zero boundary term and hence the displayed antisymmetry. Equivalently, for any constant vector $C$, the vector triple product identity yields $\int (C\cdot r)JdV=\frac{1}{2}C\times \int (J\times r)dV$. A piecewise smooth field must have no normal jump at its support boundary; otherwise its weak divergence contains a surface contribution and the argument fails.

Suppose two oriented surfaces have the same oriented boundary and $S_1\cup (-S_2)$ bounds a volume inside a domain where the vector field $\mathbf{F}$ is smooth and solenoidal. The divergence theorem gives zero flux through their closed union, proving the displayed equality. If a global vector potential $\mathbf{F}=\nabla \times \mathbf{G}$ exists on the relevant surfaces, the Stokes theorem proves the equality directly from their common boundary. The domain condition matters: a singularity inside the joining volume can obstruct the divergence argument.

A spherical poloidal-toroidal decomposition represents a solenoidal vector field with zero flux through each sphere by scalar potentials $S,T$. The second term is tangent to concentric spheres; the first carries the radial component. Requiring zero angular means for the potentials fixes their otherwise irrelevant radial additions. In a plane layer, the analogous fluctuating components are $\nabla \times \nabla \times (\varphi \hat{\mathbf{z}})$ and $\nabla \times (\psi \hat{\mathbf{z}})$, with any horizontally uniform components handled separately. These constructions enforce the divergence-free condition because the divergence of a curl vanishes.

// Target: fluid-mechanics.bigb

For an isolated solenoidal magnetic field, zero flux through every sphere removes the degree-zero spherical harmonic of $P=\mathbf{B}\cdot \mathbf{x}$. In a poloidal-toroidal decomposition, let ${\mathcal{L}}^2=-{\Delta }_{S^2}$. The spherical poloidal component obeysIf $P={\sum }_{l\ge 1,m}{p}_{lm}(r)Y_{lm}$, put $D_{lm}={\int }_0^{\infty }[r^2|p_{lm}^{\prime }{|}^2+l(l+1)|p_{lm}{|}^2]dr$. The two sides of the bound are respectively $\sum D_{lm}/[l(l+1)]$ and $\frac{1}{2}\sum D_{lm}$. Since $l(l+1)\ge 2$, the bound follows, with equality for degree-one radial scalars. Together with the radial-flow dynamo energy bound it gives $(q/\eta {)}^2\ge 2E_P/E_{\mathrm{total}}$ for a marginal or growing dynamo field.

// Target: fluid-mechanics.bigb