己烷质量分数为0.02 进料温度126 压力为120kpa 临界温度234.25 临界压力3.012 用PR方程计算pvt

好的，根据PR方程，可以得到以下计算公式：

a = 0.45724 * R^2 * Tc^2 / Pc
b = 0.07780 * R * Tc / Pc
m = 0.37464 + 1.54226 * w - 0.26992 * w^2
Tr = T / Tc
alpha = [1 + m * (1 - sqrt(Tr))]^2
A = alpha * a * P / (R^2 * T^2)
B = b * P / (R * T)
Z^3 - (1 - B) * Z^2 + (A - 2B - 3B^2) * Z - (AB - B^2 - B^3) = 0

其中，a，b分别为常数，R为气体常数，w为己烷偏离度，Tc，Pc为己烷的临界温度和临界压力，T，P分别为进料温度和压力，Tr为还原温度，A，B为PR方程中的两个参数，Z为己烷的压缩因子。

将己烷的参数代入公式中，可以得到：

a = 0.1391 kPa * m^3 / mol^2
b = 2.578E-4 m^3 / mol
m = 0.446
R = 8.314 J / mol / K
w = 0.0115
Tc = 234.25 K
Pc = 3.012 MPa
T = 126 ℃ = 399.15 K
P = 120 kPa = 1.2 bar

代入计算后，得到：

Tr = T / Tc = 1.7067
alpha = [1 + m * (1 - sqrt(Tr))]^2 = 0.9804
A = alpha * a * P / (R^2 * T^2) = 0.00374
B = b * P / (R * T) = 2.026E-5

然后，需要使用数值方法（例如牛顿迭代法）求解PR方程，得到Z的值。最后，可以使用以下公式计算密度、粘度和体积系数：

ρ = Z * M / (V * R * T)
μ = 1E-6 * 3.03E-5 * exp(520 / T) * sqrt(M) / Z
β = 1 / V

其中，M为己烷的摩尔质量，V为摩尔体积。

这里的密度单位为kg/m³，粘度单位为Pa·s，体积系数为1/K。