lib.PPR78 module¶
- class lib.EoS.Cubic.PPR78.PPR78(T, P, mezcla, **kwargs)[source]¶
Bases:
PR78Predictive Peng-Robinson cubic equation of state
\[\begin{split}\begin{array}[t]{l} P = \frac{RT}{V-b}-\frac{a}{V\left(V+b\right)+b\left(V-b\right)}\\ a = 0.45747\frac{R^2T_c^2}{P_c}\alpha\\ b = 0.0778\frac{RT_c}{P_c}\\ \alpha^{0.5} = 1 + m\left(1-Tr^{0.5}\right)\\ m = 0.37464 + 1.54226\omega-0.26992\omega^2 if \omega < 0.491\\ m = 0.379642 + 1.48503\omega - 0.164423*\omega^2 + 0.016666*\omega^3\\ k_{ij} = \frac{-\frac{1}{2}\sum_{k=1}^{Ng} \sum_{l=1}^{Ng} \left(\alpha_{ik}-\alpha_{jk}\right)\left(\alpha_{il}-\alpha_{jl}\right) A_{kl} \left(\frac{298.15}{T}\right)^{\left(\frac{B_{kl}}{A_{kl}}-1\right)} -\left(\frac{\sqrt{a_i}}{b_i}-\frac{\sqrt{a_j}}{b_j}\right)^2} {2\frac{\sqrt{a_i a_j}}{b_i b_j}} \end{array}\end{split}\]Methods
Convert UNIFAC group saved in compound database to group contribution as used in PPR78
Examples
kij between propane and n-butane at T=303.15 from Appendix A of [1], tiny variation because differences in critical properties of components
>>> mix = Mezcla(5, ids=[4, 6], caudalMolar=1, fraccionMolar=[0.5, 0.5]) >>> eq = PPR78(303.15, 101325, mix) >>> print("%0.4f" % self.kij[0][1]) 0.0029