cascade.dynamics.mascon_asteroid

cascade.dynamics.mascon_asteroid(Gconst: float, points, masses, omega) → List[Tuple[expression, expression]]

Dynamics around a rotating asteroid.

This function generates the dynamical equations for objects orbiting around an irregular rotating body. The equations are written in the body fixed rotating frame.

Specifically, the returned system of differential equations will be:

\[\begin{split}\begin{array}{ll} \left\{ \begin{array}{l} \dot{\mathbf r} = \mathbf v \\ \dot {\mathbf v} = -G \sum_{j=0}^N \frac {m_j}{|\mathbf r - \mathbf r_j|^3} (\mathbf r - \mathbf r_j) - 2 \boldsymbol\omega \times \mathbf v - \boldsymbol \omega \times\boldsymbol\omega \times \mathbf r \end{array} \right. \end{array}\end{split}\]

Parameters:

• Gconst (float) – Cavendish constant

• points (array (N,3)) – Cartesian coordinates of the mascons

• masses (array (N,)) – mass of each mascon

• omega (array (3,)) – asteroid angular velocity

Returns:

The dynamics in the same units used by the input arguments. Can be used directly to instantiate a sim.