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
36

use super::PropertyAD;
use crate::DensityInitialization::Liquid;
use crate::{FeosResult, Residual, State};
use nalgebra::DefaultAllocator;
use nalgebra::allocator::Allocator;
use num_dual::{DualNum, DualStruct, Gradients};
use quantity::{_MolarEntropy, KELVIN, MolarEntropy, PASCAL, Pressure, Temperature};

/// Liquid residual isobaric heat capacity of a pure component as function of temperature
/// and pressure.
pub struct ResidualIsobaricHeatCapacity(pub Temperature, pub Pressure);

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

impl<N: Gradients> PropertyAD<N> for ResidualIsobaricHeatCapacity
where
    DefaultAllocator: Allocator<N>,
{
    type Unit = _MolarEntropy;
    const REFERENCE: MolarEntropy = MolarEntropy::new(1.0);

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