化工热力学(第三版)课后答案完整版_朱自强
化工热力学第三版课后习题答案全
化工热力学第三版课后习题答案第一章比较简单略第二章2-1.使用下述方法计算1kmol 甲烷贮存在体积为0.1246m 3、温度为50℃的容器中产生的压力:(1)理想气体方程;(2)R-K 方程;(3)普遍化关系式。
解:甲烷的摩尔体积V =0.1246 m 3/1kmol=124.6 cm 3/mol查附录二得甲烷的临界参数:T c =190.6K P c =4.600MPa V c =99 cm 3/mol ω=0.008 (1) 理想气体方程P=RT/V=8.314×323.15/124.6×10-6=21.56MPa(2) R-K 方程22.522.560.5268.314190.60.427480.42748 3.2224.610c cR T a Pa m K mol P -⨯===⋅⋅⋅⨯53168.314190.60.086640.08664 2.985104.610c c RT b m mol P --⨯===⨯⋅⨯ ∴()0.5RT aP V b T V V b =--+()()50.5558.314323.15 3.22212.46 2.98510323.1512.461012.46 2.98510---⨯=--⨯⨯⨯+⨯=19.04MPa (3) 普遍化关系式323.15190.61.695r c T T T === 124.699 1.259r c V V V ===<2∴利用普压法计算,01Z Z Z ω=+∵ c r ZRTP P P V == ∴c r PV Z P RT =654.61012.46100.21338.314323.15cr r r PV Z P P P RT -⨯⨯⨯===⨯迭代:令Z 0=1→P r0=4.687 又Tr=1.695,查附录三得:Z 0=0.8938 Z 1=0.462301Z Z Z ω=+=0.8938+0.008×0.4623=0.8975此时,P=P c P r =4.6×4.687=21.56MPa同理,取Z 1=0.8975 依上述过程计算,直至计算出的相邻的两个Z 值相差很小,迭代结束,得Z 和P 的值。
化工热力学第三版(完全版)课后习题问题详解
化工热力学课后答案第1章 绪言一、是否题1. 封闭体系的体积为一常数。
(错)2. 封闭体系中有两个相βα,。
在尚未达到平衡时,βα,两个相都是均相敞开体系;达到平衡时,则βα,两个相都等价于均相封闭体系。
(对)3. 理想气体的焓和热容仅是温度的函数。
(对)4. 理想气体的熵和吉氏函数仅是温度的函数。
(错。
还与压力或摩尔体积有关。
)5. 封闭体系的1mol 气体进行了某一过程,其体积总是变化着的,但是初态和终态的体积相等,初态和终态的温度分别为T 1和T 2,则该过程的⎰=21T T V dT C U ∆;同样,对于初、终态压力相等的过程有⎰=21T T P dT C H ∆。
(对。
状态函数的变化仅决定于初、终态与途径无关。
)二、填空题1. 状态函数的特点是:状态函数的变化与途径无关,仅决定于初、终态 。
2. 封闭体系中,温度是T 的1mol 理想气体从(P i ,V i )等温可逆地膨胀到(P f ,V f ),则所做的功为()f i rev V V RT W ln =(以V 表示)或()i f rev P P RT W ln = (以P 表示)。
3. 封闭体系中的1mol 理想气体(已知igP C ),按下列途径由T 1、P 1和V 1可逆地变化至P 2,则A 等容过程的 W = 0 ,Q =()1121T P P R C igP⎪⎪⎭⎫ ⎝⎛--,错误!未找到引用源。
U =()1121T PP R C igP⎪⎪⎭⎫⎝⎛--,错误!未找到引用源。
H = 1121T P P C ig P ⎪⎪⎭⎫ ⎝⎛-。
B 等温过程的 W =21lnP P RT -,Q =21ln P PRT ,错误!未找到引用源。
U = 0 ,错误!未找到引用源。
H = 0 。
C 绝热过程的 W =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛--11211igPC RigPP P R V P R C ,Q = 0 ,错误!未找到引用源。
化工热力学第三版(完全版)课后习题答案
化工热力学课后答案第1章 绪言一、是否题1. 封闭体系的体积为一常数。
(错)2. 封闭体系中有两个相βα,。
在尚未达到平衡时,βα,两个相都是均相敞开体系;达到平衡时,则βα,两个相都等价于均相封闭体系。
(对)3. 理想气体的焓和热容仅是温度的函数。
(对)4. 理想气体的熵和吉氏函数仅是温度的函数。
(错。
还与压力或摩尔体积有关。
)5. 封闭体系的1mol 气体进行了某一过程,其体积总是变化着的,但是初态和终态的体积相等,初态和终态的温度分别为T 1和T 2,则该过程的⎰=21T T V dT C U ∆;同样,对于初、终态压力相等的过程有⎰=21T T P dT C H ∆。
(对。
状态函数的变化仅决定于初、终态与途径无关。
)二、填空题1. 状态函数的特点是:状态函数的变化与途径无关,仅决定于初、终态 。
2. 封闭体系中,温度是T 的1mol 理想气体从(P i ,V i )等温可逆地膨胀到(P f ,V f ),则所做的功为()f i rev V V RT W ln =(以V 表示)或()i f rev P P RT W ln = (以P 表示)。
3. 封闭体系中的1mol 理想气体(已知igP C ),按下列途径由T 1、P 1和V 1可逆地变化至P 2,则A 等容过程的 W = 0 ,Q =()1121T P P R C igP ⎪⎪⎭⎫ ⎝⎛--,U =()1121T PPR C igP ⎪⎪⎭⎫⎝⎛--,H =1121T P P C ig P ⎪⎪⎭⎫ ⎝⎛-。
B 等温过程的 W =21lnP P RT -,Q =21ln P PRT ,U = 0 ,H = 0 。
C 绝热过程的 W =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛--11211igPC RigPP P R V P R C ,Q = 0 ,U =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛-11211ig P C R igPP P R V P R C ,H =1121T P P C igP C R ig P ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛。
化工热力学第三版(完全版)课后习题答案
化工热力学课后答案第1章 绪言一、是否题1. 封闭体系的体积为一常数。
(错)2. 封闭体系中有两个相βα,。
在尚未达到平衡时,βα,两个相都是均相敞开体系;达到平衡时,则βα,两个相都等价于均相封闭体系。
(对) 3. 理想气体的焓和热容仅是温度的函数。
(对)4. 理想气体的熵和吉氏函数仅是温度的函数。
(错。
还与压力或摩尔体积有关。
)5. 封闭体系的1mol 气体进行了某一过程,其体积总是变化着的,但是初态和终态的体积相等,初态和终态的温度分别为T 1和T 2,则该过程的⎰=21T T V dT C U ∆;同样,对于初、终态压力相等的过程有⎰=21T T P dT C H ∆。
(对。
状态函数的变化仅决定于初、终态与途径无关。
) 二、填空题1. 状态函数的特点是:状态函数的变化与途径无关,仅决定于初、终态 。
2. 封闭体系中,温度是T 的1mol 理想气体从(P i ,V i )等温可逆地膨胀到(P f ,V f ),则所做的功为()f i rev V V RT W ln =(以V 表示)或()i f rev P P RT W ln = (以P 表示)。
3. 封闭体系中的1mol 理想气体(已知ig P C ),按下列途径由T 1、P 1和V 1可逆地变化至P 2,则A 等容过程的 W = 0 ,Q =()1121T P P R C igP ⎪⎪⎭⎫ ⎝⎛--,错误!未找到引用源。
U =()1121T PPR C igP ⎪⎪⎭⎫⎝⎛--,错误!未找到引用源。
H = 1121T P P C ig P ⎪⎪⎭⎫ ⎝⎛-。
B 等温过程的 W =21lnP P RT -,Q =21ln P PRT ,错误!未找到引用源。
U = 0 ,错误!未找到引用源。
H = 0 。
C 绝热过程的 W =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛--11211igPC RigPP P R V P R C ,Q = 0 ,错误!未找到引用源。
化工热力学(第三版)第3章答案
化工热力学(第三版)习题解答集朱自强、吴有庭、李勉编著前言理论联系实际是工程科学的核心。
化工热力学素以概念抽象、难懂而深深印在学生的脑海之中。
特别使他们感到困惑的是难以和实际问题进行联系。
为了学以致用,除选好教科书中的例题之外,很重要的是习题的安排。
凭借习题来加深和印证基本概念的理解和运用,补充原书中某些理论的推导,更主要的是使学生在完成习题时能在理论联系实际的锻炼上跨出重要的一步。
《化工热力学》(第三版)的习题就是用这样的指导思想来安排和编写的。
《化工热力学》自出版以来,深受国内同行和学生的关注和欢迎,但认为习题有一定的难度,希望有一本习题集问世,帮助初学者更有效地掌握基本概念,并提高分析问题和解决问题的能力。
为此我们应出版社的要求把该书第三版的习题解撰并付印,以飨读者。
在编写过程中除详尽地进行习题解答外,还对部分习题列出了不同的解题方法,便于读者进一步扩大思路,增加灵活程度;对部分有较大难度的习题前加上“*”号,如果教学时间较少,可以暂时不做,但对能力较强的学生和研究生也不妨一试。
使用本题解的学生,应该先对习题尽量多加思考,在自学和独自完成解题的基础上加以利用和印证,否则将与出版此书的初衷有悖。
参加本习题题解编写的人员是浙江大学化工系的朱自强教授、南京大学化工系的吴有庭教授、以及李勉博士等,浙江大学的林东强教授、谢荣锦老师等也对本习题编写提供了有益的帮助。
在此深表感谢。
由于编写时间仓促,有些地方考虑不周,习题题解的写作方法不善,甚至尚有解题不妥之处,希望读者能不吝赐教,提出宝贵意见,以便再版时予以修改完善。
第二章 流体的压力、体积、浓度关系:状态方程式2-1 试分别用下述方法求出400℃、4.053MPa 下甲烷气体的摩尔体积。
(1) 理想气体方程;(2) RK 方程;(3)PR 方程;(4) 维里截断式(2-7)。
其中B 用Pitzer 的普遍化关联法计算。
[解] (1) 根据理想气体状态方程,可求出甲烷气体在理想情况下的摩尔体积idV 为33168.314(400273.15)1.381104.05310id RT V m mol p --⨯+===⨯⋅⨯ (2) 用RK 方程求摩尔体积将RK 方程稍加变形,可写为0.5()()RT a V b V b p T pV V b -=+-+ (E1)其中2 2.50.427480.08664c c ccR T a p RT b p ==从附表1查得甲烷的临界温度和压力分别为c T =190.6K, c p =4.60MPa ,将它们代入a, b 表达式得2 2.56-20.560.427488.314190.6 3.2217m Pa mol K 4.6010a ⨯⨯==⋅⋅⋅⨯ 53160.086648.314190.6 2.9846104.6010b m mol --⨯⨯==⨯⋅⨯ 以理想气体状态方程求得的idV 为初值,代入式(E1)中迭代求解,第一次迭代得到1V 值为5168.314673.15 2.9846104.05310V -⨯=+⨯⨯350.563353.2217(1.38110 2.984610)673.15 4.05310 1.38110(1.38110 2.984610)-----⨯⨯-⨯-⨯⨯⨯⨯⨯⨯+⨯ 3553311.381102.984610 2.1246101.389610m mol -----=⨯+⨯-⨯=⨯⋅第二次迭代得2V 为353520.563353553313.2217(1.389610 2.984610)1.381102.984610673.154.05310 1.389610(1.389610 2.984610)1.381102.984610 2.1120101.389710V m mol ------------⨯⨯-⨯=⨯+⨯-⨯⨯⨯⨯⨯⨯+⨯=⨯+⨯-⨯=⨯⋅1V 和2V 已经相差很小,可终止迭代。
化工热力学答案课后总习题答案详解
化工热力学答案—课后总习题答案详解第二章习题解答一.问答题:2-1为什么要研究流体的"VT关系?【参考答案】:流体P-V-T关系是化工热力学的基石,是化工过程开发和设讣、安全操作和科学研究必不可少的基础数据。
(I)流体的PVT关系可以直接用于设汁。
(2)利用可测的热力学性质(T, P, V等)计算不可测的热力学性质(H, S, G.等)。
只要有了旷/T关系加上理想气体的C;;, 可以解决化工热力学的大多数问题匚以及该区域的特征:同时指岀其它重要的点、2- 2 ⅛ P-V图上指出超临界萃取技术所处的区域,而以及它们的特征。
【参考答案】:1)超临界流体区的特征是:环、P>Pco2)临界点C的数学特征:(^PM Z)/ =° (在C点)($2p/刃2)・0 (在C点)3)饱和液相线是不同压力下产生第一个气泡的那个点的连线:4)饱和汽相线是不同压力下产生第一个液滴点(或露点)那个点的连线。
5)过冷液体区的特征:给左压力下液体的温度低于该压力下的泡点温度。
6)过热蒸气区的特征:给左压力下蒸气的温度髙于该压力下的露点温度。
7)汽液共存区:在此区域温度压力保持不变,只有体积在变化。
2-3要满足什么条件,气体才能液化?【参考答案】:气体只有在低于7;条件下才能被液化。
2-4不同气体在相同温度压力下,偏离理想气体的程度是否相同?你认为哪些是决左偏离理想气体程度的最本质因素?【参考答案】:不同。
真实气体偏离理想气体程度不仅与7∖ P有关,而且与每个气体的临界特性有关,即最本质的因素是对比温度、对比压力以及偏心因子7;, /和Q。
2-5偏心因子的概念是什么?为什么要提出这个槪念?它可以直接测呈:吗?【参考答案】:偏心因子。
为两个分子间的相互作用力偏离分子中心之间的作用力的程度。
其物理意义为:一般流体与球形非极性简单流体(氮,氟、毎)在形状和极性方而的偏心度。
为了提高计算复杂分子压缩因子的准确度。
化工热力学第三版(完全版)课后习题答案解析
化工热力学课后答案第1章 绪言一、是否题1. 封闭体系的体积为一常数。
(错)2. 封闭体系中有两个相βα,。
在尚未达到平衡时,βα,两个相都是均相敞开体系;达到平衡时,则βα,两个相都等价于均相封闭体系。
(对)3. 理想气体的焓和热容仅是温度的函数。
(对)4. 理想气体的熵和吉氏函数仅是温度的函数。
(错。
还与压力或摩尔体积有关。
)5. 封闭体系的1mol 气体进行了某一过程,其体积总是变化着的,但是初态和终态的体积相等,初态和终态的温度分别为T 1和T 2,则该过程的⎰=21T T V dT C U ∆;同样,对于初、终态压力相等的过程有⎰=21T T P dT C H ∆。
(对。
状态函数的变化仅决定于初、终态与途径无关。
)二、填空题1. 状态函数的特点是:状态函数的变化与途径无关,仅决定于初、终态 。
2. 封闭体系中,温度是T 的1mol 理想气体从(P i ,V i )等温可逆地膨胀到(P f ,V f ),则所做的功为()f i rev V V RT W ln =(以V 表示)或()i f rev P P RT W ln = (以P 表示)。
3. 封闭体系中的1mol 理想气体(已知igP C ),按下列途径由T 1、P 1和V 1可逆地变化至P 2,则A 等容过程的 W = 0 ,Q =()1121T P P R C igP⎪⎪⎭⎫ ⎝⎛--,错误!未找到引用源。
U =()1121T PP R C igP⎪⎪⎭⎫⎝⎛--,错误!未找到引用源。
H = 1121T P P C ig P ⎪⎪⎭⎫ ⎝⎛-。
B 等温过程的 W =21lnP P RT -,Q =21ln P PRT ,错误!未找到引用源。
U = 0 ,错误!未找到引用源。
H = 0 。
C 绝热过程的 W =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛--11211igPC RigPP P R V P R C ,Q = 0 ,错误!未找到引用源。
化工热力学第三版课后答案完整版朱自强
第二章 流体的压力、体积、浓度关系:状态方程式2-1 试分别用下述方法求出400℃、下甲烷气体的摩尔体积;1 理想气体方程;2 RK 方程;3PR 方程;4 维里截断式2-7;其中B 用Pitzer 的普遍化关联法计算;解 1 根据理想气体状态方程,可求出甲烷气体在理想情况下的摩尔体积id V 为2 用RK 方程求摩尔体积将RK 方程稍加变形,可写为0.5()()RT a V b V b p T pV V b -=+-+ E1其中从附表1查得甲烷的临界温度和压力分别为c T =, c p =,将它们代入a, b 表达式得以理想气体状态方程求得的id V 为初值,代入式E1中迭代求解,第一次迭代得到1V 值为第二次迭代得2V 为353520.563353553313.2217(1.389610 2.984610)1.381102.984610673.154.05310 1.389610(1.389610 2.984610)1.381102.984610 2.1120101.389710V m mol ------------⨯⨯-⨯=⨯+⨯-⨯⨯⨯⨯⨯⨯+⨯=⨯+⨯-⨯=⨯⋅1V 和2V 已经相差很小,可终止迭代;故用RK 方程求得的摩尔体积近似为3用PR 方程求摩尔体积将PR 方程稍加变形,可写为()()()RT a V b V b p pV V b pb V b -=+-++- E2式中 220.45724c cR T a p α=从附表1查得甲烷的ω=;将c T 与ω代入上式 用c p 、c T 和α求a 和b,以RK 方程求得的V 值代入式E2,同时将a 和b 的值也代入该式的右边,藉此求式E2左边的V 值,得563563355353558.314673.152.68012104.053100.10864(1.39010 2.6801210)4.05310[1.39010(1.39010 2.6801210) 2.6801210(1.39010 2.6801210)]1.381102.6801210 1.8217101.3896V ------------⨯=+⨯-⨯⨯⨯-⨯⨯⨯⨯⨯⨯+⨯+⨯⨯⨯-⨯=⨯+⨯-⨯=33110m mol --⨯⋅再按上法迭代一次,V 值仍为3311.389610m mol --⨯⋅,故最后求得甲烷的摩尔体积近似为3311.39010m mol --⨯⋅; 4维里截断式求摩尔体积根据维里截断式2-711()c r c rBp p BpZ RT RT T =+=+ E3 01ccBp B B RT ω=+ E40 1.60.0830.422/r B T =- E5 1 4.20.1390.172/r B T =-E6其中已知甲烷的偏心因子ω=,故由式E4~E6可计算得到从式E3可得 因pVZ RT=,故 四种方法计算得到的甲烷气体的摩尔体积分别为31.38110-⨯、31.39010-⨯、31.39010-⨯和31.39110-⨯31m mol -⋅;其中后三种方法求得的甲烷的摩尔体积基本相等,且与第一种方法求得的值差异也小,这是由于该物系比较接近理想气体的缘故;2-2 含有丙烷的3m 的容器具有的耐压极限;出于安全考虑,规定充进容器的丙烷为127℃,压力不得超过耐压极限的一半;试问可充入容器的丙烷为多少千克解 从附表1查得丙烷的c p 、c T 和ω,分别为,和;则用普遍化压缩因子关联求该物系的压缩因子Z;根据r T 、r p 值,从附表7-2,7-3插值求得:(0)0.911Z = ,(1)0.004Z =,故丙烷的分子量为,即丙烷的摩尔质量M 为0.00441 kg;所以可充进容器的丙烷的质量m 为从计算知,可充9.81 kg 的丙烷;本题也可用合适的EOS 法和其它的普遍化方法求解;2-3 根据RK 方程、SRK 方程和PR 方程,导出其常数a 、b 与临界常数的关系式;解 1RK 方程式,0.5()RT ap V b T V V b =--+ E1利用临界点时临界等温线拐点的特征,即22()()0c c T T T T p pV V==∂∂==∂∂ E2将式E1代入式E2得到两个偏导数方程,即20.52211()0()()c c c c c RT a V b T b V V b -+-=-+E330.53311()0()()c c c c c RT a V b T b V V b --=-+ E4临界点也符合式E1,得0.5()c c c c c c RT ap V b T V V b =--+ E5式E3~E5三个方程中共有a 、b 、c p 、c T 和c V 五个常数,由于c V 的实验值误差较大,通常将其消去,用c p 和c T 来表达a 和b;解法步骤如下:令c c c c p V Z RT =临界压缩因子,即 c c c cZ RTV p =; 同理,令2 2.5a c cR T a p Ω=,b c c RT b p Ω=,a Ω和b Ω为两个待定常数;将a 、b 、c V 的表达式代入式E3~E5,且整理得222(2)1()()a cbc c b c b Z Z Z Z Ω+Ω=+Ω-Ω E622333(33)1()()a cbc b c c b c b Z Z Z Z Z Ω+Ω+Ω=+Ω-Ω E711()a c c b c bZ Z Z Ω=-+Ω-ΩE8式E6除以式E7,式E6除以式E8得3223330c b c b c b Z Z Z -Ω-Ω-Ω=E9322232320c c b c b c b b Z Z Z Z -++Ω-Ω-Ω-Ω=E10对式E8整理后,得()(1)c c b c b a c bZ Z Z Z +Ω-+ΩΩ=-ΩE11式E9减去E10,得22(13)(2)0c b b c c Z Z Z -Ω+Ω-=E12由式E12解得13c Z =,或1)b c Z Ω=此解不一定为最小正根,或1)b c Z Ω=-b Ω不能为负值,宜摒弃再将13c Z =代入式E9或式E10,得32110327b b b Ω+Ω+Ω-=E13解式E13,得最小正根为将13c Z =和0.08664b Ω=代入式E11,得0.42748a Ω=,故2 2.50.42748c cR T a p =E140.08664ccRT b p =E15式E14和式E15即为导出的a 、b 与临界常数的关系式;2 SRK 方程立方型状态方程中的a 、b 与临界常数间的通用关系式可写为 SRK 方程的α是c T 与ω的函数,而RK 方程的0.5r T α=,两者有所区别;至于a Ω与b Ω的求算方法对RK 和SRK 方程一致;因此就可顺利地写出SRK 方程中a 、b与临界常数间的关系式为220.42748c cR T a p α=⋅E160.08664ccRT b p =E173PR 方程由于PR 方程也属于立方型方程,a 、b 与临界常数间的通用关系式仍然适用,但a Ω、b Ω的值却与方程的形式有关,需要重新推导PR 方程由下式表达 因()c T T pV=∂∂=0 22()20()[()()]c c c T T c c c c c RT V b pa V Vb V V b b V b =+∂=-+=∂-++- E18 经简化,上式可写为2222222()()()4()c c c c c c c RT a V b V b V b bV V b +=-++-E19把c c c c Z RT V p =、22a c c cR T a p Ω=、b c c RT b p Ω=代入式E19中,化简得出222222()1()()4()a cbc b c b c b c b Z Z Z Z Z Ω+Ω=-Ω+Ω-Ω-ΩE20对式E18再求导,得22222322322322222222[()4()()(44124)]()()[()4()]c c c c c c c c c c T T c c c c RT a V b bV V b V b V b V bV b pV V b V b bV V b =++--+++-∂=+∂-++- 0= E21将上式化简后得出E22再将c c c c Z RT V p =、22a c c cR T a p Ω=、b c c RT b p Ω=代入式E22中,化简得出432234387263544536278(3121445)1()8208268208a c b c b c b c b c b c b c b c b c b c b c b c b c b Z Z Z Z Z Z Z Z Z Z Z Z Z Ω+Ω+Ω+Ω-Ω=-Ω+Ω+Ω+Ω-Ω-Ω+Ω-Ω+ΩE23 PR 方程的c Z =,将其分别代入式E21和E23后,就可联立解出a Ω与b Ω,得到a Ω=和b Ω=;最后得到2 2.50.45724c cR T a p =和 2-4 反应器的容积为3m ,内有45.40kg 乙醇蒸气,温度为227℃;试用下列四种方法求算反应器的压力;已知实验值为;1RK 方程;2SRK 方程;3PR 方程;4 三参数普遍化关联法;解 1用R-K 方程法计算从附表1查得乙醇的c p 和T c 分别为 和;则RK 方程参数a, b 为 再求乙醇在该状态下的摩尔体积,V 按R-K 方程求算压力,有350.5335668.314(227273.15)28.0391.22910 5.82810500.15 1.229*10(1.22910 5.82810)(3.55190.7925)10 2.75910 2.759Pa MPa-----⨯+=-⨯-⨯⨯⨯⨯+⨯=-⨯=⨯=2用SRK 方程计算从附表1查得乙醇的ω为;SRK 方程中的a 和b 分别计算如下: 在给定条件下乙醇摩尔体积为3311.22910m mol --⨯⋅,将上述有关数值代入SRK 方程,得3用PR 方程计算 将上述数值代入PR 方程,得3533553568.314500.151.22910 5.233410 1.372031.22910(1.22910 5.233410) 5.233410(1.22910 5.233410)(3.53390.83848)10 2.695p Pa MPa--------⨯=⨯-⨯-⨯⨯⨯+⨯+⨯⨯-⨯=-⨯=3用普遍化维里系数法计算根据临界常数和以RK 方程求出的p 为初值,求出对比温度和对比压力,即2.7590.43246.38r c p p p ===, 500.150.9689516.2r c T T T === 故已知乙醇的偏心因子ω=,按下式求压缩因子Z 的值, 所以因和比较接近,不需再迭代;将4种方法计算得到的结果列表比较;由上表知,所用四种方法的计算误差都不大,但以RK 方程法求得的值和实验值最为接近;其余的方法稍差;第一和第四种方法得到的是负偏差,而第二和第三种方法却是正偏差;2-5 某气体的p -V -T 关系可用RK 方程表述,当温度高于c T 时,试推导出以下两个极限斜率的关系式:10lim()T P Z p→∂∂ ;2lim()T P Zp →∞∂∂ ;两式中应包含温度T和RK 方程的常数a 和b;解 根据压缩因子的定义pVZ RT=E1将式E1在恒T 下对p 求偏导,得1()()()T T T Z V p V V p p p RT RT p RT RT V-∂∂∂=+=+∂∂∂ E2根据RK 方程 可求出()T pV∂∂, 20.522(2)()()()T p RT a V b V V b T V V b ∂+=-+∂-+ E3将E3代入E2,得120.522(2)()[]()()T Z V p RT a V b p RT RT V b T V V b -∂+=+-+∂-+ E4pRT也用RK 方程来表达,即 1.51()p a RT V b RT V V b =--+ E5将E5代入E4,得(1) 当0p →,V →∞,故 (2) 当p →∞,V b →,故1、2两种情况下得到的结果即为两个极限斜率的关系式;2-6 试分别用普遍化的RK 方程、SRK 方程和PR 方程求算异丁烷蒸气在350K 、下的压缩因子;已知实验值为;解 1 将RK 方程普遍化,可见原书中的2-20c 和2-20d,即1.51 4.9340() 11r hZ h T h =--+ E10.08664h=rrP ZT E2式E2的右边的Z 以1为初值代入进行迭代,直至得到一收敛的Z 值;由附表1查得异丁烷的c p 、c T 分别为c p = ,c T =,则3500.8576408.1r c T T T ===, 1.20.32883.65r c p P p === 以Z=1代入式E2右边,得 把1h 代入式E1右边,得再把1=0.8346Z 代入式E2,解得2h ,代入式E1,得 按此方法不断迭代,依次得3=0.7965Z , 4=0.7948Z , 5=0.7944Z5Z 和4Z 已非常接近,可终止迭代;异丁烷蒸气的压缩因子为=0.7944Z2 SRK 的普遍化形式如下见原书式2-211 4.934011FhZ h h=--+ E3 0.521[1(1)]r rF m T T =+- E4 20.480 1.5740.176m ωω=+- E50.08664rrp h ZT =E6迭代的过程为:求m 和F 值→取0Z =1→求h 值−−−−→←−−−−循环迭代求Z 值→得收敛的Z 值;查得异丁烷的偏心因子,0.176ω=,故根据式E5和式E4可得 以0Z =1代入式E6右边,得 再由式E3可得 按上述方法,依次可得2=0.7947Z ,3=0.7864Z ,4=0.7843Z ,5=0.7839Z ,6=0.7837Z6Z 和5Z 已非常接近,可终止迭代;故=0.7837Z(3) 用普遍化的PR 方程计算若要按例2-4的思路来计算,必先导出类似于式2-21的普遍化的PR 方程; 令bh V=,则 (1)bV b h h+=+,(1)bV b h h-=-,hZRTp b=将上述4式代入式2-18,并简化后,得(1)(1)(1)RTahZRTp b b b bbh h b h h h h h =-=-++-,即 211[][](1)(1)11(1)(1)hRT a a h Z h b h b RTh h h bRT h h h h h=-=-+---++-+ E7将PR 方程中的a 、b 代入式E7,则1 5.8771[]1(1)(1)r h h T h h hα=--++- E8令0.5220.5211[1(1)][1(0.37464 1.542260.26992)(1)]r r r rF k T T T T ωω=+-=++--,则1 5.8771[]1(1)(1)hZ F h h h h =--++- E9且0.0778/0.0778/0.0778/c c c c rrRT p RT p p b h V V ZRT p ZT ==== E10 通过式E9和E10就可迭代求得Z; 第一次迭代,设0Z =1,则继续迭代,依次可得Z 2=,Z 3=,Z 4=,Z 5=,Z 6=;由于前后两次迭代出的Z 值已很接近,从而得出异丁烷的Z =,与实验值相比,误差为%;由RK 和SRK 方程计算得到的异丁烷的Z 分别为和,它们与实验值的计算误差分别为%和%;可见,三种方法中,普遍化PR 方程计算结果显得更好些;2-7 试用下列三种方法计算250℃、2000Kpa 水蒸气的Z 和V ;1维里截断式2-8,已知B 和C 的实验值分别为310.1525B m kmol -=-⋅和2620.580010C m kmol --=-⨯⋅;2式2-7,其中的B 用Pitzer 普遍化关联法求出;3用水蒸气表计算;解 1用维里截断式2-8计算先求出理想气体状态时的摩尔体积,id V 维里截断式2-8为21pV B CZ RT V V==++ 2-8以id V 为初值,即0id V V =,代入上式右边,求得1V10200(1)B C V V V V =++ E1将1V 再代入式E1右边得同理,3313 2.00710V m mol --=⨯⋅;2V 和3V 很接近,停止迭代,则水蒸气的摩尔体积为3312.00710V m mol --=⨯⋅;所以2用维里截断式2-7计算 维里截断式2-7为11()c r c rBp p BpZ RT RT T =+=+ E201ccBp B B RT ω=+ E3由附表1查得水蒸气的c p 、c T 和ω分别为, 和,则2.00.090722.05r c p p p ===, 250273.150.8082647.3r c T T T +=== 根据Pitzer 的普遍化关联式,有再由式E3和式E2得 故33310.9319 2.17510 2.02710id ZRTV ZV m mol p---===⨯⨯=⨯⋅ 3用水蒸气表计算从水蒸气表附表3查得250℃,2000Kpa 时的水蒸气的比容为 由于水的摩尔质量为,故 同理 2.0080.92322.175id pV V Z RT V ==== 将三种方法计算得到的结果列表比较;计算结果表明,1、3两种方法所得的结果比较接近;2方法偏差较大,主要是忽略了第三维里系数之故;2-8 试用Magoulas 等法、Teja 等法、CG 法和Hu 等法等估算正十九烷的临界温度、临界压力原书中有误,没有计算压缩因子的要求;查阅其文献值,并与所得计算值进行比较;解 正十九烷的分子式为1940C H ,故19c N = 1用Magoulas 等法 按式2-36, 按式2-37,2用Teja 等式按式2-38, 按式2-39,3用CG 法 按式2-40, 按式2-41,4用Hu 等式 按式2-42, 按式2-43,经查阅,c T 、c p 的手册值如下表所示:从上表知,文献中的c T 、c p 手册值并不完全一致,特别c p 间的差值还有些大;由于Nikitin 等的数据和Poling B E 等专着的手册值更为接近,以Nikitin 等的数据为基准手册值,计算出上述各法的误差列于下表;由表知,对c T 、c p 的推算,分别以Magoulas 等法和Hu 等法为最好,且c p 的推算误差比c T 要大;Nikitin 等也给出了c T 和c p 的推算方程如下:据此也可推算正十九烷的c T 和c p ;误差:756754.611000.18%756-⨯= 误差:11.6011.551000.43%11.60-⨯=由Nikitin 等法估算正十九烷的T c ,其误差仅比Magoulas 等法稍差,但比其它三种方法都要优越些;相反,该法估算p c 的误差却最小,比以上四种方法都好,误差要小近半个数量级,甚至更好;由此可见经常查阅文献,与时俱进是很重要的;2-9 试用Constantinou, Gani 和O ’Connell 法估算下列化合物的偏心因子和时液体摩尔体积;1甲乙酮,2环乙烷,3丙烯酸;解 此题如何计算首先要查阅原书P34脚注中的文献4;从该文献中知晓应用何种方程、并查表此两表已在附表9和附表10中给出获得一阶和二阶的数据1i ω、1i υ和2j ω、2j υ等;1甲乙酮应注意到式2-48仅能用于正烷烃的偏心因子估算;对于甲乙酮则应从查阅的文献中得出求算方程;先估算甲乙酮的偏心因子,查得一阶计算的方程为0.50501exp() 1.15070.4085i i N ωω-=∑E1式中,i N 为要估算化合物分子中基团i 出现的次数;1i ω为i 的偏心因子一阶基团贡献值;甲乙酮可分解为3CH 、2CH 和3CH CO 三种基团,从附表9中可以查得1i ω和1i υ,并列表如下:将有关数据代入式E1,得 解得 0.376ω=;从附表1中查得甲乙酮的0.329ω=,0.3290.37610014.28%0.329-=⨯=-误差;一阶估算的误差较大,试进行二阶估算;从文献得出的计算方程为0.505012exp() 1.15070.4085i i j j N A M ωωω-=∑+∑ E2式中 1A =;j M 是在要估算的化合物分子中基团j 出现的次数;2j ω为j 的偏心因子二阶基团贡献值;经查附表10知,甲乙酮的二阶基团为32CH COCH ,其2j ω和2j υ分别为了和31m kmol -⋅;将相关1i ω和2j ω值代入式E2,得0.5050exp() 1.150710.2960210.146911 1.015221(0.20789)0.40851.458150.20789 1.25026ω-=⨯+⨯+⨯+⨯-=-=将上式简化并解得 0.314ω=,0.3290.314100 4.56%0.329-=⨯=误差;从文献查得估算298K 时的l V 估算式为120.01211l i i j j V N A M ωω-=∑+∑E3一阶估算时,0A =,将已查得的各基团的一阶饱和液体贡献值代入式E3,得从化学化工物性数据手册查得甲乙酮在20℃和40℃时的密度分别为3kg m -⋅和3kg m -⋅;内插得25℃时液体的摩尔密度为3kmol m -⋅,则可得出其摩尔体积为31m kmol -⋅;以此为文献值,进行一阶估算结果的误差计算,得二阶估算时,A=1,除1i υ外,尚需要2j υ,以上都已查得备用,依次代入式E3,得2环乙烷偏心因子的一阶估算时,环乙烷可作如下分解,得出基团,并查出基团贡献值:按式E1从附表1查得环乙烷的偏心因子为,0.2130.207100 2.82%0.213-=⨯=误差偏心因子的二阶估算时,从附表10中查得六元环的基团贡献值为,A=1,则按式E2得298K 时环乙烷的摩尔体积按式E3作一阶估算,此时A=0,则从Poling B E 等着的气体物性估算手册中查得时环乙烷的饱和液体摩尔体积为31m kmol -⋅;以此为文献值,则0.108750.11057100 1.67%0.10875-=⨯=-误差;按式E3作二阶估算时,A=1,从附表10中查得六元环的基团贡献值为31m kmol -⋅,因此对环乙烷而言,不论是ω或是l V ,二阶估算的结果都没有一阶估算的精确; 3丙烯酸丙烯酸可分解成如下的基团,并查得其基团贡献值;一阶估算ω,按式E1,从化学化工物性数据手册查得丙烯酸的ω值为,以此为文献值,进行误差计算,二阶估算ω,按式E2,A=1,一阶估算V,按式E3,A=0,l丙烯酸的密度数据来自化学化工物性数据手册,经换算,丙烯酸在25℃时的液体摩尔体积为31⋅,以此为文献值,则m kmol-二阶估算V,按式E3,A=1,l二阶估算结果显示出,ω的估算结果不如一阶的好,而V则相反,二阶估算结l果要比一阶的好;现把全部所得结果示于下表;由表的结果可以得出如下一些看法和观点:aConsfantinou, Gani 和O ’Connell 法预测估算法,对上述三种不同化合物的偏心因子和298K 饱和液体的摩尔体积都比较成功地进行了预测,误差也不算太大,在工程计算中应该有其应用价值;b 从预期来说,二阶估算的结果应该要比一阶估算的好;但从实际估算结果知,并非如此,例如环乙烷的ω和l V 两者的二阶估算结果都比一阶估算结果差;丙烯酸的ω估算,情况也与上述相同;估计出现相仿情况的场合,恐怕为数不少,说明该法应有改进的需要;2-10 估算150℃时乙硫醇的液体的摩尔体积;已知实验值为31m kmol -⋅;乙硫醇的物性参数为c T =499K 、c p =、c V =31m kmol -⋅、ω=,20℃的饱和液体密度为8393kg m -⋅;解 方法1:用Rackett 方程计算液体摩尔体积; Rackett 方程为 其中: 635.49100.207100.27398.314 4.99c c c c p V Z RT -⨯⨯⨯===⨯故0.2857(10.8480)310.207(0.2739)0.0972SL V m kmol --=⨯=⋅乙硫醇的摩尔体积为31m kmol -⋅,该值和实验值31m kmol -⋅相比,误差为%; 方法2:用Lyderson 方法计算由20℃时的饱和液体密度求出此状态的摩尔体积1V ,M 为乙硫醇的摩尔质量,则20℃时的对比温度为根据1r T 值,从图2-11的饱和液体线上查得对比度密度,1r ρ=;同理,根据此值,从图2-11的饱和液体线上查得2 2.15r ρ=;故根据Lyderson 方程,有乙硫醇的摩尔体积计算值为31m kmol -⋅,和实验值相比,误差为%; 2-11 50℃、由摩尔分数的氮和摩尔分数的乙烯组成混合气体,试用下列4种方法求算混合气体的摩尔体积;已知从实验数据, 1.40Z =实;1理想气体方程;2Amagat 定律和普遍化压缩因子图;3 虚拟临界常数法Kay 规则;4 混合物的第二维里系数法;解 1 理想气体方程法根据理想气体方程求得混合物的摩尔体积id m V 为 2 Amagat 定律和普遍化压缩因子图法 根据Amagat 定律()(/)()id id m i i i i i i m m m iiiV V y y Z RT p y Z V Z V ====∑∑∑E1从附表1查得2N 和24C H 的c p 和c T ,2N 1: c p =3.39MPa ,c T =126.2K 24C H 2: c p =5.04MPa ,c T =282.4K根据c p 、c T 值,求出2N 1和24C H 2的r T 和r p 为2N 1:150273.15 2.561126.2r T +==, 160.9717.993.39r p ==24C H 2:250273.15 1.144282.4r T +==, 260.9712.105.04r p ==从普遍化二参数压缩因子图查得相应的i Z 值2N : 1 1.49Z =;24C H :2 1.34Z =代入式E1得3 虚拟临界常数法Kay 规则法根据Kay 规则计算混合物虚拟临界常数, 故可求出混合物的对比温度和对比压力,50273.15 1.470219.8rm T +==, 60.9713.924.38rm p ==根据rm T 和rm p ,查二参数普遍化压缩因子图2-4,得 1.45m Z =,故 4混合物的第二维里系数法 根据式2-712-72e,2211112122222m B y B y y B y B =++E2 01()cij ij ij ij ij cijRT B B B p ω=+E31/2()(0)cij ci c j ij T T T ==这里KE41/31/33()2ci c j cij V V V +=E5()/2cij ci c j Z Z Z =+E6()/2ij i j ωωω=+ E7cij cij cij cijZ RT p V =E80ij B 和1ij B 用Pitzer 的普遍化关联法计算,即0 1.60.0830.422/ij rij B T =- E9 1 4.20.1390.174/ij rij B T =-E10其中 /rij cij T T T =,/rij cij p p p = E11 纯组分的第二维里系数,可按通常的方法求出,即只须用式E3、式E9和式E10,当然此时i=j;而对交叉第二维里系数,须从式E3式E11求出;先从附表1查得各组分的c p 、c T 、c V 、c Z 和ω,具体数值见后面的表1,具体的计算步骤如下: 对2N 1,根据式E11,1111273.1550/ 2.5606126.2r c T T T +===, 111160.97/17.9853.39r c p p p ===根据式E9和E10,011 1.60.4220.0830.010752.5606B =-=-, 1114.20.1740.1390.13572.5606B =-= 代入式E3,得 对24C H 2,根据式E11,2222273.1550/ 1.1443282.4r c T T T +===, 222260.97/12.0975.04r c p p p ===根据式E9和E10,022 1.60.4220.0830.25711.1443B =-=-, 122 4.20.1740.1390.041351.1443B =-= 代入式E3,得交叉第二维里系数12B 的计算如下: 根据式E4式E8, 根据式E11, 代入式E9和E10,012 1.60.4220.0830.095561.7118B =-=-, 112 4.20.1740.1390.12101.7118B =-= 代入式E3得将上述计算结果综合成表1;表1、维里方程计算混合气体的摩尔体积时的一些中间参数i-j T cKp c /MPa V c /m 3kmol -1Z c ω T r B 0 B 1 B/m 3kmol -11-2注:方框中的数值系从附表1查得,其余的分别根据式E3式E11求得;根据式E2求出m B ,得根据维里截断式2-7,求出混合物的压缩因子为若压缩因子为“负值”,意味着摩尔体积为负值;这是没有任何物理意义的,也是不合理的;说明方法4在高达的压力下是不适合的;将四种方法计算结果综合成表2;由表可知,2、3两种方法求出的结果和实验值很接近,而方法1也即理想气体方程求得的结果偏差很大,这是由于系统非理想的缘故;比较2、3两种方法,可以看出2法,也即Amagat 定律,求出的结果为最优;表2、由4种方法计算混合气体的压缩因子和摩尔体积计算方法压缩因子Z m摩尔体积V m/ m 3kmol -1误差 / %实验值计算值1 2 3 4无意义无意义2-12 以化学计量比的2N 和2H 合成氨,在25℃和下,混合气以3311.666710m s --⨯⋅的流速进入反应器;氨的转化率为15%;从反应器出来的气体经冷却和凝缩,将氨分离出后,再行循环;1计算每小时合成氨的量;2若反应器出口的条件为,150℃,求内径为2510m -⨯的出口管中气体的流速;解 先求出2N 1+2H 2混合气体的摩尔体积m V ,拟用Amagat 定律求解; 由附表1分别查得2N 和2H 的c p 、c T 为2N :1 3.39c p MPa =, 1126.2c T K = 2H :2 1.30c p MPa =, 233.2c T K =然后求2N 和2H 的r p 、r T ,2N : 130.3958.9663.39r P ==, 125273.157.237126.2r T +== 2H : 230.39514.401.300.8106r P ==+, 2298.158.96633.28r T ==+根据r P 、r T 查二参数普遍化Z 图得1 1.13Z =,2 1.22Z =因为2N 和2H 是以化学计量比输入,故10.25y =, 20.75y =根据Amagat 定律 故53161.208.314298.159.791030.39510m m Z RT V m mol p --⨯⨯===⨯⋅⨯ 已知混合气体的进口体积流量,331.666710in v m s -=⨯⋅,则混合气体的进口摩尔流速in m 为根据反应的计量关系,22332N H NH −−→+←−−总量 开始 1 3 0 4 结束 330.15-⨯ 20.15⨯ 则每小时合成氨的量可由下式计算得出,2 先求出口气体的组成;因为出口气体中223::(10.15):(330.15):(20.15)N H NH =--⨯⨯,故20.2297N y =,20.6892H y =,30.0811NH y =, 1.000i iy =∑再求出口气体的摩尔流速利用Amagat 定律求出口气体的摩尔体积m V ;先从附表查得3NH 的11.28c p MPa =,405.6c T K =,则可求出各组分的对比性质为2H : 27.8613.201.300.8106r p ==+, 150273.1510.2733.28r T +==+2N : 27.868.2183.39r p ==, 273.151503.353126.2r T +==3NH : 27.86 2.47011.28r p ==, 273.151501.043405.6r T +==根据上述对比参数,查二参数普遍化Z 图,得2 1.15H Z =,2 1.14N Z =,30.380NH Z =则 1.150.0892 1.140.22970.3800.0811 1.085m i i iZ y Z ==⨯+⨯+⨯=∑故43161.0858.314423.15 1.371027.8610m m Z RT V m mol p --⨯⨯===⨯⋅⨯ 出口管中气体的体积流速为出口管中气体的流速,μ,可按下式计算,式中:A 为管子的截面积;计算得出出口管中混合气体的流速为11.10m s -⋅; 58页第2章2-1 求温度、压力的甲烷气体摩尔体积; 解:a 理想气体方程 b 用R-K 方程① 查表求c T 、c p ;② 计算a 、b ;③ 利用迭代法计算V; c 用PR 方程步骤同b,计算结果:1331103893.1--+⋅⋅⋅=mol m V i ; d 利用维里截断式查表可计算r p 、r T 、0B 、1B 和Z 由13310391.1--⋅⋅⋅==⇒=mol m pZRTV RT pV Z 2-2 V=0.5 m 3,耐压 MPa 容器;规定丙烷在T=时,p<;求可充丙烷多少千克 解:a 用理想气体方程136948.815.400314.85.01035.10441.0--⋅⋅=⋅⋅⋅⋅==⇒=⇒=mol m RT MpV m RT M m pV nRT pV b用R-K 方程① 查表求c T 、c p ;② 计算a 、b ;③ 利用迭代法计算V; 则可充丙烷质量计算如下: c 利用维里截断式:查表可计算r p 、r T 、0B 、1B 和Z 由133610257.21035.115.400314.8916.0--⋅⋅⋅=⋅⋅⋅=⇒=mol m V RT pV Z m m 则可充丙烷质量计算如下:2-4 V=1.213 m 3,乙醇45.40 kg,T=,求压力; 解:a 理想气体状态方程 b 用R-K 方程 c 用SRK 方程计算 d 用PR 方程计算 e 用三参数普遍化关联2-6 本题以计算机编程计算提供参考,考试时不能使用aR-K 方程 S-R-K2-7 计算T=,p=2 MPa 的水蒸气的Z 和V 解:a 用维里截断式221pVCRTpV BRT p RT V V C V B RT pV Z ++=⇒++≈=采用迭代法计算V= 之后求得Z= d 利用维里截断式查表可计算r p 、r T 、0B 、1B 可得到Z=; 由13310025.2--⋅⋅⋅==⇒=mol m pZRTV RT pV Z c 水蒸气表 92页第三章 3-4利用三参数压缩因子计算方法,查图表,得到压缩因子: 3-7: 解:注意:J kPa m ⋅=⋅3310 或者 3-9解:乙腈的Antonie 方程为 160℃时,乙腈的蒸气压 2乙腈的标准沸点320℃、40℃和标准沸点时的汽化焓 117页第四章 4-1 4-2 方法一: 经计算得体积流速为:()132210132.02075.0314.32/-⋅⋅=⎪⎭⎫ ⎝⎛⋅⋅=⋅⋅=s m d u V π摩尔流速为:1015.41500000/15.593314.80132.0/-⋅⋅=⋅===s mol p RT V V V n m 根据热力学第一定律,绝热时W s = -△H,所以 方法二:根据过热蒸汽表,内插法应用可查得35kPa 、80℃的乏汽处在过热蒸汽区,其焓值h 2= kJ ·kg -1; 1500 kPa 、320℃的水蒸汽在过热蒸汽区,其焓值h 1= kJ ·kg -1; 4-6 解: 通过112T C H T pmhR+=迭代计算温度,T 2= K 4-7 解:经迭代计算参考101页例题4-3得到T 2=;∆HT1T 2T C p T ()⎛⎜⎜⎠d H 2R T 2()+H 1R-8.32725-103⨯J ⋅mol1-⋅;146页第五章5-1:b 5-2: c 5-4: a 5-5: a 5-1:解:可逆过程熵产为零,即005<∆⇒=--∆=∆-∆=∆sys sys f sys g S T S S S S ; 5-2:解:不可逆过程熵产大于零,即00505T S T S S S S sys sys f sys g ->∆⇒>--∆=∆-∆=∆;即系统熵变可小于零也可大于零; 5-4:解:不可逆绝热过程熵产大于零,即0>∆=∆-∆=∆sys f sys g S S S S ;所以流体熵变大于零; 5-5:解:不可逆过程熵产大于零,即0010010T S T S S S S sys sys f sys g >∆⇒>-∆=∆-∆=∆; 5-3:解:电阻器作为系统,温度维持100℃,即,属于放热;环境温度,属于吸热,根据孤立体系的熵变为系统熵变加环境熵变,可计算如下: 5-6:解:理想气体节流过程即是等焓变化,温度不变,而且过程绝热,所以系统的熵变等于熵产,计算如下: 所以过程不可逆; 5-7: 解:∆S g204.184⋅ln 339363⎛⎝⎫⎪⎭⋅304.184⋅ln 339323⎛ ⎝⎫⎪⎭⋅+0.345kJ ⋅K1-⋅s1-⋅;不同温度的S 值也可以直接用饱和水表查得;计算结果是; 5-12解:1循环的热效率 2 水泵功与透平功之比H 2= kJ ·kg -1,H 3= kJ ·kg -1,H 4= kJ ·kg -1,H 5= kJ ·kg -1,3 提供1 kw 电功的蒸气循环量 5-15题: 194页第六章 6-1:解:水蒸气的摩尔流量为:nm M 3600⋅16801000⋅183600⋅25.926mol ⋅s1-⋅a 通过内插法求出时对应的温度,如下b 6-3 6-6:解:理想气体经一锐孔降压过程为节流过程,0=∆H ,且0=Q ,故0=S W ,过程恒温;则绝热膨胀过程的理想功和损耗功计算如下: 6-8:解:1产品是纯氮和纯氧时,2产品是98% N 2和50% O 2的空气时,设计计算流程如下:总的功6-12:解:6-13解:由1pmh 31()2pmh 32(),可得3 使用内插法可求得66.03℃时的熵值, 1利用熵分析法计算损耗功,2利用火用分析法:或者241页第七章7-2解:假设需水m kg,则产品酒中含水产品酒中含醇所以酒的体积7-3解:7-498% N 2 50% O 2解:根据吉布斯-杜亥姆公式,恒温恒压时 则有0=∑ii i V d x ,所以所以设计的方程不合理;。
化工热力学第三版(完全版)课后习题答案(I).doc
化工热力学课后答案第1章 绪言一、是否题1. 封闭体系的体积为一常数。
(错)2. 封闭体系中有两个相βα,。
在尚未达到平衡时,βα,两个相都是均相敞开体系;达到平衡时,则βα,两个相都等价于均相封闭体系。
(对)3. 理想气体的焓和热容仅是温度的函数。
(对)4. 理想气体的熵和吉氏函数仅是温度的函数。
(错。
还与压力或摩尔体积有关。
)5. 封闭体系的1mol 气体进行了某一过程,其体积总是变化着的,但是初态和终态的体积相等,初态和终态的温度分别为T 1和T 2,则该过程的⎰=21T T V dT C U ∆;同样,对于初、终态压力相等的过程有⎰=21T T P dT C H ∆。
(对。
状态函数的变化仅决定于初、终态与途径无关。
)二、填空题1. 状态函数的特点是:状态函数的变化与途径无关,仅决定于初、终态 。
2. 封闭体系中,温度是T 的1mol 理想气体从(P i ,V i )等温可逆地膨胀到(P f ,V f ),则所做的功为()f i rev V V RT W ln =(以V 表示)或()i f rev P P RT W ln = (以P 表示)。
3. 封闭体系中的1mol 理想气体(已知igP C ),按下列途径由T 1、P 1和V 1可逆地变化至P 2,则A 等容过程的 W = 0 ,Q =()1121T P P R C igP⎪⎪⎭⎫ ⎝⎛--,错误!未找到引用源。
U =()1121T PP R C igP⎪⎪⎭⎫⎝⎛--,错误!未找到引用源。
H = 1121T P P C ig P ⎪⎪⎭⎫ ⎝⎛-。
B 等温过程的 W =21lnP P RT -,Q =21ln P PRT ,错误!未找到引用源。
U = 0 ,错误!未找到引用源。
H = 0 。
C 绝热过程的 W =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛--11211igPC RigPP P R V P R C ,Q = 0 ,错误!未找到引用源。
化工热力学第三版(完全版)课后习题答案
化工热力学课后答案第1章 绪言一、是否题1. 封闭体系的体积为一常数。
(错)2. 封闭体系中有两个相βα,。
在尚未达到平衡时,βα,两个相都是均相敞开体系;达到平衡时,则βα,两个相都等价于均相封闭体系。
(对)3. 理想气体的焓和热容仅是温度的函数。
(对)4. 理想气体的熵和吉氏函数仅是温度的函数。
(错。
还与压力或摩尔体积有关。
)5. 封闭体系的1mol 气体进行了某一过程,其体积总是变化着的,但是初态和终态的体积相等,初态和终态的温度分别为T 1和T 2,则该过程的⎰=21T T V dT C U ∆;同样,对于初、终态压力相等的过程有⎰=21T T P dT C H ∆。
(对。
状态函数的变化仅决定于初、终态与途径无关。
)二、填空题1. 状态函数的特点是:状态函数的变化与途径无关,仅决定于初、终态 。
2. 封闭体系中,温度是T 的1mol 理想气体从(P i ,V i )等温可逆地膨胀到(P f ,V f ),则所做的功为()f i rev V V RT W ln =(以V 表示)或()i f rev P P RT W ln = (以P 表示)。
3. 封闭体系中的1mol 理想气体(已知igP C ),按下列途径由T 1、P 1和V 1可逆地变化至P 2,则A 等容过程的 W = 0 ,Q =()1121T P P R C igP ⎪⎪⎭⎫ ⎝⎛--,U =()1121T PPR C igP ⎪⎪⎭⎫⎝⎛--,H =1121T P P C ig P ⎪⎪⎭⎫ ⎝⎛-。
B 等温过程的 W =21lnP P RT -,Q =21ln P PRT ,U = 0 ,H = 0 。
C 绝热过程的 W =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛--11211igPC RigPP P R V P R C ,Q = 0 ,U =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛-11211ig P C R igPP P R V P R C ,H =1121T P P C igP C R ig P ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛。
化工热力学(第三版)答案
化工热力学(第三版)第二章 流体的压力、体积、浓度关系:状态方程式2-1 试分别用下述方法求出400℃、4.053MPa 下甲烷气体的摩尔体积。
(1) 理想气体方程;(2) RK 方程;(3)PR 方程;(4) 维里截断式(2-7)。
其中B 用Pitzer 的普遍化关联法计算。
[解] (1) 根据理想气体状态方程,可求出甲烷气体在理想情况下的摩尔体积id V 为33168.314(400273.15)1.381104.05310idRT Vm molp--⨯+===⨯⋅⨯(2) 用RK 方程求摩尔体积将RK 方程稍加变形,可写为0.5()()RT a V b V b p TpV V b -=+-+ (E1)其中22.50.427480.08664ccccR T a p RT b p ==从附表1查得甲烷的临界温度和压力分别为c T =190.6K, c p =4.60MPa ,将它们代入a, b 表达式得2 2.56-20.560.427488.314190.63.2217m Pa mol K4.6010a ⨯⨯==⋅⋅⋅⨯53160.086648.314190.62.9846104.6010b m m ol --⨯⨯==⨯⋅⨯以理想气体状态方程求得的idV为初值,代入式(E1)中迭代求解,第一次迭代得到1V 值为5168.314673.15 2.9846104.05310V -⨯=+⨯⨯350.563353.2217(1.381102.984610)673.154.05310 1.38110(1.381102.984610)-----⨯⨯-⨯-⨯⨯⨯⨯⨯⨯+⨯3553311.381102.9846102.1246101.389610m m ol-----=⨯+⨯-⨯=⨯⋅第二次迭代得2V 为353520.563353553313.2217(1.3896102.984610)1.381102.984610673.154.05310 1.389610(1.389610 2.984610)1.381102.9846102.1120101.389710V m m ol------------⨯⨯-⨯=⨯+⨯-⨯⨯⨯⨯⨯⨯+⨯=⨯+⨯-⨯=⨯⋅1V 和2V 已经相差很小,可终止迭代。
化工热力学习题参考答案 第三版
化工热力学习题参考答案第三版化工热力学习题参考答案第三版化工热力学习题参考答案(第三版)第二章p582-1求温度673.15k、压力4.053mpa的甲烷气体摩尔体积。
解:(a)理想气体方程 pv rt v8.314673.154.053101.38110(b)用r-k方程①换算谋tc、pc;②排序a、b;③利用迭代法排序v。
rtv brtpat v v ba vib t vi vi b m molv0 1.38110vi1 1.389610(c)用pr方程步骤同(b),计算结果:vi1(d)利用维里截断式bprt0.422tr1.3893100.0830.1390.172tr换算可以排序pr、tr、b0、b1和z由z1.391102-2v=0.5m3,耐压2.7mpa容器。
规定丙烷在t=400.15k时,p0.0441 1.35108.314400.158.948m(b)用r-k方程①查表求tc、pc;②计算a、b;③利用迭代法计算v。
rtv brtpat v v ba vib t vi vi b mv0 2.46410vi1 2.24110则可充丙烷质量计算如下:0.04410.52.241109.838kg(c)利用维里切断式:0.0830.422tr0.1390.172tr查表可计算pr、tr、b0、b1和z由z0.9168.314400.151.35102.25710则可充丙烷质量计算如下:0.04410.52.257109.77kg2-4v=1.213m3,乙醇45.40kg,t=500.15k,谋压力。
求解:(a)理想气体状态方程45.40468.314500.153.383mpa(b)用r-k方程0.42748r tc0.08664r tca28.039rtvm bvnt vm vm b1.229m1.21345.40/46p 2.759mpa(c)用srk方程计算(d)用pr方程排序(e)用三参数普遍化关联mmvnpvrt45.4461.2130.98710.987kmol 1.229mbprtvm z bpkmolrt0.361,b0.057,0.635b0.267p2.779mpa2-6求解:(1)将rk方程普遍化,可知原书中的(2-20c)和(2-20d),即为4.9340tr0.08664pr式(e2)的右边的z以1为初值代入进行迭代,直至得到一收敛的z值。
化工热力学第三版课后习题答案
《化工热力学》(第三版)习题参考答案58页第2章2-1 求温度673.15K 、压力4.053MPa 的甲烷气体摩尔体积。
解:(a )理想气体方程133610381.110053.415.673314.8--⋅⋅⋅=⋅⋅==⇒=molm p RT V RT pV(b )用R-K 方程① 查表求c T 、c p ;② 计算a 、b ;③ 利用迭代法计算V 。
()()()133113301103896.110381.1--+--+⋅⋅⋅=⋅⋅⋅⋅⋅⋅=+⋅⋅--+=+⋅⋅--=molm V molm V b V V T b V a b p RT V b V V T a b V RT p i i i i i(c )用PR 方程步骤同(b ),计算结果:1331103893.1--+⋅⋅⋅=molm V i 。
(d )利用维里截断式2.416.101172.0139.0422.0083.0111rrrr rr rr cc T B T BT p B T p B T p RT Bp RT Bp RTpV Z -=-=⋅⋅+⋅+=⋅+=+==ω查表可计算r p 、r T 、0B 、1B 和Z 由13310391.1--⋅⋅⋅==⇒=molm pZRT V RTpV Z2-2 V=0.5 m 3,耐压2.7 MPa 容器。
规定丙烷在T=400.15K 时,p<1.35MPa 。
求可充丙烷多少千克?解:(a )用理想气体方程 136948.815.400314.85.01035.10441.0--⋅⋅=⋅⋅⋅⋅==⇒=⇒=molm RTMpV m RT Mm pV nRT pV (b )用R-K 方程① 查表求c T 、c p ;② 计算a 、b ;③ 利用迭代法计算V 。
()()()13311330110241.210464.2--+--+⋅⋅⋅=⋅⋅⋅⋅⋅⋅=+⋅⋅--+=+⋅⋅--=molm V molm V b V V T b V a b pRT V b V V T a b V RT p i i i i i则可充丙烷质量计算如下:kg M V V M n m i ⋅=⋅⋅=⋅=⋅=-+838.910241.25.00441.031(c )利用维里截断式:2.416.101172.0139.0422.0083.0111rrrr rr rr cc m T B T BT p B T p B T p RT Bp RTBp RTpV Z -=-=⋅⋅+⋅+=⋅+=+==ω查表可计算r p 、r T 、0B 、1B 和Z 由133610257.21035.115.400314.8916.0--⋅⋅⋅=⋅⋅⋅=⇒=molm V RTpV Z m m则可充丙烷质量计算如下:kg M V V M n m i ⋅=⋅⋅=⋅=⋅=-+77.910257.25.00441.0312-4 V=1.213 m 3,乙醇45.40 kg ,T=500.15K ,求压力。
化工热力学第三版(完全版)课后习题问题详解
化工热力学课后答案第1章 绪言一、是否题1. 封闭体系的体积为一常数。
(错)2. 封闭体系中有两个相βα,。
在尚未达到平衡时,βα,两个相都是均相敞开体系;达到平衡时,则βα,两个相都等价于均相封闭体系。
(对) 3. 理想气体的焓和热容仅是温度的函数。
(对)4. 理想气体的熵和吉氏函数仅是温度的函数。
(错。
还与压力或摩尔体积有关。
)5. 封闭体系的1mol 气体进行了某一过程,其体积总是变化着的,但是初态和终态的体积相等,初态和终态的温度分别为T 1和T 2,则该过程的⎰=21T T V dT C U ∆;同样,对于初、终态压力相等的过程有⎰=21T T P dT C H ∆。
(对。
状态函数的变化仅决定于初、终态与途径无关。
) 二、填空题1. 状态函数的特点是:状态函数的变化与途径无关,仅决定于初、终态 。
2. 封闭体系中,温度是T 的1mol 理想气体从(P i ,V i )等温可逆地膨胀到(P f ,V f ),则所做的功为()f i rev V V RT W ln =(以V 表示)或()i f rev P P RT W ln = (以P 表示)。
3. 封闭体系中的1mol 理想气体(已知ig P C ),按下列途径由T 1、P 1和V 1可逆地变化至P 2,则A 等容过程的 W = 0 ,Q =()1121T P P R C igP ⎪⎪⎭⎫ ⎝⎛--,错误!未找到引用源。
U =()1121T PPR C igP ⎪⎪⎭⎫⎝⎛--,错误!未找到引用源。
H = 1121T P P C ig P ⎪⎪⎭⎫ ⎝⎛-。
B 等温过程的 W =21lnP P RT -,Q =21ln P PRT ,错误!未找到引用源。
U = 0 ,错误!未找到引用源。
H = 0 。
C 绝热过程的 W =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛--11211igPC RigPP P R V P R C ,Q = 0 ,错误!未找到引用源。
化工热力学第三版(完全版)课后习题答案
临界参数 Tc=425.4K, Pc=3.797MPa, ω=0.193
修正的 Rackett 方程常数: α=0.2726, β=0.0003
ln P S 6.8146 2151.63 36.24 T
P S 0.504 MPa
由软件计算知 V sl 103.0193cm3mol 1, V sv 4757.469cm3mol 1
化工热力学课后答案
第 1 章 绪言
一、是否题
1. 封闭体系的体积为一常数。 (错)
2. 封闭体系中有两个相 , 。在尚未达到平衡时, , 两个相都是均相敞开体系;
达到平衡时,则 , 两个相都等价于均相封闭体系。 (对)
3. 理想气体的焓和热容仅是温度的函数。 (对)
4. 理想气体的熵和吉氏函数仅是温度的函数。 (错。还与压力或摩尔体积有关。 )
Wrev PdV
C
ig P
R dT
RT dV
V
a bT cT 2 R dT Rd ln V 0
T
T2 a R
T1
T
b cT dT
R ln V2 V1
0,又 V 2 V1
P1 T2 ,故 P2 T1
a ln T2 b T2 T1 T1
c
T
2 2
T12
2
R ln P2 0 P1
3. 一个 0.057m 3气瓶中贮有的 1MPa 和 294 K的高压气体通过一半开的阀门放入一个压力
33
4. 对于三混合物,展开 PR 方程常数 a的表达式, a
yi y j aii a jj (1 k ij ) =
i1 j 1
y12 a1 y 22a 2 y32 a 3 2 y1 y 2 a1 a2 1 k12 2 y2 y3 a2 a 3 1 k 23 2 y3 y1 a3 a1 1 k31 ,其
化工热力学第三版(完全版)课后习题答案
化工热力学课后答案第1章 绪言一、是否题1. 封闭体系的体积为一常数。
(错)2. 封闭体系中有两个相βα,。
在尚未达到平衡时,βα,两个相都是均相敞开体系;达到平衡时,则βα,两个相都等价于均相封闭体系。
(对)3. 理想气体的焓和热容仅是温度的函数。
(对)4. 理想气体的熵和吉氏函数仅是温度的函数。
(错。
还与压力或摩尔体积有关。
)5. 封闭体系的1mol 气体进行了某一过程,其体积总是变化着的,但是初态和终态的体积相等,初态和终态的温度分别为T 1和T 2,则该过程的⎰=21T T V dT C U ∆;同样,对于初、终态压力相等的过程有⎰=21T T P dT C H ∆。
(对。
状态函数的变化仅决定于初、终态与途径无关。
)二、填空题1. 状态函数的特点是:状态函数的变化与途径无关,仅决定于初、终态 。
2. 封闭体系中,温度是T 的1mol 理想气体从(P i ,V i )等温可逆地膨胀到(P f ,V f ),则所做的功为()f i rev V V RT W ln =(以V 表示)或()i f rev P P RT W ln = (以P 表示)。
3. 封闭体系中的1mol 理想气体(已知igP C ),按下列途径由T 1、P 1和V 1可逆地变化至P 2,则A 等容过程的 W = 0 ,Q =()1121T P P R C igP⎪⎪⎭⎫ ⎝⎛--,错误!未找到引用源。
U =()1121T PP R C igP⎪⎪⎭⎫⎝⎛--,错误!未找到引用源。
H = 1121T P P C ig P ⎪⎪⎭⎫ ⎝⎛-。
B 等温过程的 W =21lnP P RT -,Q =21ln P PRT ,错误!未找到引用源。
U = 0 ,错误!未找到引用源。
H = 0 。
C 绝热过程的 W =()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎪⎪⎭⎫ ⎝⎛--11211igPC RigPP P R V P R C ,Q = 0 ,错误!未找到引用源。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第二章 流体的压力、体积、浓度关系:状态方程式2-1 试分别用下述方法求出400℃、4.053MPa 下甲烷气体的摩尔体积。
(1) 理想气体方程;(2) RK 方程;(3)PR 方程;(4) 维里截断式(2-7)。
其中B 用Pitzer 的普遍化关联法计算。
[解] (1) 根据理想气体状态方程,可求出甲烷气体在理想情况下的摩尔体积idV 为33168.314(400273.15)1.381104.05310id RT V m mol p --⨯+===⨯⋅⨯ (2) 用RK 方程求摩尔体积将RK 方程稍加变形,可写为0.5()()RT a V b V b p T pV V b -=+-+ (E1)其中2 2.50.427480.08664c c c cR T a p RT b p ==从附表1查得甲烷的临界温度和压力分别为c T =190.6K, c p =4.60MPa ,将它们代入a, b 表达式得2 2.56-20.560.427488.314190.6 3.2217m Pa mol K 4.6010a ⨯⨯==⋅⋅⋅⨯ 53160.086648.314190.6 2.9846104.6010b m mol --⨯⨯==⨯⋅⨯ 以理想气体状态方程求得的idV 为初值,代入式(E1)中迭代求解,第一次迭代得到1V 值为5168.314673.152.9846104.05310V -⨯=+⨯⨯ 350.563353.2217(1.38110 2.984610)673.15 4.05310 1.38110(1.38110 2.984610)-----⨯⨯-⨯-⨯⨯⨯⨯⨯⨯+⨯ 3553311.381102.984610 2.1246101.389610m mol -----=⨯+⨯-⨯=⨯⋅ 第二次迭代得2V 为353520.563353553313.2217(1.389610 2.984610)1.381102.984610673.154.05310 1.389610(1.389610 2.984610)1.381102.984610 2.1120101.389710V m mol ------------⨯⨯-⨯=⨯+⨯-⨯⨯⨯⨯⨯⨯+⨯=⨯+⨯-⨯=⨯⋅1V 和2V 已经相差很小,可终止迭代。
故用RK 方程求得的摩尔体积近似为3311.39010V m mol --=⨯⋅(3)用PR 方程求摩尔体积将PR 方程稍加变形,可写为()()()RT a V b V b p pV V b pb V b -=+-++- (E2)式中 220.45724c cR T a p α=0.07780ccRT b p = 0.520.51(0.37464 1.542260.26992)(1)r T αωω=++--从附表1查得甲烷的ω=0.008。
将c T 与ω代入上式0.520.5673.151(0.37464 1.542260.0080.269920.008)(1())190.60.659747α=++⨯-⨯-= 0.435266α=用c p 、c T 和α求a 和b ,226268.314190.60.457240.4352660.108644.6010a m Pa mol -⨯=⨯=⋅⋅⨯ 53168.314190.60.077802.68012104.6010b m mol --⨯==⨯⋅⨯ 以RK 方程求得的V 值代入式(E2),同时将a 和b 的值也代入该式的右边,藉此求式(E2)左边的V 值,得563563355353558.314673.15 2.68012104.053100.10864(1.39010 2.6801210)4.05310[1.39010(1.39010 2.6801210) 2.6801210(1.39010 2.6801210)]1.38110 2.6801210 1.8217101.3896V ------------⨯=+⨯-⨯⨯⨯-⨯⨯⨯⨯⨯⨯+⨯+⨯⨯⨯-⨯=⨯+⨯-⨯=33110m mol --⨯⋅再按上法迭代一次,V 值仍为3311.389610m mol --⨯⋅,故最后求得甲烷的摩尔体积近似为3311.39010m mol --⨯⋅。
(4)维里截断式求摩尔体积根据维里截断式(2-7)11()c r c rBp p BpZ RT RT T =+=+ (E3)01ccBp B B RT ω=+(E4)0 1.60.0830.422/r B T =- (E5) 1 4.20.1390.172/r B T =-(E6)其中673.15 3.5317190.6r c T T T === 4.0530.88114.60r c p p p === 已知甲烷的偏心因子ω=0.008,故由式(E4)~(E6)可计算得到0 1.60.0830.422/3.53170.02696B =-= 1 4.20.1390.172/3.53170.1381B =-= 0.026960.0080.13810.02806ccBp RT =+⨯= 从式(E3)可得0.881110.02806 1.0073.5317Z =+⨯=因pVZ RT=,故 33311.007 1.38110 1.39110id ZRTV ZV m mol p---===⨯⨯=⨯⋅ 四种方法计算得到的甲烷气体的摩尔体积分别为31.38110-⨯、31.39010-⨯、31.39010-⨯和31.39110-⨯31m mol -⋅。
其中后三种方法求得的甲烷的摩尔体积基本相等,且与第一种方法求得的值差异也小,这是由于该物系比较接近理想气体的缘故。
2-2 含有丙烷的0.53m 的容器具有2.7Mpa 的耐压极限。
出于安全考虑,规定充进容器的丙烷为127℃,压力不得超过耐压极限的一半。
试问可充入容器的丙烷为多少千克?[解] 从附表1查得丙烷的c p 、c T 和ω,分别为4.25MPa ,369.8K 和0.152。
则127373.15 1.08369.8r c T T T +=== 2.70.3184.252r c p p p ===⨯ 用普遍化压缩因子关联求该物系的压缩因子Z 。
根据r T 、r p 值,从附表(7-2),(7-3)插值求得:(0)0.911Z = ,(1)0.004Z =,故 (0)(1)0.9110.1520.0040.912Z Z Z ω=+=+⨯=丙烷的分子量为44.1,即丙烷的摩尔质量M 为0.00441 kg 。
所以可充进容器的丙烷的质量m 为61.35100.50.04419.810.9128.314(127373.15)tpV m M ZRTkg =⋅⨯⨯⨯==⨯⨯+从计算知,可充9.81 kg 的丙烷。
本题也可用合适的EOS 法和其它的普遍化方法求解。
2-3 根据RK 方程、SRK 方程和PR 方程,导出其常数a 、b 与临界常数的关系式。
[解] (1)RK 方程式,0.5()RT ap V b T V V b =--+ (E1)利用临界点时临界等温线拐点的特征,即22()()0c c T T T T p pV V==∂∂==∂∂ (E2)将式(E1)代入式(E2)得到两个偏导数方程,即20.52211()0()()c c c c c RT a V b T b V V b -+-=-+ (E3)30.53311()0()()c c c c c RT a V b T b V V b --=-+(E4)临界点也符合式(E1),得0.5()c c c c c c RT ap V b T V V b =--+ (E5)式(E3)~(E5)三个方程中共有a 、b 、c p 、c T 和c V 五个常数,由于c V 的实验值误差较大,通常将其消去,用c p 和c T 来表达a 和b 。
解法步骤如下:令c c c c p V Z RT =(临界压缩因子),即 c c c cZ RTV p =。
同理,令2 2.5a c cR T a p Ω=,b c c RT b p Ω=,a Ω和b Ω为两个待定常数。
将a 、b 、c V 的表达式代入式(E3)~(E5),且整理得222(2)1()()a c b c c b c b Z Z Z Z Ω+Ω=+Ω-Ω(E6)22333(33)1()()a c b c b c c b c b Z Z Z Z Z Ω+Ω+Ω=+Ω-Ω (E7)11()a c c b c bZ Z Z Ω=-+Ω-Ω(E8)式(E6)除以式(E7),式(E6)除以式(E8)得3223330c b c b c b Z Z Z -Ω-Ω-Ω=(E9)322232320c c b c b c b b Z Z Z Z -++Ω-Ω-Ω-Ω=(E10)对式(E8)整理后,得()(1)c c b c b a c bZ Z Z Z +Ω-+ΩΩ=-Ω(E11)式(E9)减去(E10),得22(13)(2)0c b b c c Z Z Z -Ω+Ω-=(E12)由式(E12)解得13c Z =,或1)b c Z Ω=(此解不一定为最小正根),或1)b c Z Ω=-(b Ω不能为负值,宜摒弃)再将13c Z =代入式(E9)或式(E10),得 32110327b b b Ω+Ω+Ω-=(E13)解式(E13),得最小正根为0.08664b Ω=将13c Z =和0.08664b Ω=代入式(E11),得0.42748a Ω=,故 2 2.50.42748c cR T a p =(E14)0.08664ccRT b p =(E15)式(E14)和式(E15)即为导出的a 、b 与临界常数的关系式。
(2) SRK 方程立方型状态方程中的a 、b 与临界常数间的通用关系式可写为22c a c ccb cR T a a p RT b p αα=⋅Ω==⋅ΩSRK 方程的α是c T 与ω的函数,而RK 方程的0.5r T α=,两者有所区别。
至于a Ω与b Ω的求算方法对RK 和SRK 方程一致。
因此就可顺利地写出SRK 方程中a 、b 与临界常数间的关系式为220.42748c cR T a p α=⋅(E16)0.08664ccRT b p =(E17)(3)PR 方程由于PR 方程也属于立方型方程,a 、b 与临界常数间的通用关系式仍然适用,但a Ω、b Ω的值却与方程的形式有关,需要重新推导PR 方程由下式表达()()RT ap V b V V b b V b =--++- 因()c T T pV=∂∂=022()20()[()()]c c c T T c c c c c RT V b pa V Vb V V b b V b =+∂=-+=∂-++- (E18) 经简化,上式可写为2222222()()()4()c c c c c c c RT a V b V b V b bV V b +=-++- (E19)把c c c c Z RT V p =、22a c c cR T a p Ω=、b c c RT b p Ω=代入式(E19)中,化简得出222222()1()()4()a cbc b c b c b c b Z Z Z Z Z Ω+Ω=-Ω+Ω-Ω-Ω(E20)对式(E18)再求导,得22222322322322222222[()4()()(44124)]()()[()4()]c c c c c c c c c c T T c c c c RT a V b bV V b V b V b V bV b pV V b V b bV V b =++--+++-∂=+∂-++- 0=(E21)将上式化简后得出43223438726354453627822(3121445)()8208268208c c c c c c c c c c c c c c c RT a V bV b V b V b V b V bV b V b V b V b V b V bV b +++-=-+++--+-+ (E22)再将c c c c Z RT V p =、22a c c cR T a p Ω=、b c c RT b p Ω=代入式(E22)中,化简得出432234387263544536278(3121445)1()8208268208a c b c b c b c b c b c b c b c b c b c b c b c b c b Z Z Z Z Z Z Z Z Z Z Z Z Z Ω+Ω+Ω+Ω-Ω=-Ω+Ω+Ω+Ω-Ω-Ω+Ω-Ω+Ω(E23)PR 方程的c Z =0.3074,将其分别代入式(E21)和(E23)后,就可联立解出a Ω与b Ω,得到a Ω=0.45724和b Ω=0.0778。