Skip to content

v1.0.0 - Release: Core Astrodynamics & Adaptive Solvers

Choose a tag to compare

@mukid1805 mukid1805 released this 02 Sep 21:36
· 28 commits to main since this release

Overview

Major milestone release introducing high-performance adaptive numerical integration with local truncation error control, complete Circular Restricted Three-Body Problem (CR3BP) dynamics, and validation benchmarks.


Key Features & Scenarios

  • Adaptive Dormand-Prince Integrator: Implemented embedded RK5(4) solver with local truncation error (LTE) step-size scaling and First Same As Last (FSAL) efficiency (core/integrators.py).
  • CR3BP Cislunar Dynamics & Equilibria: Non-dimensional rotating-frame equations of motion with root-finding solvers for $L_1\text{-}L_5$ Lagrange points and continuous Jacobi constant ($C_J$) tracking (core/cr3bp.py).
  • Scenario 6 Free-Return Upgrade: Benchmarked Earth-Moon figure-8 free-return trajectory comparing classical RK4 against adaptive RK45 in both synodic and inertial frames (examples/ex06_cislunar_free_return.py).
  • Scenario 7 Earth–Mars Transfer Optimization: Integrated RK45 propagation into the porkchop-derived Lambert transfer arc, reducing deep-space cruise evaluations by >99% (examples/ex07_earth_mars_transfer.py).
  • Scenario 8 Multi-Leg Flyby Dynamics: Patched-conic Earth–Venus–Mars (EVM) gravity-assist trajectory verification updated with dual-leg RK4 vs. RK45 telemetry (examples/ex08_gravity_assist_transfer.py).

Verification

  • Full Test Suite Coverage (25/25 Passing): Added tests/test_cr3bp.py and tests/test_forces.py to validate equilibrium acceleration zeroing, analytical triangular geometry, Jacobi conservation, and individual perturbation forces.
  • Extreme Orbit Energy Conservation: Verified specific mechanical energy preservation on a high-eccentricity GTO ($e = 0.7265$) yielding a relative energy drift of $\Delta \mathcal{E}/\mathcal{E}_0 = 8.14 \times 10^{-10}$.
  • Empirical Speedup Telemetry: Demonstrated runtime speedups ranging from $19.3\times$ to $79.4\times$ across interplanetary cruise and cislunar free-return arcs while maintaining $0.00\text{ km}$ arrival miss distances.
  • Continuous Integration: Verified automated test matrices passing across Python 3.10, 3.11, and 3.12 on both Ubuntu and Windows runners (tests.yml).