2.2.2 相平衡压缩因子计算

由于PR方程、SRK方程系数较多，本书采用C语言编程的方法计算压缩因子。由PR压缩因子方程求导，得：

　（2-62）

由SRK压缩因子方程求导，得：

　（2-63）

通过求解式（2-49）中f（Z）的导数的根，求压缩因子、判断混合物的状态。

方程f（Z）为开口向上的抛物线方程。

　（2-64）

　（2-65）

为f·Δ（Z）求根判据，如果Δ<0，则方程·f（Z）无实根，f（Z）在实数区间上为增函数，·f（Z）恒大于零，f仅有一个根，用牛顿迭代法求该根。如果Δ=0，证明三次方程f（Z）仅存在一个根，该根必然大于零，小于某个较大值。如果Δ>0，方程·f（Z）有两个实根，则方程f（Z）必有三个实根，混合物为气相，则最大实根为方程的压缩因子；混合物为液相，则最小的实根为压缩因子，中间一解为无用解。

在本章的流程中，温度计算范围为-180～200℃，混合物最高压力为5.0MPa，最低压力为200kPa，因此设定该上限值为12MPa，利用分区间依据迭代法求解各根（图2-4）。

图2-4　压缩因子的求解程序框图