---
title: Jet Bundle Geometry in Variational Integrators
url: https://www.emergentmind.com/topics/jet-bundle-geometry
type: topic
---

# Jet Bundle Geometry in Variational Integrators

Variational integrators are a class of geometric numerical methods for mechanical systems whose defining property is derivation from a discretization of Hamilton’s principle of stationary action, rather than from direct discretization of the differential equations of motion. This approach yields numerical schemes that preserve key geometric invariants of the original system, such as symplecticity and momentum maps, and exhibit superior long-time energy behavior. Variational integrators have been systematically developed for finite- and infinite-dimensional systems, including general Lagrangian and Hamiltonian dynamics, Lie group systems, systems with constraints, time-dependent and stochastic systems, and applications across areas such as optimal control, plasma physics, electric circuits, and flexible multibody mechanics.

## 1. Foundations: Discrete Variational Principle and Construction

The fundamental construction is based on the discrete analog of Hamilton’s variational principle. Given a continuous Lagrangian system $(Q, L: TQ \rightarrow \mathbb{R})$, the discrete Lagrangian $L_d: Q \times Q \to \mathbb{R}$ approximates the fixed-time action integral
$$ L_d(q_k, q_{k+1}) \approx \int_0^h L(q(t), \dot{q}(t))\,dt $$
along a solution of the Euler–Lagrange equations interpolating $q(0) = q_k$, $q(h) = q_{k+1}$. The discrete action is
$$ S_d(\{q_k\}) = \sum_{k=0}^{N-1} L_d(q_k, q_{k+1}) $$
and the discrete variational principle requires stationarity with respect to interior variations, yielding the discrete Euler–Lagrange (DEL) equations:
$$ D_2 L_d(q_{k-1}, q_k) + D_1 L_d(q_k, q_{k+1}) = 0 $$
with $D_i$ denoting partial differentiation with respect to the $i$th argument. The resulting two-step map $(q_{k-1}, q_k) \mapsto (q_k, q_{k+1})$ defines the variational integrator [1102.2685].

Two central strategies for constructing high-order variational integrators are the Galerkin approach—replacing the space of interpolating curves by a finite-dimensional ansatz and then approximating the action by quadrature—and the shooting-based approach—solving the two-point boundary-value problem approximately using a one-step integrator and quadrature. The order of the resulting integrator is governed by the accuracy in matching the exact discrete Lagrangian, as quantified by variational error analysis [1102.2685, 1703.06599].

For systems evolving on Lie groups, the interpolating curves and discretizations are constructed in local charts or using group retractions (exponential or Cayley map), leading to Lie group variational integrators [1402.3327, 2502.05692].

## 2. Structure-Preservation: Symplecticity, Momentum, and Energy Behavior

A distinguishing feature of variational integrators is their exact preservation of a discrete symplectic form. The DEL update map $F_{L_d}:(q_{k-1},q_k)\mapsto(q_k,q_{k+1})$ satisfies
$$ F_{L_d}^* \Omega_{L_d} = \Omega_{L_d} $$
where $\Omega_{L_d}$ is the discrete symplectic form built from the discrete Poincaré–Cartan one-forms [1211.4534]. This discrete symplecticity leads to bounded numerical energy error over exponentially long times; the numerical energy typically exhibits bounded oscillations rather than systematic drift [1102.2685, 1211.4534].

If the discrete Lagrangian is invariant under a group action (e.g., for mechanical systems with symmetries), a discrete version of Noether’s theorem guarantees preservation of discrete momentum maps:
$$ J_d(q_k, q_{k+1}) = \langle D_2L_d(q_k, q_{k+1}), \xi_Q(q_{k+1}) \rangle $$
which remain constant along the discrete flow [1102.2685, 1402.3327].

In time-dependent and forced/dissipative systems, backward error analysis of the symplectic discrete flow in extended phase space shows that discrete energies related to the modified Hamiltonian are conserved up to exponentially small errors [2202.01471].

## 3. High-Order and Spectral Variational Integrators

By varying the interpolation trial space and the quadrature scheme, variational integrators admit systematic construction to arbitrarily high order. Spectral (Galerkin) variational integrators use polynomial or spectral approximations in each time step; under appropriate analytic assumptions, this yields geometric (exponential in degree) convergence rates in both trajectory and invariant errors [1211.4534]. For systems on Lie groups, analogous Lie group spectral integrators have been developed [1402.3327].

In classical (flat) mechanical systems, equivalent high-order variational integrators can also be recast as continuous-stage or partitioned Runge–Kutta (csPRK, SPRK) methods that match the variational structure and order-optimality [1809.06825, 1001.1408].

Taylor variational integrators exploit high-order Taylor expansions both in solving the boundary value problem (shooting formulation) and evaluating quadrature nodes—resulting in integrators that, for given computational cost, can achieve one order higher than traditional shooting integrators [1703.06599].

## 4. Constraints, Extensions, and Generalizations

Variational integrators generalize systematically to systems with holonomic, nonholonomic, or higher-order constraints through the discrete Lagrange multiplier and augmented discrete action principle [1306.0298]. For higher-order Lagrangian systems, the discrete phase space is $Q^{k+1}$, and the variational principle yields higher-order discrete Euler–Lagrange equations that are symplectic and momentum-preserving [1306.0298].

For nonvariational or degenerate systems (e.g., certain electric circuits, advection equations, or underactuated mechanical systems), the framework extends by constructing a formal (adjoint) Lagrangian, embedding the dynamics into a larger system that admits a Lagrangian structure, then applying the variational discretization to the extended system [1103.1859, 1412.2011, 1307.5665].

Non-autonomous (time-dependent) Lagrangian and Hamiltonian systems are accommodated by introducing a family of discrete Lagrangians $L_d^k$ indexed by time step, with symplecticity demonstrated in the extended (time, configuration) space [2202.01471].

## 5. Stochastic and Dissipative Generalizations

Stochastic variational integrators are constructed by discretizing stochastic action functionals (with Stratonovich integrals) and producing discrete Lagrange–d’Alembert principles. Such integrators preserve symplecticity, momentum, and yield improved long-time statistical behavior compared to non-geometric integrators [1904.06205, 1609.00463]. For systems with multiplicative noise, the construction yields stochastic symplectic Runge–Kutta methods as special cases [1609.00463]. Applications include Vlasov–Fokker–Planck equations and collisional kinetic plasma models [1904

Source: https://www.emergentmind.com/topics/jet-bundle-geometry