astrodyn 0.2.0

Gateway to the astrodyn orbital-dynamics framework — a pure-Rust, engine-agnostic port of NASA JEOD (spherical-harmonics gravity, RNP Earth rotation, atmosphere, drag/SRP, multi-body dynamics) composing the astrodyn_* physics crates into one pipeline API any host can drive
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42

//! Geostationary orbit scenario.
//!
//! ```
//! use astrodyn::recipes::Mission;
//! let sb = Mission::geo().into_builder();
//! assert_eq!(sb.bodies.len(), 1);
//! ```

use astrodyn_gravity::{GravityControl, GravityGradient};

use crate::recipes::{constants, earth, epoch, orbital_elements, vehicle};
use crate::vehicle_builder::VehicleBuilder;
use crate::SimulationBuilder;
use astrodyn_quantities::frame::{Earth, PlanetInertial};
use astrodyn_quantities::frame_descriptor::FrameUid;

/// Geostationary point-mass satellite at 42164 km, 0° inclination.
/// Earth point-mass gravity, RK4, J2000 epoch, 300 s step.
///
/// ```
/// use astrodyn::recipes::scenarios::geostationary;
/// let sb = geostationary::geo();
/// assert_eq!(sb.sources.len(), 1);
/// assert_eq!(sb.bodies.len(), 1);
/// assert_eq!(sb.dt, 300.0);
/// ```
pub fn geo() -> SimulationBuilder {
    let mut sb = SimulationBuilder::new(epoch::j2000(), 300.0);
    let _ = sb.add_source("Earth", earth::point_mass());
    let vehicle = VehicleBuilder::new()
        .vehicle_named("geo-sat")
        .from_orbital_elements(orbital_elements::geostationary(), constants::mu_ggm05c())
        .three_dof_point_mass(vehicle::unit_sphere_mass())
        .rk4()
        .gravity(GravityControl::new_spherical(
            FrameUid::of::<PlanetInertial<Earth>>(),
            GravityGradient::Skip,
        ))
        .build();
    sb.add_body(vehicle);
    sb
}