lib.PRYuLu module

class lib.EoS.Cubic.PRYuLu.PRYuLu(T, P, mezcla, **kwargs)[source]

Bases: Cubic

Peng-Robinson equation of state modified by Yu-Lu

\[\begin{split}\begin{array}[t]{l} P = \frac{RT}{V-b}-\frac{a}{V\left(V+c\right)+b\left(3V+c\right)}\\ a = \Omega_{ac}\frac{R^2T_c^2}{P_c}\alpha\\ b = \Omega_{bc}\frac{RT_c}{P_c}\\ c = \Omega_{cc}\frac{RT_c}{P_c}\\ \Omega_{ac} = 0.468630-0.0379304\omega+0.00751969\omega^2\\ \Omega_{bc} = 0.0892828-0.0340903\omega-0.00518289\omega^2\\ \Omega_{cc} = w \Omega_{bc} c = wb\\ w+3 = u = 1.70083+0.648463\omega+0.895926\omega^2\\ \alpha = 10^{M\left(A_0+A_1T_r+A_2T_r"2\right)\left(1-T_r\right)}\\ \end{array}\end{split}\]

The alpha parameter depend of acentric factor. For ω ≤ 0.49:

\[\begin{split}\begin{array}[t]{l} M = 0.406846+1.87907\omega-0.792636\omega^2+0.737519\omega^3\\ A_0 = 0.536843\\ A_1 = -0.39244\\ A_2 = 0.26507\\ \end{array}\end{split}\]

For 0.49 < ω ≤ 1:

\[\begin{split}\begin{array}[t]{l} M = 0.581981-0.171416\omega-1.84441\omega^2+1.19047\omega^3\\ A_0 = 0.79355\\ A_1 = -0.53409\\ A_2 = 0.37273\\ \end{array}\end{split}\]
_cubicDefinition(T)[source]

Definition of individual components coefficients

_GEOS(xi)[source]

Definition of parameters of generalized cubic equation of state, each child class must define in this procedure the values of mixture a, b, delta, epsilon. The returned values are not dimensionless.

Parameters:
xilist

Molar fraction of component in mixture, [-]

Returns:
parameterslist

Mixture parameters of equation, a, b, c, d

_lib(cmp, T)[source]

References