# Euler spiral

Last updated

An Euler spiral is a curve whose curvature changes linearly with its curve length (the curvature of a circular curve is equal to the reciprocal of the radius). Euler spirals are also commonly referred to as spiros, clothoids, or Cornu spirals.

## Contents

Euler spirals have applications to diffraction computations. They are also widely used as transition curves in railroad engineering/highway engineering for connecting and transitioning the geometry between a tangent and a circular curve. A similar application is also found in photonic integrated circuits. The principle of linear variation of the curvature of the transition curve between a tangent and a circular curve defines the geometry of the Euler spiral:

• Its curvature begins with zero at the straight section (the tangent) and increases linearly with its curve length.
• Where the Euler spiral meets the circular curve, its curvature becomes equal to that of the latter.

## Applications

### Track transition curve

To travel along a circular path, an object needs to be subject to a centripetal acceleration (for example: the Moon circles around the Earth because of gravity; a car turns its front wheels inward to generate a centripetal force). If a vehicle traveling on a straight path were to suddenly transition to a tangential circular path, it would require centripetal acceleration suddenly switching at the tangent point from zero to the required value; this would be difficult to achieve (think of a driver instantly moving the steering wheel from straight line to turning position, and the car actually doing it), putting mechanical stress on the vehicle's parts, and causing much discomfort (due to lateral jerk).

On early railroads this instant application of lateral force was not an issue since low speeds and wide-radius curves were employed (lateral forces on the passengers and the lateral sway was small and tolerable). As speeds of rail vehicles increased over the years, it became obvious that an easement is necessary, so that the centripetal acceleration increases linearly with the traveled distance. Given the expression of centripetal acceleration v2/r, the obvious solution is to provide an easement curve whose curvature, 1/R, increases linearly with the traveled distance. This geometry is an Euler spiral.

Unaware of the solution of the geometry by Leonhard Euler, Rankine cited the cubic curve (a polynomial curve of degree 3), which is an approximation of the Euler spiral for small angular changes in the same way that a parabola is an approximation to a circular curve.

Marie Alfred Cornu (and later some civil engineers) also solved the calculus of the Euler spiral independently. Euler spirals are now widely used in rail and highway engineering for providing a transition or an easement between a tangent and a horizontal circular curve.

### Optics

The Cornu spiral can be used to describe a diffraction pattern. [1] Consider a plane wave with phasor amplitude E0ejkz which is diffracted by a "knife edge" of height h above x = 0 on the z = 0 plane. Then the diffracted wave field can be expressed as

${\displaystyle {\textbf {E}}(x,z)=E_{0}e^{-jkz}{\frac {\mathrm {Fr} (\infty )-\mathrm {Fr} \left({\sqrt {\frac {2}{\lambda z}}}(h-x)\right)}{\mathrm {Fr} (\infty )-\mathrm {Fr} (-\infty )}}}$,

where Fr(x) is the Fresnel integral function, which forms the Cornu spiral on the complex plane.

So, to simplify the calculation of plane wave attenuation as it is diffracted from the knife-edge, one can use the diagram of a Cornu spiral by representing the quantities Fr(a) − Fr(b) as the physical distances between the points represented by Fr(a) and Fr(b) for appropriate a and b. This facilitates a rough computation of the attenuation of the plane wave by the knife edge of height h at a location (x, z) beyond the knife edge.

### Integrated optics

Bends with continuously varying radius of curvature following the Euler spiral are also used to reduce losses in photonic integrated circuits, either in singlemode waveguides, [2] [3] to smoothen the abrupt change of curvature and coupling to radiation modes, or in multimode waveguides, [4] in order to suppress coupling to higher order modes and ensure effective singlemode operation. A pioneering and very elegant application of the Euler spiral to waveguides had been made as early as 1957, [5] with a hollow metal waveguide for microwaves. There the idea was to exploit the fact that a straight metal waveguide can be physically bent to naturally take a gradual bend shape resembling an Euler spiral.

### Auto racing

Motorsport author Adam Brouillard has shown the Euler spiral's use in optimizing the racing line during the corner entry portion of a turn. [6]

### Typography and digital vector drawing

Raph Levien has released Spiro as a toolkit for curve design, especially font design, in 2007 [7] [8] under a free licence. This toolkit has been implemented quite quickly afterwards in the font design tool Fontforge and the digital vector drawing Inkscape.

### Map projection

Cutting a sphere along a spiral with width 1/N and flattening out the resulting shape yields an Euler spiral when n tends to the infinity. [9] If the sphere is the globe, this produces a map projection whose distortion tends to zero as n tends to the infinity. [10]

### Whisker shapes

Natural shapes of rat's mystacial pad vibrissae (whiskers) are well approximated by pieces of the Euler spiral. When all these pieces for a single rat are assembled together, they span an interval extending from one coiled domain of the Euler spiral to the other. [11]

## Formulation

### Symbols

 R Radius of curvature Rc Radius of circular curve at the end of the spiral θ Angle of curve from beginning of spiral (infinite R) to a particular point on the spiral. This can also be measured as the angle between the initial tangent and the tangent at the concerned point. θs Angle of full spiral curve L, s Length measured along the spiral curve from its initial position Ls, so Length of spiral curve

### Expansion of Fresnel integral

If a = 1, which is the case for normalized Euler curve, then the Cartesian coordinates are given by Fresnel integrals (or Euler integrals):

{\displaystyle {\begin{aligned}C(L)&=\int _{0}^{L}\cos \left(s^{2}\right)\,ds\\S(L)&=\int _{0}^{L}\sin \left(s^{2}\right)\,ds\end{aligned}}}

### Normalization and conclusion

For a given Euler curve with:

${\displaystyle 2RL=2R_{c}L_{s}={\frac {1}{a^{2}}}}$

or

${\displaystyle {\frac {1}{R}}={\frac {L}{R_{c}L_{s}}}=2a^{2}L}$

then

{\displaystyle {\begin{aligned}x&={\frac {1}{a}}\int _{0}^{L'}\cos \left(s^{2}\right)\,ds\\y&={\frac {1}{a}}\int _{0}^{L'}\sin \left(s^{2}\right)\,ds\end{aligned}}}

where

{\displaystyle {\begin{aligned}L'&=aL\\a&={\frac {1}{\sqrt {2R_{c}L_{s}}}}.\end{aligned}}}

The process of obtaining solution of (x, y) of an Euler spiral can thus be described as:

• Map L of the original Euler spiral by multiplying with factor a to L of the normalized Euler spiral;
• Find (x′, y′) from the Fresnel integrals; and
• Map (x′, y′) to (x, y) by scaling up (denormalize) with factor 1/a. Note that 1/a > 1.

In the normalization process,

{\displaystyle {\begin{aligned}R'_{c}&={\frac {R_{c}}{\sqrt {2R_{c}L_{s}}}}\\&={\sqrt {\frac {R_{c}}{2L_{s}}}}\\L'_{s}&={\frac {L_{s}}{\sqrt {2R_{c}L_{s}}}}\\&={\sqrt {\frac {L_{s}}{2R_{c}}}}\end{aligned}}}

Then

{\displaystyle {\begin{aligned}2R'_{c}L'_{s}&=2{\sqrt {\frac {R_{c}}{2L_{s}}}}{\sqrt {\frac {L_{s}}{2R_{c}}}}\\&={\frac {2}{2}}=1\end{aligned}}}

Generally the normalization reduces L to a small value (less than 1) and results in good converging characteristics of the Fresnel integral manageable with only a few terms (at a price of increased numerical instability of the calculation, especially for bigger θ values.).

#### Illustration

Given:

{\displaystyle {\begin{aligned}R_{c}&=300\,{\mbox{m}}\\L_{s}&=100\,{\mbox{m}}\end{aligned}}}

Then

{\displaystyle {\begin{aligned}\theta _{s}&={\frac {L_{s}}{2R_{c}}}\\&={\frac {100}{2\times 300}}\\&={\tfrac {1}{6}}\ {\mbox{radian}}\\\end{aligned}}}

and

${\displaystyle 2R_{c}L_{s}=60\,000}$

We scale down the Euler spiral by 60000, i.e. 1006 to normalized Euler spiral that has:

{\displaystyle {\begin{aligned}R'_{c}&={\tfrac {3}{\sqrt {6}}}\,{\mbox{m}}\\L'_{s}&={\tfrac {1}{\sqrt {6}}}\,{\mbox{m}}\\2R'_{c}L'_{s}&=2\times {\tfrac {3}{\sqrt {6}}}\times {\tfrac {1}{\sqrt {6}}}\\&=1\end{aligned}}}

and

{\displaystyle {\begin{aligned}\theta _{s}&={\frac {L'_{s}}{2R'_{c}}}\\&={\frac {\frac {1}{\sqrt {6}}}{2\times {\frac {3}{\sqrt {6}}}}}\\&={\tfrac {1}{6}}\ {\mbox{radian}}\\\end{aligned}}}

The two angles θs are the same. This thus confirms that the original and normalized Euler spirals are geometrically similar. The locus of the normalized curve can be determined from Fresnel Integral, while the locus of the original Euler spiral can be obtained by scaling up or denormalizing.

### Other properties of normalized Euler spirals

Normalized Euler spirals can be expressed as:

{\displaystyle {\begin{aligned}x&=\int _{0}^{L}\cos \left(s^{2}\right)\,ds\\y&=\int _{0}^{L}\sin \left(s^{2}\right)\,ds\end{aligned}}}

or expressed as power series:

{\displaystyle {\begin{aligned}x&=\left.\sum _{i=0}^{\infty }{\frac {(-1)^{i}}{(2i)!}}{\frac {s^{4i+1}}{4i+1}}\right|_{0}^{L}&&=\sum _{i=0}^{\infty }{\frac {(-1)^{i}}{(2i)!}}{\frac {L^{4i+1}}{4i+1}}\\y&=\left.\sum _{i=0}^{\infty }{\frac {(-1)^{i}}{(2i+1)!}}{\frac {s^{4i+3}}{4i+3}}\right|_{0}^{L}&&=\sum _{i=0}^{\infty }{\frac {(-1)^{i}}{(2i+1)!}}{\frac {L^{4i+3}}{4i+3}}\end{aligned}}}

The normalized Euler spiral will converge to a single point in the limit as the parameter L approaches infinity, which can be expressed as:

{\displaystyle {\begin{aligned}x^{\prime }&=\lim _{L\to \infty }\int _{0}^{L}\cos \left(s^{2}\right)\,ds&&={\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}\approx 0.6267\\y^{\prime }&=\lim _{L\to \infty }\int _{0}^{L}\sin \left(s^{2}\right)\,ds&&={\frac {1}{2}}{\sqrt {\frac {\pi }{2}}}\approx 0.6267\end{aligned}}}

Normalized Euler spirals have the following properties:

{\displaystyle {\begin{aligned}2R_{c}L_{s}&=1\\\theta _{s}&={\frac {L_{s}}{2R_{c}}}=L_{s}^{2}\end{aligned}}}

and

{\displaystyle {\begin{aligned}\theta &=\theta _{s}\cdot {\frac {L^{2}}{L_{s}^{2}}}=L^{2}\\{\frac {1}{R}}&={\frac {d\theta }{dL}}=2L\end{aligned}}}

Note that 2RcLs = 1 also means 1/Rc = 2Ls, in agreement with the last mathematical statement.

## Related Research Articles

In mathematics, the arithmetic–geometric mean (AGM) of two positive real numbers x and y is defined as follows:

Euler's formula, named after Leonhard Euler, is a mathematical formula in complex analysis that establishes the fundamental relationship between the trigonometric functions and the complex exponential function. Euler's formula states that for any real number x:

In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices which are Hermitian and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

In mathematics, the polar coordinate system is a two-dimensional coordinate system in which each point on a plane is determined by a distance from a reference point and an angle from a reference direction. The reference point is called the pole, and the ray from the pole in the reference direction is the polar axis. The distance from the pole is called the radial coordinate, radial distance or simply radius, and the angle is called the angular coordinate, polar angle, or azimuth. The radial coordinate is often denoted by r or ρ, and the angular coordinate by φ, θ, or t. Angles in polar notation are generally expressed in either degrees or radians.

The Archimedean spiral is a spiral named after the 3rd-century BC Greek mathematician Archimedes. It is the locus corresponding to the locations over time of a point moving away from a fixed point with a constant speed along a line that rotates with constant angular velocity. Equivalently, in polar coordinates (r, θ) it can be described by the equation

In mathematics, a spiral is a curve which emanates from a point, moving farther away as it revolves around the point.

A helix, plural helixes or helices, is a shape like a corkscrew or spiral staircase. It is a type of smooth space curve with tangent lines at a constant angle to a fixed axis. Helices are important in biology, as the DNA molecule is formed as two intertwined helices, and many proteins have helical substructures, known as alpha helices. The word helix comes from the Greek word ἕλιξ, "twisted, curved". A "filled-in" helix – for example, a "spiral" (helical) ramp – is called a helicoid.

A hyperbolic spiral is a plane curve, which can be described in polar coordinates by the equation

A Fermat's spiral or parabolic spiral is a plane curve named after Pierre de Fermat. Its polar coordinate representation is given by

In mathematics, a Gaussian function, often simply referred to as a Gaussian, is a function of the form

The Fresnel integralsS(x) and C(x) are two transcendental functions named after Augustin-Jean Fresnel that are used in optics and are closely related to the error function. They arise in the description of near-field Fresnel diffraction phenomena and are defined through the following integral representations:

In mathematics, theta functions are special functions of several complex variables. They are important in many areas, including the theories of Abelian varieties and moduli spaces, and of quadratic forms. They have also been applied to soliton theory. When generalized to a Grassmann algebra, they also appear in quantum field theory.

In mathematics, physics and engineering, the sinc function, denoted by sinc(x), has two forms, normalized and unnormalized.

In geometry, the area enclosed by a circle of radius r is πr2. Here the Greek letter π represents the constant ratio of the circumference of any circle to its diameter, approximately equal to 3.1416.

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

In optics, the Fresnel diffraction equation for near-field diffraction is an approximation of the Kirchhoff–Fresnel diffraction that can be applied to the propagation of waves in the near field. It is used to calculate the diffraction pattern created by waves passing through an aperture or around an object, when viewed from relatively close to the object. In contrast the diffraction pattern in the far field region is given by the Fraunhofer diffraction equation.

In mathematics, the sine is a trigonometric function of an angle. The sine of an acute angle is defined in the context of a right triangle: for the specified angle, it is the ratio of the length of the side that is opposite that angle, to the length of the longest side of the triangle. For an angle , the sine function is denoted simply as .

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 integral calculus, the Weierstrass substitution or tangent half-angle substitution is a method for evaluating integrals, which converts a rational function of trigonometric functions of into an ordinary rational function of by setting . No generality is lost by taking these to be rational functions of the sine and cosine. The general transformation formula is

## References

### Notes

1. Eugene Hecht (1998). Optics (3rd ed.). Addison-Wesley. p. 491. ISBN   978-0-201-30425-1.
2. Kohtoku, M.; et al. (7 July 2005). "New Waveguide Fabrication Techniques for Next-generation PLCs" (PDF). NTT Technical Review. 3 (7): 37–41. Retrieved 24 January 2017.
3. Li, G.; et al. (11 May 2012). "Ultralow-loss, high-density SOI optical waveguide routing for macrochip interconnects". Optics Express. 20 (11): 12035–12039. doi:. PMID   22714189.
4. Cherchi, M.; et al. (18 July 2013). "Dramatic size reduction of waveguide bends on a micron-scale silicon photonic platform". Optics Express. 21 (15): 17814–17823. arXiv:. doi:10.1364/OE.21.017814. PMID   23938654.
5. Unger, H.G. (September 1957). "Normal Mode Bends for Circular Electric Waves". The Bell System Technical Journal. 36 (5): 1292–1307. doi:10.1002/j.1538-7305.1957.tb01509.x.
6. Development, Paradigm Shift Driver; Brouillard, Adam (2016-03-18). The Perfect Corner: A Driver's Step-By-Step Guide to Finding Their Own Optimal Line Through the Physics of Racing. Paradigm Shift Motorsport Books. ISBN   9780997382426.
7. Bartholdi, Laurent; Henriques, André (2012). "Orange Peels and Fresnel Integrals". The Mathematical Intelligencer. 34 (3): 1–3. arXiv:. doi:10.1007/s00283-012-9304-1. ISSN   0343-6993.
8. Starostin, E.L.; et al. (15 January 2020). "The Euler spiral of rat whiskers". Science Advances. 6 (3): eaax5145. doi:.