Brachistochrone curve

The solution to the brachistochrone problem is not a straight line or some combination thereof but a cycloid.

In mathematics and physics, a brachistochrone curve (from Ancient Greek βράχιστος χρόνος (brákhistos khrónos) 'shortest time'), or curve of fastest descent, is the one lying on plane between a point A and a lower point B, where B is not directly below A, on which a bead slides frictionlessly under the influence of a uniform gravitational field to a given end point in the shortest time. Incidentally, for a given starting point, the brachistochrone curve is the same as the tautochrone curve. More specifically, the solution to the brachistochrone and tautochrone problem are one and the same, the cycloid.

The problem can be solved with the tools from the calculus of variations and optimal control.

The curve is independent of both the mass of the test body and the (local) strength of gravity. Only a parameter is chosen so that the curve fits the starting point A and the ending point B. If the body is given an initial velocity at A, or if friction is taken into account, then the curve that minimizes time will differ from the one described above.

Johann Bernoulli's solution

According to Fermat’s principle, the actual path between two points taken by a beam of light is one that takes the least time. In 1697 Johann Bernoulli used this principle to derive the brachistochrone curve by considering the trajectory of a beam of light in a medium where the speed of light increases following a constant vertical acceleration (that of gravity g).

By the conservation of energy, the instantaneous speed of a body v after falling a height y in a uniform gravitational field is given by:

$v={\sqrt {2gy}}$,

The speed of motion of the body along an arbitrary curve does not depend on the horizontal displacement.

Johann Bernoulli noted that the law of refraction gives a constant of the motion for a beam of light in a medium of variable density:

${\frac {\sin {\theta }}{v}}={\frac {1}{v}}{\frac {dx}{ds}}={\frac {1}{v_{m}}}$,

where v_m is the constant and $\theta$ represents the angle of the trajectory with respect to the vertical.

The equations above lead to two conclusions:

1. At the onset, the angle must be zero when the particle speed is zero. Hence, the brachistochrone curve is tangent to the vertical at the origin.
2. The speed reaches a maximum value when the trajectory becomes horizontal and the angle θ = 90°.

Assuming for simplicity that the particle (or the beam) with coordinates (x,y) departs from the point (0,0) and reaches maximum speed after falling a vertical distance D:

$v_{m}={\sqrt {2gD}}$.

Rearranging terms in the law of refraction and squaring gives:

$v_{m}^{2}dx^{2}=v^{2}ds^{2}=v^{2}(dx^{2}+dy^{2})$

which can be solved for dx in terms of dy:

$dx={\frac {v\,dy}{\sqrt {v_{m}^{2}-v^{2}}}}$.

Substituting from the expressions for v and v_m above gives:

$dx={\sqrt {\frac {y}{D-y}}}dy$

which is the differential equation of an inverted cycloid generated by a circle of diameter D.

Jakob Bernoulli's solution

Johann's brother Jakob showed how 2nd differentials can be used to obtain the condition for least time. A modernized version of the proof is as follows. If we make a negligible deviation from the path of least time, then, for the differential triangle formed by the displacement along the path and the horizontal and vertical displacements,

$ds^{2}=dx^{2}+dy^{2}$.

On differentiation with dy fixed we get,

$2ds\ d^{2}s=2dx\ d^{2}x$.

And finally rearranging terms gives,

${\frac {dx}{ds}}d^{2}x=d^{2}s=v\ d^{2}t$

where the last part is the displacement for given change in time for 2nd differentials. Now consider the changes along the two neighboring paths in the figure below for which the horizontal separation between paths along the central line is d²x (the same for both the upper and lower differential triangles). Along the old and new paths, the parts that differ are,

$d^{2}t_{1}={\frac {1}{v_{1}}}{\frac {dx_{1}}{ds_{1}}}d^{2}x$

$d^{2}t_{2}={\frac {1}{v_{2}}}{\frac {dx_{2}}{ds_{2}}}d^{2}x$

For the path of least times these times are equal so for their difference we get,

$d^{2}t_{2}-d^{2}t_{1}=0={\bigg (}{\frac {1}{v_{2}}}{\frac {dx_{2}}{ds_{2}}}-{\frac {1}{v_{1}}}{\frac {dx_{1}}{ds_{1}}}{\bigg )}d^{2}x$

And the condition for least time is,

${\frac {1}{v_{2}}}{\frac {dx_{2}}{ds_{2}}}={\frac {1}{v_{1}}}{\frac {dx_{1}}{ds_{1}}}$

History

Johann Bernoulli posed the problem of the brachistochrone to the readers of Acta Eruditorum in June, 1696. He published his solution in the journal in May of the following year, and noted that the solution is the same curve as Huygens's tautochrone curve. After deriving the differential equation for the curve by the method given above, he went on to show that it does yield a cycloid. However, his proof is marred by his use of a single constant instead of the three constants, v_m, 2g and D, above.

Bernoulli allowed six months for the solutions but none were received during this period. At the request of Leibniz, the time was publicly extended for a year and a half. On 29 January 1697 the challenge was received by Isaac Newton, who found it in his mail, in a letter from Johann Bernoulli, when he arrived home from the Royal Mint at 4 p.m., and stayed up all night to solve it and mailed the solution anonymously by the next post. Upon reading the solution, Bernoulli immediately recognized its author, exclaiming that he "recognizes a lion from his claw marks." This story gives some idea of Newton's power, since Johann Bernoulli took two weeks to solve it.

At the end, five mathematicians responded with solutions: Newton, Jakob Bernoulli (Johann's brother), Gottfried Leibniz, Ehrenfried Walther von Tschirnhaus and Guillaume de l'Hôpital. Four of the solutions (excluding l'Hôpital's) were published in the same edition of the journal as Johann Bernoulli's. In his paper Jakob Bernoulli gave a proof of the condition for least time similar to that above before showing that its solution is a cycloid. According to Newtonian scholar Tom Whiteside,

In an attempt to outdo his brother, Jakob Bernoulli created a harder version of the brachistochrone problem. In solving it, he developed new methods that were refined by Leonhard Euler into what the latter called (in 1766) the calculus of variations. Joseph-Louis Lagrange did further work that resulted in modern infinitesimal calculus.

Earlier, in 1638, Galileo had tried to solve a similar problem for the path of the fastest descent from a point to a wall in his Two New Sciences. He draws the conclusion (Third Day, Theorem 22, Prop. 36) that the arc of a circle is faster than any number of its chords,

"From the preceding it is possible to infer that the quickest path of all [lationem omnium velocissimam], from one point to another, is not the shortest path, namely, a straight line, but the arc of a circle.

...

Consequently the nearer the inscribed polygon approaches a circle the shorter is the time required for descent from A to C. What has been proven for the quadrant holds true also for smaller arcs; the reasoning is the same."

We are warned earlier in the Two New Sciences (just after Theorem 6) of possible fallacies and the need for a "higher science." In this dialogue Galileo reviews his own work. The actual solution to Galileo's problem is half a cycloid. Galileo studied the cycloid and gave it its name, but the connection between it and his problem had to wait for advances in mathematics.

See also

• Beltrami identity
• Calculus of variations
• Catenary
• Cycloid
• Tautochrone curve
• Uniformly accelerated motion

External links

• "Brachistochrone", Encyclopedia of Mathematics, EMS Press, 2001 [1994]
• Weisstein, Eric W. "Brachistochrone Problem". MathWorld.
• Courbe Brachistrocrone (in French with excellent animated examples)
• The Brachistochrone, Whistler Alley Mathematics.
• Table IV from Bernoulli's article in Acta Eruditorum 1697
• Brachistochrones by Michael Trott and Brachistochrone Problem by Okay Arik, Wolfram Demonstrations Project.
• The Brachistochrone problem at MacTutor
• Geodesics Revisited — Introduction to geodesics including two ways of derivation of the equation of geodesic with brachistochrone as a special case of a geodesic.
• Optimal control solution to the Brachistochrone problem in Python.
• The straight line, the catenary, the brachistochrone, the circle, and Fermat Unified approach to some geodesics.