Beltrami identity

Template:Short description {{#invoke:sidebar|collapsible | class = plainlist | titlestyle = padding-bottom:0.25em; | pretitle = Part of a series of articles about | title = Calculus | image = $\underset{a}{\overset{b}{\int }}{f}^{\prime }(t)dt=f(b)-f(a)$ | listtitlestyle = text-align:center; | liststyle = border-top:1px solid #aaa;padding-top:0.15em;border-bottom:1px solid #aaa; | expanded = specialized | abovestyle = padding:0.15em 0.25em 0.3em;font-weight:normal; | above =

• Fundamental theorem

• Limits
• Continuity

• Rolle's theorem
• Mean value theorem
• Inverse function theorem

| list2name = differential | list2titlestyle = display:block;margin-top:0.65em; | list2title = Template:Bigger | list2 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | heading1 = Definitions
 | content1 =

• Derivative (generalizations)
• Differential
  • infinitesimal
  • of a function
  • total

 | heading2 = Concepts
 | content2 =

• Differentiation notation
• Second derivative
• Implicit differentiation
• Logarithmic differentiation
• Related rates
• Taylor's theorem

 | heading3 = Rules and identities
 | content3 =

• Sum
• Product
• Chain
• Power
• Quotient
• L'Hôpital's rule
• Inverse
• General Leibniz
• Faà di Bruno's formula
• Reynolds

}}

| list3name = integral | list3title = Template:Bigger | list3 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

• Lists of integrals
• Integral transform
• Leibniz integral rule

| heading2 = Definitions

 | content2 =

• Antiderivative
• Integral (improper)
• Riemann integral
• Lebesgue integration
• Contour integration
• Integral of inverse functions

 | heading3 = Integration by
 | content3 =

• Parts
• Discs
• Cylindrical shells
• Substitution (trigonometric, tangent half-angle, Euler)
• Euler's formula
• Partial fractions (Heaviside's method)
• Changing order
• Reduction formulae
• Differentiating under the integral sign
• Risch algorithm

}}

| list4name = series | list4title = Template:Bigger | list4 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

• Geometric (arithmetico-geometric)
• Harmonic
• Alternating
• Power
• Binomial
• Taylor

 | heading2 = Convergence tests
 | content2 =

• Summand limit (term test)
• Ratio
• Root
• Integral
• Direct comparison
• 
  Limit comparison
• Alternating series
• Cauchy condensation
• Dirichlet
• Abel

}}

| list5name = vector | list5title = Template:Bigger | list5 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

• Gradient
• Divergence
• Curl
• Laplacian
• Directional derivative
• Identities

 | heading2 = Theorems
 | content2 =

• Gradient
• Green's
• Stokes'
• Divergence
• generalized Stokes
• Helmholtz decomposition

}}

| list6name = multivariable | list6title = Template:Bigger | list6 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | heading1 = Formalisms
 | content1 =

• Matrix
• Tensor
• Exterior
• Geometric

 | heading2 = Definitions
 | content2 =

• Partial derivative
• Multiple integral
• Line integral
• Surface integral
• Volume integral
• Jacobian
• Hessian

}}

| list7name = advanced | list7title = Template:Bigger | list7 ={{#invoke:sidebar|sidebar|child=yes

 |contentclass=hlist
 | content1 =

• Calculus on Euclidean space
• Generalized functions
• Limit of distributions

}}

| list8name = specialized | list8title = Template:Bigger | list8 =

• Fractional
• Malliavin
• Stochastic
• Variations

| list9name = miscellanea | list9title = Template:Bigger | list9 =

• Precalculus
• History
• Glossary
• List of topics
• Integration Bee
• Mathematical analysis
• Nonstandard analysis

}}

The Beltrami identity, named after Eugenio Beltrami, is a special case of the Euler–Lagrange equation in the calculus of variations.

The Euler–Lagrange equation serves to extremize action functionals of the form

$I[u]=\underset{a}{\overset{b}{\int }}L[x,u(x),u^{\prime }(x)]dx,$

where $a$ and $b$ are constants and $u^{\prime }(x)=\frac{du}{dx}$.

If $\frac{\partial L}{\partial x}=0$, then the Euler–Lagrange equation reduces to the Beltrami identity,

$L-u^{\prime }\frac{\partial L}{\partial u^{\prime }}=C,$

where Template:Math is a constant.

Derivation

By the chain rule, the derivative of Template:Math is

$\frac{dL}{dx}=\frac{\partial L}{\partial x}\frac{dx}{dx}+\frac{\partial L}{\partial u}\frac{du}{dx}+\frac{\partial L}{\partial u^{\prime }}\frac{d{u}^{\prime }}{dx}.$

Because $\frac{\partial L}{\partial x}=0$, we write

$\frac{dL}{dx}=\frac{\partial L}{\partial u}{u}^{\prime }+\frac{\partial L}{\partial u^{\prime }}{u}^{″}.$

We have an expression for $\frac{\partial L}{\partial u}$ from the Euler–Lagrange equation,

$\frac{\partial L}{\partial u}=\frac{d}{dx}\frac{\partial L}{\partial u^{\prime }}$

that we can substitute in the above expression for $\frac{dL}{dx}$ to obtain

$\frac{dL}{dx}=u^{\prime }\frac{d}{dx}\frac{\partial L}{\partial u^{\prime }}+u^{″}\frac{\partial L}{\partial u^{\prime }}.$

By the product rule, the right side is equivalent to

$\frac{dL}{dx}=\frac{d}{dx}\left(u^{\prime }\frac{\partial L}{\partial u^{\prime }}\right).$

By integrating both sides and putting both terms on one side, we get the Beltrami identity,

$L-u^{\prime }\frac{\partial L}{\partial u^{\prime }}=C.$

Applications

Solution to the brachistochrone problem

The solution to the brachistochrone problem is the cycloid.

An example of an application of the Beltrami identity is the brachistochrone problem, which involves finding the curve $y=y(x)$ that minimizes the integral

$I[y]=\underset{0}{\overset{a}{\int }}\sqrt{\frac{1+y{\text{'}}^2}{y}}dx.$

The integrand

$L(y,y^{\prime })=\sqrt{\frac{1+y{\text{'}}^2}{y}}$

does not depend explicitly on the variable of integration $x$, so the Beltrami identity applies,

$L-y^{\prime }\frac{\partial L}{\partial y^{\prime }}=C.$

Substituting for $L$ and simplifying,

$y(1+y{\text{'}}^2)=1/C^2\text{(constant)},$

which can be solved with the result put in the form of parametric equations

$x=A(\varphi -\sin \varphi )$

$y=A(1-\cos \varphi )$

with $A$ being half the above constant, $\frac{1}{2C^2}$, and $\varphi$ being a variable. These are the parametric equations for a cycloid.

Solution to the catenary problem

A chain hanging from points forms a catenary.

Consider a string with uniform density $\mu$ of length $l$ suspended from two points of equal height and at distance $D$. By the formula for arc length, $$l=\underset{S}{\int }dS=\underset{s_1}{\overset{s_2}{\int }}\sqrt{1+y{\text{'}}^2}dx,$$ where $S$ is the path of the string, and $s_1$ and $s_2$ are the boundary conditions.

The curve has to minimize its potential energy $$U=\underset{S}{\int }g\mu y\cdot dS=\underset{s_1}{\overset{s_2}{\int }}g\mu y\sqrt{1+y{\text{'}}^2}dx,$$ and is subject to the constraint $$\underset{s_1}{\overset{s_2}{\int }}\sqrt{1+y{\text{'}}^2}dx=l,$$ where $g$ is the force of gravity.

Because the independent variable $x$ does not appear in the integrand, the Beltrami identity may be used to express the path of the string as a separable first order differential equation

$$L-y\prime \frac{\partial L}{\partial y\prime }=\mu gy\sqrt{1+y{\prime }^2}+\lambda \sqrt{1+y{\prime }^2}-[\mu gy\frac{y{\prime }^2}{\sqrt{1+y{\prime }^2}}+\lambda \frac{y{\prime }^2}{\sqrt{1+y{\prime }^2}}]=C,$$ where $\lambda$ is the Lagrange multiplier.

It is possible to simplify the differential equation as such: $$\frac{g\rho y-\lambda }{\sqrt{1+y{\text{'}}^2}}=C.$$

Solving this equation gives the hyperbolic cosine, where $C_0$ is a second constant obtained from integration

$$y=\frac{C}{\mu g}\cosh [\frac{\mu g}{C}(x+C_0)]-\frac{\lambda }{\mu g}.$$

The three unknowns $C$, $C_0$, and $\lambda$ can be solved for using the constraints for the string's endpoints and arc length $l$, though a closed-form solution is often very difficult to obtain.