《流体的压力体积温度关系状态方程精选课件.ppt》由会员分享，可在线阅读，更多相关《流体的压力体积温度关系状态方程精选课件.ppt（101页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、关于流体的压力体积温度关系状态方程第一页，本课件共有101页2本章目的：本章目的：计算真实气体的计算真实气体的PVTPVT关系；关系；为后续内容（由为后续内容（由PVTPVT关系计算关系计算H H、S S，相平衡，化学相平衡，化学 平衡）打基础。平衡）打基础。本章重点：本章重点：偏心因子定义及其物理意义；偏心因子定义及其物理意义；EOSEOS的特点及使用范围。的特点及使用范围。基本热力学性质基本热力学性质P P、V V、T T、C CP P、CvCv等等热力学性质热力学性质 计算热力学性质计算热力学性质H H、S S、F F、G G、f f、E EX X、A AN N等等第二页，本课件共有10
2、1页3热力学基本性质有两大类热力学基本性质有两大类P,V,T,Cp,xv如何解决？如何解决？U,H,S,G但存在问题但存在问题：1)1)有限有限的的P-V-T数据，数据，无法无法全面全面了解流体的了解流体的P-V-T 行为；行为；2)离散离散的的P-V-T数据数据，不便于求导和积分，无法获得数据点，不便于求导和积分，无法获得数据点以外的以外的P-V-T 和和H，U，S，G数据数据。易测易测难或不能测难或不能测!从从容易容易获得的物性数据获得的物性数据（P、V、T、x）来推算来推算较难测定较难测定的数据的数据（H，U，S，G）怎么办怎么办?第三页，本课件共有101页4由基本热力学关系式应用全微分
3、的性质得由基本热力学关系式应用全微分的性质得Maxwell关系式关系式Maxwell关系式特点是将关系式特点是将难测的量难测的量用用易测的量代替。易测的量代替。如如 用用 代代;用用 代代 ;Maxwell方程方程 建立了建立了S=S(T，P)。第四页，本课件共有101页5 引言引言v如何解决？如何解决？只有建立能反映流体只有建立能反映流体P-V-T关系的解析形式关系的解析形式才能解决。才能解决。这就是状态方程这就是状态方程Equation of State(EOS)的由来。的由来。EOS反映了体系的特征，是推算实验数据之反映了体系的特征，是推算实验数据之外信息和其它物性数据不可缺少的模型。外
4、信息和其它物性数据不可缺少的模型。v流体流体P-V-T数据数据+状态方程状态方程EOS是计算热力学性是计算热力学性质最重要的模型之一。质最重要的模型之一。vEOS+CPig所有的热力学性质所有的热力学性质第五页，本课件共有101页62.1 纯物质的纯物质的PVT行为行为 2.1.1 纯流体的纯流体的P-V-T相图相图纯物质的纯物质的P-V-T立体相图立体相图 投影图投影图纯物质的纯物质的P-T图图 等容线等容线纯物质的纯物质的P-V图图 等温线等温线第六页，本课件共有101页7在常压下加热水在常压下加热水（1atm=0.10133MPa）带有活塞的汽缸保带有活塞的汽缸保持恒压持恒压液体水液体水
5、第七页，本课件共有101页8Tv12534第八页，本课件共有101页9第九页，本课件共有101页10液体和蒸汽液体和蒸汽液体液体气体气体 临界点临界点饱和液相线饱和液相线（泡点线）（泡点线）饱和汽相线饱和汽相线（露点线）（露点线）第十页，本课件共有101页11水的水的P-V-T立体相图立体相图第十一页，本课件共有101页12纯物质的纯物质的P-V-T立体相图立体相图 投影图投影图纯物质的纯物质的P-V图图纯物质的纯物质的P-T图图第十二页，本课件共有101页13临界点临界点过冷液体区过冷液体区纯物质的纯物质的P-V图图饱和液相线饱和液相线饱和汽相线饱和汽相线汽液两相平衡区汽液两相平衡区F=C-
6、P+2=1v超临界超临界流体流体区区（TTc和和PPc）过热蒸汽区过热蒸汽区恒温线恒温线第十三页，本课件共有101页14P-V图的特征、相关概念图的特征、相关概念v单相区（单相区（V，G，L）v两相共存区（两相共存区（V/L）v饱和线（饱和液体线、饱和气体线）饱和线（饱和液体线、饱和气体线）v过热蒸汽过热蒸汽v过冷液体过冷液体v等温线（等温线（T=Tc、TTc、TTc和和PPc）第十四页，本课件共有101页15固体区固体区液体区液体区气体区气体区三相点三相点F=C-P+2=0水的三相点：水的三相点：0.0098 临界点临界点汽固平衡线汽固平衡线液固平衡液固平衡线线汽液汽液平衡平衡线线纯物质的纯
7、物质的P-T图图临界等容线临界等容线超临界流体区超临界流体区（TTc和和PPc）第十五页，本课件共有101页16临界点及超临界流体临界点及超临界流体临界点：临界点：气液两相共存气液两相共存的最高温度的最高温度或最高压力。或最高压力。超临界超临界流体流体区区（TTc和和PPc）Tc临界温度；临界温度；Pc 临界压力；临界压力；Vc 临界体积临界体积临界点的数学特征：临界点的数学特征：P-V-T相图中最重要的性质之一相图中最重要的性质之一重要！重要！第十六页，本课件共有101页17v超临界流体（超临界流体（Supercritical Fluid,SCF）在在TTc和和PPc区域内，区域内，气气体、
8、液体变得不可区分，体、液体变得不可区分，形成的一形成的一种特殊状态的流体，种特殊状态的流体，称为称为超临界流体。超临界流体。多种物理化学性质介于气体和液体之间多种物理化学性质介于气体和液体之间,并兼具两者的优点。并兼具两者的优点。具有具有类似于类似于液体液体的密度、溶解能力和传热系数的密度、溶解能力和传热系数,具有具有类似类似于于气体气体的低粘度和高扩散系数。的低粘度和高扩散系数。物质的溶解度对物质的溶解度对T、P的变化很敏感的变化很敏感，特别是在临界状态附近，特别是在临界状态附近，T、P微小变化会导致溶质的溶解度发生几个数量级的突微小变化会导致溶质的溶解度发生几个数量级的突变变，超临界流体正
9、是利用了这一特性，通过对超临界流体正是利用了这一特性，通过对T、P的调控的调控来进行物质的分离。来进行物质的分离。第十七页，本课件共有101页18超临界流体特性超临界流体特性 性质气体超临界流体液体1atm,1530PcTc1atm,1530密度/(g/mL)(0.62)10-30.20.50.61.6粘度/g/(cms)(13)10-4(13)10-4(0.23)10-2扩散系数/(cm2/s)0.10.40.710-4(比液体大1个数量级)(0.23)10-5第十八页，本课件共有101页19 超临界萃取技术的工业应用超临界萃取技术的工业应用v现研究较多的超临界流体包括现研究较多的超临界流体
10、包括:CO2、H2O、NH3、甲醇、乙醇、戊烷、乙烷、乙烯等。、甲醇、乙醇、戊烷、乙烷、乙烯等。v但受溶剂来源、价格、安全性等因素限制。只有但受溶剂来源、价格、安全性等因素限制。只有CO2 应用最多应用最多.临界条件温和临界条件温和Tc=31 ；Pc=7.4MPa。萃取温度低萃取温度低(3 0 50)能保留产品的能保留产品的天然有效活性。天然有效活性。溶解能力强。溶解能力强。惰性（不污染产品）惰性（不污染产品）、价廉易得、价廉易得、选择性良好选择性良好和和产物易产物易从反应混合物中分离出来。从反应混合物中分离出来。第十九页，本课件共有101页20超临界萃取技术的工业应用超临界萃取技术的工业应用
11、v大规模超临界流体萃取兴起于大规模超临界流体萃取兴起于用超临界用超临界CO2成功地从咖啡中提取成功地从咖啡中提取咖啡因。咖啡因。用超临界用超临界戊烷戊烷流体从石油中分离重油组分流体从石油中分离重油组分。v现在用于提取油脂、现在用于提取油脂、香精、色素、药物香精、色素、药物、酶的有效、酶的有效成分成分,制造出真正的天然制造出真正的天然“绿色食品绿色食品”。v脑白金脑白金、鱼油鱼油、小麦胚芽油、沙棘油、啤酒花浸、小麦胚芽油、沙棘油、啤酒花浸膏等天然和医药保健品。膏等天然和医药保健品。第二十页，本课件共有101页21利用超临界利用超临界CO2流体技术从植物流体技术从植物/动物中动物中提取纯天然的高附
12、加值的物质提取纯天然的高附加值的物质v高附加值的天然产品高附加值的天然产品 有樱桃味有樱桃味 US$120/Ib无樱桃味无樱桃味 US$14/Ib v（从樱桃核中提取）（从樱桃核中提取）US$140/IbUS$14/Ib苯甲酮苯甲酮苯甲醇苯甲醇第二十一页，本课件共有101页22v物质特殊的物质特殊的超临界超临界性质，近年来在化性质，近年来在化学工业中得到较多应用学工业中得到较多应用,并成为并成为化化工热力学工热力学领域研究的领域研究的热点热点之一。之一。v可以肯定，超临界流体在化学工业的可以肯定，超临界流体在化学工业的应用将会越来越广泛。应用将会越来越广泛。超临界萃取技术的工业应用超临界萃取技
13、术的工业应用第二十二页，本课件共有101页232.2 流体的状态方程流体的状态方程（EOS）什么是状态方程（什么是状态方程（EOS）？v f(P,V,T)=0 v 状态方程状态方程 Equation of State(EOS)v状态方程的状态方程的类型：类型：a.由实验数据归纳而来由实验数据归纳而来（较少）（较少）（纯经验型纯经验型）b.由理论分析归纳而来由理论分析归纳而来（不多）（不多）（纯理论型纯理论型）如维里如维里方程方程c.由一定理论分析与实验数据相结合归纳而来（多数）由一定理论分析与实验数据相结合归纳而来（多数）（理论（理论+实验实验 即即半经验型半经验型）如）如R-K等等 第二十三
14、页，本课件共有101页24一切实际气体的状态方程必须符合下列条件：一切实际气体的状态方程必须符合下列条件：P0时，该方程趋向理想气体状态方程；时，该方程趋向理想气体状态方程；在临界点处在临界点处 ；由此方程计算得到的由此方程计算得到的PVT关系应与实测数据关系应与实测数据尽可能吻合。尽可能吻合。第二十四页，本课件共有101页25电解质电解质化工热力学化工热力学相平衡（汽液两相逸度相等相平衡（汽液两相逸度相等）汽相汽相液相液相状态方程状态方程EOS聚合物聚合物非电解质非电解质焓平衡焓平衡由由P-V-T得到得到H,S,G活度系数模型活度系数模型i状态方程状态方程EOS经验型经验型第二十五页，本课件
15、共有101页26vEOSEOS是计算热力学性质是计算热力学性质最重要的模型之一最重要的模型之一。为何？为何？.EOS EOS是物质是物质P-V-TP-V-T关系的解析式即关系的解析式即用一个用一个EOSEOS即可精确地代表相当广泛范围内的即可精确地代表相当广泛范围内的P P、V V、T T实验数据，借此可精确地计算所需的实验数据，借此可精确地计算所需的P P、V V、T T数据。数据。2.2.用用EOSEOS可计算不能直接从实验测定的其它热力可计算不能直接从实验测定的其它热力学性质学性质（H H，S S，G G）数据。）数据。3.3.用用EOSEOS可进行相平衡和化学反应平衡计算可进行相平衡和
16、化学反应平衡计算第二十六页，本课件共有101页27v目前已有目前已有150150多种多种EOSEOS。但没有一个。但没有一个EOSEOS能描述在工能描述在工程应用范围内程应用范围内任何气体任何气体的行为。的行为。v状态方程包含的状态方程包含的规律愈多规律愈多，方程就，方程就愈可靠愈可靠；准确性准确性越高，范围越广，模型越有价值。越高，范围越广，模型越有价值。v状态方程的准确度和方程型式的简单性是一对矛状态方程的准确度和方程型式的简单性是一对矛盾。盾。v建立建立EOSEOS的方法：或以理论法为主、或以经验法为主。的方法：或以理论法为主、或以经验法为主。实际应用以实际应用以半经验半理论半经验半理论
17、和和纯经验纯经验的的EOSEOS为主。为主。v本课程介绍各种本课程介绍各种EOSEOS的特点和应用范围，并不要求的特点和应用范围，并不要求建立！建立！第二十七页，本课件共有101页282.2.1 理想气体的状态方程理想气体的状态方程理想气体状态方程：理想气体状态方程：PV=nRT;当当nmol，PV=RT；Z=PV/RT=1v理想气体：分子间的相互作用力可忽略不计；气体分理想气体：分子间的相互作用力可忽略不计；气体分子本身的体积可忽略不计。子本身的体积可忽略不计。v理想气体理想气体EOS是是f(P,V,T)=0 中中最简单的一种形式。最简单的一种形式。理想气体理想气体EOS只只适合压力非常低的
18、气体适合压力非常低的气体,不适合真实气体。不适合真实气体。第二十八页，本课件共有101页29研究理想气体的实际意义研究理想气体的实际意义v1 在在较低压力较低压力和和较高温度较高温度下可用理想气体方程进行计算。下可用理想气体方程进行计算。例如：例如：在大气环境下的空气、氮气以及常压高温条件下在大气环境下的空气、氮气以及常压高温条件下的燃烧烟气等均可视为理想气体。的燃烧烟气等均可视为理想气体。v 2 为真实气体状态方程计算提供初始值。为真实气体状态方程计算提供初始值。v 3 判断真实气体状态方程的极限情况的正确程度，当判断真实气体状态方程的极限情况的正确程度，当v 或者或者 时，时，任何的状态方
19、程都任何的状态方程都还原为理还原为理想气体方程。想气体方程。第二十九页，本课件共有101页30 真实气体的状态方程真实气体的状态方程v真实气体偏离理想行为，真实气体偏离理想行为，理想气体状态方程理想气体状态方程不能描述真实气体的不能描述真实气体的状态，因此出现了：状态，因此出现了：van der Waals（vdW范德华）方程范德华）方程Redlich-Kwong方程方程Soave-Redlich-Kwong方程方程Peng-Robinson方程方程Virial（维里）方程（维里）方程 BWR方程、马丁方程、马丁侯方程等等侯方程等等立立方方型型多参数多参数高次型高次型v真实气体偏离理想行为，真
20、实气体偏离理想行为，理想气体状态方程理想气体状态方程不能描述真实气体不能描述真实气体的状态，因此出现了：的状态，因此出现了：van der Waals（vdW范德华）方程范德华）方程Redlich-Kwong方程方程Soave-Redlich-Kwong方程方程Peng-Robinson方程方程Virial（维里）方程（维里）方程 BWR方程、马丁方程、马丁侯方程等等侯方程等等多参数多参数高次型高次型第三十页，本课件共有101页312.2.2 Virial(维里维里)方程方程Virial方程的形式方程的形式Z为压缩因子为压缩因子当当第三十一页，本课件共有101页32vB B、C C（或（或B
21、B、C C）称作）称作VirialVirial系数，物理意义：系数，物理意义：微观上，微观上，VirialVirial系数系数反映了反映了分子间的相互作用，分子间的相互作用，第二维里系数第二维里系数B B反反映了两个映了两个分子之间的相互作用；第三维里系数分子之间的相互作用；第三维里系数C C反映了反映了三重分子的三重分子的相互作用。相互作用。宏观上，宏观上，VirialVirial系数系数仅是温度的函数。仅是温度的函数。最初的最初的最初的最初的 Virial Virial Virial Virial 方程是以经验式提出的，之后由统计力学得到证明。方程是以经验式提出的，之后由统计力学得到证明。
22、方程是以经验式提出的，之后由统计力学得到证明。方程是以经验式提出的，之后由统计力学得到证明。vVirial 系数的获取系数的获取(1)由统计力学进行理论计算：目前应用很少由统计力学进行理论计算：目前应用很少(2)由实验测定：精度较高由实验测定：精度较高(3)用普遍化关联式计算：方便，但精度不如实验测定的数据用普遍化关联式计算：方便，但精度不如实验测定的数据Virial系数系数v任何状态方程都可以通过幂级数展开，转化为任何状态方程都可以通过幂级数展开，转化为Virial方程的形式方程的形式.（如何将（如何将vdW方程方程展开展开Virial方程的形式？试试看。）方程的形式？试试看。）第三十二页，
23、本课件共有101页33v二种形式的二种形式的Virial方程是等价的，其系数之间也有相方程是等价的，其系数之间也有相互关系。互关系。如何证明？如何证明？Virial方程不同形式的关系方程不同形式的关系试试看试试看第三十三页，本课件共有101页34适用于适用于TTc,P1.5MPa蒸汽蒸汽适用于适用于TTc,1.5MPa P Tc2)T=Tc3)TTc仅有仅有一个实根一个实根，对，对应于超临界流体和应于超临界流体和气体的摩尔体积。气体的摩尔体积。三个重实根三个重实根 V=Vc第三十八页，本课件共有101页39CPVV sl V x V svP*TTc3)TTc三个不同实根，发生于两相区三个不同实
24、根，发生于两相区V大大对应于饱和汽摩尔体积对应于饱和汽摩尔体积V小小对应于饱和液摩尔体积对应于饱和液摩尔体积V中中无物理意义。无物理意义。第三十九页，本课件共有101页40范德华方程常数范德华方程常数a、b的确定的确定对于对于Van der Waals方程方程应用临界条件，即应用临界条件，即（2）v临界等温线在临界等温线在C点点的斜率等于零的斜率等于零v临界等温线在临界等温线在C点点的曲率等于零的曲率等于零把把van der Waals方程代入上述条件，即可得方程代入上述条件，即可得（3）（4）第四十页，本课件共有101页41联立求解方程（联立求解方程（3）和（）和（4）,得得(5)（6）将方
25、程（将方程（1 1）用于临界点，即与（）用于临界点，即与（5 5）式联立，）式联立，即得即得vdW常数常用常数常用形式形式将（将（6）式与）式与（5）联立，即得联立，即得由于由于VC的实验值的实验值误差大误差大，a，b要用要用Pc,Tc来表示来表示第四十一页，本课件共有101页42Zc值是状态方程优劣的标志之一值是状态方程优劣的标志之一范德华方程计算值与真实气体实验值的比较vdW计算值真实气体VC3b2.2bZC0.3750.230.29第四十二页，本课件共有101页43 van der Waals(vdW)EOS1.vdW EOS的缺点：的缺点：两项修正项过于简单，两项修正项过于简单，准确度
26、低，准确度低，不能在任何情况下都能精不能在任何情况下都能精确描述真实气体的确描述真实气体的P-V-T关系。关系。实际应用少。实际应用少。2.vdW EOS的改进的改进改进形式为改进形式为Redlich-Kwong(RK)；Soave RK(SRK)；Peng-Robinson(PR)状态方程等状态方程等但改进形式均以但改进形式均以vdW状态方程为基础状态方程为基础第四十三页，本课件共有101页442.2.3.2 Redlich-Kwong（RK）EOSv改变了方程的引力项改变了方程的引力项Patt，以使得计算的，以使得计算的V减小，试减小，试图改进方程计算图改进方程计算P-V-T的准确性；的准
27、确性；vRK方程明显优于方程明显优于vdW方程，方程，是真正实用的是真正实用的EOS。第四十四页，本课件共有101页45RK方程常数的获取方程常数的获取v用同于用同于vdW方程的方法得到常数方程的方法得到常数a,b值值，即临界等温线在临界点的条件得到：即临界等温线在临界点的条件得到：第四十五页，本课件共有101页46如何用如何用RK方程计算方程计算1）已知）已知T,V，如何求，如何求P？显压型，直接计算，很方便。显压型，直接计算，很方便。在计算时，一定要注意单位，在计算时，一定要注意单位，1atm=0.101325106Pa=0.101325MPa2）已知）已知P,T，如何求，如何求V？工程上
28、最常用的情况，因为工程上最常用的情况，因为P，T易测。用各种易测。用各种迭代法迭代法求解求解。3）已知）已知P,V，如何求，如何求T？用各种用各种迭代法迭代法求解求解。第四十六页，本课件共有101页47方程两边乘以方程两边乘以初值取初值取方法方法1直接迭代法直接迭代法2）已知）已知P,T，如何求，如何求V？求求 ：第四十七页，本课件共有101页48求求 时时：初初值值：两种迭代式只是数学式的表示方法不同，便于选迭代初值，两种迭代式只是数学式的表示方法不同，便于选迭代初值，加快收敛。加快收敛。第四十八页，本课件共有101页49求方程求方程f(x)=0 f(x)=0 的实根。的实根。解：设解：设x
29、0是一个初值是一个初值（初值的设置关系到是否收敛）初值的设置关系到是否收敛）若若|xn-xn-1|0.0001则则xn即为合适的根即为合适的根牛顿迭代法是收牛顿迭代法是收敛最快的方法！敛最快的方法！方法方法2牛顿迭代法牛顿迭代法2）已知）已知P,T，如何求，如何求V？第四十九页，本课件共有101页50已知已知P,T，如何用，如何用牛顿牛顿迭代迭代法法求求V？第五十页，本课件共有101页51小结：已知小结：已知P,T，如何求，如何求 V？1解题的解题的两种方两种方法法迭代法迭代法解析法解析法直接直接迭代法迭代法牛顿迭代法牛顿迭代法2后面讲！后面讲！第五十一页，本课件共有101页523）已知）已知
30、P,V，如何求，如何求 T？牛顿迭代法（方法牛顿迭代法（方法同已知同已知P,T，求，求 V）第五十二页，本课件共有101页53RK方程的方程的特点特点v1、Zc=1/3=0.333，RK方程计算气相体积准确性有了很大提高。方程计算气相体积准确性有了很大提高。v2、RK方程方程能较成功地用于气相能较成功地用于气相能较成功地用于气相能较成功地用于气相P-V-TP-V-T的计算，但的计算，但的计算，但的计算，但计算液相体积计算液相体积的准确性不够，的准确性不够，不能同时用于汽、液两相。不能同时用于汽、液两相。v3、RK方程用于烃类、氮、氢等非极性气体时，即使在压力为方程用于烃类、氮、氢等非极性气体时
31、，即使在压力为几十几十MPa精度都较高，误差仅在精度都较高，误差仅在2%左右左右；但对于氨、水蒸气；但对于氨、水蒸气等等极性较强的气体则精度较差极性较强的气体则精度较差，误差在，误差在10-20%。v4、对对RK方程进行修正，但同时降低了方程进行修正，但同时降低了 RK的简便性和易算性。的简便性和易算性。成功的有成功的有Soave的修正式（的修正式（SRK）等）等第五十三页，本课件共有101页542.2.3.3 Soave-Redlich-Kwong(SRK)方程方程Soave把把RK方程中的方程中的常数常数a看作是看作是温度的函数温度的函数：偏心因偏心因子子v 与与RK方程相比，方程相比，S
32、RK方程可计算方程可计算极性物质极性物质，更主要的是可计算，更主要的是可计算饱和液体密度，使之能用于混合物的汽液平衡计算，故在工业上获饱和液体密度，使之能用于混合物的汽液平衡计算，故在工业上获得了得了广泛应用。广泛应用。第五十四页，本课件共有101页55用用SRK状态方程计算的状态方程计算的甲烷丙烷的恒温压力甲烷丙烷的恒温压力组成图。组成图。计算值与实验值非计算值与实验值非常符合常符合第五十五页，本课件共有101页56 例例2-2 P17 v应用状态方程时要注意：应用状态方程时要注意：v(1)(1)式中式中P P是是绝对压力，绝对压力，SISI制单位：制单位：papav(2)(2)式中的式中的
33、T T是是绝对温度绝对温度,K,Kv(3)(3)式中的式中的V V是是单位单位体积，体积，SISI制单位：制单位：mm3 3 v(4)(4)式中的式中的n n为为1 1 molmolv(5)(5)通用气体常数通用气体常数R R的的单位必须和单位必须和P,V,TP,V,T的单位相适应的单位相适应。v注：注：R=8.314R=8.314mm3 3.pa/mol.K.pa/mol.K=8.3148.314 J/mol.KJ/mol.K常出现的问题常出现的问题：因：因量纲量纲的问题，易出现计算结果错误。的问题，易出现计算结果错误。解决方法解决方法：为避免此类错误的产生，且计算过程方便、简易起见，各物理
34、量均换为避免此类错误的产生，且计算过程方便、简易起见，各物理量均换算为算为SI制单位制单位，再代入方程式中进行计算。，再代入方程式中进行计算。第五十六页，本课件共有101页572.2.3.4 PengRobinson 方程RK 方程和RKS 方程在计算临界压缩因子Zc 和液体密度时都会出现较大的偏差，为了弥补这一明显的不足，PengRobinson 于1976 年提出了他们的方程，简称PR 方程：见P18PR方程预测液体摩尔体积的准确度较方程预测液体摩尔体积的准确度较SRK有明显的改善，有明显的改善，也可以用于极性物质。也可以用于极性物质。2.2.4 多常数状态方程多常数状态方程 自阅自阅第五
35、十七页，本课件共有101页58状态方程的选用状态方程的选用真实流体真实流体Zc=0.230.29EOSZc适合范围优缺点理想气体1压力极低的气体不适合真实气体vdW0.375同时能计算汽，液两相准确度低RK0.333计算气相体积准确性高，很实用不能同时用于汽、液两相SRK0.333能同时用于汽液两相平衡，广泛应用精度高于RK，能预测液相体积PR0.307能同时用于汽液两相平衡，广泛应用能预测液相体积VirialTTc，P5MPa的气相不能同时用于汽液两相第五十八页，本课件共有101页592.3 对应态原理对应态原理的应用的应用2.3.1 2.3.1 普遍化状态方程普遍化状态方程 虽然所有立方型
36、状态方程均可用虽然所有立方型状态方程均可用解析法求解解析法求解V V，但工程上大多但工程上大多采用简便的采用简便的迭代法迭代法.为了迭代方便，将为了迭代方便，将RKRK方程对方程对V V的的隐式隐式，变换成，变换成Z Z的的显式显式即方程即方程(a)(a)，这，这是是RKRK方程的另一种形式。方程的另一种形式。（a）(b)式中：式中：第五十九页，本课件共有101页60是是否否v计算步骤：计算步骤：第六十页，本课件共有101页61应用：应用：简单迭代即简单迭代即取取Z=1Z=1，代入（，代入（b b）式）式求出求出h h，将将h h代入（代入（a a）式）式求得一个新求得一个新Z Z，再将此，再
37、将此Z Z代入（代入（b b）式重复计算至两次得）式重复计算至两次得到的到的Z Z值满足某一精度为止。值满足某一精度为止。SRKSRK方程方程的普遍化形式为的普遍化形式为P22（2-21a）（）（2-20d）。例例 例例2-4 2-4 P22P22应注意的问题：单位一致。应注意的问题：单位一致。第六十一页，本课件共有101页62气体的对比态原理：气体的对比态原理：通过大量的实验发现，许多物质的气体当接近临界点时，都显示出相通过大量的实验发现，许多物质的气体当接近临界点时，都显示出相似的性质，因而引出了似的性质，因而引出了对比参数对比参数的概念。的概念。第六十二页，本课件共有101页63 2.3
38、.2 两参数普遍化压缩因子图两参数普遍化压缩因子图v对比态原理认为：对比态原理认为：在相同的对比状态下，所有物质表现出相同的性质。在相同的对比状态下，所有物质表现出相同的性质。即：即：组成、结构、分子大小相近的物质有相近的性质。组成、结构、分子大小相近的物质有相近的性质。在相同对比温度、对比压力下，任何气体或液体的对比体在相同对比温度、对比压力下，任何气体或液体的对比体积（或压缩因子）是相同的积（或压缩因子）是相同的。即当即当 相同时，相同时，也相同也相同。写成状态方程为：写成状态方程为：f(Pr，Tr，Vr)=0 （A）v对比（应）态原理（对比（应）态原理（CSPCorresponding
39、State Principle）是一种特别的）是一种特别的EOS。第六十三页，本课件共有101页64基于此对比态原理得出的两参数的基于此对比态原理得出的两参数的图，实际上只适用于图，实际上只适用于“简单分子简单分子”：即分子四周的力场是高度对称的，就是：即分子四周的力场是高度对称的，就是说，分子之间的作用力仅与距离有关，与方位无关。严格的讲，仅球形非极说，分子之间的作用力仅与距离有关，与方位无关。严格的讲，仅球形非极性分子（如氩、氪、氙等）属于这类简单分子。氧、氮和一氧化碳等，接近性分子（如氩、氪、氙等）属于这类简单分子。氧、氮和一氧化碳等，接近于简单分子。于简单分子。图应用于图应用于H H2
40、 2,He,Ne,He,Ne 等量子气体时，误差太大，等量子气体时，误差太大，进行经验校正：进行经验校正：(P(P和和PcPc单位为单位为MPaMPa)(T T和和TcTc单位为单位为K K)简单对比态原理得出的两参数普遍化压缩因子图简单对比态原理得出的两参数普遍化压缩因子图5050年代应用它，目前应用很少。讲义年代应用它，目前应用很少。讲义 P25P25 图图2-42-4（a a）（）（b b）（）（c c）三张图。）三张图。（当初是使用（当初是使用atm，所以校正都是加上，所以校正都是加上8）第六十四页，本课件共有101页65在低压下任何真实气体都服从在低压下任何真实气体都服从 PV=RT
41、PV=RT，将，将 代入式代入式PV=RT得：得：则则由此可看出：只有在由此可看出：只有在各种气体的临界压缩因子各种气体的临界压缩因子都相等的条件下都相等的条件下（A A）式才成立。实际上物质的）式才成立。实际上物质的Zc0.20.3之间，不是一个常数。之间，不是一个常数。可见范德华提出的对比态原理只是一个近似的关系。后人称它为可见范德华提出的对比态原理只是一个近似的关系。后人称它为简单对比态原理。简单对比态原理。对比状态原理尽管不太严密，但在实际当中很有指导意义。对比状态原理尽管不太严密，但在实际当中很有指导意义。两参数对比态原理的缺点：两参数对比态原理的缺点：第六十五页，本课件共有101页
42、661 1）两参数）两参数CSPCSP只能适合于简单的球形流体只能适合于简单的球形流体(如如 Ar,Kr,Xe)Ar,Kr,Xe)。2 2）为了提高对比态原理的精度）为了提高对比态原理的精度 ,人们引入了第三参数人们引入了第三参数.第三参数的特性：第三参数的特性：最灵敏最灵敏反映物质分子间相互作用力的物性参数，反映物质分子间相互作用力的物性参数，当分子间的作用力稍有不同，就有明显的变化。当分子间的作用力稍有不同，就有明显的变化。Lydersen Lydersen 等人于等人于19551955年提出以年提出以ZcZc作为第三参数，即：作为第三参数，即：想法很好，但想法很好，但VcVc数据不准确导
43、致数据不准确导致ZcZc不准确。不准确。3 3）比较成功的第三参数为）比较成功的第三参数为Pitzer提出的提出的偏心因子偏心因子 。第六十六页，本课件共有101页67球形分子球形分子 称称为简单为简单流体（分子流体（分子间间作用力作用力简单简单）非球形分子非球形分子 称为复杂流体（称为复杂流体（分子间作用力变得复杂了分子间作用力变得复杂了）如水分子如水分子 偏离球形偏离球形用一（物理）量用一（物理）量来表示分子偏离球形的程度，那么来表示分子偏离球形的程度，那么的物理意的物理意义是什么？义是什么？表征物质分子的偏心度，即非球形分子偏离球对称的程度。表征物质分子的偏心度，即非球形分子偏离球对称的
44、程度。其值大小是反映物质分子形状与物质极性大小的量度。其值大小是反映物质分子形状与物质极性大小的量度。如何计算此值的大小呢？如何计算此值的大小呢？因为经典热力学不研究分子结构，所以不用几何的方法来计算偏离球因为经典热力学不研究分子结构，所以不用几何的方法来计算偏离球形的程度（统计热力学这样做）我们形的程度（统计热力学这样做）我们采用宏观的物理量来表征微观分采用宏观的物理量来表征微观分子的行为！子的行为！学会这种解决问题的思路和方法，对今后的工作会有帮助。学会这种解决问题的思路和方法，对今后的工作会有帮助。2.3.2 三参数普遍化压缩因子图三参数普遍化压缩因子图第六十七页，本课件共有101页68
45、 采用宏观的物理量来采用宏观的物理量来表征微观分子的行为：表征微观分子的行为：从克克方程出发从克克方程出发认为该方程也符合对比态原理：认为该方程也符合对比态原理：克克方程的适用范围：从三相点克克方程的适用范围：从三相点临界点临界点对于临界点对于临界点 方程简化为：方程简化为：即为一个直线方程，斜率为即为一个直线方程，斜率为-a-a。第六十八页，本课件共有101页69若简单对比态原理准确，若简单对比态原理准确，T Tr r相同时，相同时，P Pr r也一定相同。因为两相区也一定相同。因为两相区F=C-P+2=1-2+2=1F=C-P+2=1-2+2=1，只有一个自由度。所以只有一个自由度。所以图
46、应为一条直线（所有物质都符合）图应为一条直线（所有物质都符合）实际上，只有氩、氪、氙等简单流体在一条直线上。其它物质不在实际上，只有氩、氪、氙等简单流体在一条直线上。其它物质不在这条直线上。用这些直线在某一点上偏离氩、氪、氙等简单流体直这条直线上。用这些直线在某一点上偏离氩、氪、氙等简单流体直线的距离（纵坐标之差）表示其分子偏离球形的程度。取线的距离（纵坐标之差）表示其分子偏离球形的程度。取TrTr0.70.7时，时，lgPlgPr rS S与简单流体该值的差值为与简单流体该值的差值为 的值。因为：的值。因为：简单流体简单流体TrTr0.70.7时，时，lgPlgPr rS S-1.0-1.0
47、，便于计算便于计算第六十九页，本课件共有101页701.01.21.41.61.8-1-1.2-1.8lgPrs1/Tr 1 2简单流体简单流体(Ar,Kr,Xe(Ar,Kr,Xe）非球形分子非球形分子1 1(正癸烷正癸烷)非球形分子非球形分子2 2 对比蒸汽压与温度的近似关系对比蒸汽压与温度的近似关系1/Tr=1.43即即Tr=0.7第七十页，本课件共有101页71 偏心因子。偏心因子。0 1,愈大，愈大，分子分子偏离球对称的偏离球对称的程度愈大程度愈大。CH4=0.008。球形分子或者说简单流体的。球形分子或者说简单流体的 为为0。各物质的各物质的 值见值见 附表附表1 P332第七十一页
48、，本课件共有101页72所以所以简单流体值简单流体值 复杂流体校正项复杂流体校正项 即即Z Z0 0 查查 图图2 26 6，图，图2 27 7 也即两参数普遍化也即两参数普遍化Z Z图图Z Z1 1 查查 图图 2 28 8，图，图2 29 9复杂流体校正项（复杂流体校正项（P26-27P26-27）查附表查附表1 1（Tc,PcTc,Pc也在此表查也在此表查）因为蒸汽压数据比临界性质因为蒸汽压数据比临界性质VcVc更易测准，所以更易测准，所以 为第三参数的为第三参数的Z Z值精值精度更高，对于烃类物质，如果不在临界区内，误差在度更高，对于烃类物质，如果不在临界区内，误差在1 1以内。所以以
49、内。所以此图目前应用较多。但对极性物质精度稍差。此图目前应用较多。但对极性物质精度稍差。Z Z0 0还可查还可查P349 P349 附表附表7 72 2；Z Z1 1 查查 P351 P351 附表附表7 73 3。第七十二页，本课件共有101页732.3.4 2.3.4 普遍化第二维里系数关系式普遍化第二维里系数关系式 大家可能注意到了：大家可能注意到了：P18P181919页的页的Z Zo o、Z Z1 1 图左上角有一块区域有的未图左上角有一块区域有的未划出来，即使划出来了线条也比较密，不好用。这是因为低压、高温划出来，即使划出来了线条也比较密，不好用。这是因为低压、高温时时Z Z值接近
50、值接近1 1，所以相对误差较大，即使划出来了也尽量不要使用。怎么，所以相对误差较大，即使划出来了也尽量不要使用。怎么办呢？采用计算的方法办呢？采用计算的方法用普遍化第二维里系数关系式计算：用普遍化第二维里系数关系式计算：定义定义 这里这里所以所以B与压力无关与压力无关求出求出T Tr r、P Pr r即可计算出即可计算出Z Z，很方便。倒着算：，很方便。倒着算：T Tr r B B0 0、B B1 1 Z Z 第七十三页，本课件共有101页74注意：注意：(1)(1)使用范围中低压高温，查普遍化使用范围中低压高温，查普遍化Z Z图困难时。（第二维里截断图困难时。（第二维里截断式只适先查图，查不