Reliable event detection for Taylor methods in astrodynamics
arXiv:2204.09948 · doi:10.1093/mnras/stac1092
Abstract
We present a novel approach for the detection of events in systems of ordinary differential equations. The new method combines the unique features of Taylor integrators with state-of-the-art polynomial root finding techniques to yield a novel algorithm ensuring strong event detection guarantees at a modest computational overhead. Detailed tests and benchmarks focused on problems in astrodynamics and celestial mechanics (such as collisional N-body systems, spacecraft dynamics around irregular bodies accounting for eclipses, computation of Poincare' sections, etc.) show how our approach is superior in both performance and detection accuracy to strategies commonly employed in modern numerical integration works. The new algorithm is available in our open source Taylor integration package heyoka.
Accepted for publication in MNRAS
References in corpus (3)
- Turbulence sets the length scale for planetesimal formation: Local 2D simulations of streaming instability and planetesimal formation
- Revisiting high-order Taylor methods for astrodynamics and celestial mechanics
- A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives