Finding extremals of Lagrangian actions
Document Type
Journal Article
Publication Title
Communications in Nonlinear Science and Numerical Simulation
Volume
116
Publisher
Elsevier
School
School of Science
RAS ID
56436
Abstract
Given a smooth m-manifold M, a smooth Lagrangian L:TM→R and endpoints x0,xT∈M, we look for an extremal x:[0,T]→M of the action ∫0TL(x(t),ẋ(t))dt satisfying x(0)=x0 and x(T)=xT. When interpolating between endpoints, this amounts to a 2-point boundary value problem for the Euler–Lagrange equation. Single or multiple shooting is one of the most popular methods to solve boundary value problems, but the efficiency of shooting and the quality of solutions depends heavily on initial guesses. In the present paper, by dividing the interval [0,T] into several sub-intervals, on which extremals can be found efficiently by shooting when good initial guesses are available from the geometry of a variational problem, we then adjust all junctions by finding zeros of vector fields associated with the velocities at junctions with Newton's method. We discuss the cases where L is the difference between kinetic energy and potential, M is a hypersurface in Euclidean space, or M is a Lie group. We make some comparisons in numerical experiments for a double pendulum, for obstacle avoidance by a moving particle on the 2-sphere, and for obstacle avoidance by a planar rigid body.
DOI
10.1016/j.cnsns.2022.106826
Access Rights
free_to_read
Comments
Noakes, L., Zhang, E. (2023). Finding extremals of Lagrangian actions. Communications in Nonlinear Science and Numerical Simulation, 116, article 106826. https://doi.org/10.1016/j.cnsns.2022.106826