冶金物理化学课件热力学-第七章

《冶金物理化学课件热力学-第七章》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第七章多相多元系平衡计算7.1几个基本问题7.2化学平衡法7.3最小自由能法（White法）为什么要进行多相多元系的平衡计算？在过去的的热力学平衡计算中，判断一个反应能否发生，我们总是计算单独这一个反应的，实际上在做这项工作的同时，忽略了同一体系中的其他组元对这个反应的影响，也就对我们的判断的正误产生了影响。最典型的例子是假如我们由热力学判断的这个反应是不能进行的，但有这个反应的一个偶合反应存在的情况下，这个反应是可以进行的（见本书8.3.1），若我们不考虑偶合反应的存在，所判断的反应发生与否就是不确定的；若考虑偶合反应的存在，就必须同时计算两个以上的反应的自由能。实际上

2、在一个多元多相热力学体系中，反应与反应之间都存在这或多或少的“偶合”，又如炼钢过程中的氧化反应，好几个元素都在和氧反应，况且都是同时，单独计算一个氧化反应，其结果是值得怀疑的。综上所述，单独研究发生在一个多元多相体系中的反应实际上意义不大，应该同时研究这个体系中全部反应的平衡结果。这在过去计算机不是很发达是，计算是很困难的，目前这项工作已经变的很容易了，所以我们要改变观念。多元多相体系的平衡研究，已经有大量的研究成果，从以下方面分类：平衡常数法（Bbinkey为代表）最小自由能原理（White为代表）从热力学原理分计算方法的特征计算方法的收敛速度计算程序a．求解非线性方程

3、组b．梯度法c．最优化法专用程序通用程序一阶法二阶法……P阶法按计算方法分7.1几个基本问题7.1.1独立反应数1）必要条件：在一个多元多相体系中，如果：①体系中各独立反应间线性无关；②描述每一独立反应的反应进程可独立进行；③由一组独立反应式可以描述整个体系的组成变化。这是描述独立反应数的必要条件。2）确定方法设一体系中存在C个组分，分别为A1、A2、…、AC。其中第j个化学反应可以表达成如下通式：(7-1)代表i组分在第j个反应中的化学计量系数，规定i组分为反应物时为负，为生成物时为正，若存在S个反应，则……(7-2)其系数矩阵为求的秩，若为r，即为该多组元多相体系的独

4、立反应数。(7-3)7.1.2反应进度（ExtentofReaction）反应进度，早期由T．de．Donder引入，IUPAC推荐。在反应热，化学平衡，反应速度中普遍采用。对反应j定义或（7-4）def即或为反应进度。其中ni为Ai的摩尔数；ni0为Ai的初始摩尔数。7.1.2反应进度（ExtentofReaction）反应进度，早期由T．de．Donder引入，IUPAC推荐。在反应热，化学平衡，反应速度中普遍采用。对反应j7.2化学平衡法7.2.1热力学原理设体系中存在C个组分，T、P下平衡时，有r个独立化学反应j=1，2，…，r（7-5）则j=1，2，…，r（7-

5、6）若体系有m元素，则存在m个物料平衡条件e=1，2，…，m。（7-7）其中ai,e为组分i中元素e的原子系数（例如：若组元i为Fe2O3，Fe元素的编号为1，O元素的编号为2。则ai1=2，ai2=3，）；be为体系中元素e的总物质的量。联立求解方程式（7-6），（7-7）（共r+m个），可求出平衡时各组元的浓度。解的存在性讨论：由热力学原理，c－m=r所以c=r+m对上述联立方程式，ni（i＝1，2，…，c）为c个变量，有c=r+m个独立方程；为活度，活度系数可以由其他方法确定；vji，Gi0，R，T，ai,e，be皆为常数。上式可求解。例7－1：对于CO、H2、H2

6、O、CO2、CH4体系，c=5，m=3（C、H、O）独立反应数r=c－m=2，如下：热力学模型若,（i=CO、H2、H2O、CO2、CH4）5个未知数，5个方程，可求解。以上已经建立了求解多元多相体系化学反应的平衡浓度的热力学模型，模型怎么求出具体的解？7.2.2Newton-Raphson法二阶非线性方程组(2-7-8)设式2-7-8的一组近似解为，则有设f1、f2对x、y的二阶偏导数存在，利用泰勒展开或(2-7-9)(2-7-10)方程组2-7-10的系数矩阵为若其非奇异，则可解出∆x、∆y，于是以x1，y1代替x0，y0，进行迭代，可得一系列近似解n=1,2,…。若

7、相邻两组近似解满足条件其中式中ε1---允许误差C----控制常数，通常取1(2)N阶非线性方程组或(2-7-11)第k次迭代后(K=0,1,2,…)设x0为方程的近似解，且方程在x0附近二阶可微，将f(x)用泰勒公式展开：Jacobi迭代矩阵：若此矩阵非奇异，则方程2-7-11有唯一解迭代求解可得一系列近似解x1，x2，…，xn。(n=0,1,2,…)误差函数为当时，各也应趋于零。当初始值充分接近于解时，N-R法迭代按平方收敛速度收敛：K为常数（3）N-R法的优缺点优点：具有良好的收敛性缺点：每次迭代都要计算Jacobi迭代