Fermat's spiral

From HandWiki
Short description: Spiral that surrounds equal area per turn
Fermat's spiral: a>0, one branch [math]\displaystyle{ r=+a\sqrt{\varphi} }[/math]
Fermat's spiral, both branches

A Fermat's spiral or parabolic spiral is a plane curve with the property that the area between any two consecutive full turns around the spiral is invariant. As a result, the distance between turns grows in inverse proportion to their distance from the spiral center, contrasting with the Archimedean spiral (for which this distance is invariant) and the logarithmic spiral (for which the distance between turns is proportional to the distance from the center). Fermat spirals are named after Pierre de Fermat.[1]

Their applications include curvature continuous blending of curves,[1] modeling plant growth and the shapes of certain spiral galaxies, and the design of variable capacitors, solar power reflector arrays, and cyclotrons.

Coordinate representation

Polar

The representation of the Fermat spiral in polar coordinates (r, φ) is given by the equation [math]\displaystyle{ r=\pm a\sqrt{\varphi} }[/math] for φ ≥ 0.

The two choices of sign give the two branches of the spiral, which meet smoothly at the origin. If the same variables were reinterpreted as Cartesian coordinates, this would be the equation of a parabola with horizontal axis, which again has two branches above and below the axis, meeting at the origin.

Cartesian

The Fermat spiral with polar equation [math]\displaystyle{ r=\pm a\sqrt{\varphi} }[/math] can be converted to the Cartesian coordinates (x, y) by using the standard conversion formulas x = r cos φ and y = r sin φ. Using the polar equation for the spiral to eliminate r from these conversions produces parametric equations for one branch of the curve:

[math]\displaystyle{ \begin{cases} x(\varphi) = + a\sqrt{\varphi} \cos(\varphi) \\ y(\varphi) = + a\sqrt{\varphi} \sin(\varphi) \end{cases} }[/math]

and the second one

[math]\displaystyle{ \begin{cases} x(\varphi) = - a\sqrt{\varphi} \cos(\varphi) \\ y(\varphi) = - a\sqrt{\varphi} \sin(\varphi) \end{cases} }[/math]

They generate the points of branches of the curve as the parameter φ ranges over the positive real numbers.

For any (x, y) generated in this way, dividing x by y cancels the aφ parts of the parametric equations, leaving the simpler equation x/y = cot φ. From this equation, substituting φ by φ = r2/a2 (a rearranged form of the polar equation for the spiral) and then substituting r by r = x2 + y2 (the conversion from Cartesian to polar) leaves an equation for the Fermat spiral in terms of only x and y: [math]\displaystyle{ \frac{x}{y}=\cot\left( \frac{x^2 + y^{2}}{a^{2}} \right). }[/math] Because the sign of a is lost when it is squared, this equation covers both branches of the curve.

Geometric properties

A Fermat's spiral divides the plane into two connected and congruent regions (diagram: black and white)

Division of the plane

A complete Fermat's spiral (both branches) is a smooth double point free curve, in contrast with the Archimedean and hyperbolic spiral. Like a line or circle or parabola, it divides the plane into two connected regions.

Definition of sector (light blue) and polar slope angle α

Polar slope

From vector calculus in polar coordinates one gets the formula

[math]\displaystyle{ \tan\alpha = \frac{r'}{r} }[/math]

for the polar slope and its angle α between the tangent of a curve and the corresponding polar circle (see diagram).

For Fermat's spiral r = aφ one gets

[math]\displaystyle{ \tan\alpha=\frac{1}{2\varphi}. }[/math]

Hence the slope angle is monotonely decreasing.

Curvature

From the formula

[math]\displaystyle{ \kappa = \frac{r^2 + 2(r')^2 - r\,r''}{\left(r^2+(r')^2\right)^\frac32} }[/math]

for the curvature of a curve with polar equation r = r(φ) and its derivatives

[math]\displaystyle{ \begin{align} r' &= \frac{a}{2\sqrt{\varphi}}=\frac{a^2}{2r}\\ r''&= -\frac{a}{4\sqrt{\varphi}^3}=-\frac{a^4}{4r^3} \end{align} }[/math]

one gets the curvature of a Fermat's spiral: [math]\displaystyle{ \kappa(r) = \frac{2r\left(4r^4+3a^4\right)}{\left(4r^4+a^4\right)^\frac32}. }[/math]

At the origin the curvature is 0. Hence the complete curve has at the origin an inflection point and the x-axis is its tangent there.

Area between arcs

The area of a sector of Fermat's spiral between two points (r(φ1), φ1) and (r(φ2), φ2) is

[math]\displaystyle{ \begin{align} \underline A&=\frac{1}{2}\int_{\varphi_1}^{\varphi_2} r(\varphi)^2\, d\varphi \\ &=\frac{1}{2}\int_{\varphi_1}^{\varphi_2} a^2 \varphi\, d\varphi \\ &=\frac{a^2}{4}\left(\varphi_2^2-\varphi_1^2\right) \\ &=\frac{a^2}{4}\left(\varphi_2+\varphi_1\right)\left(\varphi_2-\varphi_1\right). \end{align} }[/math]
Fermat's spiral: area between neighbored arcs

After raising both angles by 2π one gets

[math]\displaystyle{ \overline A= \frac{a^2}{4}\left(\varphi_2+\varphi_1 +4\pi\right)\left(\varphi_2-\varphi_1\right)=\underline A + a^2\pi\left(\varphi_2-\varphi_1\right). }[/math]

Hence the area A of the region between two neighboring arcs is [math]\displaystyle{ A=a^2\pi\left(\varphi_2-\varphi_1\right). }[/math] A only depends on the difference of the two angles, not on the angles themselves.

For the example shown in the diagram, all neighboring stripes have the same area: A1 = A2 = A3.

This property is used in electrical engineering for the construction of variable capacitors.[2]

the regions in between (white, blue, yellow) have all the same area, which is equal to the area of the drawn circle.

Special case due to Fermat

In 1636, Fermat wrote a letter [3] to Marin Mersenne which contains the following special case:

Let φ1 = 0, φ2 = 2π; then the area of the black region (see diagram) is A0 = a2π2, which is half of the area of the circle K0 with radius r(2π). The regions between neighboring curves (white, blue, yellow) have the same area A = 2a2π2. Hence:

  • The area between two arcs of the spiral after a full turn equals the area of the circle K0.

Arc length

The length of the arc of Fermat's spiral between two points (r(φi), φi) can be calculated by the integral:

[math]\displaystyle{ \begin{align} L&=\int_{\varphi_1}^{\varphi_2}\sqrt{\left(r^\prime(\varphi)\right)^2+r^2(\varphi)}\,d\varphi \\ &=\frac{a}{2}\int_{\varphi_1}^{\varphi_2}\sqrt{\frac{1}{\varphi}+4\varphi}\,d\varphi . \end{align} }[/math]

This integral leads to an elliptical integral, which can be solved numerically.

The arc length of the positive branch of the Fermat's spiral from the origin can also be defined by hypergeometric functions 2F1(a, b; c; z) and the incomplete beta function B(z; a, b):[4]

[math]\displaystyle{ \begin{align} L &= a \cdot \sqrt{\varphi} \cdot \operatorname{_{2}F_{1}}\left( -\tfrac12,\, \tfrac14;\, \tfrac54;\, -4 \cdot \varphi^{2} \right)\\ &= a \cdot\frac{1 - i}{8} \cdot \operatorname{B}\left( -4 \cdot \varphi^{2};\, \tfrac14,\, \tfrac32 \right)\\ \end{align} }[/math]

The inversion of Fermat's spiral (green) is a lituus (blue)

Circle inversion

The inversion at the unit circle has in polar coordinates the simple description (r, φ) ↦ (1/r, φ).

  • The image of Fermat's spiral r = aφ under the inversion at the unit circle is a lituus spiral with polar equation [math]\displaystyle{ r=\frac{1}{a\sqrt{\varphi}}. }[/math] When φ = 1/a2, both curves intersect at a fixed point on the unit circle.
  • The tangent (x-axis) at the inflection point (origin) of Fermat's spiral is mapped onto itself and is the asymptotic line of the lituus spiral.

The golden ratio and the golden angle

In disc phyllotaxis, as in the sunflower and daisy, the mesh of spirals occurs in Fibonacci numbers because divergence (angle of succession in a single spiral arrangement) approaches the golden ratio. The shape of the spirals depends on the growth of the elements generated sequentially. In mature-disc phyllotaxis, when all the elements are the same size, the shape of the spirals is that of Fermat spirals—ideally. That is because Fermat's spiral traverses equal annuli in equal turns. The full model proposed by H. Vogel in 1979[5] is

[math]\displaystyle{ \begin{align} r &= c \sqrt{n},\\ \theta &= n \times 137.508^\circ, \end{align} }[/math]

where θ is the angle, r is the radius or distance from the center, and n is the index number of the floret and c is a constant scaling factor. The angle 137.508° is the golden angle which is approximated by ratios of Fibonacci numbers.[6]

The pattern of florets produced by Vogel's model (central image). The other two images show the patterns for slightly different values of the angle.

The resulting spiral pattern of unit disks should be distinguished from the Doyle spirals, patterns formed by tangent disks of geometrically increasing radii placed on logarithmic spirals.

Solar plants

Fermat's spiral has also been found to be an efficient layout for the mirrors of concentrated solar power plants.[7]

See also

References

  1. 1.0 1.1 Lekkas, Anastasios M.; Dahl, Andreas R.; Breivik, Morten; Fossen, Thor I. (2013). "Continuous-Curvature Path Generation Using Fermat's Spiral". Modeling, Identification and Control 34 (4): 183–198. ISSN 1890-1328. http://www.mic-journal.no/pdf/2013/MIC-2013-4-3.pdf. 
  2. Wicke, Fritz (2013). Einführung in die höhere Mathematik. Springer-Verlag. p. 414. ISBN 978-3-662-36804-6. 
  3. Tannery, Paul, ed. "Lettre de Fermat à Mersenne du 3 juin 1636". Œuvres de Fermat. 3. p. 277. https://gallica.bnf.fr/ark:/12148/bpt6k62145354/f303.vertical. 
  4. Weisstein, Eric W.. "Fermat's Spiral". http://mathworld.wolfram.com/FermatsSpiral.html. 
  5. Vogel, H. (1979). "A better way to construct the sunflower head". Mathematical Biosciences 44 (3–4): 179–189. doi:10.1016/0025-5564(79)90080-4. 
  6. Prusinkiewicz, Przemyslaw; Lindenmayer, Aristid (1990). The Algorithmic Beauty of Plants. Springer-Verlag. pp. 101–107. ISBN 978-0-387-97297-8. https://archive.org/details/algorithmicbeaut0000prus/page/101. 
  7. Noone, Corey J.; Torrilhon, Manuel; Mitsos, Alexander (December 2011). "Heliostat Field Optimization: A New Computationally Efficient Model and Biomimetic Layout". Solar Energy 86 (2): 792–803. doi:10.1016/j.solener.2011.12.007. 

Further reading

External links