Solenoidal Fields and Zero Divergence

Zero Divergence and Source-Free Flow

A vector field can describe a velocity, a fluid flux, or another quantity moving through space. Its divergence measures local outward flow: positive divergence indicates net outward flow from a small region, while negative divergence indicates net inward flow. When the divergence vanishes, outward and inward contributions balance to first order locally.

A field with this property is called solenoidal. The condition is local and must be checked throughout the field’s domain. In particular, a formula that has zero divergence wherever it is defined may still have a singular point excluded from its domain. That distinction will matter for radial fields.

Definition: Let \(\mathbf{F}=(P,Q,R)\) be a continuously differentiable vector field on an open set \(D\subseteq\mathbb{R}^3\). Its divergence is \[ \nabla\cdot\mathbf{F} = \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z}. \] The field is solenoidal on \(D\) if \(\nabla\cdot\mathbf{F}=0\) at every point of \(D\).

Continuous differentiability is a convenient sufficient condition for the standard pointwise divergence formulas used here. The domain condition is essential: the equation must hold at every point where the field is being claimed to be solenoidal.

Why Divergence Measures Local Source Strength

The local interpretation follows by examining flux through a small rectangular box. Let a box centered at \((x,y,z)\) have side lengths \(\Delta x,\Delta y,\Delta z\), and suppose \(\mathbf{F}=(P,Q,R)\) is continuously differentiable near the box. On the two faces perpendicular to the \(x\)-axis, the outward flux is, to first order, the area \(\Delta y\Delta z\) times the difference between the \(x\)-components at the right and left faces:

\bigl(P(x+\tfrac{\Delta x}{2},y,z)-P(x-\tfrac{\Delta x}{2},y,z)\bigr)\Delta y\Delta z \approx \frac{\partial P}{\partial x}\Delta x\Delta y\Delta z.

The corresponding first-order contributions from the other pairs of faces are \(\frac{\partial Q}{\partial y}\Delta x\Delta y\Delta z\) and \(\frac{\partial R}{\partial z}\Delta x\Delta y\Delta z\). Adding them and dividing by the box volume gives the local net outward flux per unit volume, in the limit as the box shrinks:

\lim_{\Delta x,\Delta y,\Delta z\to 0} \frac{\text{net outward flux}}{\Delta x\Delta y\Delta z} = \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z} = \nabla\cdot\mathbf{F}.

This is a local, first-order interpretation. Zero divergence means the net outward flux per unit volume tends to zero at each point; it does not say that the field itself is zero. A field can circulate, or have nonzero components pointing in different directions, while having no local source strength.

Testing Linear Fields

For a linear field, divergence is particularly easy to identify. If \(\mathbf{F}(\mathbf{x})=\mathbf{A}\mathbf{x}\), where \(\mathbf{x}=(x,y,z)^{\top}\) and \(\mathbf{A}\) is a constant \(3\times3\) matrix, then the diagonal entries of \(\mathbf{A}\) are the component derivatives that appear in the divergence. Thus the divergence is the trace of the matrix: the sum of its diagonal entries. Off-diagonal terms contribute no divergence because they differentiate a component with respect to a different coordinate.

Worked Example: A Linear Field with Zero Divergence

Consider \(\mathbf{F}(x,y,z)=(2x-y,\ 3y+z,\ -5z+4x)\). Its divergence is the sum of the three indicated component derivatives:

\nabla\cdot\mathbf{F} = \frac{\partial}{\partial x}(2x-y) + \frac{\partial}{\partial y}(3y+z) + \frac{\partial}{\partial z}(-5z+4x) = 2+3-5 = 0.

Therefore \(\mathbf{F}\) is solenoidal on all of \(\mathbb{R}^3\). The terms \(-y\), \(z\), and \(4x\) do not affect this calculation: each is independent of the coordinate with respect to which its component is differentiated.

Check: the three diagonal derivatives are \(2\), \(3\), and \(-5\), whose sum is \(0\). Since the field is polynomial, its derivatives exist and are continuous everywhere, so the zero-divergence conclusion holds throughout its domain.

This calculation also gives a way to construct linear examples: choose diagonal coefficients whose sum is zero, while the off-diagonal coefficients may be chosen freely.

Worked Example: Constructing a Linear Solenoidal Field

Construct a linear field of the form \(\mathbf{F}=(ax+y,\ -2x+by,\ 3z)\) that is solenoidal, and find \(b\) when \(a=4\). The divergence condition is:

\nabla\cdot\mathbf{F} = a+b+3 = 0.

With \(a=4\), this gives \(4+b+3=0\), so \(b=-7\). One resulting field is \(\mathbf{F}=(4x+y,\ -2x-7y,\ 3z)\). Direct verification gives:

\nabla\cdot\mathbf{F} = \frac{\partial}{\partial x}(4x+y) + \frac{\partial}{\partial y}(-2x-7y) + \frac{\partial}{\partial z}(3z) = 4-7+3 = 0.

Thus this field is solenoidal everywhere. Substituting the chosen coefficients into the constraint gives \(a+b+3=4-7+3=0\), confirming that the construction satisfies the required condition.

Rotational Fields Can Also Be Solenoidal

A field that circulates around an axis need not have a source or sink. For example, the horizontal vectors in the field \((-y,x,0)\) point tangentially around the \(z\)-axis. Its components vary with position, but each component is independent of the coordinate used in its divergence derivative. Consequently, circulation and divergence measure different aspects of a field.

Worked Example: A Rotational Field

Let \(\mathbf{F}(x,y,z)=(-y,x,0)\). Computing divergence component by component gives:

\nabla\cdot\mathbf{F} = \frac{\partial(-y)}{\partial x} + \frac{\partial x}{\partial y} + \frac{\partial 0}{\partial z} = 0+0+0 = 0.

The field is therefore solenoidal throughout \(\mathbb{R}^3\), including on the \(z\)-axis, where its value is simply \(\mathbf{0}\). Its vectors are nonzero away from that axis, so zero divergence does not mean a zero field.

Check: differentiating the displayed components with respect to \(x\), \(y\), and \(z\), respectively, again yields \(0\), \(0\), and \(0\). The sum is zero at every point.

Radial Fields: A Condition on the Magnitude

A radial field in three dimensions has the form \(\mathbf{F}(\mathbf{x})=f(r)\mathbf{x}\), where \(r=\sqrt{x^2+y^2+z^2}\), and \(f\) depends only on the distance from the origin. Away from the origin, its divergence can be derived directly using the product rule. Since \(\nabla\cdot\mathbf{x}=3\) and \(\nabla f(r)=f'(r)\mathbf{x}/r\) for \(r>0\), we obtain:

\nabla\cdot\bigl(f(r)\mathbf{x}\bigr) = \nabla f(r)\cdot\mathbf{x} + f(r)\nabla\cdot\mathbf{x} = r f'(r)+3f(r).

Thus, on a region where \(r>0\), zero divergence requires \(r f'(r)+3f(r)=0\). Solving this condition reveals the radial magnitude needed to balance the increasing area of spheres around the origin.

Result: For a continuously differentiable radial field \(\mathbf{F}(\mathbf{x})=f(r)\mathbf{x}\) on a region with \(r>0\), the field is solenoidal precisely where \[ r f'(r)+3f(r)=0. \] Equivalently, \(\frac{d}{dr}(r^3f(r))=0\), so \(r^3f(r)\) is constant on each connected radial interval.

Worked Example: A Radial Field with Zero Divergence Away from the Origin

Take \(f(r)=1/r^3\) for \(r>0\). Then \(\mathbf{F}(\mathbf{x})=\mathbf{x}/r^3\), which is defined on \(\mathbb{R}^3\setminus\{\mathbf{0}\}\). Differentiate \(f\): \(f'(r)=-3/r^4\). Substituting into the radial divergence formula gives:

\nabla\cdot\mathbf{F} = r\left(-\frac{3}{r^4}\right) + 3\left(\frac{1}{r^3}\right) = -\frac{3}{r^3}+\frac{3}{r^3} = 0 \qquad (r>0).

Hence this field is solenoidal on its stated domain, which excludes the origin. The exclusion is not optional: the formula is undefined at \(r=0\), so the field cannot be claimed to have zero divergence there.

Check: the two terms in the radial formula cancel for every \(r>0\). This verifies zero divergence at every point in the domain, not at the omitted origin.

How Coordinate Formulas Fit the Same Test

The divergence formulas in cylindrical and spherical coordinates express the same local quantity using the coordinate system suited to the field. They are useful when the components have those symmetries, but their coordinate singularities require care. In cylindrical coordinates \((\rho,\phi,z)\), for a field with components \((F_\rho,F_\phi,F_z)\), the formula is:

\nabla\cdot\mathbf{F} = \frac{1}{\rho}\frac{\partial}{\partial \rho}(\rho F_\rho) + \frac{1}{\rho}\frac{\partial F_\phi}{\partial \phi} + \frac{\partial F_z}{\partial z}, \qquad \rho>0.

For a radial field with \(F_\rho=1/\rho^2\), \(F_\phi=0\), and \(F_z=0\), this formula yields \(\frac{1}{\rho}\frac{\partial}{\partial \rho}(1/\rho)= -1/\rho^3\), not zero. This is not the same radial field as \(\mathbf{x}/r^3\) in three dimensions: the latter depends on spherical distance \(r\) and has spherical radial component \(1/r^2\). Keeping the geometry and the meaning of the radial coordinate straight is essential.

In spherical coordinates \((r,\theta,\phi)\), with radial, polar, and azimuthal components \((F_r,F_\theta,F_\phi)\), the divergence is:

\nabla\cdot\mathbf{F} = \frac{1}{r^2}\frac{\partial}{\partial r}(r^2F_r) + \frac{1}{r\sin\theta}\frac{\partial}{\partial\theta}(\sin\theta F_\theta) + \frac{1}{r\sin\theta}\frac{\partial F_\phi}{\partial\phi}.

For \(\mathbf{x}/r^3\), the spherical components are \(F_r=1/r^2\), \(F_\theta=0\), and \(F_\phi=0\). The radial term becomes \(\frac{1}{r^2}\frac{\partial}{\partial r}(1)=0\); the other terms vanish as well. This confirms the result for \(r>0\). The formula does not include the origin, where spherical coordinates and the field itself are singular.

A Useful Pitfall: Ignoring the Domain

A zero result from a divergence calculation applies only where the field and the calculation are defined. It does not erase a singularity from the field’s domain. For example, it would be incorrect to write that \(\mathbf{x}/r^3\) is solenoidal on all of \(\mathbb{R}^3\): substitution into the radial formula proves zero divergence only for \(r>0\). The field is undefined at \(\mathbf{0}\).

Another common error is to conclude that a field must be zero if its divergence is zero. The rotational example disproves this: \((-y,x,0)\) is generally nonzero, but its divergence vanishes. Divergence tests local net outward flow, not the size or direction of the field by itself.

Key takeaway: To establish that a field is solenoidal, compute its divergence and verify that it vanishes at every point of the stated domain. If the formula has singularities, state those exclusions explicitly; zero divergence away from a singularity is not a claim about the singular point.

Solenoidal fields are useful models of incompressible flow and source-free flux. The local condition \(\nabla\cdot\mathbf{F}=0\) can identify balanced flow in a linear field, a circulating field, or a radial field whose strength changes with distance. The domain check determines precisely where each conclusion holds.

Check Your Understanding

Use component derivatives and state the domain whenever a field is not defined everywhere.

Compute the divergence of \(\mathbf{F}=(3x+y,\ -x+2y,\ -5z)\). Is it solenoidal on \(\mathbb{R}^3\)?
For \(\mathbf{F}=(ax,\ by,\ cz)\), what condition on \(a,b,c\) makes the field solenoidal?
Explain why \((-y,x,0)\) shows that zero divergence does not imply a zero vector field.
For \(\mathbf{F}=f(r)\mathbf{x}\) with \(r>0\), write the differential condition on \(f\) required for zero divergence.
Why can \(\mathbf{x}/r^3\) be called solenoidal on \(\mathbb{R}^3\setminus\{\mathbf{0}\}\) but not on all of \(\mathbb{R}^3\)?