化工热力学第三版课件第二章3专用课件.pptVIP

聊城大学 化学化工学院 化工热力学 第二章 流体的热力学性质 第三节 §2-3 对应态原理的应用 立方形状态方程的普适化形式 两参数普遍化压缩因子图 三参数普遍化压缩因子图 普遍化第二维里系数关联式 一、立方形状态方程的普适形式 + ?SRK(Tr,?)={1+(0.480+1.574?-0.176?2)(1-Tr1/2)}2 ++ ?PR(Tr,?)={1+(0.37464+1.54226?-0.26992?2)(1-Tr1/2)}2 0.30740 0.45724 0.07779 1-√2 1+√2 ?PR(Tr,?)++ PR 1/3 0.42748 0.08664 0 1 ?SRK(Tr,?)+ SRK 1/3 0.42748 0.08664 0 1 Tr-1/2 RK 3/8 27/64 1/8 0 0 1 Vdw ZC ? ? ? ? ?(Tr) 状态方程 三次型普适方程的数值求解 1. 气相体积或压缩因子的求解 乘以(v-b)/p 迭代计算，用 作为初值 (a-1) 将式V=ZRT/P代入式(a-1)中， 再乘以P/RT， 令 ∴ 又∵ ∴ (a-2) 用Z=1为初值进行迭代计算 对比温度、对比压力 1. 液相体积或压缩因子的求解 由式(a-1)中最右边项的分子中，解出V，表示成下式： (b-1) 迭代计算，用 作为初值 由式(a-2)中最右边项的分子中，解出Z，表示成下式： (b-2) 用 为初值进行迭代计算 用迭代计算： 查 TC、 PC 求 Tr、 Pr 求q、β 由式右边求 Z1 No Yes 气相：Z0=1 液相：Z0=β 例题：异丁烷在360K时的蒸汽压是1.541MPa，试用立方型状态方程计算：（1）饱和汽相的摩尔体积；（2）饱和液相的摩尔体积。已知饱和汽相与液相的摩尔体积实验值分别为1393cm3.mol-1和128.0cm3.mol-1 从附表 1 查得异丁烷的 [解] ∴ （b）用 VDW 方程 σ=0 ， ε=0 ，Ω=1/8 ， ψ=27/64 ， α(Tr)=1 ∴ 对于汽相 将ε、 σ、q、β的值代入得： 利用 起始值进行迭代， 直到求得收敛结果为Z=0.79297 迭代示例 用Excel迭代 对于液相 将ε、 σ、q、β的值代入得： 利用 起始值进行迭代， 直到求得收敛结果为Z=0.10970 用Excel迭代 （b）用 R-K 方程 σ=1 ， ε=0 ，Ω=0.08664 ， ψ=0.42748 ， α(Tr)=Tr-1/2 ∴ 对于汽相 将ε、 σ、q、β的值代入得： 利用 起始值进行迭代， 直到求得收敛结果为Z=0.74495 迭代示例 用Excel迭代 对于液相 将ε、 σ、q、β的值代入得： 利用 起始值进行迭代， 直到求得收敛结果为Z=0.07716 用Excel迭代 （c）用 SRK 方程 σ=1 ， ε=0 ，Ω=0.08664 ， ψ=0.42748 ? (Tr)={1+(0.480+1.574?-0.176?2)(1-Tr1/2)}2 利用 起始值进行迭代， 直到求得收敛结果为Z=0.73225 对于汽相 将ε、 σ、q、β的值代入得： 用Excel迭代 对于液相 将ε、 σ、q、β的值代入得： 利用 起始值进行迭代， 直到求得收敛结果为Z=0.074034 用Excel迭代 （d）用 PR 方程 σ= ， ε= ，Ω= ， ψ= ? (Tr)={1+(0.37464+1.54226?-0.26992?2)(1-Tr1/2)}2 对于汽相 将ε、 σ、q、β的值代入得： 利用 起始值进行迭代， 直到求得收敛结果为Z=0.71509 用Excel迭代 对于液相 将ε、 σ、q、β的值代入得： 利用 起始值进行迭代， 直到求得收敛结果为Z=0.065238 用Excel迭代 项 目 实验值 vdW RK SRK PR 饱和汽相摩尔体积Vv/cm3mol-1 1393 1540 1446 1423 1389 相对偏差/% 0.00 10.55 3.81 2.15 -0.28 饱和液相摩尔体积Vl/ cm3mol-1 128.0 213.1 149.5 144.8 126.7 相对偏差/% 0.00 66.48 16.80 13.13 -1.02 为了比较计算结果，在例2-4的情况下，运用四种立方型状态方程所计算得到的V v及V l 值列表如下： 液体的两个常见物性 定义：