4第2章PVT-3.pptVIP

2.6 混合法则 2.6.1 virial 方程的混合法则 2.6.3 BWR方程 2.6.1 virial 方程的混合法则 思考：三元物系和四元物系如何展开？ 2.6.1 virial 方程的混合法则 交叉第二Virial系数计算式: 2.6.1 virial 方程的混合法则 交叉第二Virial系数计算式: B0和B1可用对比温度Tr = T/Tcij计算见式(2-27); kij-相互作用参数，可由实验数据拟合； 一般kij=0-0.2，是温度的函数,Kij=a0+a1T,近似计算可取 Kij = 0 。 2.6.1 virial 方程的混合法则 (2-27) 2.6.2 立方型方程 2.6.2 立方型方程 二次型、LC、DDLC、UNIWALLS 26篇文献 PHV、MHV2、WS、SR型 PR方程的混合规则 混合物虚拟临界温度和真实临界温度是不同的!!! 其中r的值见表2-4。 3 2 1 2 r a C0 B0 A0 χ 2 3 3 3 γ α c b 表2-4 * * 重点回顾 2.5 多常数状态方程 Virial (维里)方程 Virial 系数获取 统计法 P-V-T数据拟合 普遍化关联(普维法) 普遍化第二维里系数关系式 (2-27) 2.5 多常数状态方程 重点回顾 对比第二维里系数 普遍化关联式的适用范围 普遍化维里系数法使用区 普遍化压缩因子法使用区 Pr Tr Vr≥ 2 B法 Vr2 Z法 4 3 2 1 0 1 2 3 4 5 6 7 8 普遍化维里系数法使用区 普遍化压缩因子法使用区 Vr≥ 2 B法 Vr2 Z法 0 1 2 3 4 5 6 7 8 Vr≥2 B法 2.5 多常数状态方程 重点回顾 两相区 Benedict-Webb-Rubin(BWR)方程（1940 ） 2.5 多常数状态方程 重点回顾 八参数方程; 适用范围广、精度高(极性、高压、VLE) 具有参数的有30多种，主要是烃（P-V-T数据和蒸气压数据关联，或用普遍化关系估算) BWRS方程 2.5.3 Martin-Hou（MH）方程1955年 2.5 多常数状态方程 MH-55 9参数方程 可用于非极性和强极性气体的计算！ MH-81 10参数方程 可用于汽-液相平衡(VLE)和液液相平衡(LLE)的计算！ 临界参数和一点的蒸气压数据可计算所有方程常数！ 重点回顾 1.已知T,V，如何求P？ 显压型[P=P(T,V)]，可直接计算，很方便; 在计算时，一定要注意单位! 2.已知T,P，如何求V？ 　工程上最常用的情况，因为T, P易测; 解析法及各种数值法(迭代法)求解。 3.已知P, V，如何求T ？ 用各种迭代法求解。 工程计算的类型: 重点回顾 *2.7 状态方程体积根的求解 1.解析法 2.数值法 二分法 牛顿法(切线法、割线法） 一般迭代法 图解法 软件法 重点回顾 [例题2-2] 用RK方程计算异丁烷（a）在420K和2Mpa时的摩尔体积（实验值是1411.2cm3.mol-1）（b）在380K时的饱和汽、液相摩尔体积，已知该温度下的蒸汽压为2.25Mpa（实验值分别为866.1cm3.mol-1和140.8cm3.mol-1）. 审题:已知T,P求V 解:RK方程→a,b 参数→临界参数(TC,PC,ω) TC=408.1K,PC=3.648MPa,ω=0.176 P267附录A-1,注意单位! C 方法1——直接迭代法1 方法2——直接迭代法2(RK方程) 方法3——牛顿迭代法 方法4--Excel电子表格计算 单变量求解体积根 直接迭代法 用Excel电子表格 直接迭代法计算 例2-2(a) V0=RT/P V=1404.54cm3/mol 实验值:1411.2cm3/mol 误差:0.5% 单变量求解 V=1404.68cm3/mol 单变量求解体积根 例题2-2 (b T=380K P=2.25MPa) 实验值:866.1 误差:5.8% 实验值:140.8 误差23.7% RK不适合计算液相! 方法5--图解法求体积根 b=80.58 Vb时,方程才有意义(参见P22图2-4) 例题2-2 图解法求体积根 V=14