lib.EoS.Lee_Kesler module

class lib.EoS.Lee_Kesler.Lee_Kesler(T, P, mezcla)[source]

Bases: EoS

Corresponding state equation of state of Lee-Kesler

\[\begin{split}\begin{array}[t]{l} Z = Z^{(0)} + \omega Z^{(1)}\\ Z = \frac{P_rV_r}{T_r} = 1 + \frac{B}{V_r} + \frac{C}{V_r^2} + \frac{D}{V_r^5} + \frac{c_4}{T_r^3V_r^2}\left(\beta+\frac{\gamma} {V_r^2}\right)\exp{\left(-\frac{\gamma}{V_r^2}\right)}\\ B = b_1 - \frac{b_2}{T_r} - \frac{b_3}{T_r^2}-\frac{b_4}{T_r^3}\\ C = c_1 - \frac{c_2}{T_r} - \frac{c_3}{T_r^3}\\ D = d_1 - \frac{d_2}{T_r}\\ \end{array}\end{split}\]

Using the mixing rules defined by Plöcker [2]

\[\begin{split}\begin{array}[t]{l} T_{CM} = \frac{1}{v_{CM}^\nu \sum_j \sum_k Z_j Z_k v_{Cjk}^\nu T_{Cjk}\\ v_{CM} = \sum_j \sum_k Z_j Z_k v_{Cjk}\\ \omega_M = \sum_j Z_j \omega_j\\ P_{CM} = \left(0.2905-0.085\omega_m\right) R \frac{T_{CM}}{v_{CM}\\ \end{array}\end{split}\]

with the critical cross parameters

\[\begin{split}\begin{array}[t]{l} T_{Cjk} = \left(T_{Cj} T_{Ck}\right)^{1/2) k_{jk}\\ v_{Cjk} = \frac{1}{8} \left(v_{Cj}^{1/3} + v_{Ck}^{1/3}\right)^3\\ \end{array}\end{split}\]

Examples

Example 1.17 from [3], Propane compressibility

>>> from lib.mezcla import Mezcla
>>> mix = Mezcla(1, ids=[8], caudalUnitarioMasico=[1.])
>>> eq = Lee_Kesler(0.8*mix.Tc, 0.4*mix.Pc, mix)
>>> '%0.4f' % (eq.Z[0])
'0.0592'
>>> eq = Lee_Kesler(0.9*mix.Tc, 0.4*mix.Pc, mix)
>>> '%0.4f' % (eq.Z[1])
'0.7520'
>>> eq = Lee_Kesler(1*mix.Tc, 0.4*mix.Pc, mix)
>>> '%0.4f' % (eq.Z[1])
'0.8437'
>>> eq = Lee_Kesler(1.1*mix.Tc, 0.4*mix.Pc, mix)
>>> '%0.4f' % (eq.Z[1])
'0.8939'
>>> eq = Lee_Kesler(1.2*mix.Tc, 0.4*mix.Pc, mix)
>>> '%0.4f' % (eq.Z[1])
'0.9253'
__init__(T, P, mezcla)[source]
_mix(zi)[source]

Mixing rules, Eq 20-25

_lib(zi=None, T=None, P=None, rho0=None)[source]
_Z(zi=None, T=None, P=None, rho0=None)[source]
_fug(xi, yi, T, P)[source]

Fugacities of component in mixture calculation

Parameters:
xilist

Molar fraction of component in liquid phase, [-]

yilist

Molar fraction of component in vapor phase, [-]

Returns:
titallist

List with liquid phase component fugacities

titavlist

List with vapour phase component fugacities

_fugacity(ref, Tr, Pr, vr)[source]

Calculate fugacity coefficient of mixture

Parameters:
refinteger

Index of parameter, 0 for simple fluid and 1 for reference fluid

Trfloat

Reduced temperature

Prfloat

Reduced pressure

vrfloat

Reduced volume

_Hexc(zi, T, P)[source]
_Hexc_lib(ref, Tr, Pr, vr)[source]

Calculate enthalpy excess of mixture

Parameters:
refinteger

Index of parameter, 0 for simple fluid and 1 for reference fluid

Trfloat

Reduced temperature

Prfloat

Reduced pressure

vrfloat

Reduced volume

lib.EoS.Lee_Kesler.Lee_Kesler_lib(Tr, Pr, fase=2, Vvo=0.0001, Vlo=5)[source]

Librería para el cálculo de la EoS de Lee-Kesler Procedure API 6B1.8 pag 518 Perry pag 2-358 fase: fase para la que se realiza el cálculo

0 - Liquido 1 - Vapor 2 - Ambas

lib.EoS.Lee_Kesler.Z_Lee_Kesler(T, P, mezcla)[source]

Factor de compresibilidad según la ecuación de estado de Lee-Kesler

lib.EoS.Lee_Kesler.Lee_Kesler_lib_Cp(Tr, Pr, fase=1)[source]

Librería para el cálculo de capacidades calorificas, usada a continuación en diferentes funciones Procedure API 7E1.6 Pag.726

lib.EoS.Lee_Kesler.Lee_Kesler_Entalpia_lib(Tr, Pr, w, fase=1)[source]

Librería para el cálculo del factor adimensional de influencia de la presión sobre la temperatura. Usado en diversos métodos a continuación eq 7B3.7-1 pag 643

lib.EoS.Lee_Kesler.Lee_Kesler_Fugacidad_lib(Tr, Pr, w, fase=1)[source]

Librería para el cálculo de la fugacidad, entropia… Procedure API 7G1.8 Pag.752

lib.EoS.Lee_Kesler.Entalpia_Lee_Kesler(self)[source]

Librería para el cálculo del factor adimensional de influencia de la presión sobre la temperatura. Usado en diversos métodos a continuación eq 7B3.7-1 pag 643

References