Integrating rotation from angular velocity

Abstract\ud The integration of the rotation from a given angular velocity is often required in practice. The present paper explores how the choice of the parametrization of rotation, when employed in conjuction with different numerical time-integration schemes, effects the accuracy and the computational efficiency. Three rotation parametrizations – the rotational vector, the Argyris tangential vector and the rotational quaternion – are combined with three different numerical time-integration schemes, including classical explicit Runge–Kutta method and the novel midpoint rule proposed here. The key result of the study is the assessment of the integration errors of various parametrization–integration method combinations. In order to assess the errors, we choose a time-dependent function corresponding to a rotational vector, and derive the related exact time-dependent angular velocity. This is then employed in the numerical solution as the data. The resulting numerically integrated approximate rotations are compared with the analytical solution. A novel global solution error norm for discrete solutions given by a set of values at chosen time-points is employed. Several characteristic angular velocity functions, resulting in small, finite and fast oscillating rotations are studied

A minimal-variable symplectic integrator on spheres

We construct a symplectic, globally defined, minimal-coordinate, equivariant integrator on products of 2-spheres. Examples of corresponding Hamiltonian systems, called spin systems, include the reduced free rigid body, the motion of point vortices on a sphere, and the classical Heisenberg spin chain, a spatial discretisation of the Landau-Lifschitz equation. The existence of such an integrator is remarkable, as the sphere is neither a vector space, nor a cotangent bundle, has no global coordinate chart, and its symplectic form is not even exact. Moreover, the formulation of the integrator is very simple, and resembles the geodesic midpoint method, although the latter is not symplectic

Energy-Entropy-Momentum integration of discrete thermo-visco-elastic dynamics.

A novel time integration scheme is presented for the numerical solution of the dynamics of discrete systems consisting of point masses and thermo-visco-elastic springs. Even considering fully coupled constitutive laws for the elements, the obtained solutions strictly preserve the two laws of thermo dynamics and the symmetries of the continuum evolution equations. Moreover, the unconditional control over the energy and the entropy growth have the effect of stabilizing the numerical solution, allowing the use of larger time steps than those suitable for comparable implicit algorithms. Proofs for these claims are provided in the article as well as numerical examples that illustrate the performance of the method

Long-time behaviour of discretizations of the simple pendulum equation

We compare the performance of several discretizations of the simple pendulum equation in a series of numerical experiments. The stress is put on the long-time behaviour. We choose for the comparison numerical schemes which preserve the qualitative features of solutions (like periodicity). All these schemes are either symplectic maps or integrable (preserving the energy integral) maps, or both. We describe and explain systematic errors (produced by any method) in numerical computations of the period and the amplitude of oscillations. We propose a new numerical scheme which is a modification of the discrete gradient method. This discretization preserves (almost exactly) the period of small oscillations for any time step.Comment: 41 pages, including 18 figures and 4 table