feos-core 0.10.1

Core traits and functionalities for the `feos` project.
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

use super::PropertyAD;
use crate::DensityInitialization::Liquid;
use crate::density_iteration::density_iteration;
use crate::{FeosResult, Residual};
use nalgebra::allocator::Allocator;
use nalgebra::{DefaultAllocator, Dim};
use num_dual::{DualNum, DualStruct};
use quantity::{_Density, Density, KELVIN, PASCAL, Pressure, Temperature};

/// Liquid density of a pure component as function of temperature and pressure.
pub struct LiquidDensity(pub Temperature, pub Pressure);

impl<'a> From<&'a [f64]> for LiquidDensity {
    fn from(value: &'a [f64]) -> Self {
        Self(value[0] * KELVIN, value[1] * PASCAL)
    }
}

impl<N: Dim> PropertyAD<N> for LiquidDensity
where
    DefaultAllocator: Allocator<N>,
{
    type Unit = _Density;
    const REFERENCE: Density = Density::new(1000.0);

    fn evaluate<E: Residual<N, D>, D: DualNum<f64, Inner = f64> + Copy>(
        &self,
        eos: &E,
    ) -> FeosResult<Density<D>> {
        let x = E::pure_molefracs();
        let t = Temperature::from_inner(&self.0);
        let p = Pressure::from_inner(&self.1);
        density_iteration(eos, t, p, &x, Some(Liquid))
    }
}