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

第三章纯流体的热力学性质计算3.1热力学关系式3.2以T、p为变量的焓变和熵变计算3.3剩余性质3.4用剩余性质计算气体热力学性质3.5液体的热力学性质3.6两相系统.第三章纯流体的热力学性质计算在工业设计过程中物质的流体热力学性质和需要能量的数值十分重要，否则难以实现即符合实际需要，又能节能的优良工业设计。简而言之，对一台膨胀机的设计，当其在绝热条件下操作，压力从p1降到p2，做多少功？从热力学第一定律在开放稳流过程的能量平衡知-Ws=△H=H2-H1.3.1.1热力学基本关系式单位摩尔数定组成均相流体的热力学基本关系式：dU=TdS–PdV(3-1)dH=TdS+VdP(3-2)dA=-PdV–SdT(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） 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示（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)p]dP(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=VdP–SdT(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），得出G／RT(或G）与变量T、P的关系，且能给出G／RT＝G(T，P)的关系式，经过某些数学处理就可以求出所有其他的热力学性质。从这个意义上说，自由焓可视为其他热力学性质的母函数，它包含了完整的热力学性质的信息。但遗憾的是目前尚不能用简便的实验方法直接测量G／RT（或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剩余性质的引入按此，可以写出剩余性质的通式MRM-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真实气体的焓与熵压缩因子定义式为Z＝PV/RT，Z与(Z/T)P之值可根据实验的PVT数据求取。积分式(2—64)至式(2—66)可用数值解，也可以用图解法求解。若用状态力程来表达Z，这些积分式就可用解析法求解。因此只要给出PVT数据或合适的状态方程，就能求出HR和SR或其他的剩余性质。正是由于剩余性质与实验数据有直接联系，因此在实际应用中，这些剩余性质是非常重要的。’对于焓和熵，式3—20可写成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）用对比 参数 转速和进给参数表a氧化沟运行参数高温蒸汽处理医疗废物pid参数自整定算法口腔医院集中消毒供应 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值，可根据图2—6与图2—7用数值或图解积分求得。在后面的积分项根据图2—8与图2—9用类似的方法求得。若将此两式中的第一项积分值分别用(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.083–0.422/Tr1.6(3-31a)dB0/dTr=0.675/Tr2.5(3-46)B1=0.139–0.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-H1H=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.446×106/4.884×106=0.71由Tr，pr查表得Z=0.86因为pV=ZRT，则V=ZRT/p=0.86×8.314×366.48/3.446×106=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ρ+Cρ2(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）在等温下积分，当p→0，ρ→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都会有类似的情况，但也有例外，在汽化过程时，等温和等压条件下，两相中的G值是不变的，即按相平衡的判据理应如此，纯物质的汽化乃是单元组分的相平衡，无疑可写成GL=GV式中，GL和GV分别代表液相和蒸汽相（也可称为气相）的摩尔自由焓。在平衡系统中要特别关注G的作用，同时也要充分理解G之所以会不变的原因。.3.6.1Clapeyron方程根据式（3-87）做如下推导。按G定义，可写出GV=HV-TSV和GL=HL-TSL在平衡条件下，上两式相减，得△GLV=△HLV-T△SLV（3088）式中，△GLV、△HLV和△SLV分别为汽化自由焓、汽化焓和汽化熵。.3.6.1Clapeyron方程此外，还应示出△SLV和△VLV的关系。对平衡态时的微分变化，借助式（3-4）可写出dGV=-SVdT+VLdp和dGL=-SLdT+VLdp上两式相减得d△GLV=-△SLVdT+△VLVdp（3-89）.3.6.1Clapeyron方程根据式（3-87），对单元系气掖平衡下的任何状态，可得d△GLV=0（3-90）结合式（3-89）和式（3-90）得（3-91）式中ps=p因△GLV=0，则由式（3-88）得△HLV=T△SLV故式（3-91）也可写成（3-92）式（3-92）即为著名的Clapeyron方程.3.6.1Clapeyron方程汽化过程的相变体积变化，△VLV可用△ZLV=ZV-ZL来表示，因为，在饱和条件下，汽化过程的压缩因子变化，△ZLV可写为（3-93）式（3-93）与式（3-92）结合可写成（3-94）式（3-92）和式（3-94）都是Clapeyron方程式的表达形式，可用来从蒸汽压和体积数据求算汽化焓，或从汽化焓与体积数据来酸蒸汽压。.3.6.1Clapeyron方程若将设为常数，切与温度无关，则式（3-94）便可积分化简后得（3-95）式中A为积分常数，B=△HLV/R△ZLV。因VV＞＞VL，故△VLV=VV，则△ZLV=ZV，若再假设蒸汽具有理想气体行为，则ZV=1。将这些条件用于式（3-94），化简后得（3-96）因此式（3-95）与式（3-96）都是近似式，是Clausius-Clapeyron方程式的不同表现形式，当然Clausius-Clapeyron方程只适用与低压情况。.