Eirik Hoel Høiseth

Elena Celledoni, Marta Farré Puiggalí, Eirik Hoel Høiseth et al.

We introduce energy-preserving integrators for nonholonomic mechanical systems. We will see that the nonholonomic dynamics is completely determined by a triple $({\mathcal D}^*, Π, \mathcal{H})$, where ${\mathcal D}^*$ is the dual of the vector bundle determined by the nonholonomic constraints, $Π$ is an almost-Poisson bracket (the nonholonomic bracket) and $\mathcal{H}:{\mathcal D}^{*}\rightarrow \mathbb{R}$ is a Hamiltonian function. For this triple, we can apply energy-preserving integrators, in particular, we show that discrete gradients can be used in the numerical integration of nonholonomic dynamics. By construction, we achieve preservation of the constraints and of the energy of the nonholonomic system. Moreover, to facilitate their applicability to complex systems which cannot be easily transformed into the aforementioned almost-Poisson form, we rewrite our integrators using just the initial information of the nonholonomic system. The derived procedures are tested on several examples: A chaotic quartic nonholonomic mechanical system, the Chaplygin sleigh system, the Suslov problem and a continuous gearbox driven by an asymmetric pendulum. Their performace is compared with other standard methods in nonholonomic dynamics, and their merits verified in practice.

Elena Celledoni, Eirik Hoel Høiseth, Nataliya Ramzina

A rigid body model for the dynamics of a marine vessel, used in simulations of offshore pipe-lay operations, gives rise to a set of ordinary differential equations with controls. The system is input-output passive. We propose passivity-preserving splitting methods for the numerical solution of a class of problems which includes this system as a special case. We prove the passivity-preservation property for the splitting methods, and we investigate stability and energy behaviour in numerical experiments. Implementation is discussed in detail for a special case where the splitting gives rise to the subsequent integration of two completely integrable flows. The equations for the attitude are reformulated on $SO(3)$ using rotation matrices rather than local parametrizations with Euler angles.

Elena Celledoni, Eirik Hoel Høiseth

In this paper we design discrete port-Hamiltonian systems systematically in two different ways, by applying discrete gradient methods and splitting methods respectively. The discrete port-Hamiltonian systems we get satisfy a discrete notion of passivity, which lets us, by choosing the input appropriately, make them globally asymptotically stable with respect to an equilibrium point. We test methods designed using the discrete gradient approach in numerical experiments, and the results are encouraging when compared to relevant existing integrators of identical order.