In Section 10.4, we defined \(\vT\) and \(\vN\) to measure the direction of travel and the direction of turning for an oriented curve in space. Both \(\vT\) and \(\vN\) are unit vectors that are used to capture direction aspects of motion. In the next couple of sections, we will look at a few different ways to measure the magnitude of various aspects of our motion and turning on a curve in space.
Given a parameterization of a curve in space, the speed will give us an idea about how fast an object moves along the curve at any given parameter value. We learned in the previous section that \(\vN\) tells us the direction of turning, but we have not yet developed a way to measure how much turning is happening at a given location. For example, consider the three curves shown in Figure 10.5.1. We would like to be able to measure how much turning is happening at the indicated point on each curve. Right now, you probably intuitively would say there is no turning in the first curve, a small amount in the second, and considerably more in the third. However, how can we assign a number to measure this amount of turning?
In this section, we focus on measuring how quickly a curve is turning. Keeping our running example in mind, we can think of this as measuring a property of the road. How quickly a curve is turning is a different measurement from how fast an object is turning as it moves along a path. That is a property of the driver to which we will return in Section 10.6.
As Chief Engineer and CEO at Steer Clear, you have done great work in transforming information from your location tracking system (LTS) into information on the position, velocity, and distance traveled by your car. In order to make the “self-driving” part of the self-driving car, you will need to compare information from the LTS (which describes how your car is moving) to GPS location data for the network of roads (which describe the paths your car should be taking). Since everyone on the road drives differently, you realize that you need to measure how quickly a stretch of road is turning in different places, which leads you to contact a friend who is a civil engineer for help understanding highway construction specifications.
Figure 10.5.2 shows a map of a section of road on your testing route for your self-driving car with points labeled \(P_0\) and \(P_1\) which are 10 meters apart (along the road). Draw a vector in the direction of travel at both \(P_0\) and \(P_1\text{.}\) Use these vectors to describe how the direction of travel is changing along the path from \(P_0\) to \(P_1\text{.}\)
Figure 10.5.3 shows a map of a section of road on your testing route for your self-driving car with points labeled \(Q_0\) and \(Q_1\) which are 10 meters apart (along the road). Draw a vector in the direction of travel at both \(Q_0\) and \(Q_1\text{.}\) Use these vectors to describe how the direction of travel is changing along the path from \(Q_0\) to \(Q_1\text{.}\)
Is the direction of travel changing faster over the path from \(P_0\) to \(P_1\) or over the path from \(Q_0\) to \(Q_1\text{.}\) Write a few sentences to explain your reasoning and connect to your arguments and plots for the first two tasks.
The Preview Activity showed how looking at the change in \(\vT\text{,}\) the direction of travel, for small steps in terms of arc length along our curve will allow us to measure the rate at which a curve is turning. This is different than measuring how quickly an object traveling along the path needs to turn in order to stay on the path. The arguments in the Preview Activity depend only on intrinsic geometric properties of the path (road) rather than how an object (car) moves along that path.
Hopefully it is not surprising that the unit speed parameterization of a curve \(\vr_1(s)\text{.}\) One way to think about this is that if two drivers follow the same path, each at a speed of exactly 1 meter per second, then for any point on the path, the two drivers both reach that point in the same amount of time. For instance, it will take both drivers exactly 15 seconds to reach the point that is 15 meters down the path. (Remember the most useful property of a unit speed parameterization is that the parameter value for \(\vr_1\) is equal to the numeric value of the distance traveled along the curve, albeit with different units.) Conceptually, you can see how using the unit speed parameterization allows us to measure geometric properties of the curve independent of how an object moves along the path.
We have already explored how \(\vN\) points in the direction of turning. Now, however, we are looking for a scalar quantity that measures how quickly a curve turns. To measure how quickly a path turns, we examine how quickly the direction of travel changes in terms of the unit speed parameterization. We apply the classic calculus approach to this measurement in order to obtain the exact value of this measurement. Suppose \(\vr_1(s)\) is the unit speed parameterization of a curve \(C\) and consider a section of \(C\) from \(\vr_1(s_0)\) to \(\vr_1(s_0+\Delta s)\text{.}\) We first consider how different the direction of travel is over this curve of length \(\Delta s\text{.}\) (Remember, the change in parameter value here corresponds to the arc length along the curve because \(\vr_1\) is a unit speed parameterization.) A plot a sample curve and two points separated by an arc length of \(\Delta s\) is shown in Figure 10.5.4.
A two-dimensional curve with points corresponding to the tips of \(\vr_1(s_0)\) and \(\vr_1(s_0+\Delta s)\) marked. Also depictred are the unit tangent vectors to the curve at these points.
We want to know how fast the direction of travel is changing at the point \(\vr_1(s_0)\text{.}\) To do this, we approximate the change in the direction of travel on an interval around this point and take the limit as the size of the interval goes to zero. We compute the difference between the directions of travel at the beginning and end of our interval by putting both \(\vT(s_0+\Delta s)\) and \(\vT(s_0)\) starting at the same location. This allows us to use vector subtraction to find \(\vT(s_0+\Delta s)-\vT(s_0)\text{,}\) which is the change in the direction of travel over the given interval. This vector subtraction is shown using the triangle of vectors in Figure 10.5.5.
The unit tangent vectors from Figure 10.5.4 drawn with their initial points at the same location. The vector \(\vT(s_0+\Delta s)-\vT(s_0)\) pointing from the tip of \(\vT(s_0)\) to the tip of \(\vT(s_0+\Delta s)\) is also shown.
So for a step of size \(\Delta s\) in arc length along our curve, we have a change in direction of \(\vT(s_0+\Delta s)-\vT(s_0)\) in the direction of travel. In order to approximate the rate of change for the direction of travel, we will need to look at the ratio
The unit normal vector \(\vN\) points in the direction of turning. We want to measure as a scalar how quickly the direction of turning is changing along our path. To do this, we compute the magnitude of this difference quotient as an approximation of rate of change for the direction of travel on the interval \([s_0,s_0+\Delta s]\text{.}\) This gives
which completes the first step of the classic calculus approach. We now consider smaller step sizes in the parameter, which corresponds to smaller changes in arc length along the curve because we have a unit speed parameterization. In doing so, we obtain a better approximation of how fast the direction of travel is changing along the path at the point \(\vr_1(s_0)\text{.}\)
Taking the limit as \(\Delta s \rightarrow 0\text{,}\) completes step three of the classic calculus approach to measure how quickly a curve is turning at the point \(P=\vr_1(s_0)\text{.}\) Calculating, we have
If \(C\) is a smooth space curve and \(s\) is the parameter for the unit speed parameterization of \(C\text{,}\) then the curvature, \(\kappa\text{,}\) of \(C\) is
Remember that the parameter \(s\) for the unit speed parameterization \(C\) is the arc length traveled along the curve. Curvature is denoted by \(\kappa\text{,}\) the lowercase Greek letter “kappa”.
In this activity we will use our previous work in finding unit speed parameterizations of lines and circles to make sense of the definition of curvature.
Recall that in Example 10.3.7 we found that the unit speed parameterization of a line through the points \((1,2,3)\) and \((2,0,5)\) to be \(\vr_1(t)= \langle 1,2,3\rangle + \frac{t}{3} \langle 1,-2,2 \rangle\text{.}\) Calculate \(\vT(s)\text{,}\)\(\frac{d\vT}{ds}\text{,}\) and \(\vecmag{ \frac{d\vT}{ds} }\) for the unit speed parameterization of this line.
Write a few sentences to explain why your result for the calculation of the curvature \(\kappa=\vecmag{ \frac{d\vT}{ds} }\) for this line makes sense. Exercise 8 asks you to do the calculations to confirm this property of an arbitrary line.
Recall that we concluded Activity 10.3.2 by noting giving a parameterization of circle in 2-space of radius \(R\) centered at the origin as \(\langle R\cos(t),R\sin(t)\text{.}\) A quick calculation can verify that this parameterization has constant speed and that speed is \(R\text{.}\) Modifying this parameterization by multiplying the parameter by \(\frac{1}{\text{speed}}=\frac{1}{R}\) gives the following parameterization:
\begin{equation*}
\vr(s) = \left\langle R \cos\left(\frac{s}{R}\right), R \sin\left(\frac{s}{R}\right)\right\rangle
\end{equation*}
Compute the speed of this parameterization to verify that this has unit speed. This will also show that \(\vr\, '(s)=\vv(s)=\vT(s)\text{.}\) (Recall that this fact is not true in general!)
Calculate \(\vT(s)\text{,}\)\(\frac{d\vT}{ds}\text{,}\) and \(\vecmag{ \frac{d\vT}{ds} }\) for the unit speed parameterization of a circle of radius \(R\) centered at the origin.
Write a few sentences to explain why your result for the calculation of \(\kappa=\vecmag{ \frac{d\vT}{ds} }\) of circle will be constant. You should also address why circles with larger radii will have have smaller curvature.
Notice that a unit speed parameterization will have the property that \(\vr_1\, '(s)=\vv(s)=\vT(s)\text{.}\) While this fact is not true in general, it does make the calculation and interpretation of \(\vT\) as simple as possible. It is also helpful to observe that a unit speed parameterization has the property that \(\vT(s) = \va(s)\text{,}\) the acceleration vector of the curve. Thus, when we have a unit speed parameterization for a path, we can think of the curvature of the path at a point as being the magnitude of the acceleration vector.
The definition of curvature relies on our ability to parameterize curves with unit speed. Since finding unit speed parametrization is frequently difficult to compute directly, the following theorem gives a way of computing the curvature of a curve in terms of any parametrization \(\vr(t)\text{.}\)
If \(\vr\) is a parameterization of a curve \(C\text{,}\) and if \(\vr\, '(t)\neq \vzero\) and \(\vr\, ''(t)\) exists, then the curvature \(\kappa\) of \(C\) can be calculated by
Let \(\vr(t)\) be some parameterization of the curve \(C\text{.}\) To calculate the curvature, \(\vT\) must have an input corresponding to the parameter of the unit speed parameterization. This means that the input of \(\vT\) must measure the arc length of \(C\) to the point of interest. We will use a function \(s(t)\text{.}\) The input to this function is the parameter of \(\vr(t)\text{,}\) and the output of \(s(t)\) is the arc length on \(C\) to the point \(\vr(t)\text{.}\) Note that \(s(t)\) is computed using Theorem 10.3.5 and that the derivative \(\frac{ds}{dt}\) will be \(\vecmag{\vr\, '(t)}\text{.}\)
We can write \(\vT\) as a function of the parameter \(t\) by the composition \(\vT(t)=\vT(s(t))\text{.}\) Taking the derivative of this composition with respect to \(t\) uses a chain rule argument:
There is also a nice conceptual interpretation of the equation above. Specifically, if we want to convert a measurement such as curvature or the derivative of \(\vT\) between the unit speed parameterization and any other parameterization, we adjust by the speed. Equation (10.5.2) shows that the adjustment between the derivative of \(\vT\) with respect to the unit speed parameter and another parameter is to divide by the speed, which may vary depending on time (or position).
We omit the proof of (10.5.1) because the argument involves technical algebraic ideas that do not offer geometric insight into what curvature measures. Exercise 6 guides you through this calculation.
Soon we will turn to an example of calculating curvature and then some activities for you to practice. Frist, however, we consider how we can use (10.5.1) to make sense of when curvature is undefined (or does not make sense). From an algebraic perspective, the curvature cannot be defined when the speed is zero. This should make sense because if you are not moving along the curve, as indicated by having speed equal to zero, then it does not make sense to measure how quickly the curve is turning. In Section 10.6, we consider why \(\vecmag{ \vv(t) \times \va(t) }\) being zero means that there will be zero curvature. This is different than the curvature being undefined!
In this example, we will calculate the curvature of the parabola given by \(y=x^2\) for \(-2 \lt x \lt 2\text{.}\) We would like to use (10.5.1) to calculate the curvature, so we must parameterize this curve in three dimensions with the \(z\)-component constantly zero. This is because the cross product is only defined for three-dimensional vectors. The vector-valued function \(\vr(t)=\langle t,t^2,0\rangle\) with \(-2\leq t\leq 2\) will parameterize the curve given by \(y=x^2\) for \(-2 \lt x \lt 2\text{.}\)
Comparing the points \((0,0)\) and \((-2,4)\) on Figure 10.5.9, you can see that the path bends more at \((0,0)\) than at \((-1,1)\text{,}\) which corresponds to the parabola having larger curvature at \((0,0)\) than at \((-1,1)\text{.}\) Using the formula we obtained above, we can see that at \((0,0)\text{,}\) when \(t=0\text{,}\) the curvature is \(\kappa(0)=2\) while at \((-1,1)\text{,}\) when \(t=-1\text{,}\) the curvature is \(\kappa(-1)=\frac{2}{(5)^{3/2}}\approx0.179\text{.}\) This calculation of curvature matches with our visual description from above.
Suppose at the instant \(t=x=0\) you were driving along the path given by \(y=x^2\) at unit speed and your steering wheel locked up preventing you from changing how you were turning. Because the steering wheel is held fixed, you will no longer stay on the path given by \(y=x^2\text{.}\) Instead, your car will drive in a circle, and we can figure out what the radius of that circle would be by using curvature.
In Activity 10.5.2, we saw that the curvature of a circle of radius \(R\) is \(\frac{1}{R}\text{,}\) so it hopefully is not too surprising that the radius of the circle in which the car would match. In fact, the radius of this circle will be \(\frac{1}{\kappa(0)}=\frac{1}{2}\text{.}\) In general, the radius of curvature is defined to be \(\frac{1}{\kappa}\) and the circle that is tangent to the curve at the point where the curvature is measured with radius \(\frac{1}{\kappa}\) is called the osculating circle. The osculating circle can be described as the circle that best approximates the curve at the point since the osculating circle goes through the point and has the same curvature as the original path.
In Figure 10.5.10, you can see a plot of the path \(y=x^2\) with the osculating circle drawn at a particular point. Adjust the slider to change the parameter value (and thus point) at which the osculating circle is drawn. Notice how the radius of the circle becomes larger as the point moves farther away from the vertex in either direction. Remember that curvature and the radius of curvature have an reciprocal relationship: large curvature means small radius and small curvature means large radius.
Consider the ellipse shown below. We have omitted the scale on the axes because we want to think of this as the general ellipse \(\frac{x^2}{a^2} + \frac{y^2}{b^2} = 1\text{,}\) which has parameterization
For the given parameterization of the ellipse, find \(\vv(t)\) and \(\va(t)\text{.}\) Also identify the parameter values that correspond to the points of the ellipse with the largest and smallest \(x\) or \(y\) coordinates.
Using your calculations for \(\vv(t)\) and \(\va(t)\) from part b in Equation (10.5.1) to find a formula for the curvature of the ellipse with the given parameterization. Use your formula to find the curvature of the ellipse at the parameter values you identified in part b as corresponding to the extreme points of the ellipse.
While it is possible to use single-variable calculus techniques to find the parameter values that maximize and minimize \(\kappa(t)\) on the interval \([0,2\pi]\text{,}\) the algebra involved masks much of the insight. Use graphing technology such as Desmos to plot the curvature function you found above, preferably with sliders that allow you to adjust the values of \(a\) and \(b\text{.}\) Write a couple of sentences about what you observe about the parameter values that maximize and minimize curvature and how they depend on \(a\) and \(b\text{.}\) Then use the general form of the curvature function from above to find the curvature at these parameter values, simplifying as much as possible.
The standard helix has parameterization \(\vr_2(t) = \cos(t) \vi + \sin(t) \vj + t \vk\text{.}\) Find the curvature of the helix using this parameterization.
Before concluding this section, it may be helpful to think about what the units are on curvature. To do this, we again consider the formula that is in terms of velocity and acceleration:
If we assume that our parameterization has distance with units of meters and input is thought of as time measured in seconds, then the units on \(\vv\) are meters per second and the units on \(\va\) are meters per second-squared. This means the numerator of \(\kappa\) has units meters per second-cubed and the denominator is measured in meters-cubed per second-cubed. Simplifying, we find that curvature has units one over meters. This aligns with the idea of the radius of curvature being \(1/\kappa\text{,}\) as it results radius of curvature having meters as its unit. That makes sense for something that we call a radius!
In this activity we will look at the measurement of curvature and determine if curvature is a property of the driver or a property of the road. This is a continuation of Activity 10.3.3. As a reminder, a measurement is a property of the driver if the value(s) of that measurement can be different for different drivers (when measured at the same location on the racetrack). A measurement is a property of the road when different drivers must have the same value(s) (when measured at the same location on the racetrack).
Is curvature a property of the driver or the road? Write a few sentences about your reasoning and be sure to address why every driver may or may not have the same measurement at the same location on the track.
What is the curvature on a straightaway (the race track is a line segment)? What would the radius of curvature be for a straightaway? Write a couple of sentences to justify your answers.
We define the curvature \(\kappa\) of a curve in 2- or 3-space to be the rate of change of the magnitude of the unit tangent vector with respect to arc length, or
The WeBWorK problems are written by many different authors. Some authors use parentheses when writing vectors, e.g., \((x(t),y(t),z(t))\) instead of angle brackets \(\langle x(t),y(t),z(t) \rangle\text{.}\) Please keep this in mind when working WeBWorK exercises.
A factory has a machine which bends wire at a rate of 6 unit(s) of curvature per second. How long does it take to bend a straight wire into a circle of radius 8?
Consider the single variable function defined by \(y = 4x^2 - x^3.\)
Find a parameterization of the form \(\vr(t) = \langle x(t), y(t) \rangle\) that traces the curve \(y = 4x^2 - x^3\) on the interval from \(x = -3\) to \(x = 3\text{.}\)
Write a definite integral which, if evaluated, gives the exact length of the given curve from \(x = -3\) to \(x = 3\text{.}\) Why is the integral difficult to evaluate exactly?
Use appropriate technology to approximate the absolute maximum and minimum of \(\kappa(t)\) on the parameter interval for your parameterization. Compare your results with the graph of \(y = 4x^2 - x^3\text{.}\) How do the absolute maximum and absolute minimum of \(\kappa(t)\) align with the original curve?
In Exercise 10.4.6.4 we showed that \(\lvert \vT(t) \rvert = 1\) implies that \(\vT(t)\) is orthogonal to \(\vT'(t)\) for every value of \(t\text{.}\) Explain what this tells us about \(\lvert \vT(t) \times \vT'(t) \rvert\) and conclude that
In this exercise we explore how to find the osculating circle for a given curve. As an example, we will use the curve defined by \(f(x) = x^2\text{.}\) Recall that this curve can be parameterized by \(x(t) = t\) and \(y(t)=t^2\text{.}\)
Use Definition 10.4.3 to find \(\vT(t)\) for our function \(f\text{.}\)
To find the center of the osculating circle, we will want to find a vector that points from a point on the curve to the center of the circle. Such a vector will be orthogonal to the tangent vector at that point. Recall that \(\vT(s) = \langle \cos(\phi(s)), \sin(\phi(s)) \rangle\text{,}\) where \(\phi\) is the angle the tangent vector to the curve makes with a horizontal vector. Use this fact to show that
Explain why this tells us that \(\frac{dT}{ds}\) is orthogonal to \(\vT\text{.}\) Let \(\vN\) be the unit vector in the direction of \(\frac{dT}{ds}\text{.}\) The vector \(\vN\) is called the principal unit normal vector and points in the direction toward which the curve is turning. The vector \(\vN\) also points toward the center of the osculating circle.
Find \(\vT\) at the point \((1,1)\) on the graph of \(f\text{.}\) Then find \(\vN\) at this same point. How do you know you have the correct direction for \(\vN\text{?}\)
Let \(P\) be a point on the curve. Recall that \(\rho = \frac{1}{\kappa}\) at point \(P\) is the radius of the osculating circle at point \(P\text{.}\) We call \(\rho\) the radius of curvature at point \(P\text{.}\) Let \(C\) be the center of the osculating circle to the curve at point \(P\text{,}\) and let \(O\) be the origin. Let \(\mathbf{\gamma}\) be the vector \(\overrightarrow{OC}\text{.}\) See Figure 10.5.11 for an illustration using an arbitrary function \(f\text{.}\)
Finally, use the previous work to find the center of the osculating circle for \(f\) at the point \((1,1)\text{.}\) Draw pictures of the curve and the osculating circle to verify your work.
In Activity 10.5.2, you verified for a specific line that the curvature is \(0\text{.}\) The generalization of this is that the unit speed parameterization of a line through \((x_0,y_0,z_0)\) with direction vector \(\langle a,b,c\rangle\) can be given by
Compute the speed of this parameterization to verify that this has unit speed. This will also show that \(\vr\, '(s)=\vv(s)=\vT(s)\) (which is not true in general).