Abstract:Variational integrators applied to degenerate Lagrangians that are linear in the velocities are two-step methods. The system of modified equations for a two-step method consists of the principal modified equation and one additional equation describing parasitic oscillations. We observe that a Lagrangian for the principal modified equation can be constructed using the same technique as in the case of non-degenerate Lagrangians. Furthermore, we construct the full system of modified equations by doubling the dimension of the discrete system in such a way that the principal modified equation of the extended system coincides with the full system of modified equations of the original system. We show that the extended discrete system is Lagrangian, which leads to a construction of a Lagrangian for the full system of modified equations.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is: to construct correction equations for variational integrators in degenerate Lagrangian systems (i.e., velocity - linear Lagrangians). Specifically, the author focuses on those Lagrangians that are linear in velocity, and the Euler - Lagrange equations corresponding to these Lagrangians are first - order differential equations. However, when discretized using variational integrators, the resulting difference equations are second - order, thus forming a two - step method. This leads to the emergence of parasitic solutions (parasitic solutions).
### Main problem decomposition:
1. **Construction of correction equations**:
- The author hopes to find a corrected Lagrangian for this special Lagrangian system such that its Euler - Lagrange equation matches the correction equation of the original system.
- For non - degenerate Lagrangian systems, there already exists a well - developed theory of correction equations, but for degenerate Lagrangian systems, the situation is more complex because the Legendre transformation is not invertible and the correction equations cannot be directly derived from the theory of symplectic integrators.
2. **Treatment of parasitic solutions**:
- Since the variational integrator is a two - step method in this case, it may introduce parasitic solutions, which will affect the long - term behavior of the numerical solution.
- The author hopes to deal with this problem by extending the discrete system method so that the extended system is still variational and a complete system of correction equations can be constructed.
3. **Specific applications and verification**:
- The paper verifies the effectiveness of the proposed method through specific examples (such as planar point - vortex dynamics and toy problems).
- For example, in the toy problem, the author shows how to apply the midpoint rule and the trapezoidal rule to discretize the Lagrangian and analyzes the behavior of parasitic solutions.
### Summary of mathematical formulas:
- Form of the Lagrangian:
\[
L(q, \dot{q})=\langle\alpha(q), \dot{q}\rangle - H(q)
\]
- Definitions of matrices \(A(q)\) and \(A_{\text{skew}}(q)\):
\[
A(q)=\alpha^{\prime}(q)=\left(\frac{\partial\alpha_{i}(q)}{\partial q_{j}}\right)_{i, j = 1,\dots,N}
\]
\[
A_{\text{skew}}(q)=A(q)^{T}-A(q)
\]
- Euler - Lagrange equation:
\[
\dot{q}=A_{\text{skew}}(q)^{-1}H^{\prime}(q)^{T}
\]
By solving these problems, the author hopes to provide an effective method to study and improve the performance of variational integrators when dealing with velocity - linear Lagrangians, especially in terms of long - term behavior and parasitic solution control.