The geometric discretisation of the Suslov problem: a case study of consistency for nonholonomic integrators

From MaRDI portal
Publication:524604

DOI10.3934/DCDS.2017182zbMATH Open1401.70022arXiv1605.02666OpenAlexW3104273999MaRDI QIDQ524604FDOQ524604


Authors: Luis C. García-Naranjo, Fernando Jiménez Edit this on Wikidata


Publication date: 3 May 2017

Published in: Discrete and Continuous Dynamical Systems (Search for Journal in Brave)

Abstract: Geometric integrators for nonholonomic systems were introduced by Cort'es and Mart'inez in [Nonholonomic integrators, Nonlinearity, 14, 2001] by proposing a discrete Lagrange-D'Alembert principle. Their approach is based on the definition of a discrete Lagrangian Ld and a discrete constraint space Dd. There is no recipe to construct these objects and the performance of the integrator is sensitive to their choice. Cort'es and Mart'inez claim that choosing Ld and Dd in a consistent manner with respect to a finite difference map is necessary to guarantee an approximation of the continuous flow within a desired order of accuracy. Although this statement is given without proof, similar versions of it have appeared recently in the literature. We evaluate the importance of the consistency condition by comparing the performance of two different geometric integrators for the nonholonomic Suslov problem, only one of which corresponds to a consistent choice of Ld and Dd. We prove that both integrators produce approximations of the same order, and, moreover, that the non-consistent discretisation outperforms the other in numerical experiments and in terms of energy preservation. Our results indicate that the consistency of a discretisation might not be the most relevant feature to consider in the construction of nonholonomic geometric integrators.


Full work available at URL: https://arxiv.org/abs/1605.02666




Recommendations




Cites Work


Cited In (3)





This page was built for publication: The geometric discretisation of the Suslov problem: a case study of consistency for nonholonomic integrators

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q524604)