ioss 0.0.3

Io 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
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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170

use sciforge::hub::prelude::constants::elements::atomic_mass;

pub struct ExosphereSpecies {
    pub name: &'static str,
    pub symbol: &'static str,
    pub molar_mass_kg_mol: f64,
    pub column_density_m2: f64,
    pub scale_height_m: f64,
}

pub struct VolcanicPlumeEndpoint {
    pub plume_height_m: f64,
    pub ejection_velocity_m_s: f64,
    pub so2_mass_rate_kg_s: f64,
    pub temperature_k: f64,
    pub optical_depth_at_vent: f64,
}

pub struct ExosphereEndpoint {
    pub body_radius_m: f64,
    pub species: Vec<ExosphereSpecies>,
}

fn so2_molar() -> f64 {
    (atomic_mass(16) + 2.0 * atomic_mass(8)) * 1e-3
}
fn so_molar() -> f64 {
    (atomic_mass(16) + atomic_mass(8)) * 1e-3
}
fn s_molar() -> f64 {
    atomic_mass(16) * 1e-3
}
fn o_molar() -> f64 {
    atomic_mass(8) * 1e-3
}
fn na_molar() -> f64 {
    atomic_mass(11) * 1e-3
}
fn k_molar() -> f64 {
    atomic_mass(19) * 1e-3
}
fn cl_molar() -> f64 {
    atomic_mass(17) * 1e-3
}
fn nacl_molar() -> f64 {
    (atomic_mass(11) + atomic_mass(17)) * 1e-3
}

impl ExosphereEndpoint {
    pub fn io() -> Self {
        Self {
            body_radius_m: crate::IO_RADIUS_M,
            species: vec![
                ExosphereSpecies {
                    name: "Sulfur dioxide",
                    symbol: "SO2",
                    molar_mass_kg_mol: so2_molar(),
                    column_density_m2: 1e20,
                    scale_height_m: 12_000.0,
                },
                ExosphereSpecies {
                    name: "Sulfur monoxide",
                    symbol: "SO",
                    molar_mass_kg_mol: so_molar(),
                    column_density_m2: 1e18,
                    scale_height_m: 15_000.0,
                },
                ExosphereSpecies {
                    name: "Atomic sulfur",
                    symbol: "S",
                    molar_mass_kg_mol: s_molar(),
                    column_density_m2: 5e17,
                    scale_height_m: 40_000.0,
                },
                ExosphereSpecies {
                    name: "Atomic oxygen",
                    symbol: "O",
                    molar_mass_kg_mol: o_molar(),
                    column_density_m2: 1e18,
                    scale_height_m: 50_000.0,
                },
                ExosphereSpecies {
                    name: "Sodium",
                    symbol: "Na",
                    molar_mass_kg_mol: na_molar(),
                    column_density_m2: 1e14,
                    scale_height_m: 100_000.0,
                },
                ExosphereSpecies {
                    name: "Potassium",
                    symbol: "K",
                    molar_mass_kg_mol: k_molar(),
                    column_density_m2: 1e12,
                    scale_height_m: 80_000.0,
                },
                ExosphereSpecies {
                    name: "Chlorine",
                    symbol: "Cl",
                    molar_mass_kg_mol: cl_molar(),
                    column_density_m2: 5e13,
                    scale_height_m: 30_000.0,
                },
                ExosphereSpecies {
                    name: "Sodium chloride",
                    symbol: "NaCl",
                    molar_mass_kg_mol: nacl_molar(),
                    column_density_m2: 1e13,
                    scale_height_m: 10_000.0,
                },
            ],
        }
    }

    pub fn density_at_altitude(&self, species_symbol: &str, altitude_m: f64) -> f64 {
        self.species
            .iter()
            .find(|s| s.symbol == species_symbol)
            .map(|s| {
                let n_surface = s.column_density_m2 / s.scale_height_m;
                n_surface * (-altitude_m / s.scale_height_m).exp()
            })
            .unwrap_or(0.0)
    }

    pub fn total_column_density(&self) -> f64 {
        self.species.iter().map(|s| s.column_density_m2).sum()
    }

    pub fn sodium_glow_intensity(&self) -> f64 {
        self.species
            .iter()
            .find(|s| s.symbol == "Na")
            .map(|s| s.column_density_m2 * 1e-14)
            .unwrap_or(0.0)
    }
}

impl VolcanicPlumeEndpoint {
    pub fn pele_type() -> Self {
        Self {
            plume_height_m: 300_000.0,
            ejection_velocity_m_s: 1_000.0,
            so2_mass_rate_kg_s: 1e4,
            temperature_k: 1_500.0,
            optical_depth_at_vent: 5.0,
        }
    }

    pub fn prometheus_type() -> Self {
        Self {
            plume_height_m: 75_000.0,
            ejection_velocity_m_s: 500.0,
            so2_mass_rate_kg_s: 1e2,
            temperature_k: 450.0,
            optical_depth_at_vent: 1.5,
        }
    }

    pub fn plume_density_at_height(&self, height_m: f64) -> f64 {
        if height_m > self.plume_height_m {
            return 0.0;
        }
        let frac = height_m / self.plume_height_m;
        self.optical_depth_at_vent * (1.0 - frac).powi(2)
    }

    pub fn thermal_emission_w_m2(&self) -> f64 {
        5.670374419e-8 * self.temperature_k.powi(4)
    }
}