5.3: Applications of Euler’s Equation
( \newcommand{\kernel}{\mathrm{null}\,}\)
Example 5.3.1: Shortest distance between two points
Consider the path lies in the x−y plane.
The infinitessimal length of arc is
ds=√dx2+dy2=[√1+(dydx)2]dx
Then the length of the arc is J=∫21ds=∫21[√1+(dydx)2]dx
The function f is
f=√1+(y′)2
Therefore
∂f∂y=0
and
∂f∂y′=y′√1+(y′)2
Inserting these into Euler’s equation (5.2.13) gives
0+ddx(y′√1+(y′)2)=0
that is
y′√1+(y′)2=constant=C
This is valid if
y′=C√1−C2=a
Therefore
y=ax+b
which is the equation of a straight line in the plane. Thus the shortest path between two points in a plane is a straight line between these points, as is intuitively obvious. This stationary value obviously is a minimum.
This trivial example of the use of Euler’s equation to determine an extremum value has given the obvious answer. It has been presented here because it provides a proof that a straight line is the shortest distance in a plane and illustrates the power of the calculus of variations to determine extremum paths.
Example 5.3.2: Brachistochrone problem
The Brachistochrone problem involves finding the path having the minimum transit time between two points. The Brachistochrone problem stimulated the development of the calculus of variations by John Bernoulli and Euler. For simplicity, take the case of frictionless motion in the x−y plane with a uniform gravitational field acting in the ˆy direction, as shown in the adjacent figure. The question is what constrained path will result in the minimum transit time between two points (x1y1) and (x2y2).
Consider that the particle of mass m starts at the origin x1=0,y1=0 with zero velocity. Since the problem conserves energy and assuming that initially E=KE+PE=0 then
12mv2−mgy=0
That is v=√2gy
The transit time is given by
t=∫x2x1dsv=∫x2x1√dx2+dy2√2gy=∫x2x1√(1+x′2)2gydy
where x′≡dxdy. Note that, in this example, the independent variable has been chosen to be y and the dependent variable is x(y).
The function f of the integral is
f=1√2g√(1+x′2)y
Factor out the constant √2g term, which does not affect the final equation, and note that
∂f∂x=0∂f∂x′=x′√y(1+(x′)2)
Therefore Euler’s equation gives
0+ddy(x′√y(1+(x′)2))=0
or
x′√y(1+(x′)2)=constant=1√2a
That is
x′2y(1+(x′)2)=12a
This may be rewritten as
x=∫y2y1ydy√2ay−y2
Change the variable to y=a(1−cosθ) gives that dy=asinθdθ, leading to the integral
x=∫a(1−cosθ)dθ
or
x=a(θ−sinθ)+constant
The parametric equations for a cycloid passing through the origin are
x=a(θ−sinθ)y=a(1−cosθ)
which is the form of the solution found. That is, the shortest time between two points is obtained by constraining the motion of the mass to follow a cycloid shape. Thus the mass first accelerates rapidly by falling down steeply and then follows the curve and coasts upward at the end. The elapsed time is obtained by inserting the above parametric relations for x and y, in terms of θ, into the transit time integral giving t=√agθ where a and θ are fixed by the end point coordinates. Thus the time to fall from starting with zero velocity at the cusp to the minimum of the cycloid is π√ag. If y2=y1=0 then x2=2πa which defines the shape of the cycloid and the minimum time is 2π√ag=√2πx2g. If the mass starts with a non-zero initial velocity, then the starting point is not at the cusp of the cycloid, but down a distance d such that the kinetic energy equals the potential energy difference from the cusp.
A modern application of the Brachistochrone problem is determination of the optimum shape of the low-friction emergency chute that passengers slide down to evacuate a burning aircraft. Bernoulli solved the problem of rapid evacuation of an aircraft two centuries before the first flight of a powered aircraft.
Example 5.3.3: Minimal travel cost
Assume that the cost of flying an aircraft at height z is e−κz per unit distance of flight-path, where κ is a positive constant. Consider that the aircraft flies in the (x,z)-plane from the point (−a,0) to the point (a,0) where z=0 corresponds to ground level, and where the z-axis points vertically upwards. Find the extremal for the problem of minimizing the total cost of the journey.
The differential arc-length element of the flight path ds can be written as
ds=√dx2+dz2=√1+z′2dx
where z′≡dzdx. Thus the cost integral to be minimized is
C=∫+a−ae−κzds=∫+a−ae−κz√1+z′2dx
The function of this integral is
f=e−κz√1+z′2
The partial differentials required for the Euler equations are
ddx∂f∂z′=z′′e−κz√1+z′2−κz′2e−κz√1+z′2−z′′z′2e−κz(1+z′2)3/2∂f∂z=−κe−κz√1+z′2
Therefore Euler’s equation equals
∂f∂z−ddx∂f∂z′=−κe−κz√1+z′2−z′′e−κz√1+z′2+κz′2e−κz√1+z′2+z′′z′2e−κz(1+z′2)3/2=0
This can be simplified by multiplying the radical to give
−κ−2κz′2−κz′4−z′′−z′′z′2+κz′2+κz′4+z′′z′2=0
Cancelling terms gives
z′′+κ(1+z′2)=0
Separating the variables leads to
arctanz′=∫dz′z′2+1=−∫κdx=−κz+c1
Integration gives
z(x)=∫x−adz=∫x−atan(c1−κx)dx=ln(cos(c1−κx))−ln(cos(c1+κa))κ+c2=ln(cos(c1−κx)cos(c1+κa))κ+c2
Using the initial condition that z(−a)=0 gives c2=0. Similarly the final condition z(a)=0 implies that c1=0. Thus Euler’s equation has determined that the optimal trajectory that minimizes the cost integral C is
z(x)=1κln(cos(κx)cos(κa))
This example is typical of problems encountered in economics.