polykin.properties.viscosity

MUV_Lucas

MUV_Lucas(
    T: float,
    P: float,
    M: float,
    Tc: float,
    Pc: float,
    Zc: float,
    dm: float,
) -> float

Calculate the viscosity of a pure gas at a given temperature and pressure using the method of Lucas.

References

• RC Reid, JM Prausniz, and BE Poling. The properties of gases & liquids 4th edition, 1986, p. 421.

PARAMETER	DESCRIPTION
T	Temperature. Unit = K. TYPE: float
P	Pressure. Unit = Pa. TYPE: float
M	Molar mass. Unit = kg/mol. TYPE: float
Tc	Critical temperature. Unit = K. TYPE: float
Pc	Critical pressure. Unit = Pa. TYPE: float
Zc	Critical compressibility factor. TYPE: float
dm	Dipole moment. Unit = debye. TYPE: float

RETURNS	DESCRIPTION
float	Gas viscosity, \(\mu\). Unit = Pa·s.

Examples:

Estimate the viscosity of ethylene at 350 K and 10 bar.

>>> from polykin.properties.viscosity import MUV_Lucas
>>> mu = MUV_Lucas(T=350., P=10e5, M=28.05e-3,
...                Tc=282.4, Pc=50.4e5, Zc=0.280, dm=0.)
>>> print(f"{mu:.2e} Pa·s")
1.20e-05 Pa·s

Source code in src/polykin/properties/viscosity/vapor.py

def MUV_Lucas(T: float,
              P: float,
              M: float,
              Tc: float,
              Pc: float,
              Zc: float,
              dm: float
              ) -> float:
    r"""Calculate the viscosity of a pure gas at a given temperature and
    pressure using the method of Lucas.

    **References**

    *   RC Reid, JM Prausniz, and BE Poling. The properties of gases & liquids
        4th edition, 1986, p. 421.

    Parameters
    ----------
    T : float
        Temperature. Unit = K.
    P : float
        Pressure. Unit = Pa.
    M : float
        Molar mass. Unit = kg/mol.
    Tc : float
        Critical temperature. Unit = K.
    Pc : float
        Critical pressure. Unit = Pa.
    Zc : float
        Critical compressibility factor.
    dm : float
        Dipole moment. Unit = debye.

    Returns
    -------
    float
        Gas viscosity, $\mu$. Unit = Pa·s.

    Examples
    --------
    Estimate the viscosity of ethylene at 350 K and 10 bar.

    >>> from polykin.properties.viscosity import MUV_Lucas
    >>> mu = MUV_Lucas(T=350., P=10e5, M=28.05e-3,
    ...                Tc=282.4, Pc=50.4e5, Zc=0.280, dm=0.)
    >>> print(f"{mu:.2e} Pa·s")
    1.20e-05 Pa·s
    """
    Tr = T/Tc
    Pr = P/Pc
    FP0 = _MUV_Lucas_FP0(Tr, Tc, Pc, Zc, dm)
    mu = _MUV_Lucas_mu(Tr, Pr, M, Tc, Pc, FP0)  # type: ignore
    return mu