astrodyn 0.1.1

Pipeline orchestration, VehicleBuilder, and recipes — single API surface for ECS adapters
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67

//! Mars-orbit scenario.
//!
//! ```
//! use astrodyn::recipes::Mission;
//! let sb = Mission::mars_orbit().into_builder();
//! assert_eq!(sb.sources.len(), 2);
//! ```

use astrodyn_gravity::{GravityControl, GravityGradient};
use astrodyn_quantities::ext::Vec3Ext;
use glam::DVec3;

use crate::recipes::{epoch, mars, sun, vehicle};
use crate::vehicle_builder::VehicleBuilder;
use crate::SimulationBuilder;

/// Dawn-class spacecraft in a Mars orbit. Mars central body
/// (point-mass + IAU rotation), Sun as a third-body source. Step size
/// 10 s, Dawn-arrival epoch.
///
/// The point-mass scenario produces a substantially less accurate
/// trajectory than the JEOD reference (which uses the MRO110B2
/// 110×110 SH model). Mission code wanting verification-grade
/// accuracy substitutes the central-body source with
/// [`crate::recipes::mars::mro110b2`] and adds a non-spherical
/// [`GravityControl`].
///
/// ```
/// use astrodyn::recipes::scenarios::mars_orbit;
/// let sb = mars_orbit::mars_orbit();
/// // Mars central + Sun third-body.
/// assert_eq!(sb.sources.len(), 2);
/// assert_eq!(sb.dt, 10.0);
/// ```
pub fn mars_orbit() -> SimulationBuilder {
    let mut sb = SimulationBuilder::new(epoch::dawn_mars_2009(), 10.0);
    let mars_idx = sb.add_source("Mars", mars::point_mass());
    let sun_idx = sb.add_source(
        "Sun",
        sun::third_body(
            DVec3::new(2.27e11, 0.0, 0.0).m_at::<astrodyn_quantities::frame::RootInertial>(),
        ),
    );
    sb = sb.sun(sun_idx);

    // Dawn initial state at Mars (Mars-centered inertial frame).
    use astrodyn_dynamics::state::TranslationalStateTyped;
    use astrodyn_quantities::frame::RootInertial;
    let trans = TranslationalStateTyped::<RootInertial> {
        position: DVec3::new(11_563_355.680_2, -14_356_668.897_7, 6_293_704.616_9)
            .m_at::<RootInertial>(),
        velocity: DVec3::new(-2_273.107_8, 2_380.132_4, -22.911).m_per_s_at::<RootInertial>(),
    };

    let vehicle = VehicleBuilder::new()
        .with_translational(trans)
        .three_dof_point_mass(vehicle::dawn_mass())
        .rk4()
        .gravity(GravityControl::new_spherical(
            mars_idx,
            GravityGradient::Skip,
        ))
        .gravity(GravityControl::new_third_body(sun_idx))
        .build();
    sb.add_body(vehicle);
    sb
}