sciforge-core 0.0.4

Shared engineering primitives: materials, fasteners, fluids, propulsion taxonomies, std components.
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

use crate::moleculars::Gas;
use sciforge_hub::prelude::physics::thermodynamics as sf_thermo;

pub struct DieselCycle<'a> {
    pub gas: &'a Gas,
    pub compression_ratio: f64,
    pub cutoff_ratio: f64,
    pub t1_k: f64,
    pub p1_pa: f64,
}

impl<'a> DieselCycle<'a> {
    pub fn new(
        gas: &'a Gas,
        compression_ratio: f64,
        cutoff_ratio: f64,
        t1_k: f64,
        p1_pa: f64,
    ) -> Self {
        Self { gas, compression_ratio, cutoff_ratio, t1_k, p1_pa }
    }

    pub fn efficiency(&self) -> f64 {
        sf_thermo::diesel_efficiency(self.compression_ratio, self.cutoff_ratio, self.gas.gamma())
    }

    pub fn t2_k(&self) -> f64 {
        self.t1_k * self.compression_ratio.powf(self.gas.gamma() - 1.0)
    }

    pub fn t3_k(&self) -> f64 {
        self.t2_k() * self.cutoff_ratio
    }

    pub fn t4_k(&self) -> f64 {
        let gamma = self.gas.gamma();
        self.t3_k() / (self.compression_ratio / self.cutoff_ratio).powf(gamma - 1.0)
    }

    pub fn p2_pa(&self) -> f64 {
        self.p1_pa * self.compression_ratio.powf(self.gas.gamma())
    }

    pub fn q_in_j_kg(&self) -> f64 {
        self.gas.cp_j_kgk_ref * (self.t3_k() - self.t2_k())
    }

    pub fn heat_rejected_j_kg(&self) -> f64 {
        self.gas.cv_j_kgk_ref * (self.t4_k() - self.t1_k)
    }

    pub fn net_work_j_kg(&self) -> f64 {
        self.q_in_j_kg() * self.efficiency()
    }

    pub fn mep_pa(&self) -> f64 {
        let r = self.gas.specific_gas_constant_j_kgk();
        let v1 = r * self.t1_k / self.p1_pa;
        self.net_work_j_kg() / (v1 * (1.0 - 1.0 / self.compression_ratio))
    }
}