Module

for

Curvature

     

The Circle of Curvature: It's a Limit!
by
John H. Mathews

The AMATYC Review, Vol. 25, No. 1, Fall 2003, pp. 57-63.

 

    Given the function  [Graphics:Images/CurvatureMod_gr_1.gif], the formula for the curvature (and radius of curvature) is stated in all calculus textbooks  

Definition (Curvature)            [Graphics:Images/CurvatureMod_gr_2.gif],  
          
Definition (Radius of Curvature)        [Graphics:Images/CurvatureMod_gr_3.gif].  

Definition (Osculating Circle)  At the point {x,f[x]} on the curve  y = f[x],  the osculating circle is tangent to the curve and has radius  r[x].   

 

Example 1.  Consider the parabola  [Graphics:Images/CurvatureMod_gr_4.gif] and the point {0,f[0]} = (0,0) on the curve.  Find the radius of curvature and the circle of curvature.
Solution 1.

 

Finding Curvature at Any Point

    For the above example the circle of curvature was easy to locate because it's center lies on the y-axis.  How do you locate the center if the point of tangency is not the origin?  To begin, we need the concepts of tangent and normal vectors.

Tangent and Normal Vectors

    Given the graph  [Graphics:Images/CurvatureMod_gr_17.gif],  a vector tangent to the curve at the point  (x,f(x))  is   [Graphics:Images/CurvatureMod_gr_18.gif].  The unit tangent vector is  [Graphics:Images/CurvatureMod_gr_19.gif],  which can be written as

Definition (Unit Tangent)        [Graphics:Images/CurvatureMod_gr_20.gif]
    
For vectors in  [Graphics:Images/CurvatureMod_gr_21.gif],  a corresponding perpendicular vector called the unit normal vector is given by

Lemma (Unit Normal)         [Graphics:Images/CurvatureMod_gr_22.gif].

Verification.

 

Example 2.  Consider the parabola  [Graphics:Images/CurvatureMod_gr_28.gif] and the point  (1,f(1)) = (1,1)  on the curve.  
Find the unit tangent and unit normal at point (1,1).
Solution 2.

 

Example 3.  Consider the parabola  [Graphics:Images/CurvatureMod_gr_34.gif] and the point  (1,f(1)) = (1,1)  on the curve.  
Find the radius of curvature and the circle of curvature.
Solution 3.

 

A new construction of the Circle of Curvature

        What determines a circle?  A center and a radius.  The formula for the radius of curvature is well established.  What idea could we use to help understand the situation.  We could use the fact that three points determine a circle and see where this leads.  

 

Example 4.  Consider the parabola  [Graphics:Images/CurvatureMod_gr_61.gif] and the point  (0,f(0)) = (0,0)  on the curve.  
Find collocation circle to go through the three points  [Graphics:Images/CurvatureMod_gr_62.gif],  [Graphics:Images/CurvatureMod_gr_63.gif],  and  [Graphics:Images/CurvatureMod_gr_64.gif], and explore the situation for h = 1,.1,.01.
Solution 4.

 

Derivation of the Radius of Curvature

    The standard derivation of the formula for radius curvature involves the rate of change of the unit tangent vector.  This new derivation starts with the collocation the collocation circle to go through the three points  [Graphics:Images/CurvatureMod_gr_103.gif],  [Graphics:Images/CurvatureMod_gr_104.gif],  and  [Graphics:Images/CurvatureMod_gr_105.gif]  on the curve  [Graphics:Images/CurvatureMod_gr_106.gif].  The limit as [Graphics:Images/CurvatureMod_gr_107.gif] is the osculating circle to the curve  [Graphics:Images/CurvatureMod_gr_108.gif] at the point [Graphics:Images/CurvatureMod_gr_109.gif].  The radius of curvature and formulas for the location of its center are simultaneously derived.  The computer algebra software Mathematica is used to assist in finding the limits.   

    Start with the equation [Graphics:Images/CurvatureMod_gr_110.gif],  of a circle.  Then write down three equations that force the collocation circle to go through the three points  [Graphics:Images/CurvatureMod_gr_111.gif],  [Graphics:Images/CurvatureMod_gr_112.gif],  and  [Graphics:Images/CurvatureMod_gr_113.gif]  on the curve  [Graphics:Images/CurvatureMod_gr_114.gif].  Enter the equations into Mathematica

[Graphics:Images/CurvatureMod_gr_115.gif]

[Graphics:Images/CurvatureMod_gr_116.gif]

 

Solve the equations for [Graphics:Images/CurvatureMod_gr_117.gif] and extract the formula for the radius of the collocation circle.  Since it depends on [Graphics:Images/CurvatureMod_gr_118.gif] we will store it as the function [Graphics:Images/CurvatureMod_gr_119.gif].    
[Graphics:Images/CurvatureMod_gr_120.gif]

 



[Graphics:Images/CurvatureMod_gr_121.gif]

 

The formula looks bewildering and one may wonder if it is of any value.  
However, we can demonstrate the power of Mathematica and see if it can take the limit.  

[Graphics:Images/CurvatureMod_gr_122.gif]


[Graphics:Images/CurvatureMod_gr_123.gif]

 

The computer gets the correct formula, but leaves out all the human understanding.  The formula for  [Graphics:Images/CurvatureMod_gr_124.gif] can be rewritten so that human insight and inspiration is involved, but the computer cannot think of the formulas we desire to see.  So human intervention must be called upon and the simplifications must be typed in by hand.   

[Graphics:Images/CurvatureMod_gr_125.gif]


[Graphics:Images/CurvatureMod_gr_126.gif]


[Graphics:Images/CurvatureMod_gr_127.gif]

Notice that when we take the limit as [Graphics:Images/CurvatureMod_gr_128.gif], the limiting value for each of the radicals in the numerator is  [Graphics:Images/CurvatureMod_gr_129.gif].  Mathematica is capable of finding them, we illustrate this with the third one.

[Graphics:Images/CurvatureMod_gr_130.gif]



[Graphics:Images/CurvatureMod_gr_131.gif]

Therefore, the limit in the numerator is [Graphics:Images/CurvatureMod_gr_132.gif].  The difference quotient in the denominator is recognized as the numerical approximation formula for the second derivative, hence the is  [Graphics:Images/CurvatureMod_gr_133.gif].  

[Graphics:Images/CurvatureMod_gr_134.gif]



[Graphics:Images/CurvatureMod_gr_135.gif]

We already knew that Mathematica knows the rules for taking limits of functions when the formula is given.  Now we know that it has the "artificial intelligence" to rearrange quantities involving an arbitrary function and can identify difference quotients for approximating derivatives and find their appropriate limits.  It takes quite a bit of trust to let Mathematica do our thinking.  From the steps we filled in, we can gain trust for the computer's way of doing it.  

 

The Osculating Circle

We now show that the limit of the collocation circle as [Graphics:Images/CurvatureMod_gr_136.gif] is the osculating circle.
Now we want to find the center  (a, b)  of the osculating circle.

The abscissa for the center of the collocation circle is

[Graphics:Images/CurvatureMod_gr_137.gif]


[Graphics:Images/CurvatureMod_gr_138.gif]

 

Take the limit as  [Graphics:Images/CurvatureMod_gr_139.gif],  to obtain the abscissa for the center of the circle of curvature.  

[Graphics:Images/CurvatureMod_gr_140.gif]


[Graphics:Images/CurvatureMod_gr_141.gif]

 

If we want to see "what's happening" in this limit, then we must rearrange the formula for  [Graphics:Images/CurvatureMod_gr_142.gif].  A little finesse permits us to write it as follows

[Graphics:Images/CurvatureMod_gr_143.gif]


[Graphics:Images/CurvatureMod_gr_144.gif]

 

The numerators involve three difference quotients, all of which tend to [Graphics:Images/CurvatureMod_gr_145.gif] when  [Graphics:Images/CurvatureMod_gr_146.gif] , and the difference quotient in the denominators tends to [Graphics:Images/CurvatureMod_gr_147.gif]when  [Graphics:Images/CurvatureMod_gr_148.gif].  

Thus we have established the formula  [Graphics:Images/CurvatureMod_gr_149.gif] for the abscissa of center of the circle of curvature.

 

The Abscissa the Easy Way

Subtract from  [Graphics:Images/CurvatureMod_gr_150.gif] the radius of curvature times [Graphics:Images/CurvatureMod_gr_151.gif].

[Graphics:Images/CurvatureMod_gr_152.gif]


[Graphics:Images/CurvatureMod_gr_153.gif]

 

The Ordinate for the Center of the Circle of Curvature

The ordinate for the center of the collocation circle is

[Graphics:Images/CurvatureMod_gr_154.gif]


[Graphics:Images/CurvatureMod_gr_155.gif]

 

If we want to see "what's happening" in this limit, then we must rearrange the formula for  [Graphics:Images/CurvatureMod_gr_156.gif].  A little finesse permits us to write it as follows

[Graphics:Images/CurvatureMod_gr_157.gif]

The limit of  [Graphics:Images/CurvatureMod_gr_158.gif] as  [Graphics:Images/CurvatureMod_gr_159.gif] is merely [Graphics:Images/CurvatureMod_gr_160.gif],  and we have already observed that the limit in the denominator is [Graphics:Images/CurvatureMod_gr_161.gif]and both quotients in the numerator tend to  [Graphics:Images/CurvatureMod_gr_162.gif] as  [Graphics:Images/CurvatureMod_gr_163.gif] .  

Thus we have established the formula  [Graphics:Images/CurvatureMod_gr_164.gif] for the ordinate of center of the circle of curvature.
Details

 

The Ordinate of the Circle the Easy Way

Add to  [Graphics:Images/CurvatureMod_gr_179.gif]  the radius of curvature times[Graphics:Images/CurvatureMod_gr_180.gif]

[Graphics:Images/CurvatureMod_gr_181.gif]


[Graphics:Images/CurvatureMod_gr_182.gif]

 

The Osculating Circle

At the point  [Graphics:Images/CurvatureMod_gr_183.gif] on the curve  [Graphics:Images/CurvatureMod_gr_184.gif], the center and radius of the osculating circle are given by the limits calculated in the preceding discussion.

[Graphics:Images/CurvatureMod_gr_185.gif]


[Graphics:Images/CurvatureMod_gr_186.gif]

 

Example 5.  Consider the parabola  [Graphics:Images/CurvatureMod_gr_187.gif] at the point  [Graphics:Images/CurvatureMod_gr_188.gif].  Draw the circle of curvature for various values  
[Graphics:Images/CurvatureMod_gr_189.gif]
Solution 5.

 

Generalizations for 2D

    In two dimensions, a curve can be expressed with the parametric equations  [Graphics:Images/CurvatureMod_gr_213.gif] and  [Graphics:Images/CurvatureMod_gr_214.gif].   Similarly, the formulas for the radius of curvature and center of curvature can be derived using limits.  At the point  [Graphics:Images/CurvatureMod_gr_215.gif]  the center and radius of the circle of convergence is
    
        [Graphics:Images/CurvatureMod_gr_216.gif]  
        
        [Graphics:Images/CurvatureMod_gr_217.gif]  

Remark.  The absolute value is necessary, otherwise the formula would only work for a curve that is positively oriented.
Details

 

The Abscissa the Easy Way

Subtract from  [Graphics:Images/CurvatureMod_gr_229.gif] the radius of curvature times [Graphics:Images/CurvatureMod_gr_230.gif].  The abscissa of the circle of curvature is

    [Graphics:Images/CurvatureMod_gr_231.gif]

Details

 

The Ordinate the Easy Way

Add to  [Graphics:Images/CurvatureMod_gr_234.gif] the radius of curvature times [Graphics:Images/CurvatureMod_gr_235.gif].  The ordinate of the circle of curvature is

    [Graphics:Images/CurvatureMod_gr_236.gif]  

Details

 

Example 6.  Consider the cardioid  [Graphics:Images/CurvatureMod_gr_239.gif],  [Graphics:Images/CurvatureMod_gr_240.gif].  Draw the circle of curvature at   [Graphics:Images/CurvatureMod_gr_241.gif].  
Solution 6.

 

Example 7.  Consider the cardioid  [Graphics:Images/CurvatureMod_gr_253.gif],  [Graphics:Images/CurvatureMod_gr_254.gif].  Draw the circle of curvature at   [Graphics:Images/CurvatureMod_gr_255.gif].  
Solution 7.

 

Reference

John Mathews, The Circle of Curvature: It's a Limit!, The AMATYC Review, Vol. 25, No. 1, Fall 2003, pp. 57-63.  

 

Research Experience for Undergraduates

Curvature  Curvature  Internet hyperlinks to web sites and a bibliography of articles.  

 

Download this Mathematica Notebook Curvature

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004