Outward Flux Through Sphere Patches

Why Spheres Make Flux Easier

A sphere has a useful geometric feature: at every point, its outward normal points directly away from its center. This makes flux especially straightforward when a vector field is radial, and it gives a clear way to keep track of orientation for other fields.

The main tasks are to describe the patch with parameters, find the oriented area vector, and integrate the field’s component in the chosen normal direction. The radius and the direction of the normal both matter: an incorrect area factor changes the size of the flux, while a reversed normal changes its sign.

Definition: The flux of a vector field $\mathbf{F}$ across an oriented parametrized surface $\mathbf{r}(u,v)$, with parameter domain $D$, is $$ \iint_D \mathbf{F}(\mathbf{r}(u,v))\cdot \bigl(\mathbf{r}_u\times\mathbf{r}_v\bigr)\,du\,dv. $$ The cross product specifies the orientation as well as the area scaling. Reversing the orientation reverses the sign of the flux.

Parametrizing a Sphere

For a sphere of radius $R>0$ centered at the origin, use the polar angle $\phi$, measured from the positive $z$-axis, and the azimuthal angle $\theta$, measured in the $xy$-plane from the positive $x$-axis. The standard parametrization is

\mathbf{r}(\phi,\theta) = \bigl(R\sin\phi\cos\theta,\; R\sin\phi\sin\theta,\; R\cos\phi\bigr), \qquad 0\leq\phi\leq\pi,\quad 0\leq\theta\leq 2\pi.

A spherical patch is described by restricting $(\phi,\theta)$ to a region $D$. For example, the upper hemisphere is given by $0\leq\phi\leq\pi/2$, with the full azimuthal range $0\leq\theta\leq2\pi$. A patch cut out by meridians and a latitude circle can be described by separate lower and upper bounds for each angle.

To find the oriented area vector, differentiate with respect to both parameters:

\mathbf{r}_\phi = \bigl(R\cos\phi\cos\theta,\; R\cos\phi\sin\theta,\; -R\sin\phi\bigr), $$ $$ \mathbf{r}_\theta = \bigl(-R\sin\phi\sin\theta,\; R\sin\phi\cos\theta,\; 0\bigr).

Taking the cross product in the order $\mathbf{r}_\phi\times\mathbf{r}_\theta$ gives

\mathbf{r}_\phi\times\mathbf{r}_\theta = R^2\sin\phi \bigl(\sin\phi\cos\theta,\; \sin\phi\sin\theta,\; \cos\phi\bigr).

The vector in parentheses is the unit vector pointing radially outward. Thus this parameter order gives the outward orientation wherever $\sin\phi>0$. The cross product vanishes at the poles, where the angular parametrization is singular; those isolated points have zero surface area and do not affect the flux integral.

Key formula: For an outward-oriented patch on the sphere of radius $R$ centered at the origin, $$ d\mathbf{S} = R^2\sin\phi\, \bigl(\sin\phi\cos\theta,\; \sin\phi\sin\theta,\; \cos\phi\bigr)\,d\phi\,d\theta. $$ The factor $R^2\sin\phi$ is the surface-area scaling; the remaining vector is the outward unit normal.

Why the Oriented Area Vector Points Outward

The position vector $\mathbf{r}(\phi,\theta)$ points from the sphere’s center to the surface point. Dividing it by its length $R$ gives the outward unit normal. The cross product calculated above satisfies

\mathbf{r}_\phi\times\mathbf{r}_\theta = R\sin\phi\,\mathbf{r}(\phi,\theta) = R^2\sin\phi\,\frac{\mathbf{r}(\phi,\theta)}{R}.

So the cross product has the outward direction and magnitude $R^2\sin\phi$. This also explains why swapping the parameter order changes the orientation: $\mathbf{r}_\theta\times\mathbf{r}_\phi=-(\mathbf{r}_\phi\times\mathbf{r}_\theta)$.

For a sphere centered at a point $\mathbf{c}=(a,b,d)$, add that center to the parametrization. Its points are $\mathbf{c}+R\mathbf{u}$, where $\mathbf{u}$ is the unit radial direction displayed above. The outward normal is still $\mathbf{u}$, now pointing away from $\mathbf{c}$.

Radial Fields

A radial field based at the origin points in the direction of the position vector. Write it as $\mathbf{F}(\mathbf{x})=g(\|\mathbf{x}\|)\mathbf{x}/\|\mathbf{x}\|$, wherever $\|\mathbf{x}\|>0$. On a sphere of radius $R$, the field has the same radial direction as the outward unit normal, so their dot product is simply $g(R)$. The flux over an outward patch with angular region $D$ is therefore

\iint_D \mathbf{F}(\mathbf{r}(\phi,\theta))\cdot (\mathbf{r}_\phi\times\mathbf{r}_\theta)\,d\phi\,d\theta = g(R)R^2\iint_D\sin\phi\,d\phi\,d\theta.

This simplification uses the fact that the radial magnitude is constant on a sphere centered at the field’s radial origin. For a field radial from a different center, use spheres centered at that same point to get the same simplification.

Worked Example: Radial Flux Through an Upper Hemisphere

Find the outward flux of $\mathbf{F}(x,y,z)=(x,y,z)$ through the upper hemisphere of radius $2$ centered at the origin.

The parameter bounds are $0\leq\phi\leq\pi/2$ and $0\leq\theta\leq2\pi$. On the sphere, $\mathbf{F}=\mathbf{r}=2\mathbf{u}$, where $\mathbf{u}$ is the outward unit radial vector. Thus $\mathbf{F}\cdot\mathbf{u}=2$, and the outward flux is

\begin{aligned} \Phi &=\int_0^{2\pi}\int_0^{\pi/2} 2(2^2\sin\phi)\,d\phi\,d\theta\\ &=8\int_0^{2\pi}d\theta\int_0^{\pi/2}\sin\phi\,d\phi\\ &=8(2\pi)\bigl[-\cos\phi\bigr]_0^{\pi/2}\\ &=16\pi(1)=16\pi. \end{aligned}

As a check, the hemisphere’s area is $2\pi R^2=8\pi$, and the field’s outward normal component is constantly $2$. Their product is $2(8\pi)=16\pi$, matching the integral.

Constant Fields on Spherical Patches

A constant vector field does not generally point in the surface-normal direction. Its normal component varies across the sphere, so it should be dotted with the radial unit vector before integrating. If $\mathbf{F}=(A,B,C)$, then its outward normal component at angles $(\phi,\theta)$ is

\mathbf{F}\cdot\mathbf{u} = A\sin\phi\cos\theta+ B\sin\phi\sin\theta+ C\cos\phi.

The flux is this expression multiplied by the area scaling $R^2\sin\phi$, then integrated over the patch’s angle bounds. When the orientation is inward, negate the outward oriented area vector, and hence negate the resulting flux.

Worked Example: Constant Field Through a Bounded Patch

Find the outward flux of $\mathbf{F}=(1,-2,3)$ through the patch on the sphere of radius $2$ defined by $0\leq\phi\leq\pi/3$ and $0\leq\theta\leq\pi/2$.

The outward unit normal is $\mathbf{u}=(\sin\phi\cos\theta,\sin\phi\sin\theta,\cos\phi)$. Its dot product with the field is

\mathbf{F}\cdot\mathbf{u} = \sin\phi\cos\theta -2\sin\phi\sin\theta +3\cos\phi.

Multiplying by the area scaling $R^2\sin\phi=4\sin\phi$, the flux becomes

\begin{aligned} \Phi &=4\int_0^{\pi/2}\int_0^{\pi/3} \left( \sin^2\phi\cos\theta -2\sin^2\phi\sin\theta +3\sin\phi\cos\phi \right)d\phi\,d\theta\\ &=4\left[ \left(\int_0^{\pi/2}(\cos\theta-2\sin\theta)\,d\theta\right) \left(\int_0^{\pi/3}\sin^2\phi\,d\phi\right) + 3\left(\int_0^{\pi/2}d\theta\right) \left(\int_0^{\pi/3}\sin\phi\cos\phi\,d\phi\right) \right]. \end{aligned}

Evaluate each factor:

\begin{aligned} \int_0^{\pi/2}(\cos\theta-2\sin\theta)\,d\theta &=[\sin\theta+2\cos\theta]_0^{\pi/2}=1-2=-1,\\ \int_0^{\pi/3}\sin^2\phi\,d\phi &=\left[\frac{\phi}{2}-\frac{\sin(2\phi)}{4}\right]_0^{\pi/3} =\frac{\pi}{6}-\frac{\sqrt{3}}{8},\\ \int_0^{\pi/3}\sin\phi\cos\phi\,d\phi &=\left[\frac{\sin^2\phi}{2}\right]_0^{\pi/3} =\frac{3}{8},\\ \int_0^{\pi/2}d\theta&=\frac{\pi}{2}. \end{aligned}

Substitution gives

\begin{aligned} \Phi &=4\left[ -\left(\frac{\pi}{6}-\frac{\sqrt{3}}{8}\right) +3\left(\frac{\pi}{2}\right)\left(\frac{3}{8}\right) \right]\\ &=4\left(-\frac{\pi}{6}+\frac{\sqrt{3}}{8}+\frac{9\pi}{16}\right) =\frac{19\pi}{12}+\frac{\sqrt{3}}{2}. \end{aligned}

The result is positive for the stated outward orientation. Reversing the orientation gives the inward flux $-\frac{19\pi}{12}-\frac{\sqrt{3}}{2}$, as required by the sign change.

A Sphere with a Different Center

The normal on a sphere is radial from the sphere’s own center, not necessarily from the origin. For a radial field based at that center, use the displacement from the center when identifying the normal.

Worked Example: Radial Field on a Shifted Hemisphere

Consider the upper hemisphere of radius $2$ centered at $\mathbf{c}=(1,-1,2)$. Find the outward flux of

\mathbf{F}(x,y,z) = \frac{(x-1,\;y+1,\;z-2)} {\sqrt{(x-1)^2+(y+1)^2+(z-2)^2}}.

On this sphere, the denominator is $2$, and the numerator is $2\mathbf{u}$, where $\mathbf{u}$ points outward from $\mathbf{c}$. Therefore $\mathbf{F}=\mathbf{u}$, so $\mathbf{F}\cdot\mathbf{u}=1$. The upper hemisphere has bounds $0\leq\phi\leq\pi/2$, $0\leq\theta\leq2\pi$, and its outward flux is

\begin{aligned} \Phi &=\int_0^{2\pi}\int_0^{\pi/2} 1\cdot 2^2\sin\phi\,d\phi\,d\theta\\ &=4(2\pi)\bigl[-\cos\phi\bigr]_0^{\pi/2} =8\pi. \end{aligned}

Checking by area gives the same value: the hemisphere has area $2\pi(2^2)=8\pi$, and the normal component of the field is $1$ everywhere on it. The center shift changes the location of the patch, but not its area or this flux calculation.

Common Orientation and Area-Element Errors

A frequent mistake is to use the outward unit normal but omit the surface-area scaling. For a sphere, integrating only $\mathbf{F}\cdot\mathbf{u}\,d\phi\,d\theta$ is not a flux integral: the angular parameters do not measure surface area uniformly. The correct factor is $R^2\sin\phi$.

For example, in the upper-hemisphere calculation above, replacing $4\sin\phi$ by $1$ would give $2\pi$, not $16\pi$. That calculation is wrong because it drops both the radius-squared factor $4$ and the angular factor $\sin\phi$. Restoring them gives

\int_0^{2\pi}\int_0^{\pi/2} 2\cdot4\sin\phi\,d\phi\,d\theta =16\pi.

A second sign error occurs when the parameter order is reversed without changing the intended orientation. If the cross product is $\mathbf{r}_\theta\times\mathbf{r}_\phi$, it points inward; either use $\mathbf{r}_\phi\times\mathbf{r}_\theta$ for outward flux or explicitly negate the inward-oriented area vector.

Using the Formula Reliably

For a spherical patch, first identify the sphere’s center and radius, then determine the angular bounds describing the patch. Compute the normal direction from the center to the surface point, and decide whether the requested orientation is outward or inward. Finally, include $R^2\sin\phi$ and evaluate the resulting integral over the full patch domain.

For a radial field, the dot product with the outward normal is constant on a centered sphere when the field’s magnitude depends only on distance from that center. For a constant field, calculate the dot product with the varying radial unit normal before integrating. In both cases, check the sign against the stated orientation and verify the area scaling is present.

Check Your Understanding

Use the outward orientation unless a question specifies otherwise.

What is the outward unit normal at the point $\mathbf{r}(\phi,\theta)$ on a sphere centered at the origin?
Why does $\mathbf{r}_\phi\times\mathbf{r}_\theta$ have magnitude $R^2\sin\phi$?
Set up, but do not evaluate, the outward flux integral of $\mathbf{F}(x,y,z)=(0,0,1)$ over the upper hemisphere of radius $3$.
How does the flux change if the orientation on a spherical patch is reversed?
For a radial field based at a sphere’s center, why is its outward normal component constant over the sphere?