参考答案化工热力学第二版马沛生著化学工业出版社课后答.docx
《参考答案化工热力学第二版马沛生著化学工业出版社课后答.docx》由会员分享,可在线阅读,更多相关《参考答案化工热力学第二版马沛生著化学工业出版社课后答.docx(37页珍藏版)》请在冰豆网上搜索。
![参考答案化工热力学第二版马沛生著化学工业出版社课后答.docx](https://file1.bdocx.com/fileroot1/2023-1/7/769337ea-4760-4796-93f7-448d3d9a8dc0/769337ea-4760-4796-93f7-448d3d9a8dc01.gif)
参考答案化工热力学第二版马沛生著化学工业出版社课后答
参考答案化工热力学第二版(马沛生著)化学工业出版社课后答
习题:
2-1.为什么要研究流体的pVT关系?
答:
在化工过程的分析、研究与设计中,流体的压力p、体积V和温度T是流体最基本的性质之一,并且是可以通过实验直接测量的。
而许多其它的热力学性质如内能U、熵S、Gibbs自由能G等都不方便直接测量,它们需要利用流体的pCVCT数据和热力学基本关系式进行推算;此外,还有一些概念如逸度等也通过pCVCT数据和热力学基本关系式进行计算。
因此,流体的pCVCT关系的研究是一项重要的基础工作。
2-2.理想气体的特征是什么?
答:
假定分子的大小如同几何点一样,分子间不存在相互作用力,由这样的分子组成的气体叫做理想气体。
严格地说,理想气体是不存在的,在极低的压力下,真实气体是非常接近理想气体的,可以当作理想气体处理,以便简化问题。
理想气体状态方程是最简单的状态方程:
pV=RT
气压的对数与对比温度的倒数呈近似直线关系,即符合:
对于不同的流体,α具有不同的值。
但Pitzer发现,简单流体(氩、氪、氙)的所有蒸气压
网
w
1
logp=α1T
r
s
r
ps
其中,p=
pc
sr
ww
s
案
数据落在了同一条直线上,而且该直线通过Tr=0.7,logpr=1这一点。
对于给定流体对
.k
答:
纯物质的偏心因子ω是根据物质的蒸气压来定义的。
实验发现,纯态流体对比饱和蒸
hd
2-3.偏心因子的概念是什么?
为什么要提出这个概念?
它可以直接测量吗?
http:
//.cn
s
课
Pitzer把这一差值定义为偏心因子ω,即
后
之差来表征。
答
比蒸气压曲线的位置,能够用在Tr=0.7的流体与氩、氪、氙(简单球形分子)的logpr值
ω=logprs1.00
(Tr=0.7)
任何流体的ω值都不是直接测量的,均由该流体的临界温度Tc、临界压力pc值及
Tr=0.7时的饱和蒸气压ps来确定。
2-4.纯物质的饱和液体的摩尔体积随着温度升高而增大,饱和蒸气的摩尔体积随着温度的升高而减小吗?
答:
正确。
由纯物质的pCV图上的饱和蒸气和饱和液体曲线可知。
2-5.同一温度下,纯物质的饱和液体与饱和蒸气的热力学性质均不同吗?
答:
同一温度下,纯物质的饱和液体与饱和蒸气的Gibbs自由能是相同的,这是纯物质气液
平衡准则。
气他的热力学性质均不同。
2-6.常用的三参数的对应状态原理有哪几种?
答:
常用的三参数对比态原理有两种,一种是以临界压缩因子Zc为第三参数;另外一种是以Pitzer提出的以偏心因子ω作为第三参数的对应状态原理。
2-7.总结纯气体和纯液体pVT计算的异同。
答:
许多pCV-T关系如RKS方程、PR方程及BWR方程既可以用于计算气体的pCVCT,又都可以用到液相区,由这些方程解出的最小体积根即为液体的摩尔体积。
当然,还有许多状态方程只能较好地说明气体的pCV-T关系,不适用于液体,当应用到液相区时会产生较大的误差。
与气体相比,液体的摩尔体积容易测定。
除临界区外,温度(特别是压力)对液体容积性质的影响不大。
除状态方程外,工程上还常常选用经验关系式和普遍化关系式等方法来估算。
2-8.简述对应状态原理。
答:
对比态原理认为,在相同的对比状态下,所有的物质表现出相同的性质。
对比态原理是从适用于pCV-T关系两参数对比态原理开始的,后来又发展了适用于许多热力学性质和传递性质的三参数和更多参数的对比态原理。
2-9.如何理解混合规则?
为什么要提出这个概念?
有哪些类型的混合规则?
和分体积定律无法准确地描述流体混合物的pCV-T关系。
如何将适用于纯物质的状态方程扩展到真实流体混合物是化工热力学中的一个热点问题。
目前广泛采用的方法是将状态方程中的常数项,表示成组成x以及纯物质参数项的函数,这种函数关系称作为混合规则。
数项,使其能准确地描述真实流体混合物的pCV-T关系,常常是计算混合热力学性质的关键。
常用的混合规则包括适用于压缩因子图的虚拟临界性质的混合规则、维里系数的混合规
2-10.在一个刚性的容器中,装入了1mol的某一纯物质,容器的体积正好等于该物质的摩尔临界体积Vc。
如果使其加热,并沿着习题图2-1的pCT图中的1→C→2的途径变化(C是临界点)。
请将该变化过程表示在pCV图上,并描述在加热过程中各点的状态和现象。
解:
由于加热过程是等容过程,1→C→2是一条V=VC的等容线,所以在pCV图可以表示为如图的形式。
点1表示容器中所装的是该物质的汽液混合物(由饱和蒸汽和饱和液体组成)。
沿1-2线,是表示等容加热过程。
随着过程的进行,容器中的饱和液体体积与饱和蒸汽体积的相对比例有所变化,但由图可知变化不是很大。
到了临界点C点时,汽液相界面逐渐消失。
继续加热,容器中一直是均相的超临界流体。
在整个过程中,容器内的压力是不断增加的。
课
则以及适用于立方型状态方程的混合规则。
后
答
案
对于不同的状态方程,有不同的混合规则。
寻找适当的混合规则,计算状态方程中的常
网
w
ww
答:
对于真实流体,由于组分的非理想性及由于混合引起的非理想性,使得理想的分压定律
.k
hd
http:
//.cn
C=9.0×103m6kmol2,试计算:
(1)SO2在431K、10×105Pa下的摩尔体积;
时所作的功。
解:
(1)三项维里方程为:
Z=
(A)
网
w
ww
pVBC
=1++2RTVV
将p=10×105Pa,T=431K,B=0.159mkmol,C=9.0×10mkmol代入式(A)并整理得:
3
.k
1
hd
(2)在封闭系统内,将1kmolSO2由10×105Pa恒温(431K)可逆压缩到75×105Pa
0.279V3V2+0.159V9×106=0
答
迭代求解,初值为:
V=
后
案
RT
=3.5m3kmol1p
3
1
(2)压缩功W=pdV
课
迭代结果为:
V=3.39mkmol
由(A)式得:
p=RT
V2
BC1
+2+3,则:
VVV
BC1
W=RT∫+2+3dV
V1VVVV211C11=RTlnBVV2V2V2V11212
当p=75×105Pa时,用
(1)同样的方法解出:
aw
.com
3
6
2
2-11.已知SO2在431K下,第二、第三Virial系数分别为:
B=0.159mkmol,
31
(B)
V=0.212m3kmol1
将V1=3.39mkmol,V2=0.212mkmol代入式(B)解出:
3
1
3
1
W=77×105Jkmol1
2-12.试计算一个125cm3的刚性容器,在50℃和18.745MPa的条件下能贮存甲烷多少克(实验值为17g)?
分别用理想气体方程和RK方程计算(RK方程可以用软件计算)。
解:
由附录三查得甲烷的临界参数为:
Tc=190.56K,pc=4.599MPa,ω=0.011
(1)利用理想气体状态方程pV=RT得:
V=
RT8.314×(273.15+50)==1.433×104m3mol1=143.3cm3mol1
6
p18.745×10
(2)RK方程
式中:
b=0.08664RTc/pc=
bp2.985×105×18.745×106
B=0.2083
RT8.314×323.15
按照式(2-16a)Z=
课
ap3.2207×18.745×106
A=22.5=0.*****.5
RT8.314×323.15后
和式(2-16b)
迭代计算,取初值Z=1,迭代过程和结果见下表。
迭代次数
答
1Ah1h
2.2342=1hB1+h1h1h+
h=
案
a=0.*****RT
2
2.5c
0.*****×(8.314)×(190.56)/pc==3.2207Pam6K0.5mol-2
6
4.599×10
2
2.5
0.08664×8.314×190.56531
2.985×10mmol
4.599×106
bB0.2083==VZZ
Z
h
网w
ww
p=
RTa
0.5
VbTV(V+b)
.k
hd
V125
m=M=16×=13.95g
V143.3
aw
.com
01234
10.*****.*****.*****.8823
0.*****.*****.*****.*****.2361
V=
ZRT0.8823×8.314×323.154331
==1.265×10m/mol=126.5cmmolp18.745×106
V125
=16×=15.81gV126.5
m=M
可见,用RK方程计算更接近实验值。
2-13.欲在一个7810cm3的钢瓶中装入1kg的丙烷,且在253.2℃下工作,若钢瓶的安全工作压力为10MPa,问是否安全?
解:
查得丙烷的临界性质为:
Tc=369.83K,pc=4.248MPa,ω=0.152
使用RK方程:
首先用下式计算a,b:
a=0.*****RT
2
2.5
c
后
8.3142×369.832.560.5-2
/pc=0.*****×=18.296PamKmol
4.248×106
答
案
网w
p=
RTa
0.5
VbTV(V+b)
ww
V总7810×106
V===343.63×106m3mol1
n22.727
代入RK方程得:
p=9.870MPa非常接近于10MPa,故有一定危险。
2-14.试用RKS方程计算异丁烷在300K,3.704×105Pa时的饱和蒸气的摩尔体积。
已知实验值为V=6.081×10mmol。
解:
由附录三查得异丁烷的临界参数为:
Tc=407.8K,pc=3.640MPa,ω=0.177
3
3
1
Tr=T/Tc=300/407.8=0.7357
m=0.480+1.574ω0.176ω2=0.480+1.574×0.1770.176×0.1772=0.7531
课
b=0.08664RTc/pc=0.08664×
8.314×369.83531
=6.2771×10mmol
4.248×106
.k
n=
m1000==22.727molM44
hd
aw
.com
α(T)=1+m(1Tr0.5)=1+0.7531(10.*****.5)=1.2258
2
2[][
2
]
a(T)=aα(T)=0.4278RT
2c
22
(8.314)×(407.8)62
/pcα(T)=0.*****××1.2258=1.6548(Pam)/mol6
3.640×10
b=0.08664RTc/pc=0.08664×8.314×407.8/(3.640×106)=8.0700×105m3/molap1.6548×3.704×105
A=22=0.***-*****
RT8.314×300bp8.0700×105×3.704×105
B=0.01198
RT8.314×300
按照式(2-16a)Z=
1Ah1h
8.2245=1hB1+h1h1+h
和式(2-16b)
h=
迭代计算,取初值Z=1,迭代过程和结果见下表。
迭代次数0*****
0.*****.*****.*****.*****.9061
ww
2
2.5
1
.k
Z
网
w案课
(6.0316.1015)×102/6.031×102=1.2%
误差
2-15.试分别用RK方程及RKS方程计算在273K、1000×105Pa下,氮的压缩因子值,已知实验值为Z=2.0685。
解:
由附录三查得氮的临界参数为:
Tc=126.10K,pc=3.394MPa,ω=0.040
(1)RK方程
后
V=
ZRT0.9061×8.314×30023
==6.1015×10m/mol6p3.704×10
a=0.*****RT
2
b=0.08664RTc/pc=
答
2.5
c
0.*****×(8.314)×(126.10)60.5-2
/pc==1.5546PamKmol
3.394×106
0.08664×8.314×126.10531
2.6763×10mmol
3.394×106
hd
h
0.0*****.0*****.0*****.0*****.0*****.01322
bB0.01198==VZZ
aw
.com
ap1.5546×100×106
A=22.51.*****.5
RT8.314×273bp2.6763×105×1000×105
=1.1791B=
8.314×273RT
按照式(2-16a)Z=
1Ah1h
1.5489=1hB1+h1h1+h
和式(2-16b)
h=
bB1.1791
==VZZ
迭代计算,取初值Z=2,迭代过程和结果见下表。
01234……..
21.8622.*****.*****.8823
迭代不收敛,采用RK方程解三次方程得:
pV4.422×105×1000×105
==1.9485Z=
8.314×273RT
RKS方程
Tr=T/Tc=273/126.1=2.1649
m=0.480+1.574ω0.176ω2=0.480+1.574×0.0400.176×0.0402=0.5427
课
α(T)=1+m(1Tr0.5)=1+0.5427(12.*****.5)=0.5538
2
2[
后
答
案
][
2
网
w
]
2c
V=0.***-*****m3/mol
a(T)=aα(T)=0.4278RT
22
(8.314)×(126.1)6
()/mol/pcα(T)=0.*****××0.5538=0.0*****Pam6
ww
3.394×10
b=0.08664RTc/pc=0.08664×8.314×126.1/(3.394×106)=2.6763×105m3/molap0.0*****×1000×105
A=22==1.*****
RT8.314×273.k
hd
aw
.com
0.*****.*****.*****.*****.2361
迭代次数Zh
bp2.6763×105×1000×105
B==1.1791
RT8.314×273
按照式(2-16a)Z=
1Ah1h
1.2621=1hB1+h1h1+h
和式(2-16b)同样迭代不收敛
h=
bB1.1791
==VZZ
采用RKS方程解三次方程得:
V=0.***-*****m3/mol
2-16.试用下列各种方法计算水蒸气在107.9×105Pa、593K下的比容,并与水蒸气表查出的数据(V=0.01687mkg)进行比较。
(1)理想气体定律
(2)维里方程(3)普遍化RK方程
3
1
(1)理想气体定律
误差=
(2)维里方程
Tr=
p107.9×105
pr===0.4896
pc22.055×10
使用普遍化的第二维里系数:
课
T593==0.916Tc647.13
后
0.016870.02538
×100%=50.5%
0.01687
答
案
V=
RT8.314×59363131
==4.569×10mmol=0.02538mkgp107.9×105
B(0)=0.0830.422/Tr1.6=0.083
B
(1)=0.1390.172/Tr4.2=0.139
网
w
解:
从附录三中查得水的临界参数为:
Tc=647.13K,pc=22.055MPa,ω=0.345
0.422
=0.4026Tr1.6
0.172
=0.*****.2
Tr
ww
.k
hd
http:
//.cn
pV4.512×105×1000×105
Z===1.9881
RT8.314×273
Bpc
=B(0)+ωB
(1)=0.4026+0.345×(0.1096)=0.4404RTc
BpcprBp0.489
Z=1+=1+=1+×(0.4404)=0.76490.916RTRTcTrV=
ZRT0.7649×8.314×593
==3.495×106m3mol1=0.01942m3kg1
5
p107.9×100.016870.01942
×100%=15.1%
0.01687
a1
1hbTr1.5bprZTr
h1+h
误差=
(3)普遍化R-K方程
Z=
h=
将对比温度和对比压力值代入并整理的:
h=
联立上述两式迭代求解得:
Z=0.7335
误差=
水是极性较强的物质
2-17.试分别用
(1)vanderWaals方程;
(2)RK方程;(3)RKS方程计算273.15K时将CO2压缩到体积为550.1cmmol所需要的压力。
实验值为3.090MPa。
解:
从附录三中查得CO2的临界参数为:
Tc=304.19K,pc=7.382MPa,ω=0.228
(1)vanderWaals方程
3
1
课
后
0.016870.01862
×100%=10.4%
0.01687
答
V=
ZRT0.7335×8.314×59363131
==3.3515×10mmol=0.01862mkgp107.9×105
p=
式中:
RTa
2
VbV
2
2
22c
27×(8.314)×(304.19)32
a=27RT/64pc==0.3655Pammol
64×7.382×106
案
网w
bpr0.04625
=ZTrZ
ww
Z=
a1
1hbTr1.51h15.628=
++1h1h1h
.k
hd
aw
.com
(2-38a)
(2-38b)
b=RTc/8pc=
则
8.314×304.19531
=4.282×10mmol
8×7.382×106
:
p=
RTa8.314×273.150.3655
2=-=3.269×106Pa=3.269MPa662
VbV550.1×10-42.82×*****.1×1063.0903.269
×100%=5.79%
3.090
误差%=
(2)RK方程
p=
RTa
0.5
VbTV(V+b)
式中:
0.*****×(8.314)×(304.19)6.4599Pam6K0.5mol-2a=0.*****R2Tc2.5/pc=6
7.382×10
2
2.5
(3)RKS方程
答
案
=3.138×106Pa=3.138MPa
3.0903.138
误差%=×100%=1.55%
3.090
后
课
式中,
a(T)=aα(T)=0.4278R2Tc2/pcα(T)α(T)=1+m(1Tr0.5)
2
网w
p=
RTa
0.5
VbTV(V+b)
8.314×273.156.4599
=-
550.1×106-29.68×106273.150.5×550.1×106×550.1+29.68×106
RTa(T)
VbVV+b[
ww
p=
而,m=0.480+1.574ω0.176ω=0.480+1.574×0.228-0.176×(0.228)=0.8297
2
.k
]
2
b=0.08664RTc/pc=
0.08664×8.314×304.19
2.968×105m3mol1
6
7.382×10
hd
0.5
2273.150.5
则,α(T)=1+m(1Tr)=1+0.8297×1=1.089304.19
[]
aw
.com
2
0.*****×(8.314)×(304.19)a(T)=aα(T)=0.*****RT/pcα(T)=×1.0896
7.382×10
=0.*****Pam3mol1
2
2
22c
b=0.08664RTc/pc=p=
0.08664×8.314×304.19
2.968×105m3mol1
6
7.382×10
RTa(T)8.314×273.150.*****
-
VbVV+b550.1×106-29.68×106550.1×106×550.1+29.68×106
3.0903.099
×100%=0.291%
3.090
=3.099×106Pa=3.099MPa
误差%=
最好。
RK方程的计算精度还可以。
2-18.一个体积为0.3m3的封闭储槽内贮乙烷,温度为290K、压力为25×105Pa,若将乙烷加热到479K,试估算压力将变为多少?
因此:
Tr1=T1/Tc=290/305.32=0.95pr1=p1/pc=2.5/48.72=0.513
故使用图2-11,应该使用普遍化第二维里系数计算
B(0)=0.0830.422/Tr1.6=0.083
网w
0.950.172
0.422
案
B
(1)=0.1390.172/Tr4.2=0.139
prBp
=1+(B(0)+ωB
(1))TRTr
0.95答
V=
n=
V总0.3==392.2(mol)5V76.5×10
加热后,采用RK方程进行计算。
其中:
T=479K,摩尔体积仍然为V=76.5×10mmol,首先计算:
5
3
1
课
ZRT0.7935×8.314×290531
()==76.5×10mmol5
p25×10
后
Z=1+
0.513
()=