Euler spiral

Formulation

Symbols

Derivation

The graph on the right illustrates an Euler spiral used as an easement (transition) curve between two given curves, in this case a straight line (the negative x axis) and a circle. The spiral starts at the origin in the positive x direction and gradually turns anticlockwise to osculate the circle.

The spiral is a small segment of the above double-end Euler spiral in the first quadrant.

From the definition of the curvature, $${\displaystyle \frac{1}{R}=\frac{d\theta }{ds}\propto s}$$ i.e., $${\displaystyle \begin{array}{cc}Rs=\text{constant}& =R_c{s}_o\\ \frac{d\theta }{ds}& =\frac{s}{R_c{s}_o}\end{array}}$$ We write in the format, $${\displaystyle \frac{d\theta }{ds}=2a^{2}s}$$ where $${\displaystyle 2a^2=\frac{1}{R_c{s}_o}}$$ or $${\displaystyle a=\frac{1}{\sqrt{2R_c{s}_o}}}$$ thus $${\displaystyle \theta =(as{)}^2}$$ Now $${\displaystyle x={\int }_0^L\cos \theta ds={\int }_0^L\cos \left[{\left(as\right)}^2\right]ds}$$ If $${\displaystyle s^{\prime }=as}$$ Then $${\displaystyle ds=\frac{d{s}^{\prime }}{a}}$$ Thus $${\displaystyle \begin{array}{cc}x& =\frac{1}{a}{\int }_0^{L^{\prime }}\cos \left(s^2\right)ds\\ y& ={\int }_0^L\sin \theta ds\\ & ={\int }_0^L\sin \left[{\left(as\right)}^2\right]ds\\ & =\frac{1}{a}{\int }_0^{L^{\prime }}\sin \left(s^2\right)ds\end{array}}$$

Relation to 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{array}{cc}C(L)& ={\int }_0^L\cos \left(s^2\right)ds\\ S(L)& ={\int }_0^L\sin \left(s^2\right)ds\end{array}}$$

Relation to superspirals

Euler's spiral is a type of superspiral that has the property of a monotonic curvature function.

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 smoothly with the traveled distance. Given the expression of centripetal acceleration ⁠v²/r⁠, the obvious solution is to provide an easement curve whose curvature, ⁠1/R⁠, increases linearly with the traveled distance.

Optics

In optics, the term Cornu spiral is used. The Cornu spiral can be used to describe a diffraction pattern. Consider a plane wave with phasor amplitude E₀e^−jkz 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 \mathbf{E}(x,z)=E_0{e}^{-jkz}\frac{\mathrm{F}\mathrm{r}(\mathrm{\infty })-\mathrm{F}\mathrm{r}\left(\sqrt{\frac{2}{\lambda z}}(h-x)\right)}{\mathrm{F}\mathrm{r}(\mathrm{\infty })-\mathrm{F}\mathrm{r}(-\mathrm{\infty })},}$$ where ${\displaystyle Fr(x)=C(x)+iS(x)}$ is the Fresnel integral function 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, to smoothen the abrupt change of curvature and suppress coupling to radiation modes, or in multimode waveguides, 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, 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.

Feynman's path integral

In the path integral formulation of quantum mechanics, the probability amplitude for propagation between two points can be visualized by connecting action phase arrows for each time step between the two points. The arrows spiral around each endpoint forming what is termed a Cornu 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.

Typography and digital vector drawing

Raph Levien has released Spiro as a toolkit for curve design, especially font design, in 2007 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. If the sphere is the globe, this produces a map projection whose distortion tends to zero as n tends to the infinity.

Whisker shapes

Natural shapes of rats' whiskers are well approximated by segments of Euler spirals; for a single rat all of the whiskers can be approximated as segments of the same spiral. The two parameters of the Cesàro equation for an Euler spiral segment might give insight into the keratinization mechanism of whisker growth.

Normalization

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{array}{cc}x& =\frac{1}{a}{\int }_0^{L^{\prime }}\cos \left(s^2\right)ds\\ y& =\frac{1}{a}{\int }_0^{L^{\prime }}\sin \left(s^2\right)ds\end{array}}$$ where $${\displaystyle \begin{array}{cc}{L}^{\prime }& =aL\\ a& =\frac{1}{\sqrt{2R_c{L}_s}}.\end{array}}$$

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{array}{cc}{R}_c^{\prime }& =\frac{R_c}{\sqrt{2R_c{L}_s}}=\sqrt{\frac{R_c}{2L_s}}\\ L_s^{\prime }& =\frac{L_s}{\sqrt{2R_c{L}_s}}=\sqrt{\frac{L_s}{2R_c}}\end{array}}$$ Then $${\displaystyle 2R_c^{\prime }{L}_s^{\prime }=2\sqrt{\frac{R_c}{2L_s}}\sqrt{\frac{L_s}{2R_c}}=\frac{2}{2}=1}$$

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{array}{cc}{R}_c& =300\mathrm{m}\\ L_s& =100\mathrm{m}\end{array}}$$ Then $${\displaystyle {\theta }_s=\frac{L_s}{2R_c}=\frac{100}{2\times 300}=\frac{1}{6}\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{n}}$$ and $${\displaystyle 2R_c{L}_s=60000}$$

We scale down the Euler spiral by √60000, i.e. 100√6 to normalized Euler spiral that has: $${\displaystyle \begin{array}{cc}{R}_c^{\prime }& =\frac{3}{\sqrt{6}}\mathrm{m}\\ L_s^{\prime }& =\frac{1}{\sqrt{6}}\mathrm{m}\\ 2R_c^{\prime }{L}_s^{\prime }& =2\times \frac{3}{\sqrt{6}}\times \frac{1}{\sqrt{6}}\\ & =1\end{array}}$$ and $${\displaystyle {\theta }_s=\frac{L_s^{\prime }}{2R_c^{\prime }}=\frac{\frac{1}{\sqrt{6}}}{2\times \frac{3}{\sqrt{6}}}=\frac{1}{6}\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{n}}$$

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{array}{cc}x& ={\int }_0^L\cos \left(s^2\right)ds\\ y& ={\int }_0^L\sin \left(s^2\right)ds\end{array}}$$ or expressed as power series: $${\displaystyle \begin{array}{cccc}x& ={\sum _{i=0}^{\mathrm{\infty }}\frac{(-1{)}^i}{(2i)!}\frac{s^{4i+1}}{4i+1}|}_0^L& & =\sum _{i=0}^{\mathrm{\infty }}\frac{(-1{)}^i}{(2i)!}\frac{L^{4i+1}}{4i+1}\\ y& ={\sum _{i=0}^{\mathrm{\infty }}\frac{(-1{)}^i}{(2i+1)!}\frac{s^{4i+3}}{4i+3}|}_0^L& & =\sum _{i=0}^{\mathrm{\infty }}\frac{(-1{)}^i}{(2i+1)!}\frac{L^{4i+3}}{4i+3}\end{array}}$$

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{array}{cccc}{x}^{\mathrm{\prime }}& =\underset{L\to \mathrm{\infty }}{lim}{\int }_0^L\cos \left(s^2\right)ds& & =\frac{1}{2}\sqrt{\frac{\pi }{2}}\approx 0.6267\\ y^{\mathrm{\prime }}& =\underset{L\to \mathrm{\infty }}{lim}{\int }_0^L\sin \left(s^2\right)ds& & =\frac{1}{2}\sqrt{\frac{\pi }{2}}\approx 0.6267\end{array}}$$

Normalized Euler spirals have the following properties: $${\displaystyle \begin{array}{cc}2R_c{L}_s& =1\\ {\theta }_s& =\frac{L_s}{2R_c}=L_s^2\end{array}}$$ and $${\displaystyle \begin{array}{cc}\theta & ={\theta }_s\cdot \frac{L^2}{L_s^2}=L^2\\ \frac{1}{R}& =\frac{d\theta }{dL}=2L\end{array}}$$

Note that 2R_cL_s = 1 also means ⁠1/R_c⁠ = 2L_s, in agreement with the last mathematical statement.