How can we measure how much of a vector field flows through a surface in space?
How can we calculate the amount of a vector field that flows through common surfaces, such as the graph of a function \(z=f(x,y)\text{?}\)
Section 11.6 showed how we can use vector valued functions of two variables to give a parameterization of a surface in space. For instance, the function \(\vr(s,t)=\langle 2\cos(t)\sin(s), 2\sin(t)\sin(s),2\cos(s)\rangle\) with domain \(0\leq t\leq 2 \pi\) and \(0\leq s\leq \pi\) parameterizes a sphere of radius \(2\) centered at the origin. Section 11.6 also gives examples of how to write parameterizations based on other geometric relationships like when one coordinate can be written as a function of the other two. In Subsection 11.6.2, we set up a Riemann sum based on a parameterization that would measure the surface area of our curved surfaces in space.
In Figure 12.9.1, we plot a surface using a parametrization \(\vr(s,t)=\langle{f(s,t),g(s,t),h(s,t)}\rangle\text{.}\) The magenta curves represent curves where \(s\) varies and \(t\) is held constant, while the yellow curves represent curves where \(t\) varies and \(s\) is held constant. The vector in magenta is \(\vr_s=\frac{\partial \vr}{\partial s}=\langle{f_s,g_s,h_s}\rangle\) which measures the direction and magnitude of change in the coordinates of the surface when only \(s\) is varied. Similarly, the vector in yellow is \(\vr_t=\frac{\partial \vr}{\partial t}=\langle{f_t,g_t,h_t}\rangle\) which measures the direction and magnitude of change in the coordinates of the surface when only \(t\) is varied. We also plot the parallelogram that is formed by \(\vr_s\) and \(\vr_t\text{,}\) which is tangent to the surface. The area of this parallelogram offers an approximation for the surface area of a patch of the surface.
From Section 9.4, we also know that \(\vr_s\times \vr_t\) (plotted in green) will be orthogonal to both \(\vr_s\) and \(\vr_t\) and its magnitude will be given by the area of the parallelogram.
Preview Activity12.9.1.
In this preview activity, we will explore the parameterizations of a few familiar surfaces and confirm some of the geometric properties described in the introduction above.
(a)
Use the ideas from Section 11.6 to give a parameterization \(\vr(s,t)\) of the following surface. Be sure to specify the bounds on each of your parameters.
\(S_1\text{:}\) the right circular cylinder centered on the \(x\)-axis of radius 2 when \(0\leq x\leq 5\)
(b)
Draw a graph of \(S_1\) from the previous part. Label the points that correspond to \((s,t)\) points of \((0,0)\text{,}\)\((0,1)\text{,}\)\((1,0)\text{,}\) and \((2,3)\text{.}\)
(c)
For the parameterization from part a, calculate \(\vr_s\text{,}\)\(\vr_t\text{,}\) and \(\vr_s \times \vr_t\text{.}\)
(d)
For the parameterization from part a, find the value for \(\vr_s\text{,}\)\(\vr_t\text{,}\) and \(\vr_s \times \vr_t\) at the \((s,t)\) points of \((0,0)\text{,}\)\((0,1)\text{,}\)\((1,0)\text{,}\) and \((2,3)\text{.}\)
(e)
Draw your vector results from d on your graph and confirm the geometric properties described in the introduction to this section. Namely, \(\vr_s\) and \(\vr_t\) should be tangent to the surface, while \(\vr_s \times \vr_t\) should be orthogonal to the surface (in addition to \(\vr_s\) and \(\vr_t\)).
(f)
Repeat steps (a) through (e) for the following two surfaces:
(i)
\(S_2\text{:}\) the sphere centered at the origin of radius 3
(ii)
\(S_3\text{:}\) the first octant portion of the plane \(x+2y+3z=6\)
As we saw in Section 11.6, we can set up a Riemann sum of the areas for the parallelograms in Figure 12.9.1 to approximate the surface area of the region plotted by our parametrization. Equation (11.6.2) shows that we can compute the exact surface area by taking a limit of a Riemann sum which will correspond to integrating the magnitude of \(\vr_s \times \vr_t\) over the appropriate parameter bounds. What if we wanted to measure a quantity other than the surface area? Our focus in this section we will be the exploration of a specific case of this question: How can we measure the amount of a three dimensional vector field that flows through a particular section of a surface? The geometric tools we have reviewed in this section, especially the vector \(\vr_s \times \vr_t\text{,}\) will be valuable.
Subsection12.9.1The Idea of the Flux of a Vector Field through a Surface
In this section we will look at how to measure the amount of a vector field that flows through a surface in space. Figure 12.9.2 illustrates a plot that demonstrates this idea. Our definition of divergence in Section 12.6 looked at measuring the amount of vector field flowing out of a small region on a 2D plane. In this subsection, we will set up a precise measurement of this same measurement but over a region of a curved surface in 3D.
As with understanding line integrals of vector-valued functions in Section 12.2, we don’t care about the output of the vector field at points away from the surface. We would really would like to examine the output vectors for the points on our surface. To do this, we will look at Figure 12.9.3, which plots the output of our vector field at an array of points on our surface.
The central question we would like to consider is “How can we measure the amount of a three dimensional vector field that flows through a particular section of a curved surface?”, so we only need to consider the amount of the vector field that flows through the surface. Any portion of our vector field that flows along (or tangent) to the surface will not contribute to the amount that goes through the surface. In Figure 12.9.4, we have split the vector field for points on our surface into two components. One component, plotted in green, is orthogonal to the surface. The component that is tangent to the surface is plotted in magenta.
In order to measure the amount of the vector field that moves through the plotted section of the surface, we must find the accumulation of the lengths of the green vectors in Figure 12.9.4. Notice that some of the green vectors are moving through the surface in a direction opposite of others. In other words, we will need to pay attention to the direction in which these vectors move through our surface and not just the magnitude of the green vectors.
If we have a parameterization of the surface, then the vector \(\vr_s \times \vr_t\) varies smoothly across our surface and gives a consistent way to describe which direction we choose as “through” the surface. If we define a positive flow through our surface as being consistent with the yellow vector in Figure 12.9.4, then there is more positive flow (in terms of both magnitude and area) than negative flow through the surface. Thus, the net flow of the vector field through this surface is positive.
Activity12.9.2.Visualizing flux through a surface.
In this activity, you will compare the net flow of different vector fields through our sample surface. In Figure 12.9.5 you can select between five different vector fields. Once you select a vector field, the vector field for a set of points on the surface will be plotted in blue. Each blue vector will also be split into its normal component (in green) and its tangential component (in magenta). The yellow vector defines the direction for positive flow through the surface.
Look at each vector field and order the vector fields from greatest net flow through the surface to least net flow through the surface. Remember that a negative net flow through the surface should be lower in your rankings than any positive net flow.
Subsection12.9.2The Details of Measuring the Flux of a Vector Field through a Surface
Now that we have developed a conceptual understanding of what we are trying to measure, we can set up the corresponding Riemann sum to measure the flux of a vector field through a section of a surface. Let \(Q\) be the section of our surface and suppose that \(Q\) is parameterized by \(\vr(s,t)\) with \(a\leq s\leq b\) and \(c \leq t \leq d\text{.}\) The domain of \(\vr\) is a region of the \(st\)-plane, which we call \(D\text{,}\) and the range of \(\vr\) is \(Q\text{.}\)
As with most problems in integral calculus, we slice our region of interest into smaller pieces. Specifically, we slice \(a\leq s\leq b\) into \(n\) equally-sized subintervals with endpoints \(s_0,s_1,\ldots,s_n\) and \(c \leq t \leq d\) into \(m\) equally-sized subintervals with endpoints \(t_0,t_1,\ldots,t_n\text{.}\) This divides \(D\) into \(nm\) rectangles of size \(\Delta{s}=\frac{b-a}{n}\) by \(\Delta{t}=\frac{d-c}{m}\text{.}\) We index these rectangles as \(D_{i,j}\text{.}\) Every \(D_{i,j}\) has area (in the \(st\)-plane) \(\Delta{s}\Delta{t}\text{.}\) The partition of \(D\) into the rectangles \(D_{i,j}\) also partitions \(Q\) into \(nm\) corresponding pieces which we call \(Q_{i,j}=\vr(D_{i,j})\text{.}\) From Section 11.6 (specifically (11.6.1)) the surface area of \(Q_{i,j}\) is approximated by \(S_{i,j}=\vecmag{(\vr_s \times \vr_t)(s_i,t_j)}\Delta{s}\Delta{t}\text{.}\)
We want to measure the total flow of the vector field \(\vF\) through \(Q\text{,}\) which we will approximate on each \(Q_{i,j}\) and then sum to obtain the total flow. In other words, the flux of \(\vF\) through \(Q\) is
where \(\vecmag{\vF_{\perp Q_{i,j}}}\) is the length of the component of \(\vF\) orthogonal to \(Q_{i,j}\text{.}\)
For each \(Q_{i,j}\text{,}\) we approximate the surface \(Q\) by the tangent plane to \(Q\) at a corner of that partition element. This corresponds to using the planar elements in Figure 12.9.6, which have surface area \(S_{i,j}\text{.}\) The vector \(\vw_{i,j}=(\vr_s \times \vr_t)(s_i,t_j)\) can be used to measure the orthogonal direction (and thus define which direction we mean by positive flow through \(Q\)) on the \(i,j\) partition element. This means that
Taking the limit as \(n,m\rightarrow\infty\) gives the following result.
Conceptually, the above limit shows how the argument in Subsection 11.6.2 generalizes when measuring the flux of a vector field through a surface rather than just the surface area. The key difference between the Riemann sum above and Riemann sum used to set up Surface area is that we are using the dot product of the vector field with the normal vector to the surface (given by \(\vw_{i,j}=(\vr_s \times \vr_t)(s_i,t_j)\)) as the scalar value in our sum. This realization gives the following theorem, which states that the flux through a curved surface in space can be computed using a parameterization of the surface and a double integral over a flat region in the plane of parameter values.
Theorem12.9.7.
Let a smooth surface \(Q\) be parametrized by \(\vr(s,t)\) over a domain \(D\text{.}\) The total flux of a smooth vector field \(\vF\) through \(Q\) is given by
In Figure 12.9.6, you can change the number of sections in the partition and see the geometric result of refining the partition. In the next example, we will look at how Theorem 12.9.7 is used on a part of a cone and make sense of the vector \(\vw(t,s)=\vr_s \times \vr_t\) and the scalar \(\vF \cdot (\vr_s \times \vr_t)\text{.}\)
Remark12.9.8.
In some Webwork exercises and other sources, you will see flux integrals specified with the following notations:
While we will use the notation specified in Theorem 12.9.7 primarily, the other notations come from combining the normal vector (calculated from the parameterization of the surface) and the area element into a single vector valued area element. In other words, \((\vr_s \times \vr_t)\, dA = d\vec{A} = d\vec{s}\text{.}\)
Example12.9.9.
In this example we will compute the flux of the vector field \(\vF= \langle xz-2, y+x, 1 \rangle \) through the surface of the cone given by \(x^2+y^2=z^2\) with \(1\leq z \leq 3\text{.}\) So that we can use Theorem 12.9.7, we first parameterize our surface and calculate \(\vw(t,s)=\vr_s \times \vr_t\) based on that parameterization.
In cylindrical coordinates, our cone is described by \(r=z\) which suggests that we can use \(r\) and \(\theta\) as \(s\) and \(t\text{.}\) Consider the parameterization \(\vr(s,t)=\langle s \cos(t), s \sin(t), s \rangle \) with bounds \(1\leq s \leq 3\) and \(0 \leq t \lt 2 \pi\text{.}\) Notice that \(s\) is acting as both \(r\) and \(z\text{.}\) Calculating the partial derivatives of \(\vr\) gives the following:
\begin{equation*}
\vw=\langle-s \cos(t), -s \sin(t), s \rangle \text{.}
\end{equation*}
Notice that \(\vw\) is very close to \(\vr\text{.}\) If we think of \(\vr\) as giving \(\langle x,y,z\rangle \text{,}\) then \(\vw = \langle -x, -y, z \rangle \text{.}\) The figure below shows the relationship between vectors \(\vr\) and \(\vw\) for a point on the surface.
In Figure 12.9.10, the normal vector (as calculated from the parameterization) points inside the cone. This is the direction of positive flow when measuring the flux of the vector field. Our parameterization will also transform our vector field into a function of \(s\) and \(t\text{.}\) Specifically,
\begin{equation*}
\vF(s,t)= \langle (s \cos(t))(s) -2, s \sin(t) + s \cos(t) , 1\rangle
\end{equation*}
Applying Theorem 12.9.7 to our parameterization and vector field allows us to compute \(I = \iint_D \vF \cdot (\vr_s \times \vr_t)\, dA\) as an interated integral
\begin{align*}
I \amp = \int_0^{2\pi}\int_1^3 \langle (s \cos(t))(s) -2, s \sin(t) + s \cos(t) , 1\rangle \cdot \langle -s \cos(t), -s \sin(t), s \rangle \, ds \, dt \\
\amp =\int_0^{2\pi}\int_1^3 ((s \cos(t))(s) -2) (-s \cos(t)) +(-s \sin(t)) (s \sin(t) + s \cos(t)) +s \\
\amp = -\frac{62}{3}\pi \text{.}
\end{align*}
Note that the result of this flux integral is negative, which means more of the vector field is going out of the cone than is coming into the cone.
In Figure 12.9.11 you can see much more of the vector field is flowing outside of the cone than inside, which matches with our result of calculating the flux as negative.
The next activity asks you to carefully go through the process of calculating the flux of some vector fields through a cylindrical surface.
Activity12.9.3.Checking the Visualization for Flux.
(a)
Figure 12.9.12 shows a plot of the vector field \(\vF=\langle{y,z,2+\sin(x)}\rangle\) and a right circular cylinder of radius \(2\) and height \(3\) (with open top and bottom). Consider the vector field going into the cylinder (toward the \(z\)-axis) as corresponding to positive flux.
(i)
Reasoning graphically, do you think the flux of \(\vF\) throught the cylinder will be positive, negative, or zero? Write a few sentences jutifying your answer.
(ii)
Parameterize the right circular cylinder of radius \(2\text{,}\) centered on the \(z\)-axis, for \(0\leq z \leq 3\text{.}\) Be sure to give bounds on your parameters.
(iii)
Based on your parameterization, compute \(\vr_s\text{,}\)\(\vr_t\text{,}\) and \(\vr_s \times \vr_t\text{.}\) Confirm that these vectors are either orthogonal or tangent to the right circular cylinder. Is your orthogonal vector pointing in the direction of positive flux or negative flux?
(iv)
Use your parametrization to write \(\vF\) as a function of \(s\) and \(t\text{.}\)
Hint.
The \(x\)-coordinate is given by the first component of \(\vr\text{.}\)
(v)
Compute the flux of \(\vF\) through the parameterized portion of the right circular cylinder.
(vi)
Does your computed value for the flux match your prediction from earlier?
(vii)
Use Figure 12.9.13 to make an argument about why the flux of \(\vF=\langle{y,z,2+\sin(x)}\rangle\) through the right circular cylinder is zero.
(b)
Write a couple of sentences to explalin how the results of the flux calculations would be different if we used the vector field \(\vF=\left\langle{y,z,\cos(xy)+\frac{9}{z^2+6.2}}\right\rangle\) and the same right circular cylinder.
(c)
write a couple of sentences to explain how the results of the flux calculations would be different if we used the vector field \(\vF=\langle{y,-x,3}\rangle\) and the same right circular cylinder.
In the exercises for this section, we will look at some computational ideas to help us more efficiently compute the value of a flux integral. In many cases, the surface we are looking at the flux through can be written with one coordinate as a function of the others. For simplicity, we consider \(z=f(x,y)\text{.}\) Additionally, there will be exercises that guide you through common surfaces like spheres and cylinder surfaces.
Remark12.9.14.
Note that throughout this section, we have implicitly assumed that we can parametrize the surface \(S\) in such a way that \(\vr_s\times \vr_t\) gives a well-defined normal vector. Technically, this means that the surface be orientable. Most “reasonable” surfaces are orientable. However, there are surfaces that are not orientable. Perhaps the most famous is formed by taking a long, narrow piece of paper, giving one end a half twist, and then gluing the ends together. Try doing this yourself, but before you twist and glue (or tape), poke a tiny hole through the paper on the line halfway between the long edges of your strip of paper and circle your hole. After gluing, place a pencil with its eraser end on your dot and the tip pointing away. Think of this as a potential normal vector. Keep the eraser on the paper, and follow the middle of your surface around until the first time the eraser is again on the dot. Is your pencil still pointing the same direction relative to the surface that it was before?
Subsection12.9.3Summary
A flux integral of a vector field, \(\vF\text{,}\) on a surface in space, \(S\text{,}\) measures how much of \(\vF\) goes through \(S_1\text{.}\)
Let the smooth surface, \(S\text{,}\) be parametrized by \(\vr(s,t)\) over a domain \(D\text{.}\) The total flux of a smooth vector field \(\vF\) through \(S\) is given by
Compute the flux of the vector field \(\vec{v}= -4\,\mathit{\vec i}-5\,\mathit{\vec j}-4\,\mathit{\vec k}\) through the rectangular region shown below, assuming it is oriented as shown and that \(a = 2\) and \(b = 4\text{.}\)
flux =
2.
Calculate the flux of the vector field \(\vec{F}(x,y,z) = 3 y \vec{j}\) through a square of side length \(7\) in the plane \(y = 6\text{.}\) The square is centered on the y-axis, has sides parallel to the axes, and is oriented in the positive y-direction.
Flux =
Hint.
\(y = 6\text{.}\)
3.
(a) Set up a double integral for calculating the flux of the vector field \(\vec{F}(x,y,z) = z^2 \vec{k}\) through the upper hemisphere of the sphere \(x^2+y^2+z^2 = 9\text{,}\) oriented away from the origin. If necessary, enter \(\rho\) as rho,\(\theta\) as theta, and \(\phi\) as phi.
Compute the flux of the vector field \(\vec F = 5 x^2y^2z\vec k\) through the surface \(S\) which is the cone \(\sqrt{x^2+y^2}=z\text{,}\) with \(0\le z\le R\text{,}\) oriented downward.
(a) Parameterize the cone using cylindrical coordinates (write \(\theta\) as theta).
\(x(r,\theta) =\)
\(y(r,\theta) =\)
\(z(r,\theta) =\)
with \(\le r\le\)
and \(\le \theta \le\)
(b) With this parameterization, what is \(d\vec A\text{?}\)
\(d\vec A\) =
(c) Find the flux of \(\vec F\) through \(S\text{.}\)
flux =
5.
Compute the flux of the vector field \(\vec F = z\vec i + 6 x\vec j\) through the parameterized surface \(S\) oriented upward and given, for \(0\le s\le 3, \ \ 2\le t\le 3\text{,}\) by
\begin{equation*}
x=s^2, \quad y=2s+t^2, \quad z=5 t.
\end{equation*}
flux =
6.
Are the following statements true or false?
If \(S_1\) is a rectangle with area 1 and \(S_2\) is a rectangle with area 2, then \(\displaystyle 2 \iint_{S_1} \vec{F} \cdot d\vec{A}
= \iint_{S_2} \vec{F} \cdot d\vec{A}\)
The area vector of a flat, oriented surface is parallel to the surface.
If \(S\) is the unit sphere centered at the origin, oriented outward and \(\vec{F} = x \vec{i} + y \vec{j} + z\vec{k} = \vec{r}\text{,}\) then the flux integral \(\displaystyle \iint_S \vec{F} \cdot d\vec{A}\) is positive.
If \(S\) is the unit sphere centered at the origin, oriented outward and the flux integral \(\displaystyle \iint_S \vec{F} \cdot d\vec{A}\) is zero, then \(\vec{F}(x,y,z)\) is perpendicular to \(\vec{r} = \langle x,y,z \rangle\) at every point of \(S\text{.}\)
If \(S\) is the unit sphere centered at the origin, oriented outward and the flux integral \(\displaystyle \iint_S \vec{F} \cdot d\vec{A}\) is zero, then \(\vec{F} = \vec{0}\text{.}\)
If \(S\) is an open-ended circular cylinder centered about the z-axis, oriented away from the z-axis, and \(\vec{F} = \langle x,y,0 \rangle\text{,}\) then the flux of \(\vec{F}\) through \(S\) is positive.
If \(S\) is an open-ended circular cylinder centered about the z-axis, oriented away from the z-axis, and \(\vec{F} = \langle 3,-2,6 \rangle\text{,}\) then the flux of \(\vec{F}\) through \(S\) is zero.
If \(S\) is the cube bounded by the six planes \(x = \pm 3\text{,}\)\(y = \pm 3\text{,}\)\(z = \pm 3\text{,}\) oriented outward, and \(\vec{F} = 2 \vec{i} - \vec{k}\text{,}\) then the flux of \(\vec{F}\) through \(S\) is zero.
7.
Let \(\mathcal{S}\) be the square in the \(xy\)-plane shown in the figure below, oriented with the normal pointing in the positive \(z\)-direction. Estimate
F = \(x\boldsymbol{i}+y\boldsymbol{j}+5z\boldsymbol{k}\) and \(\sigma\) is the portion of the cone \(z = \sqrt{x^2+y^2}\) between the planes \(z = 5\) and \(z = 10\) oriented with normal vectors pointing upwards. Find the flux of the flow field F across \(\sigma\text{:}\)
9.Flux Through Surfaces of the Form \(z=f(x,y)\).
In this exercise, we will look at how to use a parametrization of a surface that can be described as \(z=f(x,y)\) to efficiently calculate flux integrals.
(a)
Suppose that \(S\) is a surface given by \(z=f(x,y)\text{.}\) Find a parametrization \(\vr(s,t)\) of \(S\text{.}\)
Hint.
Use \(s=x\) and \(t=y\text{.}\)
(b)
Show that the vector orthogonal to the surface \(S\) has the form
For each of the three surfaces given below, compute \(\vr_s \times \vr_t\text{,}\) graph the surface, and compute \(\vr_s \times \vr_t\) for four different points of your choosing. You should make sure your vectors \(\vr_s \times \vr_t\) are orthogonal to your surface.
(i)
\(z=x^2+y^2\)
(ii)
\(x+2y+z=-4\)
(iii)
\(z=x^2-y^2\)
(d)
For each of the three surfaces in part c, use your calculations and Theorem 12.9.7 to compute the flux of each of the following vector fields through the part of the surface corresponding to the region \(D\) in the \(xy\)-plane.
(i)
\(\vF=\langle{x,y,z}\rangle\) with \(D\) given by \(0\leq x,y\leq 2\)
(ii)
\(\vF=\langle{-y,x,1}\rangle\) with \(D\) as the triangular region of the \(xy\)-plane with vertices \((0,0)\text{,}\)\((1,0)\text{,}\) and \((1,1)\)
(iii)
\(\vF=\langle{z,y-x,(y-x)^2-z^2}\rangle\) with \(D\) given by \(0\leq x,y\leq 2\)
10.Calculating the Flux through a Sphere.
For this activity, let \(S_R\) be the sphere of radius \(R\) centered at the origin.
(a)
Parametrize \(S_R\) using spherical coordinates. Give your parametrization as \(\vr(s,t)\text{,}\) and be sure to state the bounds of your parametrization.
Hint.
Use \(s=\theta\) and \(t=\phi\text{.}\)
(b)
Use your parametrization of \(S_R\) to compute \(\vr_s \times \vr_t\text{.}\)
(c)
Your result for \(\vr_s \times \vr_t\) should be a scalar expression times \(\vr(s,t)\text{.}\) Explain why the outward pointing orthogonal vector on the sphere is a multiple of \(\vr(s,t)\) and what that scalar expression means.
(d)
Use your parametrization of \(S_2\) and the results of part b to calculate the flux through \(S_2\) for each of the three following vector fields.
(i)
\(\vF_1=\langle{x,y,z}\rangle\)
(ii)
\(\vF_2=\langle{-y,x,-1}\rangle\)
(iii)
\(\vF_3=\langle{x-y,y+x,z-1}\rangle\)
(e)
Use computer software to plot each of the vector fields from part d and interpret the results of your flux integral calculations.
(f)
If we used the sphere of radius 4 instead of \(S_2\text{,}\) explain how each of the flux integrals from part d would change. You do not need to calculate these new flux integrals, but rather explain if the result would be different and how the result would be different.
Subsection12.9.5Notes to Instructors and Dependencies
This section relies on parameterized surfaces, which was first introduced in Section 11.6. While some of the activities in this section may be too much for a single student to do in a class setting, we suggest that different cases of the vector fields and surfaces can be split for small group work. Students can then present answers to the larger group.