第二章作业.doc

课堂作业13.立方型状态方程中a、b的量纲是什么？解析：质量、长度、时间和温度的量纲符号分别为M（mass）、L（length）、T（time）、 。Vander waals 方程中a的单位是：Kg*m5*s-2*mol-2，b的单位：m3*mol-1. a的量纲是M*L5*T-2,b的量纲是L3.第一次作业：1、21. 因为图2-1中蒸发线上的自由度为1，所以处在蒸发线上的1mol物质，只要温度给定，它的压力和体积就确定了。这一说法对吗？如果是错误的，您将如何改正呢？解析：自由度指的是能够维持系统原有相数而独立改变的变量（T、P、表示想组成的某些物质的相对含量Xi）。答：不对。处在蒸发线上的1mol物质，只要温度给定，它的压力是确定的，但是体积不能确定。2. 1mol理想气体在常压、25时的体积为多少？解析：理想气体状态方程:Pv=nRT，R=8.314 Jmol-1k-1=8.314 m3Pamol-1k-1，常压是101.325KPa解：v= nRT/P=（1*8.314*298.15）/101325=0.02446m3=24.46L第二次作业：男生：4（RK方程）、5；女生：4（SRK方程）、6；4.1mol丙烷放在2L的容器中，用R-K方程分别求100和6时容器内的压力。已知丙烷在6时的饱和蒸汽压为0.57MPa。 解析：R-K方程：p=RTv-b-aT12v（v+b）a=0.42748R2Tc2.5Pcb=0.08664RTcPc SRK方程：p=RTv-b-a（T）v（v+b） aT=0.42747R2Tc2Pca(Tr) b=0.08664RTcPca(Tr)=1+m(1-Tr0.5)2 m=0.480+1.574w-0.176w2Tr=T/Tc丙烷Tc=369.8k，Pc=4.25MPa=4.25*103kPa，w=0.153解：R-K方程：T=100时，T369.8k，容器中为单一气相 a=18.284，b=6.268*10-5，P=8.3145*373.152*10-3-18.284373.150.5*2*10-3*（2*10-3-6.268*10-5）=1.372 MPa T=6时，需要判断容器中是否为单一气相。P=8.3145*279.152*10-3-18.284279.150.5*2*10-3*（2*10-3-6.268*10-5）=0.933MPa0.57MPa因此，容器中丙烷已经部分液化，此时，容器中的压力P=0.57MPa。 （还要判断容器中是否为单一液相） SRK方程：T=100时，T369.8k，容器中为单一气相Tr=373.15/369.8=1.0090，m=0.7167，a(Tr)=0.9936,a(T)=0.94479,b=6.268*10-5P=8.3145*373.152*10-3-6.268*10-5-0.944792*10-3*（2*10-3+6.268*10-5）=1.3725MPa T=6，需要判断容器中是否为单一气相。Tr=279.15/369.8=0.7549，m=0.7167，a(Tr)=1.1968,a(T)=1.1380,b=6.268*10-5 P=8.3145*279.152*10-3-6.268*10-5-1.13802*10-3*（2*10-3+6.268*10-5）=0.922MPa0.57MPa，（还要判断容器中是否为单一液相）因此，容器中的丙烷已经部分液化，此时，容器中的压力位0.57MPa。5.推导式（2-24）、式（2-25）.解析：R-K方程：p=RTv-b-aT12v（v+b） a=0.42748R2Tc2.5Pc b=0.08664RTcPc 对于状态方程参数a、b的确定利用了临界等温线在临界点出现水平拐点的特征： (pv)T=TcV=Vc=(2pv2)T=TcV=Vc=0证明：pvT=TcV=Vc=-RTcvc-b2+aTc122vc+bvc2vc+b2=0 RTc(vc-b)2=aTc122vc+bvc2(vc+b)2 (1)2pv2T=TcV=Vc=2RTvc-b3+aTc122vc2vc+b2-22vc+b2v3vc+b3=0 RTvc-b3=aTc123vc2+3bvc+b2vc3vc+b3 (2) 将（1）式/（2）式，得b3+3vcb2+3vc2b-vc3=0,令t=b/Vc，则上式：t3+3t2+3t-1=0配方得：（t+1）3=2 t=32-1=0.25992 解得b=0.2599vc有pcvc=13 RTc,代入R-K方程中得a=1.2824RTc3/2vc=0.4274R2Tc2.5/Pc，因此，有a=0.42748R2Tc2.5Pcb=0.08664RTcPc6.证明式（2-33）中参数A、B均是无因次的。解析：SRK方程：p=RTv-b-a（T）v（v+b） aT=0.42747R2Tc2Pca(Tr) b=0.08664RTcPca(Tr)=1+m(1-Tr0.5)2 m=0.480+1.574w-0.176w2Tr=T/Tcv 单位是m3/molb的单位是m3/mol;RTv-b 单位是J*mol-1*K-1*Km3*mol-1=(m3*Pa*mol-1*K-1)*Km3*mol-1=Pa=Kg*m*s-2m2=Kg*m-1*s-2因此，aTvv+b 单位是Kg*m-1*s-2a（T）的单位是（Kg*m-1*s-2）*（m3*mol-1）2=Kg*m5*s-2*mol-2A=apR2T2 的单位是：Kg*m5*s-2*mol-2*（Kg*m-1*s-2）(m3*Pa*mol-1*K-1)2*K2=Kg*m5*s-2*mol-2*（Kg*m-1*s-2）(m3*（Kg*m-1*s-2）*mol-1*K-1)2*K2=1B=bpRT的单位是：(m3*mol-1)*Pa(m3*Pa*mol-1*K-1)*K=1因此，参数A、B均是无因次的。7.分别推导van der Waals、R-K、SRK、P-R状态方程的第二、第三维里系数表达式，分别用这些表达式计算甲烷在0时的第二、第三维里系数，并计算在该温度下、9MPa时的压缩因子。您认为哪个压缩因子的值更为可靠？解析:关于第二、第三维里系数的定义见课本P9， B=11!(Z)T =0 C=12!(2Z2)T =0 Vander waals 方程：P=RTv-b-av2R-K方程：p=RTv-b-aT12v（v+b）SRK方程：p=RTv-b-a（T）v（v+b）P-R方程：p=RTv-b-a(T)vv+b+b(v-b)Z=PvRT,=1v，因此，将上述方程左右两边均乘以vRT，以Vander waals方程为例，Z=PvRT=vv-b-av RT=11-b-aRT Z=b1-b2-aRT B=11!(Z)T =0=lim0(b1-b2-aRT)=b-aRT 2Z2=b1-b2-aRT=2b21-b3 C=12!(2Z2)T =0=12*lim0(2b21-b3)=b2第二类维里方程：Z = 1+BP/(RT)+(C-B2)P2/(RT)2 （见课本P10）R-K方程：B=b-aRT32 ,C=b2+abRT32 ;Z = 1+BP/(RT)+(C-B2)P2/(RT)2SRK方程: B = b-a(T)/(RT), C = b2+a(T)b/(RT), Z = 1+BP/(RT)+(C-B2)P2/(RT)2;PR方程： B = b-a(T)/(RT), C = b2+2a(T)b/(RT), Z = 1+BP/(RT)+(C-B2)P2/(RT)2 van der waals方程B = b-a/(RT) 代入数值后B = -5.81810-5C = b2 代入数值后 C = 1.85010-9Z = 1+BP/(RT)+(C-B2)P2/(RT)2 代入数值后 Z = 0.7453RK方程： B = b-a/(RT3/2) 代入数值后 B = -5.58010-5C = b2+ab/(RT3/2) 代入数值后 C = 3.44110-9Z = 1+BP/(RT)+(C-B2)P2/(RT)2 代入数值后 Z = 0.7840SRK方程： B = b-a(T)/(RT) 代入数值后 B = -5.35510-5C = b2+a(T)b/(RT) 代入数值后 C = 3.37510-9Z = 1+BP/(RT)+(C-B2)P2/(RT)2 代入数值后 Z = 0.7958PR方程： B = b-a(T)/(RT) 代入数值后 B = -6.65910-5C = b2+2a(T)b/(RT) 代入数值后 C = 5.716610-9Z = 1+BP/(RT)+(C-B2)P2/(RT)2 代入数值后 Z = 0.7562RK方程的压缩因子的值更可靠，SRK、P-R主要用于计算相平衡，R-K方程主要考虑了气体的压缩因子。8.利用附录中提供的立方方程根的求解子程序，选用适当的状态方程，计算乙烷在-40、-20、0、10、30时的饱和蒸气压、饱和液体体积、饱和蒸汽体积。已知乙烷在-40、-20、0、10、30时的饱和蒸气压实测值分别为7.748atm、14.154 atm、23.708 atm、29.930 atm、45.974 atm。请将计算值和实验值进行比较。解：9. 对所学的状态方程进行分类，并对比各方程的优缺点、适用范围。类型特点vdW型状态方程（以立方型状态方程为代表）立方型状态方程因其简单性和可靠性在工程上被认为是最为实用的状态方程。它适合于只需准确计算部分物性的场合，在过程模拟中普遍采用。对立方型状态方程的修正主要集中于温度函数a(T)和方程函数形式P(v)。前者旨在提高方程计算饱和蒸气压的精度(对相平衡计算具有重要影响)，后者旨在提高计算液相密度的精度。Virial型和BWR型状态方程截断的Virial方程形式简单，文献中拥有大量的第二维里系数的数据，而且混合规则理论上严格，但仅可用于气相；多参数的BWR和BWRS方程其主要优点是可对纯组分性质(相平衡和体积性质)作准确计算。对应状态原理型状态方程主要优点是可以利用高精度参考流体的状态方程。当所研究的流体的性质于参考流体相差不大时，计算的精度较高。11. 用普遍化第二维里系数法计算0，9MPa时甲烷的压缩因子，并和习题7结果做比较，分析一下谁更可靠？解：类似书25页例题，24页公式普遍化第二维里系数表达式为：查得甲烷物性参数Tc190.4K，Pc4.6MPa，偏心因子0.011。求得0的对比温度为1.4346，则有： 课本24页，第二维里系数法只能在较低的对比压力（Pr1-2）下采用，RK方程更可靠些，如：例2-6。12. 推导van der waals流体和RK流体的焓和熵的微分表达式（只能含PVT热容等可测变量）解析：从van der waals方程和RK方程自身出发，虽然书上29页给出了H(T,P)的dH和dS的全微分表达式，但是VTP的计算存在很大的困难，所以就尝试用H(T,V)来解决本题；另外，考虑到状态方程的表达式问题，对于解内能U、焓H、熵S的问题一般用T、V为独立变量来解决比较方便。解： H=U+Pv dH=dU+d（Pv） (1)U=U(T,v) dU=（UT）vdT+（Uv）Tdv=cvdT+（Uv）Tdv (2)dU=TdS-Pdv （Uv）T=TSvT-P (3)Maxwell 关系式：SvT=PTv （见28页） (4)将(4) (3) (2) 代入(1)，得： dH=cvdT+TPTvdv+vdP (5)dH=dU+Pdv+vdP=TdS- Pdv+Pdv+vdP= TdS+vdP dS=dH-vdPT=CvTdT+PTvdv (6)1) van der waals方程 dH=cvdT+TPTvdv+vdP=cvdT+RTv-bdv+vdP dS=CvTdT+PTvdv=CvTdT+Rv-bdv.2) RK方程： dH=cvdT+TPTvdv+vdP=cvdT+(RTv-b+a2T12v(v+b)dv+vdP dS=CvTdT+PTvdv=CvTdT+(Rv-b+a2T32v(v+b)dv13. 将以下偏导数用仅含p、v、T、Cp、Cv这些可测变量的关系式表示（1）解：由麦克斯韦关系式（见28页）直接有SPT=-VTP （2）PVS解析：如题所示，应该写成S(P,V)的微分表达式，令dS=0，则求出上式。直接求S(P,V)的微分表达式困难，一般的思路就是从S(T,V)、S(T,P)入手,先写出S(T,V)、S(T,P)的微分表达式解：12题中有dS=CVTdT+PTVdV TdSCv=dT+TCVPTVdV (1) dH=HTPdT+HPTdP=CPdT+HPTdP (2) dH = TdS + VdP HPT=TSPT+V=-TVTP+V (3) 将（3）代入（2）,得：dH = CPdT+(-TVTP+V)dP dS=dH-VdPT=CPdTT-VTPdPTdSCP=dT-TCPVTPdP (4)(注:dS有两种表达方式S(T,V)和S(T,P)，在下面的多个题中有用到)(1)-(4),（求出S（P,V）的微分表达式）得：TCV-TCPdS=TCVPTVdV+TCPVTPdP （5）固定S不变，除dv得：TCVPTVdV+TCPVTPdP=0 （6）所以有：PVS=-CPPTVCVVTP（3）解析：跟（2）中的解题思路一样，将H（T,P）的微分表达式写出来，然后令dH=0，求出所求。解：dH = CPdT+(V-TVTP)dP (1)（推导见书29页） 固定H，两边同除dP：CPTPH+V-TVTP=0 (2) 有：TPH=TVTP-VCP（4）解析：跟（2）（3）中的解题思路一样，将S（T,P）的微分表达式写出来，然后令dS=0，求出所求。解：dS=dH-VdPT=CPdTT-VTPdP （推导过程见课本29页） 固定S，同除dp：CPTTPS-vTP=0 TPS=TCPvTP（5）解析：跟（2）（3）（4）中的解题思路一样，将U（T,V）的微分表达式写出来，然后令dU=0，求出所求。解： U=U(T,V) dU=UTVdT+UVTdV=CVdT+UVTdV (1) （看到类似的这样UVT不能直接看出结果的，就联系相应的热力学基本方程）dU=TdS-PdV UVT=TSVT-P (2) Maxwell表达式： SVT=PTV (3)将（2）、（3）代入（1），得：dU=CVdT+TPTV-PdV (4)固定U，同除dT：CV+TPTV-PVTU=0 VTU=CVP-TPTV （6）解析：与上面的（2）（3）（4）（5）题不一样，解：UPS=H-PVPS=HPS-PVPS=HPS-V+PVPS (1) dH=TdS+VdP HPS=V (2) 在（2）题中有PVS=-CPPTVCVVTP VPS=-CVVTPCPPTV (3) 将（2）、（3）代入（1）中，有：UPS=PCVVTPCPPTV14. 推导衡压热容和衡容热容之差（Cp-Cv）和PVT之间的关系，并具体针对vdWs方程推导其计算式。解：dH=TdS+VdP dU=TdS-PdV CP=HTP=TSTP CV=UTV=TSTV CP-CV=TSTP-STV （两种思路：STP 和 STV各自表达成PVT的关系；另一种是（STP-STV）整体表达成PVT的关系）令S=f(T,V), 有：dS=STVdT+SVTdV STP=STV+SVTVTP CP-CV=TSVTVTP有Maxwell 关系式：SVT=PTV CP-CV=TPTVVTPvander Waals 方程：P=RTv-b-av2 PTV=Rv-b 固定P，同除dT，得： VTP=Rv-bRT(v-b)2- 2av3CP-CV=TPTVVTP=TRv-bRv-bRT(v-b)2- 2av3= R2Tv3v3RT-2a(v-b)215. 分别推导vdWs方程，RK方程的剩余焓和剩余熵表达式。并用RK方程计算CO2气体从25、0.5MPa压缩到100，1.5MPa过程中摩尔焓熵的变化解： HR=v-P+TPTvdv+Z-1RT SR=v-Rv+PTvdv+RlnZVander waals 方程：P=RTv-b-av2R-K方程：p=RTv-b-aT12v（v+b）1) Vander waals 方程：P=RTv-b-av2PTv=Rv-b HR=v-P+TPTvdv+Z-1RT=v-P+RTv-bdv+Z-1RT HR=-av+Z-1RT SR=v-Rv+PTvdv+RlnZ=v-Rv+Rv-bdv+RlnZ=Rlnv-bv+RlnZ2) R-K方程：P=RTv-b-aT12v（v+b）PTv=Rv-b+12aT32v(v+b) HR=v-P+TPTvdv+Z-1RT HR=v-P+RTv-b+12aT12v(v+b)dv+Z-1RT=v32aT12v(v+b)dv+Z-1RTHR=32abT12lnvv+b+Z-1RT SR=v-Rv+PTvdv+RlnZ=Rlnv-bv+12abT32lnvv+b+RlnZ3) CO2的 Tc=304.1K,Pc=7.38MPa=7380kPaa=0.