Algorithm for numerical integration of the rigid-body equations of motion

Abstract

A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in terms of either principal axes or quaternions. Due to specific features of the algorithm, orthonormality and unit norms of the orientational variables are integrals of motion, despite an approximate character of the produced trajectories. It is shown that the method presented appears to be the most efficient among all known algorithms of such a kind.

4 pages, 1 figure

References in corpus (3)

• On the numerical integration of motion for rigid polyatomics: The modified quaternion approach
• Numerical integration of the equations of motion for rigid polyatomics: The matrix method
• On the reaction field for interaction site models of polar systems

Cited by in corpus (12)

• Modeling a spheroidal microswimmer and cooperative swimming in thin films
• On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
• Algorithm for molecular dynamics simulations of spin liquids
• Enhanced rotational motion of spherical squirmer in polymer solutions
• Computational probes of molecular motion in the Lewis and Whanstrom model for ortho-terphenyl
• Oedometric Test, Bauer's Law and the Micro-Macro Connection for a Dry Sand
• Alignment-Induced Self-Organization of Autonomously Steering Microswimmers: Turbulence, Vortices, and Jets
• Rotating hematite cube chains
• Stability of interlocked self-propelled dumbbell clusters
• A new approach to the integration of rotational motion in systems with interacting rigid bodies
• Macroscopic emulation of microscopic magnetic particle systems
• A new leapfrog integrator of rotational motion. The revised angular-momentum approach