oxiphysics-materials 0.1.0

Material properties and material library for the OxiPhysics engine
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

//! # MieGruneisenEos - Trait Implementations
//!
//! This module contains trait implementations for `MieGruneisenEos`.
//!
//! ## Implemented Traits
//!
//! - `EquationOfState`
//!
//! 🤖 Generated with [SplitRS](https://github.com/cool-japan/splitrs)

use super::functions::EquationOfState;
#[allow(unused_imports)]
use super::functions::*;
use super::types::MieGruneisenEos;

impl EquationOfState for MieGruneisenEos {
    fn pressure(&self, density: f64) -> f64 {
        self.pressure_from_energy(density, 0.0)
    }
    fn sound_speed(&self, density: f64) -> f64 {
        let mu = density / self.rho0 - 1.0;
        let denom = (1.0 - self.s * mu).powi(2);
        if denom.abs() < f64::EPSILON {
            return self.c0;
        }
        let ph = self.hugoniot_pressure(density);
        let dp_drho =
            self.c0 * self.c0 * (1.0 + 2.0 * self.s * mu) / (denom * denom) + self.gamma_0 * ph;
        dp_drho.max(0.0).sqrt()
    }
    fn density_from_pressure(&self, pressure: f64) -> f64 {
        let mut lo = self.rho0 * 0.5;
        let mut hi = self.rho0 * 10.0;
        for _ in 0..60 {
            let mid = 0.5 * (lo + hi);
            if self.pressure(mid) < pressure {
                lo = mid;
            } else {
                hi = mid;
            }
        }
        0.5 * (lo + hi)
    }
}