Tutorials › Calculus III › Divergence as Local Source Strength

Home  /  Tutorials  /  Curl and Divergence

Curl and Divergence · Tutorial 902 of 999

Divergence as Local Source Strength

Learn why divergence measures local net outflow per unit volume and how to interpret its positive, negative, and zero values.

Advanced 10 min read

What You'll Learn

  • Define local source strength using net outward flux through a shrinking box
  • Derive the limiting flux-to-volume ratio from the Cartesian divergence formula
  • Interpret positive divergence as local net outflow and negative divergence as local net inflow
  • Explain precisely what zero divergence says—and what it does not say
  • Compute and check flux ratios for fields with positive, negative, and zero divergence

From Flux to Local Source Strength

A vector field can describe the flow of a fluid, the movement of heat, or another quantity moving through space. Flux measures how much of that field crosses a surface. For a closed surface, the net outward flux compares what leaves the enclosed region with what enters it.

Divergence gives a pointwise way to describe this balance. Its interpretation as source strength comes from looking at the net outward flux through a small box around a point and dividing by the box’s volume. As the box shrinks, this ratio approaches the divergence at that point.

Central idea: Divergence at a point is the limiting net outward flux per unit volume near that point. Positive divergence indicates local net outflow; negative divergence indicates local net inflow.

Definition and the Flux-to-Volume Ratio

Let \(\mathbf{F}=\langle P,Q,R\rangle\) be a vector field whose component functions have continuous first partial derivatives in a neighborhood of \(\mathbf{p}=(a,b,c)\). Its divergence is the scalar field

$$ \nabla\cdot\mathbf{F} = \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z}. $$

For \(h>0\), take the box centered at \(\mathbf{p}\) whose sides are parallel to the coordinate axes and whose half-length is \(h\). Its side length is \(2h\), and its volume is \(8h^3\). Write \(\Phi(h)\) for the net outward flux of \(\mathbf{F}\) through the six faces of this box.

Definition: The local source-strength interpretation of divergence at \(\mathbf{p}\) is the limit $$ \lim_{h\to 0^+}\frac{\Phi(h)}{8h^3} = \nabla\cdot\mathbf{F}(\mathbf{p}). $$ The ratio is net outward flux divided by volume, not the flux through any one face.

The continuity assumption on the first partial derivatives ensures that the face-flux expressions have the limiting behavior needed for this formula. The ratio can also be negative: outward flux is counted as positive, while inward flux contributes negatively.

Why the Limit Equals Divergence

We can derive the formula directly from the flux across the box’s faces. On the two faces perpendicular to the \(x\)-axis, the outward unit normals are \(\langle 1,0,0\rangle\) and \(\langle-1,0,0\rangle\). Thus, the combined flux through those faces is the integral over their shared \(yz\)-range of the difference between the \(x\)-components on the two faces. The \(y\)- and \(z\)-face pairs have corresponding expressions.

1
Pair the \(x\)-faces.
Their combined outward flux is the integral of \(P(a+h,y,z)-P(a-h,y,z)\) over the square \([b-h,b+h]\times[c-h,c+h]\).
2
Use the fundamental theorem of calculus.
The difference in the integrand equals \(\int_{-h}^{h} P_x(a+t,y,z)\,dt\), where \(P_x\) denotes the partial derivative of \(P\) with respect to \(x\).
3
Repeat for the other face pairs.
The combined \(y\)-face flux uses \(Q_y\), and the combined \(z\)-face flux uses \(R_z\).
4
Divide by the volume and take the limit.
The resulting ratio is the average of \(\nabla\cdot\mathbf{F}\) over the box. Continuity makes this average approach its value at the center.

In detail, applying the fundamental theorem of calculus to the paired face fluxes gives

$$ \begin{aligned} \Phi(h) &= \int_{b-h}^{b+h}\int_{c-h}^{c+h} \bigl(P(a+h,y,z)-P(a-h,y,z)\bigr)\,dz\,dy\\ &\quad+ \int_{a-h}^{a+h}\int_{c-h}^{c+h} \bigl(Q(x,b+h,z)-Q(x,b-h,z)\bigr)\,dz\,dx\\ &\quad+ \int_{a-h}^{a+h}\int_{b-h}^{b+h} \bigl(R(x,y,c+h)-R(x,y,c-h)\bigr)\,dy\,dx\\ &= \iiint_{[a-h,a+h]\times[b-h,b+h]\times[c-h,c+h]} \left(P_x+Q_y+R_z\right)\,dV. \end{aligned} $$

Here \(dV\) denotes volume integration. Each face pair becomes a triple integral: for instance, the first difference is integrated in \(x\) from \(a-h\) to \(a+h\), as well as over the full \(y\)- and \(z\)-ranges. Adding the three results gives the displayed integral of the divergence.

Dividing by the box’s volume \(8h^3\) shows that \(\Phi(h)/(8h^3)\) is the average value of the continuous function \(\nabla\cdot\mathbf{F}\) over the box. As \(h\) approaches zero, every point in the box approaches \(\mathbf{p}\). By continuity, the values of the divergence throughout the box approach \(\nabla\cdot\mathbf{F}(\mathbf{p})\), so their average approaches the same value. This proves the local source-strength formula.

What the derivation requires: The field must have continuous first partial derivatives near the center, and the box faces must be counted with outward normals. The argument establishes the limit for these axis-aligned shrinking boxes.

Worked Example: Positive Divergence

Worked Example: Positive Divergence

Let \(\mathbf{F}=\langle x,y,z\rangle\) and center a box of half-length \(h\) at \(\mathbf{p}=(1,-1,2)\). Compute the divergence and the net outward flux per unit volume.

First, the three required partial derivatives are

$$ P_x=1,\qquad Q_y=1,\qquad R_z=1, \qquad\text{so}\qquad \nabla\cdot\mathbf{F}=1+1+1=3. $$

On the positive and negative \(x\)-faces, the \(x\)-component is respectively \(1+h\) and \(1-h\). Their outward fluxes are \((1+h)4h^2\) and \(-(1-h)4h^2\). Together they contribute \(8h^3\). The corresponding \(y\)-face contributions are \((-1+h)4h^2\) and \(-(-1-h)4h^2\), also totaling \(8h^3\). The \(z\)-face contributions are \((2+h)4h^2\) and \(-(2-h)4h^2\), again totaling \(8h^3\).

Therefore, the total flux and the flux-to-volume ratio are

$$ \Phi(h)=8h^3+8h^3+8h^3=24h^3, \qquad \frac{\Phi(h)}{8h^3} = \frac{24h^3}{8h^3} =3. $$

The ratio equals the divergence for every \(h>0\), so its limit is \(3\). This checks the calculation against \(\nabla\cdot\mathbf{F}(1,-1,2)=3\). The positive sign means the box has net outward flux.

Worked Example: Negative Divergence

Worked Example: Negative Divergence

Let \(\mathbf{F}=\langle-x,-2y,-3z\rangle\), and center the box at the origin. Find the divergence and the flux-to-volume ratio.

The partial derivatives are \(P_x=-1\), \(Q_y=-2\), and \(R_z=-3\). Hence

$$ \nabla\cdot\mathbf{F}=-1-2-3=-6. $$

For the two \(x\)-faces, the component values are \(-h\) on the positive face and \(h\) on the negative face. Including their outward normals, their combined flux is \((-h)4h^2-(h)4h^2=-8h^3\). The \(y\)-faces contribute \((-2h)4h^2-(2h)4h^2=-16h^3\). The \(z\)-faces contribute \((-3h)4h^2-(3h)4h^2=-24h^3\).

Adding the signed contributions gives

$$ \Phi(h)=-8h^3-16h^3-24h^3=-48h^3, \qquad \frac{\Phi(h)}{8h^3} = \frac{-48h^3}{8h^3} =-6. $$

The ratio is exactly the divergence at the origin, confirming the result. Its negative sign indicates net inward flux: more of the field enters the box than leaves it.

Worked Example: Zero Divergence and Flux Balance

Worked Example: Zero Divergence and Flux Balance

Let \(\mathbf{F}=\langle x,-y,0\rangle\), with a box centered at the origin. Calculate the divergence and compare the contributions from the \(x\)- and \(y\)-faces.

The derivatives are \(P_x=1\), \(Q_y=-1\), and \(R_z=0\), so

$$ \nabla\cdot\mathbf{F}=1-1+0=0. $$

The \(x\)-face contribution is \(h\cdot4h^2-(-h)\cdot4h^2=8h^3\). The \(y\)-face contribution is \((-h)\cdot4h^2-h\cdot4h^2=-8h^3\). The \(z\)-component is zero on both \(z\)-faces, so their combined flux is zero. Thus

$$ \Phi(h)=8h^3-8h^3+0=0, \qquad \frac{\Phi(h)}{8h^3} = \frac{0}{8h^3}=0. $$

The flux-to-volume ratio agrees with the divergence. The field is not zero: the outward flux through the \(x\)-faces is balanced by inward flux through the \(y\)-faces. Zero divergence expresses balance in the net flux, not an absence of flow.

Interpreting the Sign Carefully

The three examples illustrate the standard local interpretation:

Divergence at the pointLocal flux interpretation
PositiveNet outward flux per unit volume is positive; the point behaves locally like a source.
NegativeNet outward flux per unit volume is negative; the point behaves locally like a sink.
ZeroThe limiting net outward flux per unit volume is zero; local outflow and inflow balance in this limit.

These descriptions concern the local flux balance, not whether the vector field itself vanishes. A field can have zero divergence and still be nonzero. Also, divergence equal to zero at a point does not, by itself, say that the net flux through every finite box around that point is exactly zero. It says that the flux divided by volume tends to zero as the box shrinks.

Common Pitfall: Ignoring the Sign of Inward Flux

A common error is to add the magnitudes of all face fluxes instead of adding their signed outward fluxes. In the zero-divergence example, each of the two \(x\)-faces has outward flux \(4h^3\), so the \(x\)-pair contributes \(8h^3\). Each \(y\)-face has flux of magnitude \(4h^3\) as well, but these fluxes are inward relative to the outward normals. Their signed contribution is \(-8h^3\), not \(+8h^3\).

$$ \text{Incorrect: }\quad \frac{8h^3+8h^3}{8h^3}=2. \qquad \text{Correct: }\quad \frac{8h^3-8h^3}{8h^3}=0. $$

The incorrect calculation counts incoming flow as outgoing. The correction is to use the outward normal on each face and retain the sign of the dot product with the field.

Why This Interpretation Matters

The definition of divergence is a local derivative calculation, while the flux-to-volume limit gives that calculation a geometric meaning. It explains why divergence is useful when studying sources, sinks, and local balance in a vector field. The box argument also shows why all three component derivatives enter: the net flux is assembled from three pairs of faces, one for each coordinate direction.

When interpreting a divergence value, keep the distinction clear: it is a signed net flux per unit volume in the shrinking-box limit. It is not the field’s magnitude, and it does not describe the flux through just one face. Keeping those points in view makes the signs and the source-strength interpretation precise.

Check Your Understanding

Use the shrinking-box interpretation to answer these questions.

  1. For \(\mathbf{F}=\langle 2x,y,-z\rangle\), find the divergence and state whether the local flux balance is outward or inward.
  2. Why must inward flux be assigned a negative sign when calculating net outward flux?
  3. If the divergence at a point is zero, what does the flux-to-volume limit say?
  4. Can a field with zero divergence be nonzero? Give a reason based on an example from this tutorial.
  5. Which regularity condition on the field is used in deriving the shrinking-box limit?

Previous

Next

Back to Tutorials