Abstract
We introduce SPIRAL, our newly developed third-order integration algorithm, which resolves many of the issues with existing ones and outperforms them in all tested metrics with only minor increases in computational cost. The algorithm does not require quaternion normalization for each timestep and is compatible with both leapfrog and non-leapfrog schemes, making it easy to adapt to most existing simulation frameworks. The integration of rotational motion is a key part of many applications, such as embedded systems, sensors, aircraft control systems, molecular dynamics (MD), and discrete element method simulations (DEM). However, current methods for integrating the rotational motion of particles in simulations are not optimal, particularly for non-spherical particles. This study compares various integration algorithms using a range of metrics, including accuracy against an analytical solution, performance, numerical stability, and energy conservation in a challenging undamped system.
| Original language | English |
|---|---|
| Title of host publication | 9th International Conference on Discrete Element Methods (DEM9): Book of Abstracts |
| Publisher | FAU Erlangen-Nürnberg |
| Pages | 248 |
| Number of pages | 1 |
| DOIs | |
| Publication status | Published - 16 Oct 2023 |
Keywords
- Discrete element method (DEM)
- non-spherical particles
- simulation accuracy
- numerical integration algorithm for rotation
- molecular dynamics (MD)
Access to Document
- Developing SPIRAL: A stable particle rotation integration algorithm
Copyright 2023 the authors. This is an open access abstract published under a Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium, provided the author and source are cited.