What is meant by rotation of a vector field in a plane?
How can a two-dimensional measurement of rotation be generalized to work in three dimensions?
How can the rotational strength of a vector field be measured?
In many physical settings, it is useful to measure the rotational strength of a vector field at a local scale. For instance, the vector field on the left of Figure 4.5.1 shows a vector field which flows in a counterclockwise fashion around the origin. The vector field on the right of Figure 4.5.1 shows a vector field which does not have a rotational aspect to its flow.
These global ideas of rotation are nice, but are not always visually apparent or may only appear in some regions. For instance, the vector field on the left of Figure 4.5.2 seems to have rotation around at least a couple of different points. In each quadrant, the vector field on the right of Figure 4.5.2 has different rotational patterns.
Preview Activity4.5.1.
We would like to understand and measure rotation of a vector field near a particular point. In order to investigate this concept, we will look at some two-dimensional vector fields and think about whether the vector field shown will rotate a small pinwheel or spinner placed at a particular location. The sort of spinner we imagine is illustrated in Figure 4.5.3. It consists of a central axis with a four-bladed paddle placed at one end of the axis. We imagine that the spinner is anchored at a point and the vector field, perhaps thought of as a fluid flow or wind velocity vector field, pushes against the blades of the spinner’s paddle.
(a)
To begin our investigation of the rotation of a spinner in a vector field, we will look at the vector field \(\vF\) in Figure 4.5.4. As you think about these questions, draw an “X” at each of the points about which you are asked and consider the vector field as being the pattern of a wind blowing across the plane.
(i)
Draw a vector on the top right blade of your spinner that represents how the wind will push on that blade. Next, draw a vector on each of the other blades of your spinner that represents how the wind will push on that blade.
(ii)
Use your vector representations from the previous part to describe if a small spinner placed at the origin would spin clockwise, counterclockwise, or not spin at all.
(iii)
Now we would like to look at the rotational strength of \(\vF\) at the point \((1,1)\text{.}\) As before, draw a spinner at this point and draw vectors on each of the blades to represent how the wind will push on that blade. It is important at this stage to draw the relative lengths of the vectors on each blade to scale so you can see which blades have a larger force due to the wind.
(iv)
Use your vector representations from the previous part to describe if a small spinner placed at at \((1,1)\) would spin clockwise, counterclockwise, or not spin at all. You should pay attention to which blades will have a stronger force (in comparison to the blade on the other side).
(v)
Would a small spinner placed at \((-2,1)\) spin clockwise, counterclockwise, or not spin at all? Draw a representation of a spinner if necessary to demonstrate your ideas.
(vi)
Are there any points on the \(xy\)-plane where the spinner would not turn counterclockwise?
(b)
For this part, look at the vector field \(\vG\) in Figure 4.5.5.
(i)
Would a small spinner placed at the origin spin clockwise, counterclockwise, or not spin at all?
(ii)
Would a small spinner placed at \((1,1)\) spin clockwise, counterclockwise, or not spin at all?
(iii)
Would a small spinner placed at \((-2,1)\) spin clockwise, counterclockwise, or not spin at all?
(iv)
Are there any points in the \(xy\)-plane where the spinner would turn differently than your descriptions in the previous part of this problem?
(c)
For this part, we consider the vector field \(\vH\) shown in Figure 4.5.6.
(i)
Would a small spinner placed at the origin spin clockwise, counterclockwise, or not spin at all?
(ii)
Would a small spinner placed at \((1,1)\) spin clockwise, counterclockwise, or not spin at all?
(iii)
Would a small spinner placed at \((-2,1)\) spin clockwise, counterclockwise, or not spin at all?
(iv)
Are there any points on the \(xy\)-plane where the spinner would turn differently than your descriptions in the previous part of this problem?
(d)
For this part, look at the vector field in Figure 4.5.7. You may find it useful to recenter the plot on your point by changing \(x_0\) and \(y_0\text{.}\) You may also zoom in on your point by using the slider labeled “zoom”.
(i)
Would a small spinner placed at the origin spin clockwise, counterclockwise, or not spin at all?
(ii)
Would a small spinner placed at \((0,2)\) spin clockwise, counterclockwise, or not spin at all?
(iii)
Would a small spinner placed at \((-2,1)\) spin clockwise, counterclockwise, or not spin at all?
(iv)
For each of the following spinner behaviors, give a point in the \(xy\)-plane (different than the ones used in the previous parts) where a spinner would:
Turn clockwise.
Turn counterclockwise.
Not turn at all.
Subsection4.5.1Measuring the Circulation Density of Vector Field in \(\R^2\)
In this subsection, we will develop the measurement of the circulation density for a two-dimensional vector field. We will use this measurement to generalize to a notion of rotational strength in higher dimensional cases in the next subsection.
We will start by measuring the circulation of a vector field on a path around the point \((a,b)\) and use this measurement to define circulation density. Specifically, we will measure the circulation of a vector field as we move around a square centered at \((a,b)\text{.}\) Using this measurement, we will calculate the circulation density by dividing our measurement by the area enclosed. This will allow us to compare our measurement across regions of different sizes. By taking the limit of this circulation density as the square’s side length goes to zero, we will have the circulation density at the point \((a,b)\text{.}\) Just as in our discussion of Subsection 4.9.1, we will look at a two dimensional setting first, then examine how our argument can be generalized to higher dimensions.
Let’s start by measuring the circulation around a square with side length \(2h\) centered at a point \((a,b)\text{.}\) Namely, we will look at the line integral of our vector field as we move along the square curve shown in Figure 4.5.8.
We can parametrize the top edge of the square using the parametrization
In order to measure the circulation, we want to know how much of the vector field \(\vF = \langle F_1,F_2\rangle\) is parallel to the direction of travel. Thus, along each of the sides of the box, we only need to look at one of the components of the vector field. The only contribution to the circulation on the top and bottom comes from the horizontal component \(F_1\text{.}\) Similarly, the vertical component \(F_2\) is all that contributes on the right and left sides. This means that we can simplify our line integrals considerably once we apply Theorem 4.3.2. For instance, the circulation on the top can be written as
Similarly, we can simplify the appropriate line integrals on the other sides (ordered below as top, left, bottom, right) to get the total circulation around the square to be
Note that the sign on the top and the left sides ends up being switched because the direction of travel on these sides is in the negative coordinate direction. It is important to note that by convention, we talk about “positive” rotation being counterclockwise and “negative” rotation being clockwise when looking at these planar graphs. This will be discussed more fully in Subsection 4.5.2.
As with our development of divergence, we apply the Mean Value Theorem to the four integrals above to simplify future calculations. Specifically, the Mean Value Theorem applied to the first integral tells us that there is a value to get a value \(a_{\text{top}}^*\) in the interval \((a-h,a+h)\) such that
Here we have \(a_{\text{top}}^*,a_{\text{bottom}}^* \in (a-h,a+h)\) and \(b_{\text{left}}^*,b_{\text{right}}^* \in (b-h,b+h)\text{.}\) We can regroup these terms to get the following expression for the total circulation around our square:
A larger square is likely to have a larger total circulation since there is more distance to accumulate how much the vector field is moving in the same direction as our path. This is noticeable in the formula above for total circulation, as the side length \(2h\) of the square appears in the formula. In order to compare our rotational or circulation ideas over different sizes of squares, we will now look at the circulation density (strength of rotation per unit area), which can be computed directly from our total circulation measurement. We can take the limit of the circulation density as we shrink the square to the central point \((a,b)\) and get a measurement rotational strength of our vector field at a point.
We will compute our circulation density as
\begin{equation*}
\text{Circulation Density}=\lim_{h\rightarrow 0} \frac{\text{Total Circulation}}{\text{Area of Square}}\text{.}
\end{equation*}
Since the area of the square is \((2h)^2 = 4h^2\text{,}\) this gives us that the circulation density is
Recall the central difference method of estimating derivatives from Section 1.5.2 1
activecalculus.org/single/sec-1-5-units.html#VjW
and notice that as \(h\to 0\text{,}\)\(a_{\text{top}}^*,a_{\text{bottom}}^*\) must go to \(a\) and \(b_{\text{left}}^*,b_{\text{right}}^*\) must go to \(b\text{.}\) Therefore, our circulation density measurement at a point is
\begin{equation*}
\text{Circulation Density at the point }(a,b)= \frac{\partial F_2}{\partial x} (a,b)-\frac{\partial F_1}{\partial y} (a,b)
\end{equation*}
You may have noticed that this argument had a similar flow to the development of divergence in Section 4.9. However, the result here is different because the property of the vector field we were trying to measure was different. Because divergence was developed to measure the change in the strength of the vector field (without regard to the direction), divergence is computed using the partial derivatives of the horizontal component of the vector field with respect to the horizontal variable and the vertical component of the vector field with respect to the vertical variable. The circulation density we just developed measures how the direction of the vector field is changing, and thus uses partial derivatives of the components with respect to the transverse variable. In other words, the rotational aspects of the vector field depend on how the horizontal component of the vector field changes when we move vertically (and vice versa).
We could have used shapes other than a square in the development of the circulation density, but the square allows for calculations that are easier to understand as the area of the shape decreases.
Activity4.5.2.Matching Visual Measurements to Algebraic Calculations.
Remember that the circulation density of a vector field \(\vF=\langle F_1,F_2\rangle\) at a point \((a,b)\) is calculated as
\begin{equation*}
\text{Circulation Density at the point }(a,b)= \frac{\partial F_2}{\partial x} (a,b)-\frac{\partial F_1}{\partial y} (a,b)\text{.}
\end{equation*}
The circulation density will be positive when a small spinner placed at \((a,b)\) will spin counterclockwise, negative when the spinner will move clockwise, and zero when the spinner does not rotate.
For each of the vector fields given below, answer the following questions:
What is the formula for the circulation density of the vector field?
For what points will you have a positive circulation density?
For what points will you have a negative circulation density?
For what points will you have a zero circulation density?
Are your answers to the above questions consistent with what you found in a part of Preview Activity 4.5.1? If so, which part and why?
(a)
\(\langle -y,x\rangle\)
(b)
\(\langle x,y\rangle\)
(c)
\(\langle 1-x,xy\rangle\)
Subsection4.5.2Measuring Rotation in Three Dimensions
The previous subsection showed how we can measure the circulation density, or strength of rotation, at a point for a two-dimensional vector field. In this subsection, we will look at how this two-dimensional measurement can be used to define circulation density for three dimensional vector fields.
We first need to consider how we want to represent and measure rotation in three dimensions. When we developed the circulation density on the \(xy\)-plane, we measured the rotational strength at a point \((a,b)\) with an axis of rotation coming out of the \(xy\)-plane. This corresponds to the axis of rotation being given by the blue vector in Figure 4.5.9. Remember that “positive” rotation corresponds to counterclockwise rotation in the \(xy\)-plane. We can generalize this idea to think of rotation in three dimensions as being represented by a vector where
the direction of the vector represents the axis of rotation and
the magnitude of the vector represents the strength of the rotation.
By convention, we consider positive “rotation” to correspond to counterclockwise rotation when the vector field is viewed looking from the terminal point of the vector to its initial point. In Figure 4.5.10, you can see that each vector corresponds to the rotation displayed and is consistent with the conventions described above.
If you look carefully at Figure 4.5.10, you can see that the vector shown in blue will be parallel to \(\vk\) and will represent rotations on planes of the form \(z=c\text{.}\) Similarly, the vector shown in yellow will be parallel to \(\vi\) and will represent rotations on planes of the form \(x=a\text{.}\) When looking down the blue and yellow vectors (from the terminal to the initial point), you can see the positive coordinate axes as pointing to the right and up 2
When placing the \(x\)-axis horizontally for the blue vector and placing the \(y\)-axis horizontally for the yellow vector
(as we would expect on a two-dimensional plot). In contrast, when you look down the magenta vector (from the terminal to the initial point), the positive coordinate axes (the \(x\) and \(z\) axes) point to the left and up. 3
When placing the \(x\)-axis horizontally.
In fact, when we look down the purple vector the \(xz\)-plane is flipped. A positive rotation on a plane of the form \(y=b\) will correspond to a rotation vector that is in the direction of \(-\vj\text{.}\) This is a consequence of the right-handed coordinate system and our right-handed idea of rotation, as reflected in the relationships amongst the vectors \(\vi,\vj\,vk\) recalled below:
With this idea of using a vector to represent the rotation of some element, we have enabled all of the tools that vectors allow, especially the ideas of projection and linear combinations. Remember that the projection of a vector \(\vv\) onto a vector \(\vu\) will give the vector component of \(\vv\) that is parallel to \(\vu\text{.}\) In other words, the projection of \(\vv\) onto \(\vu\) is how much of \(\vv\) is parallel to \(\vu\text{.}\) The projection formula will allow us to take a vector representation of rotation and determine how much rotation happens around a given axis of rotation.
Recall that we have used a set of coordinate vectors (namely \(\vi\text{,}\)\(\vj\text{,}\) and \(\vk\)) to write other vectors as a linear combination of our base vectors. In other words, a vector \(\langle a,b,c\rangle\) in \(\R^3\) can be written in the form \(a\vi+b\vj+c\vk\text{.}\)
The linear combination of vectors will allow us to build the total rotation vector by measuring the rotational strength in the direction of three coordinate direction vectors. Just as in Figure 4.5.11, we will build our total rotation vector as a sum of \(\vi\text{,}\)\(\vj\text{,}\) and \(\vk\text{.}\) A rotation vector in the direction of \(\vk\) at a point \((a,b,c)\) will measure the circulation density at our point restricted to the trace plane \(z=c\text{.}\) Similarly the rotation vectors in the directions \(\vi\) and \(\vj\) will measure the circulation density on the trace planes \(x=a\) and \(y=b\text{,}\) respectively.
Subsection4.5.3Circulation Density in Three Dimensions
The previous subsection discussed how we can view the amount of rotation of a three-dimensional vector field in a trace plane parallel to the \(xy\)-, \(xz\)-, or \(yz\)-plane and recalled that vectors in \(\R^3\) can be written as a linear combination of the vectors \(\vi\text{,}\)\(\vj\text{,}\) and \(\vk\text{.}\) We use our definition of curl for a two-dimensional vector field to measure the amount of rotation in the appropriate trace planes, which leads to the following definition.
Definition4.5.12.
Let \(\vF=\langle F_1,F_2,F_3\rangle\) be a three-dimensional vector field. The curl of \(\vF\) is given by
The curl is exactly the rotational description described at the end of Subsection 4.5.2. When evaluated at a point \((a,b,c)\text{,}\) the first component of the curl will measure the circulation density of the vector field restricted to the plane \(x=a\text{.}\) Similarly, the second and third components of the curl will measure the circulation density of the vector field restricted to the planes \(y=b\) and \(z=c\text{,}\) respectively.
Activity4.5.3.Estimating Curl in Three Dimensions.
(a)
Look at the vector field \(\vF\) plotted in Figure 4.5.13. You can adjust the size of the region around \((1,1,-2)\) over which the vector field is plotted using the “Zoom” slider. The “Density” slider allows you to adjust the number of vectors plotted. Can you see any rotation in the three dimensional vector field plot at the point \((1,1,-2)\text{?}\) Can you tell how a spinner placed at \((1,1,-2)\) would rotate? Along which axis? Try to state your answer as a vector representing the rotational strength of the vector field at \((1,1,-2)\text{.}\)
(b)
You likely found it difficult to decide how you thought a spinner might rotate in this new, three-dimensional setting. Let’s look at the vector field in the plane \(z=-2\text{,}\) as displayed in Figure 4.5.14. Do you think a spinner placed on the red point would rotate clockwise, counterclockwise, or not rotate? If the spinner will rotate, you should think about what the axis of rotation would be and whether the rotation should be positive or negative. Summarize your result as a vector representing the rotational strength of \(\vF\) in the plane \(z=-2\text{.}\)
(c)
Next we will look at the vector field in the plane \(y=1\text{,}\) as displayed in Figure 4.5.15. Do you think a spinner placed on the red point would rotate clockwise, counterclockwise, or not rotate? If the spinner will rotate, you should think about what the axis of rotation would be and whether the rotation should be positive or negative. Summarize your result as a vector representing the rotational strength of your vector field in the plane \(y=1\text{.}\) Do you think the rotation in this figure is stronger or weaker than in Figure 4.5.14?
(d)
Finally, consider the trace of \(\vF\) in the plane \(x=1\text{,}\) as displayed in Figure 4.5.16. Do you think a spinner placed on the red point would rotate clockwise, counterclockwise, or not rotate? If the spinner will rotate, you should think about what the axis of rotation would be and whether the rotation should be positive or negative. Summarize your result as a vector representing the rotational strength of your vector field in the plane \(x=1\text{.}\)
(e)
Summarize your prediction to what you think the three-dimensional rotational strength of the vector field will be at the point \((1,1,-2)\) in the form of three-dimensional vector.
Hint.
Add your results from the previous three steps.
(f)
Compute the curl of \(\vF=\langle
x-y,y+2z,x^2\rangle\text{.}\) Specifically, what is \(\curl(\vF)\) at the point \((1,1,-2)\text{?}\)
(g)
Compare the result of the curl calculation in part f to your prediction from part e. You likely found it difficult to estimate the magnitude, so your answer there may be incorrect. Hopefully, you did get the signs of the components and their relative strengths (i.e., which is biggest) correct. If you did not, go back and review the previous parts and explain why the calculated components match with the rotational strength for each of the three figures.
Subsection4.5.4Interpretation and Usage of Curl
The curl of a vector field is itself a vector field in that evaluating \(\curl(\vF)\) at a point gives a vector. As we saw earlier in this section, the vector output of \(\curl(\vF)\) represents the rotational strength of the vector field \(\vF\) as a linear combination of rotational strengths (or circulation densities) from two-dimensional planar descriptions. From our description of vectors as a representation of rotations, We can think of rotation in three dimensions as being represented by a vector where
the direction of the vector represents the axis of rotation and
the magnitude of the vector represents the strength of the rotation.
By convention, we consider positive “rotation” to correspond to counterclockwise rotation when the vector field is viewed looking from the end of the vector to the base.
With these conventions, the output vector of the curl evaluated at a point \(P\text{,}\) written as \(\curl(\vF)(P)\text{,}\) will have the following properties:
The direction of \(\curl(\vF)(P)\) will be the axis of rotation that has the strongest positive rotational strength.
The magnitude of \(\curl(\vF)(P)\) will be the circulation density of the vector field in the plane through \(P\) with normal vector \(\curl(\vF)(P)\text{.}\)
The following theorem allows us to use \(\curl(\vF)\) to measure the rotational strength of \(\vF\) to around an arbitrary axis. This will be particularly useful to us in later sections.
Theorem4.5.17.
The rotational strength of a vector field \(\vF\) around an axis given by a unit vector \(\vv\) at a point \(P=(a,b,c)\) can be computed by \(\left(\curl(\vF)(a,b,c)\right)\cdot \vv\text{.}\)
Proof.
This follows from the formula for the component of \(\curl(\vF)(a,b,c)\) along \(\vv\) and the fact that \(\vv\) was selected to be a unit vector.
Activity4.5.4.
In this activity, we will work on calculating curl algebraically and interpreting it.
(a)
Consider the vector field \(\vF=\langle x,y,z\rangle\text{.}\) If we plot this vector field in any plane through the origin, we will see the vector field shown in Figure 4.5.18. This two-dimensional vector field has no rotation. The projection of \(\curl(\vF)(0,0,0)\) onto any direction therefore must give the zero vector. The only vector that has a zero projection in every direction is the zero vector. Verify this geometric argument for the curl of \(\vF\) by doing the calculations necessary to show that \(\curl(\langle{x,y,z}\rangle)=\vzero\text{.}\)
(b)
Now consider the vector field \(\vF=\langle{-y,x,0}\rangle\text{,}\) which is shown in Figure 4.5.19. The figure includes both an interactive three-dimensional plot of \(\vF\) as well as a two-dimensional plot of \(\vF\) in the yellow plane, which can be adjusted using the “Angle of Plane” slider.
(i)
Based on the figure, would you estimate \(\curl(\vF)\) is positive, negative, or zero at the origin? Does your answer change if you pick a different point?
(ii)
Calculate \(\curl(\vF)\) algebraically. Does the curl of this vector field vary depending on the point at which the curl is measured?
(iii)
Explain why the rotational strength of the vector field on a titled plane through the origin will be given by \(2\cos(\theta)\text{,}\) where \(\theta\) is the angle between \(\vk\) and the normal vector of the tilted plane.
(c)
In Figure 4.5.20 you can plot a vector field in a region around a point of your choosing in order to look at the rotational properties of the vector field. The check box in Figure 4.5.20 will show the curl vector at the base point specified so you can make sense of your vector field and its curl.
Use the figure to estimate the direction of \(\curl(\langle
x-z,x^2+z,x+\sin(y)\rangle)\) at the point \((1,2,-1)\text{.}\) Confirm your estimate by calculating the curl at this point algebraically. Is there any point at which the direction of greatest rotational strength of this vector field has negative \(x\)-component? If there is, find such a point. If not, explain why not.
Subsection4.5.5Summary
The circulation of a vector field on a closed path measures how strong the vector field moves in the direction of travel for the path. The circulation density of a two-dimensional vector field \(\vF=\langle{F_1,F_2}\rangle\) is given by
\begin{equation*}
\text{Circulation Density at the point }(a,b)=
\frac{\partial F_2}{\partial x} (a,b)-\frac{\partial
F_1}{\partial y} (a,b)\text{.}
\end{equation*}
The circulation density can be visualized in terms of how fast a very small spinner anchored to \((a,b)\) will rotate as a result of the force of \(\vF\text{.}\)
By considering the circulation density of a three-dimensional vector field \(\vF = \langle
F_1,F_2,F_2\rangle\) in planes parallel to the coordinate planes, we can calculate the curl of \(\vF\) as
For a three-dimensional vector field \(\vF\text{,}\) the vector \(\curl(\vF)(a,b,c)\) points in the direction of greatest rotational strength at the point \((a,b,c)\text{.}\) The magnitude of this vector measures the strength of the rotation. The strength of rotation around an axis determined by a unit vector \(\vv\) is found by calculating \(\left(\curl(\vF)(a,b,c)\right)\cdot \vv\)