lib.EoS.BWRSoave module

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

Bases: EoS

Soave modification of Benedict-Webb-Rubin equation of state as give in [1]

\[Z = 1 + B\rho + D\rho^4 + E\rho^2\left(1+F\rho^2\right) \exp(-F\rho^2)\]
\[\begin{split}\begin{array}[t]{l} \beta = B \left(\frac{P_c}{RT_c}\right) = \beta_c + 0.422\left(1-\frac{1}{T_r^{1.6}}\right) + 0.234 \omega \left(1-\frac{1}{T_r^3}\right) \\ \delta = D \left(\frac{P_c}{RT_c}\right)^4 = \delta_c \left[ 1 + d_1\left(\frac{1}{T_r}-1\right) + d_2\left(\frac{1}{T_r}-1\right)^2\right] \\ \epsilon = E \left(\frac{P_c}{RT_c}\right)^2 = \epsilon_c + e_1\left(\frac{1}{T_r}-1\right) + e_2\left(\frac{1}{T_r}-1\right)^2 + e_3\left(\frac{1}{T_r}-1\right)^3 \\ \phi = F \left(\frac{P_c}{RT_c}\right)^2 = f Z_c^2 \end{array}\end{split}\]

where:

\[\begin{split}\begin{array}[t]{l} d_1 = 0.4912 + 0.6478\omega\\ d_2 = 0.3 + 0.3619\omega\\ e_1 = 0.0841 + 0.1318\omega + 0.0018\omega^2\\ e_2 = 0.075 + 0.2408\omega + 0.014\omega^2\\ e_3 = -0.0065 + 0.1798\omega + 0.0078\omega^2\\ f = 0.77\\ \beta_c = bZ_c\\ \delta_c = dZ_c^4\\ \epsilon_c= eZ_c^2\\ e = \frac{2-5Z_c}{\left(1+f+3f^2-2f^3\right)\exp(-f)}\\ d = \frac{1}{3} \left[1-2Z_c-e\left(1+f-2f^2\right)\exp(-f)\right]\\ b = Z_c - 1 - d - e\left(1+f\right)\exp(-f) \end{array}\end{split}\]

This equation use critical parameters as parameters including Zc and acentric factor.

The mixture are treated as fictitious pure compounds with the critical parameters calculated by this mixing rules:

\[\begin{split}\begin{array}[t]{l} T_{cm} = \frac{S_1}{\left(\sqrt{S_2}+\sqrt{S_3}\right)^2} \\ P_{cm} = \frac{T_{cm}}{S_3} \\ m_m = \sqrt{\frac{S_2}{S_3}} \\ Z_{cm} = \frac{\sum_ix_iZ_{ci}T{_ci}/P_{ci}}{\sum_ix_i T_{ci}/P_{ci}}\\ S_1 = \sum_i \sum_j x_ix_j\left(1-k_{ij}\right)\frac{T_{ci}T_{cj}} {\sqrt{P_{ci}P_{cj}}} \left(1+m_i\right)\left(1+m_j\right) \\ S_2 = \sum_i \sum_j x_ix_j\left(1-k_{ij}\right)\frac{T_{ci}T_{cj}} {\sqrt{P_{ci}P_{cj}}} m_i m_j \\ S_3 = \sum_i x_i \frac{T_{ci}}{P_{ci}} \\ \end{array}\end{split}\]
__init__(T, P, mezcla)[source]
_lib(T, **kw)[source]

Library for compound specified parameters

_mixture(zi)[source]

Mixing rules

_Z(zi=None, T=None, P=None, crit=None, phase=0)[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(Z, x, **kw)[source]

References