enceladuss 0.0.3

Enceladus celestial simulation crate for the MilkyWay SolarSystem workspace
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

#[derive(Debug, Clone, Copy, PartialEq)]
pub struct IceRegolith {
    pub thickness_m: f64,
    pub density_kg_m3: f64,
    pub ice_fraction: f64,
}

impl IceRegolith {
    pub fn porosity(&self) -> f64 {
        (1.0 - self.density_kg_m3 / 917.0).clamp(0.0, 1.0)
    }

    pub fn thermal_inertia(&self) -> f64 {
        self.density_kg_m3 * (150.0 + 400.0 * self.ice_fraction)
    }
}

pub fn south_polar_regolith() -> IceRegolith {
    IceRegolith {
        thickness_m: 5.0,
        density_kg_m3: 500.0,
        ice_fraction: 0.95,
    }
}

pub fn cratered_north_regolith() -> IceRegolith {
    IceRegolith {
        thickness_m: 12.0,
        density_kg_m3: 600.0,
        ice_fraction: 0.90,
    }
}