C1-continuous space-time discretization based on Hamilton's law of varying action

We develop a class of C1-continuous time integration methods that are applicable to conservative problems in elastodynamics. These methods are based on Hamilton's law of varying action. From the action of the continuous system we derive a spatially and temporally weak form of the governing equilibrium equations. This expression is first discretized in space, considering standard finite elements. The resulting system is then discretized in time, approximating the displacement by piecewise cubic Hermite shape functions. Within the time domain we thus achieve C1-continuity for the displacement field and C0-continuity for the velocity field. From the discrete virtual action we finally construct a class of one-step schemes. These methods are examined both analytically and numerically. Here, we study both linear and nonlinear systems as well as inherently continuous and discrete structures. In the numerical examples we focus on one-dimensional applications. The provided theory, however, is general and valid also for problems in 2D or 3D. We show that the most favorable candidate -- denoted as p2-scheme -- converges with order four. Thus, especially if high accuracy of the numerical solution is required, this scheme can be more efficient than methods of lower order. It further exhibits, for linear simple problems, properties similar to variational integrators, such as symplecticity. While it remains to be investigated whether symplecticity holds for arbitrary systems, all our numerical results show an excellent long-term energy behavior.Comment: slightly condensed the manuscript, added references, numerical results unchange

Stochastic Variational Integrators

This paper presents a continuous and discrete Lagrangian theory for stochastic Hamiltonian systems on manifolds. The main result is to derive stochastic governing equations for such systems from a critical point of a stochastic action. Using this result the paper derives Langevin-type equations for constrained mechanical systems and implements a stochastic analog of Lagrangian reduction. These are easy consequences of the fact that the stochastic action is intrinsically defined. Stochastic variational integrators (SVIs) are developed using a discretized stochastic variational principle. The paper shows that the discrete flow of an SVI is a.s. symplectic and in the presence of symmetry a.s. momentum-map preserving. A first-order mean-square convergent SVI for mechanical systems on Lie groups is introduced. As an application of the theory, SVIs are exhibited for multiple, randomly forced and torqued rigid-bodies interacting via a potential.Comment: 21 pages, 8 figure

Practical use of variational principles for modeling water waves

This paper describes a method for deriving approximate equations for irrotational water waves. The method is based on a 'relaxed' variational principle, i.e., on a Lagrangian involving as many variables as possible. This formulation is particularly suitable for the construction of approximate water wave models, since it allows more freedom while preserving the variational structure. The advantages of this relaxed formulation are illustrated with various examples in shallow and deep waters, as well as arbitrary depths. Using subordinate constraints (e.g., irrotationality or free surface impermeability) in various combinations, several model equations are derived, some being well-known, other being new. The models obtained are studied analytically and exact travelling wave solutions are constructed when possible.Comment: 30 pages, 1 figure, 62 references. Other author's papers can be downloaded at http://www.denys-dutykh.com

Hamilton-Pontryagin Integrators on Lie Groups: Introduction and Structure-Preserving Properties

In this paper structure-preserving time-integrators for rigid body-type mechanical systems are derived from a discrete Hamilton-Pontryagin variational principle. From this principle one can derive a novel class of variational partitioned Runge-Kutta methods on Lie groups. Included among these integrators are generalizations of symplectic Euler and St\"{o}rmer-Verlet integrators from flat spaces to Lie groups. Because of their variational design, these integrators preserve a discrete momentum map (in the presence of symmetry) and a symplectic form. In a companion paper, we perform a numerical analysis of these methods and report on numerical experiments on the rigid body and chaotic dynamics of an underwater vehicle. The numerics reveal that these variational integrators possess structure-preserving properties that methods designed to preserve momentum (using the coadjoint action of the Lie group) and energy (for example, by projection) lack.Comment: 26 pages, 4 figure

An introduction to Lie group integrators -- basics, new developments and applications

We give a short and elementary introduction to Lie group methods. A selection of applications of Lie group integrators are discussed. Finally, a family of symplectic integrators on cotangent bundles of Lie groups is presented and the notion of discrete gradient methods is generalised to Lie groups