C2.2第二节状态方程选编.pptVIP

例:将1kmol甲烷压缩储存于容积为0.125m3,温度为323.16K的钢瓶内.问此时甲烷产生的压力多大?其实验值为1.875x107Pa. 解:1）理想气体 2） RK方程 查附录得TC=190.6K，PC=4.600MPa=4.6x106Pa R=8.314x103 m3.Pa/kmol.K 误差高达14.67%!!! 误差仅为1.216%!!! 例 异丁烷是取代氟利昂的环保制冷剂，用于冰箱、冷柜、冷饮机。现需要将1kmol、300K、0.3704MPa的异丁烷装入容器，请问需设计多大的容器？试用理想气体方程、RK、SRK和PR方程分别计算，并与实际值进行比较（实际值为V=6.081m3/kmol ） 解 从附录查得异丁烷的临界参数为 Tc＝408.1K Pc＝3.648MPa ω ＝0.176 ( 1 ) RK方程 实验值为V=6.081m3/kmol,误差为：0.97% ( 2 ) SRK方程 实验值为V=6.081m3/kmol,误差为：0.38% （3）PR方程： 相同方法可求得：V=6.0685m3/kmol 实验值为V=6.081m3/kmol,误差为：-0.16% （4）各种状态方程误差比较： EOS方程 误差 理想气体方程 10.74% RK方程 0.99% SRK方程 0.38% PR方程 -0.16% 注：实验测定的误差为0.5%左右，因此SRK、 PR方程的误差已小于实验误差，非常了不起！ 应用Excel和现成软件计算EOS 迭代法用手工计算来完成是相当繁琐的，可以自己编程序，或用诸如Mathcad或Maple的软件包求解，当然必须给定初值或求解范围。 当然最简单的是利用网上免费状态方程计算软件（http://www.cheng.cam.ac.uk/~pjb10/thermo/pure.html）来求解。 应用状态方程时要注意： (1)式中的V是摩尔体积! SI制单位：[m3/mol] (2) 式中的T是绝对温度, [K] (3)式中P是绝对压力，SI制单位：[pa] (4)式中的n为1 mol (5)通用气体常数R的单位必须 和P,V,T的单位相适应。 表2-1 通用气体常数R值(P.6) 注：R=8.314[m3.pa/mol.K] =8.314 [J/mol.K] 建议各物理量均换算为SI制单位，再代入方程式中进行计算。 Virial (维里)方程 Virial方程的形式[Onnes（昂尼斯）提出] 立方型状态方程属于半经验半理论的状态方程，而在状态方程中有严格理论基础的当属Virial方程。 Z为压缩因子 当 B、C…（或B’、C’…）称作Virial系数，是有物理意义的。 微观上，Virial系数反映了分子间的相互作用，第二维里系数B反映了两个分子之间的相互作用；第三维里系数C反映了三重分子的相互作用。 宏观上，Virial系数仅是温度的函数。 Virial 系数的获取 ( 1 ) 由统计力学进行理论计算：目前应用很少 ( 2 ) 由实验测定：精度较高 ( 3 ) 用普遍化关联式计算：方便，但精度不如实验测定的数据 2. Virial系数 二种形式的Virial方程是等价的，其系数之间也有相互关系。 3. Virial方程不同形式的关系 适用于TTc,P1.5MPa蒸汽 适用于TTc, 1.5MPa P 5MPa蒸汽 4.实际中常用Virial截断式 许多气体有B;但C较少;D更少，所以只能用维里截断式。 截取的项数越少，精度也就越低，所适用的压力越低。 1.两项维里截断式 2.三项维里截断式 将 代入两项维里截断式 最常用的Virial截断式： 得到 维里截断式的优缺点： 对于更高的压力维里截断式不适用，需用立方型状态方程。 只能计算气体，不能计算液相体积。 尽管如此，Virial方程的价值已超出PVT的应用,能描述气体的粘度、声速和热容。 它是其它多参数状态方程如B-W-R方程、M-H方程的基础。 最常用的Virial截断式 ?EOS 形式 Zc 适合范围 优缺点 理想气体 ? 1 压力极低的气体 不适合真实气体 vdW ? 0.375 同时能计算汽，液两相 准确度低 RK ? 0.333 计算气相体积准确性高，很实用 不能同时用于汽、液两相 SRK 同RK 0.333 能同时用于汽液两相平衡，广泛应用 精度高于RK，能预测液相体积? PR ? ?0.307 能同时用于汽液两相平衡，广泛应用 能预测液相体积? Virial ? TTc， P 5MPa的气相 不能同时用于汽液两相 状态方程小结 PV=RT 真实流体Zc =0.23~0.29 从§2.1可知