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

第二章流体的p-V-T关系习题：2－1．为什么要研究流体的pVT关系？与设计中，流体的压力p、体积V和温度T是流体最基本的性热力学性质如内能U、熵S、Gibbsp–V–T数据和热力学基本关系式进行p–V–T数据和热力学基本关系基础工作。答：在化工过程的分析、研究质之一，并且是可以通过实验直接测量的。而许多其它的自由能G等都不方便直接测量，它们需要利用流体的推算；此外，还有一些概念如逸度等也通过式进行计算。因此，流体的p–V–T关系的研究是一项重要的2－2．理想气体的特征是什么？答：假定分子的大小如同几何点一样，分子间不存在相互作用力，由这样的分子组成的气体严格地说，理想气体是不存在的，在极低的压力下，真实气体是非常接近理想气体的，可以当作理想气体处理，以便简化问题。理想气体状态方程是最简单的状态方程：pVRT2－3．偏心因子的概念是什么？为什么要提出这个概念？它可以直接测量吗？蒸气压来定义的。实验发现，纯态流体对比饱和蒸倒数呈近似直线关系，即符合：答：纯物质的偏心因子是根据物质的气压的对数与对比温度的1其中，pspsprlogps1Trrc对于不同的流体，具有不同的值。但Pitzer发现，简单流体（氩、氪、氙）的所有蒸气压Tlogps1这一点。对于给定流体数据落在了同一条直线上，而且该直线通过=0.7，对rrTlogps值r比蒸气压曲线的位置，能够用在=0.7的流体与氩、氪、氙（简单球形分子）的r之差来表征。Pitzer把这一差值定义为偏心因子，即logps1.00r(T0.7)r任何流体的值都不是直接测量的Tp值及c，均由该流体的临界温度、临界压力cTr=0.7时的饱和蒸气压ps来确定。2－4．纯物质的升高而减小吗？答：正确。由纯物质的p–V图上的饱和蒸气和饱和液体曲线可知。2－5．同一温度下，纯物质的饱和液体与饱和蒸气的热力学性质均不同吗？答：同一温度下，纯物质的饱和液体与饱和蒸气的Gibbs自由能是相同的，这是纯物质气液饱和液体的摩尔体积随着温度升高而增大，饱和蒸气的摩尔体积随着温度的1第二章流体的p-V-T关系平衡准则。气他的热力学性质均不同。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是一条VV方对应状态原理。出相同的性质。适用于许多热的混合规则？理想的分压定律扩展热力学性质的关则包括变化，并描述在加热过程中各点的状态和现象。的等容线，所以在p–V图可以表示C为如图的形式成）。沿1－2线，是表示等汽体积的相对比例有所变化，但由图可知变化不是很大。到了临界点C点时，汽液相界面逐渐消失。继续加热，容器中一直是均相的超临界流体。在整个过程中，容器内的压力是不断增加的。。点1表示容器中所装的是该物质的汽液混合物（由饱和蒸汽和饱和液体组容加热过程。随着过程的进行，容器中的饱和液体体积与饱和蒸2第二章流体的p-V-T关系·22·C·12－11．已知SO在431K下，第二、第三Virial系数分别为：B0.159m3kmol1，2C9.0103m6kmol2，试计算：（1）SO在431K、10×10Pa下的摩尔体积；52（2）在封闭系统内，将1kmolSO由10×10Pa恒温（431K）可逆压缩到75×10Pa552时所作的功。解：（1）三项维里方程为：BCZpV1（A）RTVV2将p=10×105Pa，T=431K，B0.159m3kmol1，C9.0103m6kmol2代入式（A）并整理得：0.279V3V20.159V91060迭代求解，初值为：VRT3.5m3kmol1p迭代结果为：V3.39m3kmol1（2）压缩功WpdV由（A）式得：pRT1，则：BCVV2V31BCWRTdVV2VV2V1V3（B）C1RTlnV2B111V2V2VV21V2121当p=75×10Pa时，用（1）同样的方法解出：53第二章流体的p-V-T关系V0.212m3kmol1V3.39m3kmol1V0.212m3kmol1代入式（B）解出：将1，2W77105Jkmol12－12．试计算一个125cm3的刚性容器，在50℃和18.745MPa的条件下能贮存甲烷多少克（实验值为17g）？分别用理想气体方程和RK方程计算（RK方程可以用软件计算）。T由附录三查得甲烷的临界参数为：＝190.56K，＝4.599MPa，＝0.011cpc解：（1）利用理想气体状态方程pVRT得：VRT8.314273.15501.433104m3mol1143.3cm3mol118.745106pmMV总1612513.95gV143.3（2）RK方程RTVbT0.5V(Vb)ap式中：0.427488.314190.5622.5＝3.2207Pam6K0.5mol-2a0.42748R2T2.5/p＝4.599106ccb0.08664RT/p＝0.086648.314190.56＝2.985105m3mol14.599106ccap＝3.220718.745106A＝0.4653R2T8.314323.1522.52.5Bbp＝2.98510518.745106＝0.2083RT8.314323.151Ah1h1hZ＝2.2342按照式（2-16a）1hB1h1hbB0.2083h2－16b）VZ和式（Z迭代计算，取初值Z＝1，迭代过程和结果见下表。迭代次数Zh4第二章流体的p-V-T关系0123410.20830.23730.23600.23610.23610.87790.88260.88230.8823VZRT0.88238.314323.151.265104m3/mol＝126.5cm3mol118.745106pmMV总1612515.81gV126.5可见，用RK方程计算更接近实验值。2－13．欲在一个7810cm3的钢瓶中装入1kg的丙烷，且在253.2℃下工作，若钢瓶的安全工作压力为10MPa，问是否安全？解：查得丙烷的临界性质为：cT＝369.83K，p＝4.248MPa，＝0.152cm1000n22.727mol44M7810106343.63106m3mol1VV总n22.727RTVbT0.5V(Vb)a使用RK方程：p首先用下式计算a，b：a0.42748R2T2.5/p0.427488.314369.8322.518.296Pam6K0.5mol-24.24810cc6b0.08664RT/p0.086648.314369.836.2771105m3mol14.248106cc代入RK方程得：p9.870MPa非常接近于10MPa，故有一定危险。2－14．试用RKS方程计算异丁烷在300K，3.704×105Pa时的饱和蒸气的摩尔体积。已知实验值为V6.081103m3mol1。由附录三查得异丁烷的临界参数为：T＝407.8K，p＝3.640MPa，＝0.177cc解：TT/T300/407.80.7357rcm0.4801.5740.17620.4801.5740.1770.1760.17720.75315第二章流体的p-V-T关系1.22582(T)1m(1T)10.753110.735720.50.5r8.3142407.823.640106aTaT0.4278R2T2/pT＝0.427481.2258＝1.6548Pam/mol26ccb0.08664RT/p＝0.086648.314407.8/3.640108.0700105m3/mol6ccap＝1.65483.704105A＝0.09853R2T8.314300222Bbp＝8.07001053.704105＝0.01198RT8.3143001Ah1hZ按照式（2-16a）1hB1h1h8.22451h＝bB0.01198h2－16b）VZ和式（Z迭代计算，取初值Z＝1，迭代过程和结果见下表。迭代次数Zh01234510.011980.013100.013210.013220.013220.013220.91480.90700.90620.90610.9061VZRT0.90618.3143006.1015102m3/mol3.704106p6.0316.1015102/6.0311021.2%误差2－15．试分别用RK方程及RKS方程计算在273K、1000×105Pa下，氮的压缩因子值，已知实验值为Z＝2.0685。T临界参数为：＝126.10K，＝3.394MPa，＝0.040cp解：由附录三查得氮的c（1）RK方程0.427488.314126.1022.5＝1.5546Pam6K0.5mol-2a0.42748R2T2.5/p＝3.394106cc0.086648.314126.10＝2.6763105m3mol1b0.08664RT/p＝3.394106cc6第二章流体的p-V-T关系ap1.5546100106＝1.8264A＝R2T2.58.31427322.5Bbp＝2.67631051000105＝1.1791RT8.3142731Ah1h＝1.54891h按照式（2-16a）Z1hB1h1hbB1.1791h2－16b）VZ和式（Z迭代计算，取初值Z＝2，迭代过程和结果见下表。迭代次数Zh020.5895511.8622.12601.69260.88230.63320.55460.69660.2361234……..迭代不收敛，采用RK方程解三次方程得：V=0.00004422m3/molZpV4.42210510001051.9485RT8.314273RKS方程TT/T273/126.12.1649rcm0.4801.5740.17620.4801.5740.0400.1760.04020.54270.55382(T)1m(1T)10.542712.164920.50.5r8.314126.122aTaT0.4278R2T2/pT＝0.427480.5538＝0.076667Pam/mol63.394106ccb0.08664RT/p＝0.086648.314126.1/3.394102.6763105m3/mol6ccap0.0766671000105＝1.48828.314273A＝R2T2227第二章流体的p-V-T关系bp2.67631051000105＝1.1791B＝RT8.3142731Ah1h＝1.26211hZ按照式（2-16a）1hB1h1hbB1.17912－16b）h和式（VZZ同样迭代不收敛采用RKS方程解三次方程得：V=0.00004512m3/molZpV4.51210510001051.9881RT8.3142732－16．试用下列各种方法计算水蒸气在107.9×105Pa、593K下的比容，并与水蒸气表查出的数据（V0.01687m3kg1）进行比较。（1）理想气体定律（2）维里方程（3）普遍化RK方程解：从附录三中查得水的临界参数为：cT＝647.13K，p＝22.055MPa，＝0.345c（1）理想气体定律VRT8.314593107.91054.569106m3mol10.02538m3kg1p0.016870.02538100%50.5%误差=0.01687（2）维里方程593TT0.916T647.13rcp107.9105pr0.48922.055106pc使用普遍化的第二维里系数：B(0)0.0830.422/T1.60.0830.4220.4026T1.6rrB(1)0.1390.172/T4.20.1390.1720.1096T4.2rr8第二章流体的p-V-T关系BpRTcB(0)B(1)0.40260.3450.10960.4404cBppc0.44040.76490.4890.916Z1Bp1RT1rRTTcrVZRT0.76498.3145933.495106m3mol10.01942m3kg1107.9105p0.016870.01942100%15.1%=误差0.01687（3）普遍化R-K方程a1h1hZ1hT（2－38a）（2－38b）1.5brhpbZTrr将对比温度和对比压力值代入并整理的：a1h1h1h115.628Z1hT1h1.5brhr0.04625pbZTZr联立上述两式迭代求解得：Z=0.7335VZRT0.73358.3145933.3515106m3mol10.01862m3kg1107.9105p0.016870.01862100%10.4%=误差0.01687水是极性较强的物质2－17．试分别用（1）vanderWaals方程；（2）RK方程；（3）RKS方程计算273.15K时将CO2压缩到体积为550.1cm3mol1所需要的压力。实验值为3.090MPa。解：从附录三中查得CO2的临界参数为：T＝304.19K，p＝7.382MPa，＝0.228cc（1）vanderWaals方程RTVbVa2pa27R2T2/64p＝278.314304.1922＝0.3655Pam3mol2式中：647.382106cc9第二章流体的p-V-T关系bRT/8p8.314304.194.282105m3mol187.382106cc则：RTVba8.314273.15550.1106－42.821060.3655－＝3.269106Pa＝3.2Mp＝550.11026V23.0903.269100％＝5.79％误差％＝3.090（2）RK方程RTVbT0.5V(Vb)ap式中：a0.42748R2T2.5/p＝0.427488.314304.1922.5＝6.4599Pam6K0.5mol-27.382106ccb0.08664RT/p＝0.086648.314304.19＝2.968105m3mol17.382106ccRTVbapT0.5V(Vb)8.314273.156.4599＝－550.1106－29.68106＝3.138106Pa＝3.138MPa3.0903.138100％＝1.55％273.150.5550.1106550.1＋29.68106误差％＝3.090（3）RKS方程RTaTpVbVVbaTaT0.4278RT2/pT式中，2cc(T)1m(1T)20.5r而，m0.4801.5740.1762＝0.4801.5740.228－0.1760.2282＝0.8297＝10.8297122＝1.089273.150.5304.19则，TmT0.5)()1(1r10第二章流体的p-V-T关系7.3821060.427488.314304.1921.0892aTaT0.42748R2T2/pT＝cc＝0.40335Pam3mol1b0.08664RT/p＝0.086648.314304.19＝2.968105m3mol17.382106ccRTaT8.314273.150.40335p＝－550.1106－29.68106550.1106550.1＋29.68106VbVVb＝3.099106Pa＝3.099MPa3.0903.099100％＝0.291％误差％＝3.090比较几种方程的计算结果，可见，vanderWaals方程的计算误差最大，RKS方程的计算精度RK方程的计算精度还可以。2－18．一个体积为0.3m3的封闭储槽内贮乙烷，温度为290K、压力为5烷加热到479K，试估算压力将变为多少？最好。25×10Pa，若将乙解：乙烷的临界参数和偏心因子为c：T＝305.32K，p＝4.872MPa，＝0.099c因此：TT/T290/305.320.95cpp/p2.5/48.720.513cr11r11故使用图2-11，应该使用普遍化第二维里系数计算0.422B(0)0.0830.422/T1.60.083r0.3750.0740.951.60.172B(1)0.1390.172/T4.20.139r0.954.2r0.5130.95Z1Bp1BB1p0.3750.0990.0740.793501RTTrVZRT0.79358.31429076.5105mmol3125105p392.2mol0.3nV总V76.5105加热后，采用RK方程进行计算。其中：T=479K，摩尔体积仍然为V76.5105m3mol1，首先计算：a0.42748R2T2.5/p＝0.427488.314305.3222.5＝9.879Pam6K0.5mol-24.872106cc11第二章流体的p-V-T关系b0.08664RT/p＝0.086648.314305.32＝4.514105mmol314.872106cc代入RK方程：RTVbaT0.5V(Vb)p＝8.31447969.879－765.010－45.141064790.5765.0106765.045.14106＋＝4.804106Pa＝4.804MPa2－19．如果希望将22.7kg的乙烯在294K时装入0.085m3的钢瓶中，问压力应为多少？解：从附录三查得乙烯的临界参数为：cT＝282.34K，p＝5.041MPa，＝0.085cm22.7103810.7molnM280.085104.8106mmol31n810.7VV摩尔体积总采用RK方程进行计算。首先计算：a0.42748R2T2.5/p＝0.427488.314282.3422.5＝7.851Pam6K0.5mol-25.041106ccb0.08664RT/p＝0.086648.314282.34＝4.034105mmol315.041106cc代入RK方程：RTVbaT0.5V(Vb)p＝8.31429467.851－104.810－40.341062940.5104.8106104.840.34106＋＝7.817106Pa＝7.817MPa2-20（由于较简单省略了，忽略不计了）2－21．用Pitzer的普遍化关系式计算甲烷在323.16K时产生的压力。已知甲烷的摩尔体积mmol13为1.25×10－4，压力的实验值为1.875×107Pa。解：从附录三查得甲烷的临界参数为：cT＝190.56K，p＝4.599MPa，＝0.011cTT/T323.16/190.561.696；但是p不能直接计算，需要试差计算rcr2.14910pZRT/VZ8.314323.16/1.251047Z12第二章流体的p-V-T关系并且ppp4.599106pcrr4.599106p2.149107Z因此，结合上两式得：0.214pr（A）（B）rZZ(0)Z(1)Pitzer的普遍化关系式为：根据（A）、（B）两式进行迭代，过程为：p（1）设Z值，然后代入（A）式求出；rTpZZ和（2）根据和值查（2－9）和（2－10）得到01；rr1值代入（ZZ和（3）将查图得到的0B）式求得Z值；（4）比较Z的计算值与实验值，如果相差较大，则代入（A）式重新计算，直到迭代收敛。p依据上述迭代结果为：＝4.06时，Z＝0.877rppp4.599106p＝4.5991064.06＝1.867107Pa则：crr1.8751.867107/1.8751070.43%误差：2－22．试用RK方程计算二氧化碳和丙烷的等分子混合物在151℃和13.78MPa下的摩尔体积。解：计算所需的数据列表如下：组元T/Kp10/PaV106/(m3mol-1)Zca/(Pam)63K0.5mol-2)b/(mmol-15cccCO2（1）304.2C3H8（2）369.873.8242.4854.7294.02000.2746.4600.27718.2911.122.968×10-56.271×10-512335.4140.4由（2－51a）和（2－51b）得：ay2a2yyay2a＝0.526.46020.50.511.120.5218.29M1111212222＝11.75Pam6K0.5mol-2bybyb＝0.52.9681050.56.271105＝4.415105mmol3-1M1122ap＝11.7513.78106A＝0.6322MRT8.314424.1522.522.513第二章流体的p-V-T关系Bbp＝4.41510513.78106＝0.1725MRT8.314424.151Ah1hZ1hB1h1h3.60841h＝按照式（2-16a）（A）（B）bB0.1725h2－16b）V和式（ZZm联立求解方程（A）、）（B）进行迭代计算得：迭代次数Zh01234510.17250.25460.28310.28810.28870.28870.67760.60930.59870.59760.5975因此：Z＝0.5975，h＝0.2887ZRT0.59758.314424.151.52104m3molV混合物得摩尔体积为：13.78106p2－23．混合工质的性质是人们有兴趣的研究课题。试用RKS状态方程计算由R12（CClF）22和R22（CHClF）组成的等摩尔混合工质气体在400K和1.0MPa，2.0MPa，3.0MPa，4.0MPa2和5.0MPa时的摩尔体积。可以认为该二元混合物的相互作用参数k＝0（建议自编软件计12算）。计算中所使用的临界参数如下表组元（i）Tcp/MPac/KR22（1）R12（2）369.23854.9754.2240.2150.176解：计算过程是先计算两个纯组分的RKS常数，再由混合规则获得混合物的RKS常数a,bM后，可以进行迭代计算，也可以求解三次方程的体积根。M建议大家自编程序进行计算。所得的结果列于下表：T/K400组成yy0.512RKS方程常数组分（1）：a＝0.7568b＝5.346×10-5组分（2）：a＝1.007b＝6.565×10-514第二章流体的p-V-T关系混合物a＝0.8774b＝5.956×10-5p/MPaVcal/（cm12345·mol）3114.01442.3877.0585.5399.33-12－24．试用下列方法计算由30％（摩尔％）的氮（1）和70％正丁烷（2）所组成的二元混合物，在462K、69×10Pa下的摩尔体积。5（1）使用Pitzer三参数压缩因子关联式（2）使用RK方程，其中参数项为：b0.086640RTcipicia0.427480R2T2.5cijpijcij（3）使用三项维里方程，维里系数实验值为B14106，B265106，，1122B9.5106，（B的单位为m3mol1）。C1.3109，C3.02510911122212C4.95109，C7.27109（C的单位为m6mol2）。已知氮及正丁烷的临界112122参数和偏心因子为N2T＝126.10K，p＝3.394MPa，＝0.040ccT＝425.12K，p＝3.796MPa，＝0.199ccnC4H10解：（1）根据Kay规则求出混合物的虚拟临界参数TyT0.3126.100.7425.12335.41Kpciciipyp0.33.3940.73.7963.675MPapciciiy0.30.0400.70.1990.1513iiiT4621.377335.41虚拟对比条件为：TprTpcp6.91.8773.675pprppc查图2－9和2－10得：Z00.77,Z10.19则：ZZ(0)Z(1)0.770.15130.190.798715第二章流体的p-V-T关系VZRT0.79878.3144624.446104mmol316.9106p（2）RK方程1ah1hZ1hbRT1.5hbbpVZRTay2a2yyay2a2211112122bybyb11220.42748R2T2.5c12a12pc12组元T/Kcp105/PaV106/(m3mol-1)Zca/(Pam)63K0.5mol-2)b/(mmol-1cc112212126.1033.94425.1237.96231.5334.3790.12550.2921.5550.27429.010.2837.0122.676×10-58.067×10-5158.5a0.321.55520.30.77.0120.7229.0117.30Pam6K0.5mol-2Mb0.32.6761050.78.0671056.450105m3mol-1M1ahbRT1.5117.301h1h6.4501058.314462h1hZ1h1.51h1h1h3.25bp6.450105691050.1159hbVZRTZ8.314462Z进行试差迭代得：Z0.746,h=0.156VZRT0.7468.3144624.15104mmol316.9106p（3）三项的维里方程为：ZpV1BCRTVV216第二章流体的p-V-T关系By2B2yyBy2B0.32(14)20.30.7(9.5)0.72(265)106M11112122221.326104(m3mol-1)CMyyyCy3C3yyC3yy2Cy3C2ijkijk111112112121222222ijk0.33(1300)30.320.7(4950)30.30.7272700.73(30250)10121.455103(m6mol-2)将以上结果代入三项维里方程得：69105V1.3261041.455103VV28.3144621试差求解得：V4.25104m3mol12－25．一压缩机，每小时处理454kg甲烷及乙烷的等摩尔混合物。气体在50×105Pa、422K积流率为多少cm3h1下离开压缩机，试问离开压缩机的气体体？解：混合物的分子量为M05M0.5M0.516.040.530.0723.06C2H6CH4454混合物的流n率为：19.7kmolh-123.06利用Kay规则求虚拟临界常数：TyT0.5190.560.5305.32247.94Kpciciipyp0.54.5990.54.8724.736MPapciciiT4221.702247.94T虚拟对比条件为：Tprpcp5ppr1.056p4.736pc用图2-11判断，应该使用维里方程，现将所需数据列于下表，其中第三行数据按照（2-48a）~（2-48e）式计算。ijijTcijV/(m3kmol-1)Zcijcijp/K/MPacij112212190.56305.32241.214.5994.8724.7010.098600.14550.12050.2860.2790.28250.0110.0990.055采用二阶舍项的virial方程计算混合物的性质，需要计算混合物的交互第二virial系数，计算结果见下表，17第二章流体的p-V-T关系ijB(0)B(1)B/(m3kmol-1)ij112212-0.0353-0.1680.133－0.011650.0948－0.082870.1226－0.03528-0.0894由式（2－46）得：By2B2yyBy2B0.52(0.01165)20.50.5(0.03528)0.52(0.08287)M11112122220.04127(m3kmol-1)VRTB8.3141034220.041270.6604103(m3kmol-1)5.0106p体积流率nV19.70.660413.01mh312－26．H2和N2的混合物，按合成氨反应的化学计量比，加入到反应器中N3H2NH322mh1。其中15％的N2混合物进反应器的压力为600×10Pa，温度为298K，流率为635转化为NH3，离开反应器的气体被分离后，未反应的气体循环使用，试计算：（1）每小时生成多少公斤NH3?（2）若反应器出口物流（含NH的混合物）的压力为550×10Pa、温度为451K，53试问在内径D＝0.05m管内的流速为多少？1）这是一个二元混合物系pVT的计算问题。使用RK方程进行计算解：（1ahZ1hbRT1.51hhbbpVZRTay2a2yyay2a2211112122bybyb11220.42748R2T2.5c12a12pc12组元p105/PaV106/(m3mol-1)ZT/Kca/(Pam)63K0.5mol-2)b/(mmol-1ccc11（N2）126.1033.9490.10.2921.5552.676×10-518第二章流体的p-V-T关系22（H2）33.1813.1321.0364.20.3050.14270.2990.47271.820×10-51264.6876.45a0.2521.55520.250.750.47270.7520.14270.3547Pam6K0.5mol-2Mb0.252.6761050.751.8201052.034105m3mol-1M1ah11h1h2.0341058.3142980.3547h1hZ1hbRT1.51.51h1h1h0.4077bp2.0341056001050.4926hbVZRTZ8.314298Z进行试差迭代得：Z1.424h=0.346VZRT1.4248.31429858.8106mmol3160106p1.020105molh16106n摩尔流率58.8N的摩尔流率为：n0.251.020105＝2.551104molh12N2生成的NH量为：2.5511040.152＝7.653103molh－13（2）这是一个三元混合物系pVT的计算问题。继续使用RK方程进行计算反应器出口物流组成：以入口1molN2为基准N2：1—0.15=0.85H2：3—3×0.15＝2.55NH3：0.15×2＝0.30则总物质的量为：0.85+0.30+2.55＝3.75各物质的摩尔分率为：y0.850.2303.75N2y2.550.6893.75H2y0.300.0813.75NH3以NH3作为第三组元，补充数据如下：组元p105/PaV106/(m3mol-1)ZT/Kca/(Pam)63K0.5mol-2)b/(mmol-1ccc19第二章流体的p-V-T关系33（NH）405.65112.7872.581.068.30.2428.6830.2673.6660.2741.1092.591×10-531323226.1762.0116.0138.6ay2ay2ay2a2yya2yya2yya23M10.2321.5550.68920.14270.08128.68320.230.6890.472711222333121213132320.230.0813.66620.0810.6891.1090.617Pam6K0.5mol-2bybybyb0.232.6760.6891.8200.0812.591105M1122332.079105mmol3-11ah1h1h10.6172.0791058.314451h1hZ1hbRT1.51.510.3727h1h1hbp2.0791055501050.3050ZhbVZRTZ8.314451进行试差迭代得：Z1.250h=0.244VZRT1.2508.31445185.22106mmol3155106p所以以进口N2为1mol作基准入口总物质的量为：1＋3+0＝4mol出口总物质的量为：1×（1－0.15）＋3×（1－0.15）+1×0.15×2＝3.7mol率为：反应物摩尔流率×3.7/4＝1.022×10mmol产品的摩尔流5×3/4=9.45×1043185.221069.451048.053mh1产品的体积流率为：34103.4mh1.140ms118.0538.053速率uVA3.14/40.052D242－27．测得天然气（摩尔组成为CH84%、N9%、CH7%）在压力9.27MPa、温度37.8℃4226下的气体流速。平均时速为25m3h1。试用下述方法计算在标准状况下的（1）理想气体方程；（2）虚拟临界参数；（3）Dalton定律和普遍化压缩因子图；20第二章流体的p-V-T关系（4）Amagat定律和普遍化压缩因子图。解：（1）按理想气体状态方程；标准状况下气体流速9.272527337.82730.1013pvTv（273K，0.1013MPa）＝112010m3h12Tp12（2）虚拟临界参数法首先使用Kay规则求出虚拟的临界温度和临界压力，计算结果列表如下：组分摩尔/％T/Kcp/MPacyT/Kcyp/MPac甲烷氮气乙烷合计0.840.090.071.00190.56126.10305.324.5993.3944.872160.0711.353.8630.3050.3414.51021.37192.79虚拟临界温度为192.79K，压力为4.510MPa，混合物的平均压缩因子可由下列对比温度和对比压力求出：T37.82731.61，p9.272.055192.794.510rr查两参数普遍化压缩因子图得：Zm＝0.89将压缩因子代入方程pVZRT得：2.481104mmolVZRT0.898.31437.8273319.27106p25nvV1.008105molh1100.8kmolh12.481104在标准状态下，压缩因子Z＝1，因此体积流率可以得到：8.3142732258.5m3h1vnVnRT100.81030.1013106p（3）Dalton定律和普遍化压缩因子查普遍化压缩因子图时，各物质的压力使用分压组分ppyZyZiiTTppTrpriiicc甲烷氮气乙烷1.632.461.0287.7870.8340.6491.6930.2460.1330.900.980.960.7560.08820.067221第二章流体的p-V-T关系合计0.9114将压缩因子代入方程pVZRT得：2.541104mmolVZRT0.91148.31437.8273319.27106pv25nV2.5411049.839104molh198.39kmolh1在标准状态下，压缩因子Z＝1，因此体积流率可以得到：98.3910308..13011432731062204.5m3h1vnVnRTp（4）Amagat定律和普遍化压缩因子先查得各物质的压缩因子，再使用分体积定律进行计算组分ZyZTTpTrpriipcc甲烷氮气乙烷合计1.632.461.0282.0162.7311.9030.880.990.320.7390.08910.02240.8507ZRT0.85078.31437.8273VyVy2.371104m3mol1i9.27106piiiv25nV2.3711049.839104molh1105.44kmolh1在标准状态下，压缩因子Z＝1，因此体积流率可以得到：105.4410308..13011432731062362.5m3h1vnVnRTp2－28．试分别用下述方法计算CO2（1）和丙烷（2）以3.5：6.5的摩尔比混合的混合物在400K和13.78MPa下的摩尔体积。k（1）RK方程，采用Prausnitz建议的混合规则（令＝0.1）ij（2）Pitzer的普遍化压缩因子关系数。1）RK方程由附录三查得CO2（1）和丙烷（－48e）以及（2－13a）、（2－13b）进行计算，得解：（2）的临界参数值，并把这些值代入方程（2－48a）～（2出的结果如下：ijijT/KcijV/(m3kmol-1)Zcijcijp/MPacij22第二章流体的p-V-T关系11304.2369.8335.47.3824.2485.4720.09400.20000.14040.2740.2770.27550.2280.1520.1902212并且组元a/(Pam6K0.5mol-2)b/(m3mol-1)CO2（1）6.460C3H8（2）18.292.968×10-56.271×10-51211.12由（2－51a）和（2－51b）得：ay2a2yyay2a＝0.3526.46020.350.6511.120.65218.29M1111212222＝13.58Pam6K0.5mol-2bybyb＝0.352.9681050.656.271105＝5.115105mmol3-1M1122ap＝13.5813.78106A＝0.8460MRT8.31440022.522.5Bbp＝5.11510513.78106＝0.2119MRT8.3144