化工热力学马沛生第二版2-3章习题答案

第二章流体的p-V-T关系1习题：2－1．为什么要研究流体的pVT关系？答：在化工过程的分析、研究与设计中，流体的压力p、体积V和温度T是流体最基本的性质之一，并且是可以通过实验直接测量的。而许多其它的热力学性质如内能U、熵S、Gibbs自由能G等都不方便直接测量，它们需要利用流体的p–V–T数据和热力学基本关系式进行推算；此外，还有一些概念如逸度等也通过p–V–T数据和热力学基本关系式进行计算。因此，流体的p–V–T关系的研究是一项重要的基础工作。2－2．理想气体的特征是什么？答：假定分子的大小如同几何点一样，分子间不存在相互作用力，由这样的分子组成的气体叫做理想气体。严格地说，理想气体是不存在的，在极低的压力下，真实气体是非常接近理想气体的，可以当作理想气体处理，以便简化问题。理想气体状态方程是最简单的状态方程：RTpV2－3．偏心因子的概念是什么？为什么要提出这个概念？它可以直接测量吗？答：纯物质的偏心因子是根据物质的蒸气压来定义的。实验发现，纯态流体对比饱和蒸气压的对数与对比温度的倒数呈近似直线关系，即符合：rsrTp11log其中，cssrppp对于不同的流体，具有不同的值。但Pitzer发现，简单流体（氩、氪、氙）的所有蒸气压数据落在了同一条直线上，而且该直线通过rT=0.7，1logsrp这一点。对于给定流体对比蒸气压曲线的位置，能够用在rT=0.7的流体与氩、氪、氙（简单球形分子）的srplog值之差来表征。Pitzer把这一差值定义为偏心因子，即)7.0(00.1logrsrTp任何流体的值都不是直接测量的，均由该流体的临界温度cT、临界压力cp值及rT=0.7时的饱和蒸气压sp来确定。2－4．纯物质的饱和液体的摩尔体积随着温度升高而增大，饱和蒸气的摩尔体积随着温度的升高而减小吗？答：正确。由纯物质的p–V图上的饱和蒸气和饱和液体曲线可知。2－5．同一温度下，纯物质的饱和液体与饱和蒸气的热力学性质均不同吗？答：同一温度下，纯物质的饱和液体与饱和蒸气的Gibbs自由能是相同的，这是纯物质气液第二章流体的p-V-T关系2平衡准则。气他的热力学性质均不同。2－6．常用的三参数的对应状态原理有哪几种？答：常用的三参数对比态原理有两种，一种是以临界压缩因子Zc为第三参数；另外一种是以Pitzer提出的以偏心因子作为第三参数的对应状态原理。2－7．总结纯气体和纯液体pVT计算的异同。答：许多p–V-T关系如RKS方程、PR方程及BWR方程既可以用于计算气体的p–V–T，又都可以用到液相区，由这些方程解出的最小体积根即为液体的摩尔体积。当然，还有许多状态方程只能较好地说明气体的p–V-T关系，不适用于液体，当应用到液相区时会产生较大的误差。与气体相比，液体的摩尔体积容易测定。除临界区外，温度（特别是压力）对液体容积性质的影响不大。除状态方程外，工程上还常常选用经验关系式和普遍化关系式等方法来估算。2－8．简述对应状态原理。答：对比态原理认为，在相同的对比状态下，所有的物质表现出相同的性质。对比态原理是从适用于p–V-T关系两参数对比态原理开始的，后来又发展了适用于许多热力学性质和传递性质的三参数和更多参数的对比态原理。2－9．如何理解混合规则？为什么要提出这个概念？有哪些类型的混合规则？答：对于真实流体，由于组分的非理想性及由于混合引起的非理想性，使得理想的分压定律和分体积定律无法准确地描述流体混合物的p–V-T关系。如何将适用于纯物质的状态方程扩展到真实流体混合物是化工热力学中的一个热点问题。目前广泛采用的方法是将状态方程中的常数项，表示成组成x以及纯物质参数项的函数，这种函数关系称作为混合规则。对于不同的状态方程，有不同的混合规则。寻找适当的混合规则，计算状态方程中的常数项，使其能准确地描述真实流体混合物的p–V-T关系，常常是计算混合热力学性质的关键。常用的混合规则包括适用于压缩因子图的虚拟临界性质的混合规则、维里系数的混合规则以及适用于立方型状态方程的混合规则。2－10．在一个刚性的容器中，装入了1mol的某一纯物质，容器的体积正好等于该物质的摩尔临界体积Vc。如果使其加热，并沿着习题图2－1的p–T图中的1→C→2的途径变化（C是临界点）。请将该变化过程表示在p–V图上，并描述在加热过程中各点的状态和现象。解：由于加热过程是等容过程，1→C→2是一条CVV的等容线，所以在p–V图可以表示为如图的形式。点1表示容器中所装的是该物质的汽液混合物（由饱和蒸汽和饱和液体组成）。沿1－2线，是表示等容加热过程。随着过程的进行，容器中的饱和液体体积与饱和蒸汽体积的相对比例有所变化，但由图可知变化不是很大。到了临界点C点时，汽液相界面逐渐消失。继续加热，容器中一直是均相的超临界流体。在整个过程中，容器内的压力是不断增加的。第二章流体的p-V-T关系3C121···212－11．已知SO2在431K下，第二、第三Virial系数分别为：13kmolm159.0B，263kmolm100.9C，试计算：（1）SO2在431K、10×105Pa下的摩尔体积；（2）在封闭系统内，将1kmolSO2由10×105Pa恒温（431K）可逆压缩到75×105Pa时所作的功。解：（1）三项维里方程为：21VCVBRTpVZ（A）将p=10×105Pa，T=431K，13kmolm159.0B，263kmolm100.9C代入式（A）并整理得：0109159.0279.0623VVV迭代求解，初值为：13kmolm5.3pRTV迭代结果为：13kmolm39.3V（2）压缩功VpWd由（A）式得：321VCVBVRTp，则：212212123211211lnd121VVCVVBVVRTVVCVBVRTWVV（B）当p=75×105Pa时，用（1）同样的方法解出：第二章流体的p-V-T关系413kmolm212.0V将131kmolm39.3V，132kmolm212.0V代入式（B）解出：15kmolJ1077W2－12．试计算一个125cm3的刚性容器，在50℃和18.745MPa的条件下能贮存甲烷多少克（实验值为17g）？分别用理想气体方程和RK方程计算（RK方程可以用软件计算）。解：由附录三查得甲烷的临界参数为：cT＝190.56K，cp＝4.599MPa，＝0.011（1）利用理想气体状态方程RTpV得：131346molcm3.143molm10433.110745.185015.273314.8pRTVg95.133.14312516总VVMm（2）RK方程)(5.0bVVTabVRTp式中：2-0.5665.225.22molKmPa2207.310599.456.190314.842748.0/42748.0＝＝ccpTRa1356molm10985.210599.456.190314.808664.0/08664.0＝＝ccpRTb4653.015.323314.810745.182207.35.2265.22＝＝TRapA2083.015.323314.810745.1810985.265＝＝RTbpB按照式（2-16a）hhhhhBAhZ12342.211111＝和式（2－16b）ZZBVbh2083.0迭代计算，取初值Z＝1，迭代过程和结果见下表。迭代次数Zh第二章流体的p-V-T关系5010.208310.87790.237320.88260.236030.88230.236140.88230.236113346mol126.5cm/molm10265.110745.1815.323314.88823.0＝pZRTVg81.155.12612516总VVMm可见，用RK方程计算更接近实验值。2－13．欲在一个7810cm3的钢瓶中装入1kg的丙烷，且在253.2℃下工作，若钢瓶的安全工作压力为10MPa，问是否安全？解：查得丙烷的临界性质为：cT＝369.83K，cp＝4.248MPa，＝0.152mol727.22441000Mmn1366molm1063.343727.22107810nVV总使用RK方程：)(5.0bVVTabVRTp首先用下式计算a，b：2-0.5665.225.22molKmPa296.1810248.483.369314.842748.0/42748.0ccpTRa1356102771.610248.483.369314.808664.0/08664.0molmpRTbcc代入RK方程得：MPap870.9非常接近于10MPa，故有一定危险。2－14．试用RKS方程计算异丁烷在300K，3.704×105Pa时的饱和蒸气的摩尔体积。已知实验值为133molm10081.6V。解：由附录三查得异丁烷的临界参数为：cT＝407.8K，cp＝3.640MPa，＝0.1777357.08.407/300/crTTT7531.0177.0176.0177.0574.1480.0176.0574.1480.022m第二章流体的p-V-T关系62258.17357.017531.01)1(1)(25.025.0rTmT2662222mol/mPa6548.12258.110640.38.407314.842748.0/4278.0＝＝TpTRTaTacc/molm100700.810640.3/8.407314.808664.0/08664.0356＝ccpRTb09853.0300314.810704.36548.122522＝＝TRapA01198.0300314.810704.3100700.855＝＝RTbpB按照式（2-16a）hhhhhBAhZ12245.811111＝和式（2－16b）ZZBVbh01198.0迭代计算，取初值Z＝1，迭代过程和结果见下表。迭代次数Zh010.0119810.91480.0131020.90700.0132130.90620.0132240.90610.0132250.90610.01322/molm101015.610704.3300314.89061.0326pZRTV误差%2.110031.6/101015.6031.6222－15．试分别用RK方程及RKS方程计算在273K、1000×105Pa下，氮的压缩因子值，已知实验值为Z＝2.0685。解：由附录三查得氮的临界参数为：cT＝126.10K，cp＝3.394MPa，＝0.040（1）RK方程2-0.5665.225.22molKmPa5546.110394.310.126314.842748.0/42748.0＝＝ccpTRa1356molm106763.210394.310.126314.808664.0/08664.0＝＝ccpRTb第二章流体的p-V-T关系78264.1273314.8101005546.15.2265.22＝＝TRapA1791.1273314.8101000106763.255＝＝RTbpB按照式（2-16a）