Osculating circle

Last updated September 20, 2023

An osculating circle is a circle that best approximates the curvature of a curve at a specific point. It is tangent to the curve at that point and has the same curvature as the curve at that point.^[2] The osculating circle provides a way to understand the local behavior of a curve and is commonly used in differential geometry and calculus.

More formally, in differential geometry of curves, the osculating circle of a sufficiently smooth plane curve at a given point p on the curve has been traditionally defined as the circle passing through p and a pair of additional points on the curve infinitesimally close to p. Its center lies on the inner normal line, and its curvature defines the curvature of the given curve at that point. This circle, which is the one among all tangent circles at the given point that approaches the curve most tightly, was named circulus osculans (Latin for "kissing circle") by Leibniz.

The center and radius of the osculating circle at a given point are called center of curvature and radius of curvature of the curve at that point. A geometric construction was described by Isaac Newton in his Principia :

There being given, in any places, the velocity with which a body describes a given figure, by means of forces directed to some common centre: to find that centre.
— Isaac Newton, Principia; PROPOSITION V. PROBLEM I.

Nontechnical description

Imagine a car moving along a curved road on a vast flat plane. Suddenly, at one point along the road, the steering wheel locks in its present position. Thereafter, the car moves in a circle that "kisses" the road at the point of locking. The curvature of the circle is equal to that of the road at that point. That circle is the osculating circle of the road curve at that point.

Mathematical description

Let $γ (s)$ be a regular parametric plane curve, where $s$ is the arc length (the natural parameter). This determines the unit tangent vector $T (s)$ , the unit normal vector $N (s)$ , the signed curvature $k (s)$ and the radius of curvature $R (s)$ at each point for which $s$ is composed:

T(s)=\gamma '(s),\quad T'(s)=k(s)N(s),\quad R(s)={\frac {1}{\left|k(s)\right|}}.

Suppose that P is a point on γ where $k \neq 0$ . The corresponding center of curvature is the point Q at distance R along N, in the same direction if k is positive and in the opposite direction if k is negative. The circle with center at Q and with radius R is called the osculating circle to the curve γ at the point P.

If C is a regular space curve then the osculating circle is defined in a similar way, using the principal normal vector N. It lies in the osculating plane , the plane spanned by the tangent and principal normal vectors T and N at the point P.

The plane curve can also be given in a different regular parametrization

\gamma (t)={\begin{bmatrix}x_{1}(t)\\x_{2}(t)\end{bmatrix}}

where regular means that $\gamma '(t)\neq 0$ for all $t$ . Then the formulas for the signed curvature k(t), the normal unit vector N(t), the radius of curvature R(t), and the center Q(t) of the osculating circle are

k(t)={\frac {x_{1}'(t)\,x_{2}''(t)-x_{1}''(t)\,x_{2}'(t)}{\left(x_{1}'\left(t\right)^{2}+x_{2}'\left(t\right)^{2}\right)^{{3}/{2}}}},\qquad N(t)={\frac {1}{\|\gamma '(t)\|}}{\begin{bmatrix}-x_{2}'(t)\\x_{1}'(t)\end{bmatrix}}

R(t)=\left|{\frac {\left(x_{1}'\left(t\right)^{2}+x_{2}'\left(t\right)^{2}\right)^{{3}/{2}}}{x_{1}'(t)\,x_{2}''(t)-x_{1}''(t)\,x_{2}'(t)}}\right|\qquad {\text{and}}\qquad Q(t)=\gamma (t)+{\frac {1}{k(t)\|\gamma '(t)\|}}{\begin{bmatrix}-x_{2}'(t)\\x_{1}'(t)\end{bmatrix}}\,.

Cartesian coordinates

We can obtain the center of the osculating circle in Cartesian coordinates if we substitute $t = x$ and $y = f (x)$ for some function f. If we do the calculations the results for the X and Y coordinates of the center of the osculating circle are:

x_{c}=x-f'{\frac {1+f'^{2}}{f''}}\quad {\text{and}}\quad y_{c}=f+{\frac {1+f'^{2}}{f''}}

Direct geometrical derivation

Consider three points ${\textstyle P_{0}}$ , ${\textstyle P_{1}}$ and ${\textstyle P_{2}}$ , where ${\textstyle P_{i}=(x_{i},y_{i})}$ . To find the center of the circle that passes through these points, we have first to find the segment bisectors of ${\textstyle P_{0}P_{1}}$ and ${\textstyle P_{1}P_{2}}$ and then the point ${\textstyle C}$ where these lines cross. Therefore, the coordinates of ${\textstyle C}$ are obtained through solving a linear system of two equations:

\left(\delta x_{i}\right)x_{c}+\left(\delta y_{i}\right)y_{c}={\tfrac {1}{2}}\left(\delta ^{2}x_{i}+\delta ^{2}y_{i}\right)\quad i=1,2

where ${\textstyle \delta u_{i}=u_{i}-u_{i-1}}$ , ${\textstyle \delta ^{2}u_{i}=u_{i}^{2}-u_{i-1}^{2}}$ for ${\textstyle u=x,y}$ .

Consider now the curve ${\textstyle P=P(\tau )}$ and set ${\textstyle P_{0}=P(\tau -d\tau )}$ , ${\textstyle P_{1}=P(\tau )}$ and ${\textstyle P_{2}=P(\tau +d\tau )}$ . To the second order in ${\textstyle d\tau }$ , we have

{\begin{aligned}\delta u_{1}=&{\dot {u}}d\tau -{\frac {1}{2}}{\ddot {u}}\,d\tau ^{2}\\\delta ^{2}u_{1}=&2u{\dot {u}}\,d\tau -d\tau ^{2}\left({\dot {u}}^{2}+u{\ddot {u}}\right)\end{aligned}}

and a similar expression for ${\textstyle \delta u_{2}}$ and ${\textstyle \delta ^{2}u_{2}}$ where the sign of ${\textstyle d\tau ^{2}}$ is reversed. Developing the equation for ${\textstyle x_{c},y_{c}}$ and grouping the terms in ${\textstyle d\tau }$ and ${\textstyle d\tau ^{2}}$ , we obtain

{\begin{aligned}{\dot {x}}(x_{c}-x)+{\dot {y}}(y_{c}-y)&=0\\{\ddot {x}}(x_{c}-x)+{\ddot {y}}(y_{c}-y)&={\dot {x}}^{2}+{\dot {y}}^{2}\end{aligned}}

Denoting ${\textstyle \mathbf {r} ={\overrightarrow {P_{1}C}}}$ , the first equation means that ${\textstyle \mathbf {r} }$ is orthogonal to the unit tangent vector at ${\textstyle P_{1}}$ :

\mathbf {r} \cdot \mathbf {t} =0

The second relation means that

\mathbf {k} \cdot \mathbf {r} =1

where

\mathbf {k} ={\frac {d\mathbf {t} }{ds}}={\frac {1}{{\dot {x}}^{2}+{\dot {y}}^{2}}}{\begin{bmatrix}{\ddot {x}}\\{\ddot {y}}\end{bmatrix}}

is the curvature vector. In plane geometry, ${\textstyle \mathbf {k} }$ is orthogonal to ${\textstyle \mathbf {t} }$ because

\mathbf {t} \cdot \mathbf {k} =\mathbf {t} {\frac {d\mathbf {t} }{ds}}={\frac {1}{2}}{\frac {d}{ds}}(\mathbf {t} \cdot \mathbf {t} )={\frac {1}{2}}{\frac {d}{ds}}(1)=0

Therefore ${\textstyle \mathbf {k} \cdot \mathbf {r} =kr}$ and the radius of the osculating circle is precisely the inverse of the curvature.

Solving the equation for the coordinates of ${\textstyle C}$ , we find

{\begin{aligned}x_{c}-x=&{\frac {{\dot {y}}\left({\dot {x}}^{2}+{\dot {y}}^{2}\right)}{{\dot {y}}{\ddot {x}}-{\dot {x}}{\ddot {y}}}}\\y_{c}-y=&{\frac {-{\dot {x}}\left({\dot {x}}^{2}+{\dot {y}}^{2}\right)}{{\dot {y}}{\ddot {x}}-{\dot {x}}{\ddot {y}}}}\end{aligned}}

Osculating circle as a minimization problem

Consider a curve ${\textstyle C}$ defined intrinsically by the equation

f(x,y)=0

which we can envision as the section of the surface ${\textstyle z=f(x,y)}$ by the plane ${\textstyle z=0}$ .

The normal ${\textstyle \mathbf {n} }$ to the curve at a point ${\textstyle P_{0}=(x_{0},y_{0})}$ is the gradient at this point

\mathbf {n} =(f_{x},f_{y})

Therefore, the centers of the tangent circles ${\textstyle B_{\alpha }}$ are given by

X_{c}=x_{0}-\alpha f_{x}\,\,;\,\,Y_{c}=y_{0}-\alpha f_{y}

where ${\textstyle \alpha }$ is parameter. For a given ${\textstyle \alpha ,}$ the radius ${\textstyle R}$ of ${\textstyle B_{\alpha }}$ is

R^{2}=\alpha ^{2}(f_{x}^{2}+f_{y}^{2})

We wish to find, among all possible circles ${\textstyle B_{\alpha }}$ , the one that matches best the curve.

The coordinates of a point ${\textstyle P_{1}\in B_{\alpha }}$ can be written as

x_{1}=X_{c}+R\cos \theta \,\,;\,\,y_{1}=Y_{c}+R\sin \theta

where for ${\textstyle \theta =\theta _{0}}$ , ${\textstyle P_{1}=P_{0}}$ , i.e.

R\cos \theta _{0}=\alpha f_{x}\,\,;\,\,R\sin \theta _{0}=\alpha f_{y}

Consider now a point ${\textstyle P_{1}\in B_{\alpha }}$ close to ${\textstyle P_{0}}$ , where its "angle" is ${\textstyle \theta _{1}=\theta _{0}+d\theta }$ . Developing the trigonometric functions to the second order in ${\textstyle d\theta }$ and using the above relations, coordinates of $P_{1}$ are

{\begin{aligned}x_{1}=&x_{0}-\alpha f_{y}d\theta -{\tfrac {1}{2}}\alpha f_{x}\left(d\theta \right)^{2}\\y_{1}=&y_{0}+\alpha f_{x}d\theta -{\tfrac {1}{2}}\alpha f_{y}\left(d\theta \right)^{2}\end{aligned}}

We can now evaluate the function ${\textstyle f}$ at the point ${\textstyle P_{1}}$ and its variation $f(x_{1},y_{1})-f(x_{0},y_{0})$ . The variation is zero to the first order in ${\textstyle d\theta }$ by construction (to the first order in ${\textstyle \theta }$ , ${\textstyle P_{1}}$ is on the tangent line to the curve ${\textstyle C}$ ). The variation proportional to $(d\theta )^{2}$ is

df=-{\frac {1}{2}}\alpha \left(f_{x}^{2}+f_{y}^{2}\right)+{\frac {1}{2}}\alpha ^{2}\left(f_{y}^{2}f_{xx}+f_{x}^{2}f_{yy}-f_{x}f_{y}f_{xy}\right)

and this variation is zero if we choose

\alpha ={\frac {f_{x}^{2}+f_{y}^{2}}{f_{y}^{2}f_{xx}+f_{x}^{2}f_{yy}-f_{x}f_{y}f_{xy}}}

Therefore the radius of the osculating circle is

R=\left|{\frac {\left(f_{x}^{2}+f_{y}^{2}\right)^{3/2}}{\left(f_{y}^{2}f_{xx}+f_{x}^{2}f_{yy}-f_{x}f_{y}f_{xy}\right)}}\right|

For an explicit function $f(x,y)=y-g(x)$ , we find the results of the preceding section.

Properties

For a curve C given by a sufficiently smooth parametric equations (twice continuously differentiable), the osculating circle may be obtained by a limiting procedure: it is the limit of the circles passing through three distinct points on C as these points approach P.^[3] This is entirely analogous to the construction of the tangent to a curve as a limit of the secant lines through pairs of distinct points on C approaching P.

The osculating circle S to a plane curve C at a regular point P can be characterized by the following properties:

The circle S passes through P.
The circle S and the curve C have the common tangent line at P, and therefore the common normal line.
Close to P, the distance between the points of the curve C and the circle S in the normal direction decays as the cube or a higher power of the distance to P in the tangential direction.

This is usually expressed as "the curve and its osculating circle have the second or higher order contact" at P. Loosely speaking, the vector functions representing C and S agree together with their first and second derivatives at P.

If the derivative of the curvature with respect to s is nonzero at P then the osculating circle crosses the curve C at P. Points P at which the derivative of the curvature is zero are called vertices. If P is a vertex then C and its osculating circle have contact of order at least three. If, moreover, the curvature has a non-zero local maximum or minimum at P then the osculating circle touches the curve C at P but does not cross it.

The curve C may be obtained as the envelope of the one-parameter family of its osculating circles. Their centers, i.e. the centers of curvature, form another curve, called the evolute of C. Vertices of C correspond to singular points on its evolute.

Within any arc of a curve C within which the curvature is monotonic (that is, away from any vertex of the curve), the osculating circles are all disjoint and nested within each other. This result is known as the Tait-Kneser theorem.^[1]

Examples

Parabola

For the parabola

\gamma (t)={\begin{bmatrix}t\\t^{2}\end{bmatrix}}

the radius of curvature is

R(t)=\left|{\frac {\left(1+4t^{2}\right)^{3/2}}{2}}\right|

At the vertex $\gamma (0)={\begin{bmatrix}0\\0\end{bmatrix}}$ the radius of curvature equals $R (0) = 0.5$ (see figure). The parabola has fourth order contact with its osculating circle there. For large $t$ the radius of curvature increases ~ $t 3$ , that is, the curve straightens more and more.

Lissajous curve

A Lissajous curve with ratio of frequencies (3:2) can be parametrized as follows

\gamma (t)={\begin{bmatrix}\cos(3t)\\\sin(2t)\end{bmatrix}}.

It has signed curvature $k (t)$ , normal unit vector $N (t)$ and radius of curvature $R (t)$ given by

k(t)={\frac {6\cos(t)(8(\cos t)^{4}-10(\cos t)^{2}+5)}{\left(232(\cos t)^{4}-97(\cos t)^{2}+13-144(\cos t)^{6}\right)^{3/2}}}\,,

N(t)={\frac {1}{\|\gamma '(t)\|}}\cdot {\begin{bmatrix}-2\cos(2t)\\-3\sin(3t)\end{bmatrix}}

and

R(t)=\left|{\frac {\left(232\cos ^{4}(t)-97\cos ^{2}(t)+13-144\cos ^{6}(t)\right)^{3/2}}{6\cos(t)\left(8\cos ^{4}(t)-10\cos ^{2}(t)+5\right)}}\right|.

See the figure for an animation. There the "acceleration vector" is the second derivative ${\textstyle {\frac {d^{2}\gamma }{ds^{2}}}}$ with respect to the arc length $s$ .

Cycloid

A cycloid with radius $r$ can be parametrized as follows:

\gamma (t)={\begin{bmatrix}r\left(t-\sin t\right)\\r\left(1-\cos t\right)\end{bmatrix}}

Its curvature is given by the following formula:^[4]

\kappa (t)=-{\frac {\left|\csc \left({\frac {t}{2}}\right)\right|}{4r}}

which gives:

R(t)={\frac {4r}{\left|\csc \left({\frac {t}{2}}\right)\right|}}

Notes

1 2 Ghys, Étienne; Tabachnikov, Sergei; Timorin, Vladlen (2013). "Osculating curves: around the Tait-Kneser theorem". The Mathematical Intelligencer. 35 (1): 61–66. arXiv: 1207.5662 . doi:10.1007/s00283-012-9336-6. MR 3041992. S2CID 18183204.
↑ "12.4 Arc Length and Curvature" . Retrieved 2023-09-19.
↑ Actually, point P plus two additional points, one on either side of P will do. See Lamb (on line): Horace Lamb (1897). An Elementary Course of Infinitesimal Calculus. University Press. p. 406. osculating circle.
↑ Weisstein, Eric W. "Cycloid". MathWorld .

External links

Related Research Articles

A centripetal force is a force that makes a body follow a curved path. The direction of the centripetal force is always orthogonal to the motion of the body and towards the fixed point of the instantaneous center of curvature of the path. Isaac Newton described it as "a force by which bodies are drawn or impelled, or in any way tend, towards a point as to a centre". In the theory of Newtonian mechanics, gravity provides the centripetal force causing astronomical orbits.

<span class="mw-page-title-main">Ellipse</span> Plane curve: conic section

In mathematics, an ellipse is a plane curve surrounding two focal points, such that for all points on the curve, the sum of the two distances to the focal points is a constant. It generalizes a circle, which is the special type of ellipse in which the two focal points are the same. The elongation of an ellipse is measured by its eccentricity $, a number ranging from to .$

<span class="mw-page-title-main">Curvature</span> Mathematical measure of how much a curve or surface deviates from flatness

In mathematics, curvature is any of several strongly related concepts in geometry. Intuitively, the curvature is the amount by which a curve deviates from being a straight line, or a surface deviates from being a plane.

Kinematics is a subfield of physics, developed in classical mechanics, that describes the motion of points, bodies (objects), and systems of bodies without considering the forces that cause them to move. Kinematics, as a field of study, is often referred to as the "geometry of motion" and is occasionally seen as a branch of mathematics. A kinematics problem begins by describing the geometry of the system and declaring the initial conditions of any known values of position, velocity and/or acceleration of points within the system. Then, using arguments from geometry, the position, velocity and acceleration of any unknown parts of the system can be determined. The study of how forces act on bodies falls within kinetics, not kinematics. For further details, see analytical dynamics.

<span class="mw-page-title-main">Ellipsoid</span> Quadric surface that looks like a deformed sphere

An ellipsoid is a surface that can be obtained from a sphere by deforming it by means of directional scalings, or more generally, of an affine transformation.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space $under the operation of composition.$

Unit quaternions, known as versors, provide a convenient mathematical notation for representing spatial orientations and rotations of elements in three dimensional space. Specifically, they encode information about an axis-angle rotation about an arbitrary axis. Rotation and orientation quaternions have applications in computer graphics, computer vision, robotics, navigation, molecular dynamics, flight dynamics, orbital mechanics of satellites, and crystallographic texture analysis.

In the mathematical field of differential geometry, a metric tensor is an additional structure on a manifold $M$ that allows defining distances and angles, just as the inner product on a Euclidean space allows defining distances and angles there. More precisely, a metric tensor at a point $p$ of $M$ is a bilinear form defined on the tangent space at $p$ , and a metric tensor on $M$ consists of a metric tensor at each point $p$ of $M$ that varies smoothly with $p$ .

In physics, circular motion is a movement of an object along the circumference of a circle or rotation along a circular arc. It can be uniform, with a constant rate of rotation and constant tangential speed, or non-uniform with a changing rate of rotation. The rotation around a fixed axis of a three-dimensional body involves the circular motion of its parts. The equations of motion describe the movement of the center of mass of a body, which remains at a constant distance from the axis of rotation. In circular motion, the distance between the body and a fixed point on its surface remains the same, i.e., the body is assumed rigid.

In linear algebra, a rotation matrix is a transformation matrix that is used to perform a rotation in Euclidean space. For example, using the convention below, the matrix

<span class="mw-page-title-main">Envelope (mathematics)</span> Family of curves in geometry

In geometry, an envelope of a planar family of curves is a curve that is tangent to each member of the family at some point, and these points of tangency together form the whole envelope. Classically, a point on the envelope can be thought of as the intersection of two "infinitesimally adjacent" curves, meaning the limit of intersections of nearby curves. This idea can be generalized to an envelope of surfaces in space, and so on to higher dimensions.

Arc length is the distance between two points along a section of a curve.

In calculus, the Leibniz integral rule for differentiation under the integral sign states that for an integral of the form

A parametric surface is a surface in the Euclidean space $which is defined by a parametric equation with two parameters :\mathbb {R} ^{2}\to \mathbb {R} ^{3}} . Parametric representation is a very general way to specify a surface, as well as implicit representation. Surfaces that occur in two of the main theorems of vector calculus, Stokes' theorem and the divergence theorem, are frequently given in a parametric form. The curvature and arc length of curves on the surface, surface area, differential geometric invariants such as the first and second fundamental forms, Gaussian, mean, and principal curvatures can all be computed from a given parametrization.$

In geometry, various formalisms exist to express a rotation in three dimensions as a mathematical transformation. In physics, this concept is applied to classical mechanics where rotational kinematics is the science of quantitative description of a purely rotational motion. The orientation of an object at a given instant is described with the same tools, as it is defined as an imaginary rotation from a reference placement in space, rather than an actually observed rotation from a previous placement in space.

In classical mechanics, holonomic constraints are relations between the position variables that can be expressed in the following form:

In the differential geometry of surfaces, a Darboux frame is a natural moving frame constructed on a surface. It is the analog of the Frenet–Serret frame as applied to surface geometry. A Darboux frame exists at any non-umbilic point of a surface embedded in Euclidean space. It is named after French mathematician Jean Gaston Darboux.

In celestial mechanics, a Kepler orbit is the motion of one body relative to another, as an ellipse, parabola, or hyperbola, which forms a two-dimensional orbital plane in three-dimensional space. A Kepler orbit can also form a straight line. It considers only the point-like gravitational attraction of two bodies, neglecting perturbations due to gravitational interactions with other objects, atmospheric drag, solar radiation pressure, a non-spherical central body, and so on. It is thus said to be a solution of a special case of the two-body problem, known as the Kepler problem. As a theory in classical mechanics, it also does not take into account the effects of general relativity. Keplerian orbits can be parametrized into six orbital elements in various ways.

In differential geometry, the radius of curvature, $R$ , is the reciprocal of the curvature. For a curve, it equals the radius of the circular arc which best approximates the curve at that point. For surfaces, the radius of curvature is the radius of a circle that best fits a normal section or combinations thereof.

A proper reference frame in the theory of relativity is a particular form of accelerated reference frame, that is, a reference frame in which an accelerated observer can be considered as being at rest. It can describe phenomena in curved spacetime, as well as in "flat" Minkowski spacetime in which the spacetime curvature caused by the energy–momentum tensor can be disregarded. Since this article considers only flat spacetime—and uses the definition that special relativity is the theory of flat spacetime while general relativity is a theory of gravitation in terms of curved spacetime—it is consequently concerned with accelerated frames in special relativity.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[gtt-1] 1 2 Ghys, Étienne; Tabachnikov, Sergei; Timorin, Vladlen (2013). "Osculating curves: around the Tait-Kneser theorem". The Mathematical Intelligencer. 35 (1): 61–66. arXiv: 1207.5662 . doi:10.1007/s00283-012-9336-6. MR 3041992. S2CID 18183204.

[2] "12.4 Arc Length and Curvature" . Retrieved 2023-09-19.

[Lamb-3] Actually, point P plus two additional points, one on either side of P will do. See Lamb (on line): Horace Lamb (1897). An Elementary Course of Infinitesimal Calculus. University Press. p. 406. osculating circle.

[4] Weisstein, Eric W. "Cycloid". MathWorld .

[2]

[3]

[1]

[4]