Lagrangians and Euler-Lagrange Equations in Partial Differential Equations
📂Partial Differential Equations Lagrangians and Euler-Lagrange Equations in Partial Differential Equations Definition Lagrangian
Let’s assume a smooth function L : R n × R n → R L : \mathbb{R}^n \times \mathbb{R}^n \to \mathbb{R} L : R n × R n → R is given. This is called the Lagrangian and is denoted as follows.
L = L ( v , x ) = L ( v 1 , … , v n , x 1 , … , x n ) v , x ∈ R n D v L = ( L v 1 , … , L v n ) , D x L = ( L x 1 , … , L x n )
L = L(v,x)=L(v_{1}, \dots, v_{n}, x_{1}, \dots, x_{n}) \quad v,x\in \mathbb{R}^{n}
\\ D_{v}L = (L_{v_{1}}, \dots, L_{v_{n}}), \quad D_{x}L = (L_{x_{1}}, \dots, L_{x_{n}})
L = L ( v , x ) = L ( v 1 , … , v n , x 1 , … , x n ) v , x ∈ R n D v L = ( L v 1 , … , L v n ) , D x L = ( L x 1 , … , L x n )
The reason for using variables v , x v, x v , x is because, in physics, each variable actually signifies velocity and position.
Action , Admissible class
For two fixed points x , y ∈ R n x,y \in \mathbb{R}^{n} x , y ∈ R n and time t > 0 t>0 t > 0 , the following defined functional I I I is called action .
I [ w ( ⋅ ) ] : = ∫ 0 t L ( w ˙ ( s ) , w ( s ) ) d s ( ˙ = d d s )
I[ \mathbf{w}(\cdot)] := \int_{0}^tL(\dot{\mathbf{w}}(s), \mathbf{w}(s) ) ds \quad \left( \dot{ }=\dfrac{d}{ds}\right)
I [ w ( ⋅ )] := ∫ 0 t L ( w ˙ ( s ) , w ( s )) d s ( ˙ = d s d )
In this case, the function w ( ⋅ ) = ( w 1 ( ⋅ ) , ⋯ , w n ( ⋅ ) ) \mathbf{w}(\cdot)=\big( w^1(\cdot), \cdots, w^n(\cdot) \big) w ( ⋅ ) = ( w 1 ( ⋅ ) , ⋯ , w n ( ⋅ ) ) is an element of the admissible class defined as follows A \mathcal{A} A .
A : = { w ( ⋅ ) = ∈ C 2 ( [ 0 , t ] ; R n ) ∣ w ( 0 ) = y , w ( t ) = x }
\mathcal{A} := \left\{ \mathbf{w}(\cdot) = \in C^2 \big([0,t];\mathbb{R}^n \big) \ \big| \ \mathbf{w}(0)=y, \mathbf{w}(t)=x\right\}
A := { w ( ⋅ ) =∈ C 2 ( [ 0 , t ] ; R n ) w ( 0 ) = y , w ( t ) = x }
In other words, A \mathcal{A} A means the set of all paths that are twice continuously differentiable , starting from position y y y to x x x as time flows from 0 0 0 to t t t .
Description The goal of the calculus of variations is to find x ∈ A \mathbf{x} \in \mathcal{A} x ∈ A that makes the integral value of action I I I minimum. This x \mathbf{x} x is called the minimizer of I I I .
I [ x ( ⋅ ) ] = inf w ( ⋅ ) ∈ A I [ w ( ⋅ ) ]
I[ \mathbf{x} (\cdot) ] = \inf_{\mathbf{w}(\cdot)\in \mathcal{A}} I[\mathbf{w}(\cdot)]
I [ x ( ⋅ )] = w ( ⋅ ) ∈ A inf I [ w ( ⋅ )]
Such x \mathbf{x} x is sought because the path that minimizes the action of the Lagrangian is actually the path along which an object moves . In essence, this is for understanding the motion of an object, which fundamentally equates to solving F = m a F=ma F = ma . In classical mechanics, the Lagrangian specifically refers to the difference between kinetic energy and potential energy.
Regarding the determination of minimizers, there’s the following theorem.
Theorem Assume x ( ⋅ ) ∈ A \mathbf{x}(\cdot) \in \mathcal{A} x ( ⋅ ) ∈ A is the minimizer of action I I I . Then, x ( ⋅ ) \mathbf{x}(\cdot) x ( ⋅ ) satisfies the following equation.
− d d s [ D v L ( x ˙ ( s ) , x ( s ) ) ] + D x L ( x ˙ ( s ) , x ( s ) ) = 0 ( 0 ≤ s ≤ t )
-\dfrac{d}{ds} \Big[ D_{v}L\big( \dot{\mathbf{x}}(s), \mathbf{x}(s) \big) \Big] + D_{x}L\big( \dot{\mathbf{x}}(s), \mathbf{x}(s)\big)=0 \quad (0 \le s \le t)
− d s d [ D v L ( x ˙ ( s ) , x ( s ) ) ] + D x L ( x ˙ ( s ) , x ( s ) ) = 0 ( 0 ≤ s ≤ t )
This equation is called the Euler-Lagrange equations .
It’s important to note that although a minimizer satisfies the Euler-Lagrange equations, satisfying the Euler-Lagrange equations does not necessarily mean it’s a minimizer. Similar to how differentiating at the minimum value gives 0 0 0 , but not every point where differentiation gives 0 0 0 is a minimum. In this sense, x ( ⋅ ) ∈ A \mathbf{x}(\cdot) \in \mathcal{A} x ( ⋅ ) ∈ A satisfying the Euler-Lagrange equations is called a critical point of I I I . Therefore, while a minimizer is a critical point, not all critical points are minimizers.
Proof Assume x ∈ A \mathbf{x} \in \mathcal{A} x ∈ A is the minimizer of action I I I .
Step 1.
Let’s presume the function y : [ 0 , t ] → R n , y ( ⋅ ) = ( y 1 ( ⋅ ) , ⋯ , y n ( ⋅ ) ) \mathbf{y} : [0,t] \to \mathbb{R}^{n}, \mathbf{y}(\cdot) = (y^1(\cdot), \cdots, y^n(\cdot) ) y : [ 0 , t ] → R n , y ( ⋅ ) = ( y 1 ( ⋅ ) , ⋯ , y n ( ⋅ )) satisfies the following equation as a smooth function.
y ( 0 ) = y ( t ) = 0
\begin{equation}
\mathbf{y}(0)=\mathbf{y}(t)=\mathbf{0}
\label{eq1}
\end{equation}
y ( 0 ) = y ( t ) = 0
And for any τ ∈ R \tau \in \mathbb{R} τ ∈ R , define w ( ⋅ ) \mathbf{w}(\cdot) w ( ⋅ ) as follows.
w ( ⋅ ) : = x ( ⋅ ) + τ y ( ⋅ ) ∈ A
\mathbf{w}(\cdot) : = \mathbf{x}(\cdot) + \tau \mathbf{y}(\cdot) \in \mathcal{A}
w ( ⋅ ) := x ( ⋅ ) + τ y ( ⋅ ) ∈ A
Then w \mathbf{w} w represents a path that has the same start and end points as x \mathbf{x} x but differs by τ y ( ⋅ ) \tau \mathbf{y}(\cdot) τ y ( ⋅ ) in between. Furthermore, since x ( ⋅ ) \mathbf{x}(\cdot) x ( ⋅ ) is the minimizer of I I I , the following equation holds.
I [ x ( ⋅ ) ] ≤ I [ w ( ⋅ ) ] = I [ x ( ⋅ ) + τ y ( ⋅ ) ] = : i ( τ )
I[\mathbf{x}(\cdot)] \le I[\mathbf{w}(\cdot)]=I[\mathbf{x}(\cdot) + \tau \mathbf{y}(\cdot)] =: i(\tau)
I [ x ( ⋅ )] ≤ I [ w ( ⋅ )] = I [ x ( ⋅ ) + τ y ( ⋅ )] =: i ( τ )
Also, the function i i i has a minimum value at τ = 0 \tau=0 τ = 0 by definition of a minimizer. Therefore, if the derivative of i i i exists, it is i ′ ( 0 ) = 0 i^{\prime}(0)=0 i ′ ( 0 ) = 0 .
Step 2.
As defined above, i i i is as follows.
i ( τ ) = ∫ 0 t L ( x ˙ ( s ) + τ y ˙ ( s ) , x ( s ) + τ y ( s ) ) d s
i(\tau) = \int_{0} ^t L\big( \dot{\mathbf{x}}(s) + \tau \dot{\mathbf{y}}(s), \mathbf{x}(s)+ \tau \mathbf{y}(s) \big)ds
i ( τ ) = ∫ 0 t L ( x ˙ ( s ) + τ y ˙ ( s ) , x ( s ) + τ y ( s ) ) d s
Since L L L , y \mathbf{y} y are smooth functions, by differentiating i i i , we get the following.
i ′ ( τ ) = ∫ 0 t ∑ i = 1 n [ L v i ( x ˙ + τ y ˙ , x + τ y ) y ˙ i + L x i ( x ˙ + τ y ˙ , x + τ y ) y i ] d s
i^{\prime}(\tau) = \int_{0}^t \sum_{i=1}^{n} \left[ L_{v_{i}} ( \dot{\mathbf{x}} + \tau \dot{\mathbf{y}}, \mathbf{x}+ \tau \mathbf{y} )\dot{y}^i + L_{x_{i}} ( \dot{\mathbf{x}} + \tau \dot{\mathbf{y}}, \mathbf{x}+ \tau \mathbf{y} ) y^i \right] ds
i ′ ( τ ) = ∫ 0 t i = 1 ∑ n [ L v i ( x ˙ + τ y ˙ , x + τ y ) y ˙ i + L x i ( x ˙ + τ y ˙ , x + τ y ) y i ] d s
Substituting τ = 0 \tau=0 τ = 0 , we obtain the following result from Step 1.
0 = i ′ ( 0 ) = ∫ 0 t ∑ i = 1 n [ L v i ( x ˙ , x ) y ˙ i + L x i ( x ˙ , x ) y i ] d s
0=i^{\prime}(0) = \int_{0}^t \sum_{i=1}^{n} \left[ L_{v_{i}} ( \dot{\mathbf{x}} , \mathbf{x})\dot{y}^i + L_{x_{i}} ( \dot{\mathbf{x}} , \mathbf{x} )y^i \right] ds
0 = i ′ ( 0 ) = ∫ 0 t i = 1 ∑ n [ L v i ( x ˙ , x ) y ˙ i + L x i ( x ˙ , x ) y i ] d s
Applying partial integration to each term of L v i y ˙ i L_{v_{i}}\dot{y}^i L v i y ˙ i , and by the assumption y ( 0 ) = y ( t ) = 0 \mathbf{y}(0)=\mathbf{y}(t)=\mathbf{0} y ( 0 ) = y ( t ) = 0 , we get the following.
∫ 0 t L v i y ˙ i d s = L v i y i ] 0 t − ∫ 0 t d d s L v i y i = ∫ 0 t − d d s L v i y i
\int_{0}^t L_{v_{i}}\dot{y}^i ds = \left. L_{v_{i}}y^i \right]_{0}^t- \int_{0}^t \dfrac{d}{ds}L_{v_{i}}y^i=\int_{0}^t-\dfrac{d}{ds} L_{v_{i}}y^i
∫ 0 t L v i y ˙ i d s = L v i y i ] 0 t − ∫ 0 t d s d L v i y i = ∫ 0 t − d s d L v i y i
Therefore, the following equation holds.
0 = i ′ ( 0 ) = ∑ i = 1 n ∫ 0 t [ − d d s L v i ( x ˙ , x ) + L x i ( x ˙ , x ) ] y i d s
0=i^{\prime}(0) = \sum_{i=1}^{n} \int_{0}^t \left[ -\dfrac{d}{ds} L_{v_{i}} ( \dot{\mathbf{x}} , \mathbf{x})\ + L_{x_{i}} ( \dot{\mathbf{x}} , \mathbf{x} ) \right]y^i ds
0 = i ′ ( 0 ) = i = 1 ∑ n ∫ 0 t [ − d s d L v i ( x ˙ , x ) + L x i ( x ˙ , x ) ] y i d s
The result of Step 2. satisfies for all smooth functions y : [ 0 , t ] → R n \mathbf{y} : [0,t] \to \mathbb{R}^n y : [ 0 , t ] → R n that fulfill ( eq1 ) \eqref{eq1} ( eq1 ) . Therefore, the value inside the brackets must be 0 0 0 . Thus, the following is true.
− d d s L v i ( x ˙ , x ) + L x i ( x ˙ , x ) = 0
-\dfrac{d}{ds} L_{v_{i}}( \dot{\mathbf{x}}, \mathbf{x} ) +L_{x_{i}}( \dot{\mathbf{x}}, \mathbf{x}) =0
− d s d L v i ( x ˙ , x ) + L x i ( x ˙ , x ) = 0
■
This result leads to the Hamiltonian equations .
See Also