SPIRAL: an efficient algorithm for the integration of the equation of rotational motion. (English) Zbl 07826294
Summary: We introduce Spiral, a third-order integration algorithm for the rotational motion of extended bodies. It requires only one force calculation per time step, does not require quaternion normalization at each time step, and can be formulated for both leapfrog and synchronous integration schemes, making it compatible with many particle simulation codes. The stability and precision of Spiral exceed those of state-of-the-art algorithms currently used in popular DEM codes such as Yade, MercuryDPM, LIGGGHTS, PFC, and more, at only slightly higher computational cost. Also, beyond DEM, we see potential applications in all numerical simulations that involve the 3D rotation of extended bodies.
Keywords:
numerical integration of the equation of rotational motion; particle simulation; discrete element method (DEM); molecular dynamics (MD)Software:
MuJoCo; NeuralPDE.jl; LAMMPS; AirSim; PFC3D; TRAC-IK; DifferentialEquations.jl; GranOO; MercuryDPM; MUSEN; Gromacs; Blaze-DEMGPUReferences:
| [1] | Abe, S.; Place, D.; Mora, P., A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics. Pure Appl. Geophys., 2265-2277 (2004) |
| [2] | Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E., GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX, 19-25 (2015) |
| [3] | Altair 2023 ‘EDEM - discrete element method (DEM) software’ (2023) |
| [4] | Andre, D.; Charles, J.; Iordanoff, I., 3D Discrete Element Workbench for Highly Dynamic Thermo-Mechanical Analysis: GranOO. Discrete Element Model and Simulation of Continuous Materials Behavior (2015), Wiley |
| [5] | Banchmark, Benchmarktools.jl (2023) |
| [6] | Beeson, P.; Ames, B., TRAC-IK: an open-source library for improved solving of generic inverse kinematics, 928-935 |
| [7] | Boyce, W.; DiPrima, R.; Meade, D., Elementary Differential Equations and Boundary Value Problems (2017), Wiley |
| [8] | Buss, S. R., Accurate and efficient simulation of rigid-body rotations. J. Comput. Phys., 377-406 (2000) · Zbl 1035.70500 |
| [9] | Chen, C. W.; Kouh, J. S.; Tsai, J. F., Modeling and simulation of an AUV simulator with guidance system. IEEE J. Ocean. Eng., 2, 211-225 (2013) |
| [10] | Dosta, M.; Skorych, V., MUSEN: an open-source framework for GPU-accelerated DEM simulations. SoftwareX (2020) |
| [11] | Fincham, D., Leapfrog rotational algorithms. Mol. Simul., 165-178 (1992) |
| [12] | Goldstein, H., Classical Mechanics (2002), Pearson Education |
| [13] | Govender, N.; Wilke, D.; Kok, S., Blaze-DEMGPU: modular high performance DEM framework for the GPU architecture. SoftwareX, 62-66 (2016) |
| [14] | Grassia, F. S., Practical parameterization of rotations using the exponential map. J. Graph. Tools, 29-48 (1998) |
| [15] | Hamilton, W. R., LXXVIII. On quaternions; or on a new system of imaginaries in algebra. Lond. Edinb. Dublin Philos. Mag. J. Sci., 489-495 (1844) |
| [16] | Hockney, R.; Eastwood, J., Computer Simulation Using Particles (1988), CRC Press · Zbl 0662.76002 |
| [17] | (2008a) PFC3D — Particle Flow Code in 3 Dimensions, Version 4.0, User’s Manual (2008), Itasca: Itasca Minneapolis |
| [18] | (2021) PFC — Particle Flow Code in 2 and 3 Dimensions, Version 7.0, Documentation Set of Version 7.00.132 (2021), Itasca: Itasca Minneapolis |
| [19] | Jia, Y. B., Quaternions: Com S 477/577 notes (2017), Iowa State University, Technical report |
| [20] | Johnson, M. P., Exploiting Quaternions to Support Expressive Interactive Character Motion (2003), Massachusetts Institute of Technology, Ph.D. thesis |
| [21] | Johnson, S. M.; Williams, J. R.; Cook, B. K., Quaternion-based rigid body rotation integration algorithms for use in particle methods. Int. J. Numer. Methods Eng., 1303-1313 (2008) · Zbl 1159.70300 |
| [22] | Kloss, C.; Goniva, C.; Hager, A.; Amberger, S.; Pirker, S., Models, algorithms and validation for opensource DEM and CFD-DEM. Prog. Comput. Fluid Dyn., 2-3, 140-152 (2012) |
| [23] | Kobilarov, M.; Crane, K.; Desbrun, M., Lie group integrators for animation and control of vehicles. ACM Trans. Graph. (2009) |
| [24] | Kuipers, J., Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality (1999), Princeton University Press, Princeton paperbacks · Zbl 1053.70001 |
| [25] | Kutta, W., Beitrag zur näherungsweisen integration totaler differentialgleichungen. Z. Math. Phys., 435-453 (1901) · JFM 32.0316.02 |
| [26] | Longuski, J. M.; Tsiotras, P., Analytical solutions for a spinning rigid body subject to time-varying body-fixed torques, Part I: constant axial torque. J. Appl. Mech., 970-975 (1993) |
| [27] | Manchester, Z. R.; Peck, M. A., Quaternion variational integrators for spacecraft dynamics. J. Guid. Control Dyn., 69-76 (2016) |
| [28] | Neto, N.; Bellucci, L., A new algorithm for rigid body molecular dynamics. Chem. Phys., 259-268 (2006) |
| [29] | Omelyan, I. P., Algorithm for numerical integration of the rigid-body equations of motion. Phys. Rev. E, 1169-1172 (1998) |
| [30] | Ostanin, I.; Angelidakis, V.; Plath, T.; Pourandi, S.; Thornton, A.; Weinhart, T., Rigid clumps in the mercurydpm particle dynamics code. Comput. Phys. Commun. (2024) |
| [31] | Pons, A.; Cirak, F., Quaternion variational integration for inertial maneuvering in a biomimetic unmanned aerial vehicle. ASME Lett. Dyn. Syst. Control (2023) |
| [32] | Rackauckas, C.; Nie, Q., Differentialequations.jl – a performant and feature-rich ecosystem for solving differential equations in Julia. J. Open Res. Softw. (2017) |
| [33] | Romano, M., Exact analytic solution for the rotation of a rigid body having spherical ellipsoid of inertia and subjected to a constant torque. Celest. Mech. Dyn. Astron., 181-189 (2008) · Zbl 1335.70019 |
| [34] | Roots, Roots.jl (2023) |
| [35] | Rucker, C., Integrating rotations using nonunit quaternions. IEEE Robot. Autom. Lett., 2979-2986 (2018) |
| [36] | Sabatini, A. M., Quaternion-based strap-down integration method for applications of inertial sensing to gait analysis. Med. Biol. Eng. Comput., 94-101 (2005) |
| [37] | Seelen, L. J.H.; Padding, J. T.; Kuipers, J. A.M., Improved quaternion-based integration scheme for rigid body motion. Acta Mech., 3381-3389 (2016) · Zbl 1380.70016 |
| [38] | Shah, S.; Dey, D.; Lovett, C.; Kapoor, A., AirSim: high-fidelity visual and physical simulation for autonomous vehicles, 621-635 |
| [39] | Shu, C. W.; Osher, S., Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys., 439-471 (1988) · Zbl 0653.65072 |
| [40] | Smilauer, V.; Angelidakis, V.; Catalano, E.; Caulk, R.; Chareyre, B.; Chèvremont, W.; Dorofeenko, S.; Duriez, J.; Dyck, N.; Elias, J.; Er, B.; Eulitz, A.; Gladky, A.; Guo, N.; Jakob, C.; Kneib, F.; Kozicki, J.; Marzougui, D.; Maurin, R.; Modenese, C.; Pekmezi, G.; Scholtès, L.; Sibille, L.; Stransky, J.; Sweijen, T.; Thoeni, K.; Yuan, C., Yade documentation The Yade Project (2021) |
| [41] | Su, J.; Schroeder, C.; Fedkiw, R., Energy stability and fracture for frame rate rigid body simulations, 155-164 |
| [42] | Syamlal, M.; Rogers, W.; O‘Brien, T. J., MFiX documentation theory guide (1993), USDOE Morgantown Energy Technology Center (METC): USDOE Morgantown Energy Technology Center (METC) WV (United States), Technical report |
| [43] | Tadano, S.; Takeda, R.; Miyagawa, H., Three dimensional gait analysis using wearable acceleration and gyro sensors based on quaternion calculations. Sensors, 9321-9343 (2013) |
| [44] | Tayebi, A.; McGilvray, S., Attitude stabilization of a VTOL quadrotor aircraft. IEEE Trans. Control Syst. Technol., 562-571 (2006) |
| [45] | Thompson, A. P.; Aktulga, H. M.; Berger, R.; Bolintineanu, D. S.; Brown, W. M.; Crozier, P. S.; in ’t Veld, P. J.; Kohlmeyer, A.; Moore, S. G.; Nguyen, T. D.; Shan, R.; Stevens, M. J.; Tranchida, J.; Trott, C.; Plimpton, S. J., LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comput. Phys. Commun. (2022) · Zbl 1516.74108 |
| [46] | Thornton, C., Granular Dynamics, Contact Mechanics and Particle System Simulations: A DEM Study. Particle Technology Series (2015), Springer International Publishing |
| [47] | Todorov, E.; Erez, T.; Tassa, Y., MuJoCo: a physics engine for model-based control, 5026-5033 |
| [48] | Tsiotras, P.; Longuski, J. M., Analytic solution of Euler’s equations of motion for an asymmetric rigid body. J. Appl. Mech., 149-155 (1996) · Zbl 0874.70008 |
| [49] | (2023), unreal ‘Unity’ |
| [50] | (2023), unreal ‘Unreal engine’ |
| [51] | Valenti, R. G.; Dryanovski, I.; Xiao, J., A linear Kalman filter for MARG orientation estimation using the algebraic quaternion algorithm. IEEE Trans. Instrum. Meas., 467-481 (2016) |
| [52] | Verlet, L., Computer “experiments” on classical fluids. I. thermodynamical properties of Lennard-Jones molecules. Phys. Rev., 98-103 (1967) |
| [53] | vulkan, Vulkan (2023) |
| [54] | Weinhart, T.; Orefice, L.; Post, M.; van Schrojenstein Lantman, M. P.; Denissen, I. F.; Tunuguntla, D. R.; Tsang, J.; Cheng, H.; Shaheen, M. Y.; Shi, H.; Rapino, P.; Grannonio, E.; Losacco, N.; Barbosa, J.; Jing, L.; Alvarez Naranjo, J. E.; Roy, S.; den Otter, W. K.; Thornton, A. R., Fast, flexible particle simulations — an introduction to MercuryDPM. Comput. Phys. Commun. (2020) · Zbl 07678501 |
| [55] | woodem, WooDEM (2023) |
| [56] | Wu, Y.; Hu, X.; Hu, D.; Li, T.; Lian, J., Strapdown inertial navigation system algorithms based on dual quaternions. IEEE Trans. Aerosp. Electron. Syst., 110-132 (2005) |
| [57] | Zhao, F.; van Wachem, B. G.M., A novel quaternion integration approach for describing the behaviour of non-spherical particles. Acta Mech., 3091-3109 (2013) · Zbl 1356.70009 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
