GravDyn: A Python Framework for Gravitational Modeling of Irregular Celestial Bodies Using Polyhedral, mascon, and Series Expansion Methods

Abstract

GravDyn is an open-source Python package for computing gravitational potentials and accelerations around irregular celestial bodies. It implements three approaches within the same workflow: the constant-density polyhedral method, a layered mascon approximation based on tetrahedral decomposition, and the Potential Series Expansion Method (PSEM). The package handles shape-model preprocessing, builds the selected gravity representation, and evaluates potentials and accelerations efficiently. The three implemented methods serve different regimes: the polyhedral model provides a reference solution near the surface, the mascon model supports layered internal density structures, and PSEM offers fast evaluation outside the Brillouin sphere once the polynomial coefficients have been generated. Validation tests against the polyhedral solution show that the mascon and PSEM models reproduce the reference potential with small relative errors while reducing evaluation costs. GravDyn is intended for mission design and studies of spacecraft motion, orbital stability, and gravitational modeling around asteroids and other small bodies.