Struct EquationOfState 

pub struct EquationOfState {
    pub r_gas: f64,
    pub temperature: f64,
}

Classical cubic equations of state for real gases.

Provides van der Waals, Peng-Robinson, and Redlich-Kwong equations.

Fields

r_gas: f64

Universal gas constant R.

temperature: f64

Temperature T.

Implementations

impl EquationOfState

pub fn new(r_gas: f64, temperature: f64) -> Self

Creates a new EOS calculator.

r_gas is usually 8.314 J/(mol·K).

pub fn van_der_waals_pressure(&self, molar_volume: f64, a: f64, b: f64) -> f64

Van der Waals pressure: P = RT/(V-b) - a/V².

Arguments

• molar_volume - Molar volume V (m³/mol).
• a - Attractive parameter a (Pa·m⁶/mol²).
• b - Repulsive parameter b (m³/mol).

pub fn vdw_critical_constants(a: f64, b: f64, r_gas: f64) -> (f64, f64, f64)

Van der Waals critical constants from parameters a and b.

Returns (T_c, P_c, V_c).

pub fn peng_robinson_pressure( &self, molar_volume: f64, a_c: f64, b: f64, kappa: f64, tc: f64, ) -> f64

Peng-Robinson pressure.

P = RT/(V-b) - a(T) / (V(V+b) + b(V-b))

where a(T) = a_c · α(T), α(T) = [1 + κ(1 - √(T/T_c))]².

pub fn redlich_kwong_pressure(&self, molar_volume: f64, a: f64, b: f64) -> f64

Redlich-Kwong pressure.

P = RT/(V-b) - a / (T^0.5 V(V+b))

pub fn compressibility_factor(&self, pressure: f64, molar_volume: f64) -> f64

Compressibility factor Z = PV/(RT).

pub fn boyle_temperature_vdw(a: f64, b: f64, r_gas: f64) -> f64

Boyle temperature for van der Waals gas: T_B = a / (Rb).

Trait Implementations

impl Clone for EquationOfState

fn clone(&self) -> EquationOfState

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for EquationOfState

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.