第三章 纯流体的热力学性质计算ppt课件

.,第三章纯流体的热力学性质计算,3.1热力学关系式3.2以T、p为变量的焓变和熵变计算3.3剩余性质3.4用剩余性质计算气体热力学性质3.5液体的热力学性质3.6两相系统,.,第三章纯流体的热力学性质计算,在工业设计过程中物质的流体热力学性质和需要能量的数值十分重要，否则难以实现即符合实际需要，又能节能的优良工业设计。简而言之，对一台膨胀机的设计，当其在绝热条件下操作，压力从p1降到p2，做多少功？从热力学第一定律在开放稳流过程的能量平衡知-Ws=H=H2-H1,.,3.1.1热力学基本关系式,单位摩尔数定组成均相流体的热力学基本关系式：dU=TdSPdV(3-1)dH=TdS+VdP(3-2)dA=-PdVSdT(3-3)dG=VdP-SdT(3-4),.,3.1.2马克斯韦尔（Maxwell）式,（3-6a）(6-6b)(3-6c)(3-6d),.,3.2以T、p为变量的焓变和熵变计算,先考虑在恒压的条件下焓对温度求偏导（H/T）p=Cp(3-7)dH=TdS+VdP在恒压下除以dT可得（H/T）p=T(S/T)p结合式（3-7）可得(S/T)p=Cp/T(3-8)恒温下熵对压力求偏导，可直接由式（3-6d）表示（S/P）T=-(V/T)p（3-6d）,.,3.2焓和熵的计算式,式（3-2）在恒温下除以dP,得（H/P）T=T(S/P)T+V结合式（3-6d），上式变成（H/P）T=V-T(V/T)p(3-9),.,3.2焓和熵的计算式,焓和熵与温度、压力的函数关系为H=H(T,P),S=S(Y,P)其微分式为dH=(H/T)pdT+(H/P)TdPdS=(S/T)pdT+(S/P)TdP,.,3.2焓和熵的计算式,将式（3-7）、式（3-8）和式（3-9）代入上面两式得dH=CpdT+V-T(V/T)pdP(3-10)dS=Cp(dT/T)-(V/T)pdP(3-11)对于理想气体PV*=RT,则(V*/T)p=R/P代入式（2-49）与式（2-50）得dH*=C*pdT(3-12)dS*=C*p(dT/T)-(R/P)dP(3-13),.,3.3剩余性质,3.3.1自由焓可作为母函数式dG=VdPSdT(3-4)表示了G是P和T的函数：G=G（P，T）由于变量P和T可以直接测量或者控制，因此，自由焓G是最有潜在应用价值的热力学性质。另一个基本性质关系式是根据下述恒等式得出来的：将式（3-4）的dG代入上式的右边，同时右边的G用G=H-TS取代，整理后得：(3-14),.,3.3.1自由焓可作为母函数,根据全微分是偏微分之和的原理，由式（3-14）直接可以得到(3-15)(3-16)G/RT是T、P的函数，只要用简单求导的方法就可以得到V/RT与H/RT之值。其余的性质通过定义式均可求得，如：,.,3.3.1自由焓可作为母函数,虽然已导出式（3-14），得出GRT(或G）与变量T、P的关系，且能给出GRTG(T，P)的关系式，经过某些数学处理就可以求出所有其他的热力学性质。从这个意义上说，自由焓可视为其他热力学性质的母函数，它包含了完整的热力学性质的信息。但遗憾的是目前尚不能用简便的实验方法直接测量GRT（或G）的值，因此，根据自由焓直接导出的各种方程式尚无实用价值。但作为其他性质的一个母函数，可以从自由焓引出一个有关的相近的性质，而此性质却又是易于求值的，这样，就可以解决诸性质的计算问题。,.,3.3.2剩余性质的引入,定义剩余自由焓GR为GR=G-G*(3-17)式中GR与G*分别为在相同温度和压力下真实气体与理想气体的自由焓。剩余体积VR=V-V*（3-18）VR=V-RT/P由于V=ZRT/P，剩余体积和Z就可以进行关联VR=（RT/P）（Z-1）（3-19）,.,3.3.2剩余性质的引入,按此，可以写出剩余性质的通式MRM-M*（3-20）式中：M为摩尔性质。式（3-14）也适用于理想气体，即（3-21）式（3-14）减去上式，得(3-22),.,3.3.2剩余性质的引入,（3-22）式可应用于定组成流体，是各剩余性质之间的基本关系式。根据此式可得（3-23）（3-24）根据自由焓的定义式G=H-TS，对于理想气体，有G*=H*-TS*两着之差即为GR=HR-TSR由此式可给出熵（3-25）,.,3.3.2剩余性质的引入,从上述推导过程可见，自由焓不仅是其他性质的母函数，而且还可以将其与实验数据进行关联，由式（3-22）可得积分之，压力从零到任意压力P，则有零压态即为理想气体态。因此，当压力零时GR/RT也为零。结合式（3-19），上式写成(3-26),.,3.3.2剩余性质的引入,遵循式（3-24）的关系，将上式对温度求偏导，则有（3-27）由式（3-26）、式（3-27）代入式（2-25），可得（3-28）HR/RT和SR/RT都为无因次项，HR和SR的单位取决于通用气体常数的选择。至此已顺利地推出以p、T为变量的VR、HR、SR和GR的表达式。,.,3.3.3剩余性质与偏离性质的异同,MD=M（T，p）-M*（T，p0，ig）（3-29）MR=M（T，p）-M*（T，p，ig）（3-29）,.,3.4用剩余性质计算气体热力学性质,3.4.1真实气体的焓与熵压缩因子定义式为ZPV/RT，Z与(Z/T)P之值可根据实验的PVT数据求取。积分式(264)至式(266)可用数值解，也可以用图解法求解。若用状态力程来表达Z，这些积分式就可用解析法求解。因此只要给出PVT数据或合适的状态方程，就能求出HR和SR或其他的剩余性质。正是由于剩余性质与实验数据有直接联系，因此在实际应用中，这些剩余性质是非常重要的。对于焓和熵，式320可写成H=H*+HR(3-30)S=S*+SR(3-31)可见，H和S之值可根据相应的理想气体性质与剩余性质两者相加求得。,.,3.4.1真实气体的焓与熵,将式（3-12）和式（3-13）积分，则可求得理想气体的H*和S*。表达式的积分下限为参考态T0、P0，上限为物系所处的状态T，P：(3-12a)(3-13a)式中Cp*为理想气体的热容。将上面两式代入式（3-30）、式（3-31）。得：(3-32)(3-33),.,3.4.1真实气体的焓与熵,为了方便，可将以上两式写成：（3-32a）（3-33a）式中C*pmh和C*pms分别为理想气体求焓变和熵变需用的平均等压热容，其值可分别用下述两式求得(3-34a)(3-34b)C*p仅是温度的函数，其函数式通常为,.,3.4.1真实气体的焓与熵,式中、或a、b、c均为与物质有关的常数。为方便计，上述两式可合并成式（3-35）Cp*/R=A+BT+CT2+DT-2（3-35）式中A、B、C和D仍为与物质有关的常数。一些常用的有机和无机气体的A、B、C、D的值列于附表2。Cp*的单位要与通用气体常数的单位一致。,.,3.4.1真实气体的焓与熵,式（3-34a）与式（3-34b）用于一般的焓变和熵变计算时，当温度从T1至T2时可写成：(3-34c)(3-34d)将式（3-35）分别代入上面两式，积分后得（3-34e）,.,3.4.1真实气体的焓与熵,式中Tam为算术平均温度（3-34f）式中Tim为对数平均温度,.,3.4.1真实气体的焓与熵,根据热力学第一和第二定律导出的热力学性质方程式不能求出焓、熵的绝对值，只能求出其相对值，理想气体参考态（温度T0，压力P0）是计算焓熵的其始点，或者说参考态的焓H0*、熵S0*是计算焓、熵的基准。参考态（T0，p0）是根据计算的方便性而随意确定的，同样，H0*和S0*也是任意指定的。式（3-32）和式（3-33）的计算仅需要理想气体方程和其热容以及真实气体的p-V-T数据。在给定的p、T下求得V、H和S后，则其他的性质便可根据定义式求出。由式（3-32）和式（3-33）可见，理想气体是计算真实气体性质的基础。,.,例3-1,.,3.4.2用普遍化关联计算剩余性质,式（3-27）和式（3-28）用对比参数P=PcPrT=TcTrdp=PcdPrdT=TcdTr(3-36)(3-37)上面两式涉及的变量只有Z、Tr、Pr，因此对于任何给定的Tr和Pr值，根据普遍化因子Z的数据，就可以从上述两式求出HR/RTc和SR/R的值。对于Z的关联，基于式（2-28a）Z=Z0+Z1（2-28a）,.,3.4.2用普遍化关联计算剩余性质,恒Pr下，对Tr求偏导，得：式（3-36）和式（3-37）中的Z与用上面两式代入，得（3-38）（3-39）,.,3.4.2用普遍化关联计算剩余性质,上述两式右边第一积分项，对不同的Tr和Pr值，可根据图26与图27用数值或图解积分求得。在后面的积分项根据图28与图29用类似的方法求得。若将此两式中的第一项积分值分别用(HR)0/RTc和(SR)0/R表示，第二项积分值相应地用(HR)1/RTc和(SR)1/R表示，则可写成HR/RTc=(HR)0+(HR)1/RTc(3-40)SR/R=(SR)0+(SR)1/R(3-41),.,3.4.2用普遍化关联计算剩余性质,在不同的Tr和Pr下求得的(HR)0、(HR)1、(SR)0、(SR)1值已标绘在图3-1至3-8中。这些图都是上述诸项在各恒定Tr下对Pr标绘的。将这些图与式3-40和式3-41联合应用就可以求出HR和SR之值。显然，此计算方法是以皮策提出的三参数对应状态原理为基础的。,.,3.4.2用普遍化关联计算剩余性质,与普遍化压缩因子关联相似，由于复杂的函数关系使得（HR）0/RTC、（HR）1/RTC、（SR）0/R和（SR）1/R无法用简单的方程式表达。但是在低压下普遍化维里系数对Z的关联方法对于剩余性质仍然适用。B0和B1与Z的关联式为式（2-7）和式（2-29），将两式合并，可得（2-30）式,.,3.4.2用普遍化关联计算剩余性质,（2-30）由此式可导出将以上两式代入式（3-36）、（3-37），可得（3-42）（3-43）,.,3.4.2用普遍化关联计算剩余性质,由于B0和B1仅是温度的函数，在恒温下积分可得（3-44）（3-45）B0和B1与温度的关系由式（2-31a）和（2-31b）给出。它们对Tr求导亦较简便，下面四式可直接应用于（3-44）和（3-45）的计算。,.,3.4.2用普遍化关联计算剩余性质,B0=0.0830.422/Tr1.6(3-31a)dB0/dTr=0.675/Tr2.5(3-46)B1=0.1390.172/Tr4.2(3-31b)dB1/dTr=0.722/Tr5.2(3-47),.,3.4.2用普遍化关联计算剩余性质,根据式（3-32a）和式（3-33a），HR、SR普遍化关联和理想气体的热容联合应用可以求出任何温度和压力下的焓和值熵。设某物系从状态1变到状态2，用式（3-32a）写出这两个状态的焓值H2=H*0+C*pmh(T2-T0)+HR2H1=H*0+C*pmh(T1-T0)+HR1过程的焓变为上述两式之差，即H=H2-H1H=C*pmh(T2-T1)+HR2-HR1(3-48)同样，也有S=C*pmsln(T2/T1)-Rln(P2/P1)+SR2-SR1(3-49),.,3.4.2用普遍化关联计算剩余性质,上面两式右边诸项可以与物系从初态到终态的计算途径进行联系，（见图3-9）。实际过程是1到2（用虚线表示），可以设想用三步计算途径来实现。1到1*表示在T1和P1下由真实气体转化为理想气体，这是虚拟的，其焓变和熵变为H*1-H1=-HR1S*1-S1=-SR11*到2*是理想气体从状态1*（T1，P1）到达状态2*（T2，P2），此过程的焓变和熵变为H*=H*2-H*1=C*pmh(T2-T1)(3-50)S*=S*2-S*1=C*pmhln(T2/T1)-Rln(P2/P1)(3-51)最后2*到2，在T2、P2下由理想气体回到真实气体，这也是虚拟的过程，其焓变和熵变为：H2-H*2=HR2S2-S*2=SR2将三个过程的焓变和熵变相加，即为式（3-48）和式（3-49）,.,3.4.2用普遍化关联计算剩余性质,.,习题1,工程设计中需要乙烷在3446kPa和93.33下的体积数据，已查到的文献值为0.02527m3/kg，试计算下列诸方法进行核算：（1）两参数压缩因子法；（2）三参数压缩因子法,.,习题1,解：查表得乙烷的特性常数为：Tc=305.4K，pc=4.884MPa,=0.098，M=30.070（1）两参数压缩因子法Tr=T/Tc=366.48/305.4=1.20Pr=p/pc=3.446106/4.884106=0.71由Tr，pr查表得Z=0.86因为pV=ZRT，则V=ZRT/p=0.868.314366.48/3.446106=0.00076（m3/mol）=0.025m3/kg,.,3.4.3用状态方程计算剩余性质,以上对剩余性质的求得要用查图的方法，虽然简便，但图的精准性欠佳，当然希望能用数值估算来得出式（3-26）式（3-28）中的积分，这样就会想到利用EoS来完成，但所用的EoS必须是以V为显函数的方程，即以p和T为独立变量。这类方程形式在立方型EoS中并不常见，但容易想到维里方程能完成此项任务。,.,3.4.3.1维里方程式,由第2章中可知，可用第二维里系数B来表达压缩因子，即Z-1=Bp/RT（2-7）将（2-7）代入式（3-26）后，可得：（3-52）将式（3-52）对T微分，并结合（3-24），则（3-53）将式（3-52）和式（3-53）代入式（3-25）得（3-54）,.,3.4.3.1维里方程式,式（3-52）（3-54）乃是用二项截断维里方程导出的剩余性质的计算方程，或者也可以说为用第二维里方程系数描述的剩余性质的函数表达式，适用于较低压力条件下来计算剩余性质。,.,3.4.3.1维里方程式,VanNess和Abbott给出了从三项截断式导出的HR和SR的表达式（3-55）（3-56）按式（3-25），则GR的表达式为（3-57）,.,3.4.3.1维里方程式,式中，C为第三维里系数，为密度。由于C的因如，式（3-55）式（3-57）可在较高的压力下应用。比较式（3-52）（3-54）与式（3-55）（3-57），当T和p给定后，并已知B和T的函数关系，则式（3-52）（3-54）就直接可用来计算剩余性质，但在式（3-55）（3-57）的运用中则还需要知道C和T的函数关系外，尚要先来估算出Z和的值，一般Z和的关系可用（2-8）式来表示，即Z=1+B+C2(2-8)若已知T和的值，则在估算Z或时将会首先用到上式,.,3.4.3.2立方型状态方程式,第二章中介绍的vdW、RK、SRK和PK方程无一不是以V、T为独立变量的方程，用云此类以p为显函数的表达式是有其原因的，用个常见的例子加以说明，当在100和101325Pa下的水会有两个完全不同的平衡态：饱和水蒸汽和液态水。在恒定压力和温度下，体积显函数表达式只能给出一个体积值，但同样的p、T条件下，能适用气液平衡两相的压力显函数表达式中却能给出两个体积值。这时用EoS来计算相关的相平衡问题，至关重要。但是将以p为显函数表达式来计算剩余性质，如用式（3-26）（3-28）会有困难，必须要将式（3-26）（3-28）转换成以T、V（）为独立变量的函数式。,.,3.4.3.2立方型状态方程式,从实用角度看，用T、p为独立变量似乎更为方便，可以通过以下方式进行转换（3-58）在等温下进行微分：并用RT=p/Z代入上式，改写成（3-59）,（恒T）,.,3.4.3.2立方型状态方程式,将（3-59）式代入式（3-26），得（3-60）式（3-60）在等温下积分，当p0，0。Smith、VanNess和Abbott已知给出了HR/RT的T与函数式，即（3-61）再从式（3-25）结合式（3-60）和式（3-61）可得出SR/R的表达式（3-62）式（3-60）式（3-62）给出了以、T为独立变量的剩余性质表达式。,.,3.4.3.2立方型状态方程式,若用立方型EoS计算剩余性质，必须坚持各种具体的EoS结合式（3-60）（3-62）进行，但具体的立方型EoS形式众多，可以得出许多不同的剩余性质的函数表达式，不仅表现形式不够集中，查阅也不方便，因此在Smith等的名著中采用了通式来表达立方型EoS，具体形式为：（3-63a）,.,3.4.3.2立方型状态方程式,针对不同的EoS，式中的参数会有所不同，这在表3-3中列出。若式（3-63a）乘以（V-b）/RT后。则可写成：(3-63b)式（3-63）可以用解析法求解，但一般用迭代法更为常见。,.,3.4.3.2立方型状态方程式,表3-3,.,3.4.3.2立方型状态方程式,从式（3-60）（3-62）知，要用立方型EoS来计算剩余性质，需要有相关的表达式。若将式（3-63a）除以RT，并用1/代替V，则得：(3-64)式中(3-65),.,3.4.3.2立方型状态方程式,有了（3-64），就可相应的写出（Z/T）的表达式，从而按式（3-60）（3-62）推导出用立方型EoS通式来计算剩余性质方程。(3-66)(3-67)(3-68),.,3.4.3.2立方型状态方程式,式中。除q由式（3-65）表达外，其余的参数为I可分为两种情况当(3-70)当=(3-71)自此，可按式（3-66）（3-68）用立方型EoS来计算流体的剩余性质。,.,例3-4,.,3.5液体的热力学性质,3.5.1以T和p为变量表达焓变和熵变流体体积膨胀系数的定义为(1/V)(V/T)p(3-72)式（3-6d）和式（3-9）用表达时可消去(V/T)p相，即(S/P)T=-V(3-73)(H/P)T=(1-T)V(3-74),.,3.5.1以T和p为变量表达焓变和熵变,对于液体，当状态远离临界点时，其体积较小，也较小。通常压力对液体的焓、熵和内能等影响较小。对于不可压缩流体，值可视为零，此时(S/P)T也为零，即熵和压力无关。但由式（3-73）和式（3-74）知，压力的影响，特别是在高压下，还是应该计及的。式（3-10）和式（3-11）中的（V/T）p用代入，可得：,.,3.5.1以T和p为变量表达焓变和熵变,式（3-10）和式（3-11）中的(V/T)p用代入，可得dH=CpdT+V(1-T)dP(2-77)dS=Cp(dT/T)-VdP(2-78)人们通常用上述两式的积分式求出液体的焓变和熵变。由于液体的和V是P的弱函数，积分时可取适当的平均值简化计算。,.,3.5.2以T和V为变量表达内能、熵的变化,对内能和熵的变化，采用T和V为独立变量会更方便。根据式（3-1）可写出和根据Cv的定义可将上面左边的方程写成（3-77）将式（3-6c）代入上面右边方程得：（3-78）,.,3.5.2以T和V为变量表达内能、熵的变化,因U=U（T，V），S=S（T，V）故：（3-79）（3-80）将Cv的定义式和式（3-78）代入式（3-79），得（3-81）将式（3-77）和式（3-6c）代入式（3-80），得（3-82）,.,3.5.2以T和V为变量表达内能、熵的变化,式（3-81）和式（3-82）就是以T、V为变量来表达内能、熵的变化。为了进一步计算，再引进另外一个热力学参数，等问压缩率k，其定义式为：（3-83）由于，将与的定义式代入上式，得（3-84）,.,3.5.2以T和V为变量表达内能、熵的变化,将式（3-84）代入式（3-82），得（3-85）（3-86）若拥有Cv、和的数据就能用式（3-85）和式（3-86）计算内能、熵的变化。,.,3.6两相系统,图2-1示出了纯物质的p-T图，图中的曲线2C代表气液的相界面。每当越过此曲线时，意味着在等温和等压条件下发生了相变。导致广度热力学性质的摩尔值（或比值）发生剧烈的变化，这在图2-2中看得更为清晰，即饱和蒸汽的摩尔体积值比饱和液体体积值大得多，对于其他的广度性质，如H和S都会有类似的情况，但也有例外，在汽化过程时，等温和等压条件下，两