Construction and analysis of higher order variational integrators for dynamical systems with holonomic constraints

Wenger T, Ober-Blöbaum S, Leyendecker S (2017)

Publication Language: English

Publication Type: Journal article

Publication year: 2017

Journal

Advances in Computational Mathematics Springer Verlag (Germany)

Book Volume: 43

Pages Range: 1163-1195

Journal Issue: 5

DOI: 10.1007/s10444-017-9520-5

Abstract

In this work, variational integrators of higher order for dynamical systems with holonomic constraints are constructed and analyzed. The construction is based on approximating the configuration and the Lagrange multiplier via different polynomials. The splitting of the augmented Lagrangian in two parts enables the use of different quadrature formulas to approximate the integral of each part. Conditions are derived that ensure the linear independence of the higher order constrained discrete Euler-Lagrange equations and stiff accuracy. Time reversibility is investigated for the discrete flow on configuration level only as for the flow on configuration and momentum level. The fulfillment of the hidden constraints plays an important role for the time reversibility of the presented integrators. The order of convergence is investigated numerically. Order reduction of the momentum and the Lagrange multiplier compared to the order of the configuration occurs in general, but can be avoided by fulfilling the hidden constraints in a simple post processing step. Regarding efficiency versus accuracy a numerical analysis yields that higher orders increase the accuracy of the discrete solution substantially while the computational costs decrease. A comparison to the constrained Galerkin methods in Marsden and West (Acta Numerica 10, 357–514 2001) and the symplectic SPARK integrators of Jay (SIAM Journal on Numerical Analysis 45(5), 1814–1842 2007) reveals that the approach presented here is more general and thus allows for more flexibility in the design of the integrator.

Authors with CRIS profile

Theresa Wenger Institute of Applied Dynamics Sigrid Leyendecker Institute of Applied Dynamics

Related research project(s)

Variational integrators of higher order for constrained systems and variational integrators of mixed order for dynamical systems with multiple time scales Jan. 1, 2016 - March 31, 2022

Involved external institutions

University of Oxford GB United Kingdom (GB)

How to cite

APA:

Wenger, T., Ober-Blöbaum, S., & Leyendecker, S. (2017). Construction and analysis of higher order variational integrators for dynamical systems with holonomic constraints. Advances in Computational Mathematics, 43(5), 1163-1195. https://doi.org/10.1007/s10444-017-9520-5

MLA:

Wenger, Theresa, Sina Ober-Blöbaum, and Sigrid Leyendecker. "Construction and analysis of higher order variational integrators for dynamical systems with holonomic constraints." Advances in Computational Mathematics 43.5 (2017): 1163-1195.

BibTeX: Download