Independent study guide

Surface flux with oriented area vectors and closing surfaces

Build an oriented surface-area vector, check divergence-theorem conditions and calculate an open-surface flux with a signed closing-surface correction.

Nablio · 35 minPublished

Return to Humber Polytechnic CALC 1500 course guide. Tutoring coverage for this course must be confirmed before booking and payment.

Stuck on a step? Text us Email us Request a lesson

A flux integral measures the signed contribution of a vector field through a surface. Before integrating, decide which surface is included and which side its normal points toward. A correct divergence calculation can still answer the wrong question if a closing surface has been included but never removed.

This lesson builds a direct surface integral, explains when the divergence theorem applies, and calculates the flux through an open hemisphere in two ways. You should be comfortable with dot products, cross products, partial derivatives and double integrals. All coordinates and field values in the examples are dimensionless.

Flux and orientation

For an oriented surface SS, write

ΦS=∬SF⋅n dS.\Phi_S=\iint_S\mathbf F\cdot\mathbf n\,dS.

Here n\mathbf n is a unit normal, and dSdS is a positive surface-area element. The dot product selects the component of the field perpendicular to the surface. A positive contribution points with the chosen normal; a negative contribution points against it. A tangent field contributes zero at that point.

Choosing the opposite normal reverses the sign of the flux. It does not change the surface area. A zero integral can mean that positive and negative contributions cancel; it does not establish that the field vanishes.

An open surface needs a stated orientation, such as upward or toward the positive x direction. “Outward” is meaningful when a reference solid has been specified. For a closed boundary, outward means away from the enclosed solid, including into any cavity.

A disk includes its circular edge but is still an open surface in this calculus terminology: it has a boundary curve and does not enclose a solid by itself. Here “closed surface” means a complete enclosing boundary, not merely a closed set of points.

Build the area vector from two parameters

Let r(u,v)\mathbf r(u,v) describe SS once over a parameter region DD, apart from shared edges or coordinate seams. Its partial derivatives ru\mathbf r_u and rv\mathbf r_v are tangent to the surface. Their cross product is normal to it:

A=ru×rv,dS=∥A∥ du dv.\begin{gathered}\mathbf A=\mathbf r_u\times\mathbf r_v,\\dS=\|\mathbf A\|\,du\,dv.\end{gathered}

For a regular patch, A≠0\mathbf A\ne\mathbf0. The cross-product magnitude converts parameter area into surface area. It generally varies across the patch.

If A\mathbf A points in the required direction, then

n=A∥A∥,n dS=A du dv,\begin{gathered}\mathbf n=\frac{\mathbf A}{\|\mathbf A\|},\\\mathbf n\,dS=\mathbf A\,du\,dv,\end{gathered}

so

ΦS=∬DF(r(u,v))⋅(ru×rv) du dv.\begin{gathered}\Phi_S=\iint_D\begin{aligned}&\mathbf F(\mathbf r(u,v))\cdot{}\\& (\mathbf r_u\times\mathbf r_v)\end{aligned}\,du\,dv.\end{gathered}

For each connected regular patch with a continuous chosen orientation, check the sign at an easy point. The relative sign cannot change on that patch without losing regularity or continuity. Check disconnected patches separately. For an upward orientation, the z component of the area vector must be positive. For an outward sphere, its dot product with a radial outward vector must be positive. If the cross product points the wrong way, negate it or reverse the cross-product order, but do not do both.

There are two consistent ways to calculate: use a unit normal together with dSdS, or use the unnormalized cross product with du dvdu\,dv. Normalizing the cross product while retaining only du dvdu\,dv loses its area scale. Multiplying the unnormalized cross product by a second ∥A∥\|\mathbf A\| counts that scale twice.

A single chart need not work everywhere. Spherical coordinates, for example, degenerate at a pole. One can use regular patches and a limit; the pole and coordinate seam do not change the ordinary surface integrals used here. Do not use a parameterization that covers a positive-area part twice.

A tilted planar patch above a rectangular xy projection. Two tangent arrows lie in the patch; an upward normal arrow is perpendicular to both. The projected rectangle and tilted patch are labelled separately.
Here xx and yy are the surface parameters, and rx\mathbf r_x and ry\mathbf r_y are the tangent vectors shown. Their cross product A=rx×ry\mathbf A=\mathbf r_x\times\mathbf r_y points upward and is perpendicular to both. Its length gives the surface-area factor: dS=∥A∥ dx dydS=\|\mathbf A\|\,dx\,dy. The unit normal is n=A/∥A∥\mathbf n=\mathbf A/\|\mathbf A\|. The field vector is not shown.

Use the graph formula carefully

For z=g(x,y)z=g(x,y) over a region DD in the xy plane, choose

r(x,y)=(x,y,g(x,y)).\mathbf r(x,y)=(x,y,g(x,y)).

Then

rx=(1,0,gx),ry=(0,1,gy),rx×ry=(−gx,−gy,1).\begin{gathered}\mathbf r_x=(1,0,g_x),\\\mathbf r_y=(0,1,g_y),\\\mathbf r_x\times\mathbf r_y=(-g_x,-g_y,1).\end{gathered}

Its last component is positive, so this is the upward area-vector density. The downward one is (gx,gy,−1)(g_x,g_y,-1).

For F=(P,Q,R)\mathbf F=(P,Q,R), the upward flux is

ΦS=∬D ⁣[−P(x,y,g)gx−Q(x,y,g)gy+R(x,y,g)] dx dy.\Phi_S=\iint_D\!\left[\begin{aligned}&-P(x,y,g)g_x\\&-Q(x,y,g)g_y\\&+R(x,y,g)\end{aligned}\right]\,dx\,dy.

Evaluate every field component on the surface before simplifying. The corresponding scalar area is ∬D1+gx2+gy2 dx dy\iint_D\sqrt{1+g_x^2+g_y^2}\,dx\,dy; it is a different integral.

Worked patch

Take z=x+yz=x+y, 0≤x≤10\le x\le1, 0≤y≤10\le y\le1, upward orientation, and F=(1,0,2)\mathbf F=(1,0,2). The area-vector density is (−1,−1,1)(-1,-1,1), so

ΦS=∫01∫01 ⁣(1,0,2)⋅(−1,−1,1) dy dx=1.\begin{gathered}\Phi_S=\int_0^1\int_0^1\!\begin{aligned}&(1,0,2)\cdot{}\\&(-1,-1,1)\end{aligned}\,dy\,dx\\=1.\end{gathered}

Check with the other convention: the unit normal is (−1,−1,1)/3(-1,-1,1)/\sqrt3 and the patch area is 3\sqrt3. Their factors cancel, giving (1/3)3=1(1/\sqrt3)\sqrt3=1. The downward flux is −1-1.

If you change the projection coordinates to polar coordinates, the planar area element also changes: dx dy=ρ dρ dθdx\,dy=\rho\,d\rho\,d\theta. This factor belongs to the change of variables in the projection plane. It is not an extra surface-area factor to attach after a full polar cross product has already been computed.

Try question 1 on a tilted patch, then question 2 on a curved graph.

Check the divergence theorem conditions

For the usual introductory form, let EE be a bounded solid with a closed, piecewise-smooth boundary ∂E\partial E. Assume F\mathbf F has continuous first partial derivatives on an open set containing the whole solid and its boundary. With outward normals,

∬∂EF⋅n dS=∭E∇⋅F dV.\begin{gathered}\iint_{\partial E}\mathbf F\cdot\mathbf n\,dS\\=\iiint_E\nabla\cdot\mathbf F\,dV.\end{gathered}

Check these three things before using the equality:

  1. The surface on the left contains the entire boundary of the solid on the right.
  2. Every boundary piece uses the outward normal for that same solid.
  3. The field satisfies the differentiability condition throughout a neighbourhood of the solid, including its interior.

A polynomial field meets the third condition everywhere. A field with a denominator that vanishes inside the solid needs a separate domain check, even when it is smooth on the outer surface. Computing divergence only where a formula is defined does not repair a missing hypothesis.

If you need a closed-boundary starting example, return to the MAST 219 box calculation. It explicitly adds all six faces. Question 6 shows why a singularity inside a sphere matters.

Close an open surface and subtract the signed contribution

Suppose an open surface SS and an added surface BB together form ∂E\partial E. They may meet along an edge, but their interiors must not overlap. Orient both outward from EE. Then

ΦS+ΦB=∭E∇⋅F dV,ΦS=∭E∇⋅F dV−ΦB.\begin{gathered}\Phi_S+\Phi_B=\iiint_E\nabla\cdot\mathbf F\,dV,\\\Phi_S=\iiint_E\nabla\cdot\mathbf F\,dV-\Phi_B.\end{gathered}

The quantity subtracted is the signed flux through the added surface. Subtracting a negative cap contribution increases the result. If several pieces are needed to close the surface, subtract their sum. If the requested orientation of SS is inward, first obtain its outward flux and then reverse the sign.

Worked hemisphere with a nonzero closing disk

Let SS be the curved upper unit hemisphere

x2+y2+z2=1,z≥0,\begin{gathered}x^2+y^2+z^2=1,\\z\ge0,\end{gathered}

oriented away from the origin. The field is F=(x,y,z+1)\mathbf F=(x,y,z+1). The equatorial circle is the edge of SS; the disk it bounds is not included in the requested flux.

Add the disk B: z=0, x2+y2≤1B:\ z=0,\ x^2+y^2\le1. Together they bound the upper half-ball EE. The field is polynomial, the boundary is piecewise smooth, and ∇⋅F=3\nabla\cdot\mathbf F=3. Therefore the total outward flux is

ΦS+ΦB=3Vol⁡(E)=3(2π3)=2π.\begin{gathered}\Phi_S+\Phi_B=3\operatorname{Vol}(E)\\=3\left(\frac{2\pi}{3}\right)\\=2\pi.\end{gathered}

On the disk, the outward normal points down because the solid lies above it. Thus nB=(0,0,−1)\mathbf n_B=(0,0,-1). At z=0z=0, the field is (x,y,1)(x,y,1), giving F⋅nB=−1\mathbf F\cdot\mathbf n_B=-1. Consequently,

ΦB=−π,ΦS=2π−(−π)=3π.\begin{gathered}\Phi_B=-\pi,\\\Phi_S=2\pi-(-\pi)\\=3\pi.\end{gathered}

Reporting 2π2\pi would give the flux through the hemisphere and disk together. Reporting π\pi would subtract the disk's unsigned magnitude instead of its signed flux.

The curved upper unit hemisphere is shown first without its base. A second view closes it with a disk at z equals zero. Radial arrows point outward on the hemisphere and the disk normal points downward, away from the half-ball.
The theorem includes the added disk. Its signed flux is −π-\pi, so the hemisphere flux is 2π−(−π)=3π2\pi-(-\pi)=3\pi.

Verify the hemisphere directly

Use colatitude φ\varphi, measured down from the positive z axis, and azimuth θ\theta:

r(φ,θ)=(sin⁡φcos⁡θsin⁡φsin⁡θcos⁡φ),0≤φ≤π2,0≤θ≤2π.\begin{gathered}\mathbf r(\varphi,\theta)=\begin{pmatrix}\sin\varphi\cos\theta\\\sin\varphi\sin\theta\\\cos\varphi\end{pmatrix},\\0\le\varphi\le\frac\pi2,\\0\le\theta\le2\pi.\end{gathered}

The outward unit normal on the unit sphere is r\mathbf r. In this parameter order,

rφ×rθ=sin⁡φ r.\mathbf r_\varphi\times\mathbf r_\theta =\sin\varphi\,\mathbf r.

Its radial dot product is positive for 0<φ≤π/20<\varphi\le\pi/2, and dS=sin⁡φ dφ dθdS=\sin\varphi\,d\varphi\,d\theta. Since

F⋅n=x2+y2+z(z+1)=1+z=1+cos⁡φ,\begin{gathered}\mathbf F\cdot\mathbf n =x^2+y^2+z(z+1)\\=1+z\\=1+\cos\varphi,\end{gathered}

we obtain

ΦS=∫02π ⁣∫0π/2 ⁣(1+cos⁡φ)⋅sin⁡φ dφ dθ=2π(1+12)=3π.\begin{gathered}\Phi_S=\int_0^{2\pi}\!\int_0^{\pi/2}\!\begin{aligned}&(1+\cos\varphi)\\&{}\cdot\sin\varphi\end{aligned}\,d\varphi\,d\theta\\=2\pi\left(1+\frac12\right)\\=3\pi.\end{gathered}

The direct calculation and the closing-disk method agree. Both account for the hemisphere only and use the same orientation. For a changed boundary above the equator, try question 4.

Recognize which integral the question asks for

A scalar surface integral such as ∬S1 dS\iint_S 1\,dS measures area and does not use an orientation. The flux ∬SF⋅n dS\iint_S\mathbf F\cdot\mathbf n\,dS does use one. A line integral ∫CF⋅dr\int_C\mathbf F\cdot d\mathbf r follows an oriented curve and measures a different quantity.

The divergence theorem relates flux of F\mathbf F across a closed boundary to divergence inside a solid. Stokes' theorem relates circulation along a boundary curve to flux of ∇×F\nabla\times\mathbf F through a compatible oriented surface. An open surface can be used in Stokes, but that does not turn an arbitrary flux of F\mathbf F into a circulation calculation.

If the original question follows a curve, use the existing line-integral setup lesson. For the cylindrical side in question 3, a direct parameterization is short and also makes a two-cap check possible.

Repair the first unsupported step

  • Wrong sign: write the requested normal direction and test your cross product at one point.
  • Wrong magnitude: identify whether your normal is unit length and where the area factor appears.
  • Wrong domain: distinguish the parameter region, the projection region, the surface and the enclosed solid.
  • Missing cap: list every boundary piece and keep its signed contribution until the final subtraction.
  • Invalid theorem: check the field's domain on the whole solid before integrating its divergence.

Try the six progressive flux questions. Bring a sketch and the first step you cannot justify to a tutor. Work through the correction, then set up a changed problem without the solution.

Course connections and references

This lesson supports selected surface-integral and Gauss/divergence topics in Concordia MAST 219 and Concordia ENGR 233. Choose it when your current chapter covers surfaces and flux. It is independent practice and makes no claim about a particular assessment or its timing.

Curriculum scope was checked on 11 October 2026 against Concordia's mathematics calendar and engineering calendar. The OpenStax surface-integral section and divergence-theorem section provide further reference. The worked examples here were written for Nablio.

Useful for

Work on the step you’re missing.

Tell us the subject and your preferred time. Nablio confirms the details by text or email.

Request a lesson