热分析动力学.doc

热分析动力学 一、 基本方程 对于常见的固相反应来说，其反应方程可以表示为 （1） 其反应速度可以用两种不同形式的方程表示 微分形式 （2） 和 积分形式 （3） 式中αt时物质A已反应的分数； t时间； k反应速率常数； fα反应机理函数的微分形式； Gα反应机理函数的积分形式。 由于f（α）和G（α）分别为机理函数的微分形式和积分形式，它们之间的关系为 （4） k与反应温度T（绝对温度）之间的关系可用著名的Arrhenius方程表示 （5） 式中A表观指前因子； E表观活化能； R通用气体常数。 方程（2）～（5）是在等温条件下出来的，将这些方程应用于非等温条件时，有如下关系式 （6） 即 式中T0DSC曲线偏离基线的始点温度（K）； β加热速率（Kmin-1）。 于是可以分别得到 非均相体系在等温与非等温条件下的两个常用动力学方程式 等温 （7） 非等温 （8） 动力学研究的目的就在于求解出能描述某反应的上述方程中的“动力学三因子” E、A 和fα 对于反应过程的DSC曲线如图所示。在DSC分析中，α值等于Ht/H0，这里Ht为物质A′在某时刻的反应热，相当于DSC曲线下的部分面积，H0为反应完成后物质A′的总放热量，相当于DSC曲线下的总面积。 二、 微分法 2．1 Achar、Brindley和Sharp法 对方程进行变换得方程 （9） 对该两边直接取对数有 （10） 由式（11）可以看出，方程两边成线性关系。 通过试探不同的反应机理函数、不同温度T时的分解百分数，进行线性回归分析，就可以试解出相应的反应活化能E、指前因子A和机理函数fα. 2．2 Kissinger法 Kissinger在动力学方程时，假设反应机理函数为，相应的动力学方程表示为 （11） 该方程描绘了一条相应的热分析曲线，对方程（12）两边微分，得 （12） 在热分析曲线的峰顶处，其一阶导数为零，即边界条件为 TTp （13） （14） 将上述边界条件代入（13）式有 （15） Kissinger研究后认为与β无关，其值近似等于1，因此，从方程（16）可变换为 （16） 对方程（15）两边取对数，得方程（18），也即Kissinger方程 ，i1，2，，4 （17） 方程（18）表明，与成线性关系，将二者作图可以得到一条直线，从直线斜率求Ek，从截距求Ak，其线性相关性一般在0.9以上。 2．3 两点法 Kissinger法是在有假定条件下得到的简化方程。如果我们不作任何假设，只是利用数学的方法进行，可以得到两点法。 由方程（2）、（5）知 （18） 方程（19）两边对T微分，得 （19） 当TTp时，反应速率达到最大，ααp，从边界条件有 我们得到第一个方程 （ 20） 方程（20）两边对T微分，得 （21） 这相当于对DSC曲线求二阶导，为的是求DSC曲线的拐点。在DSC曲线的拐点处，我们有边界条件 将该条件代入方程（22），从而得到第二个方程 0 （22） 联立方程（21）和（22），即得到只与反应温度T、机理函数fα有关的方程如下 式中 通过解方程就可求出非等温反应动力学参数E和A的值。 在该方法中，只需要知道升温速率β，拐点的温度Ti、分解百分数αi，峰顶的温度Tm、分解百分数αm，就可以试算不同的fα，以求解出对应于该fα时的活化能E值、指前因子A值。 三 积分法 对于积分法， 则由方程（8）求积分得 （23） 式中 对P（u）的不同处理，构成了一系列的积分法方程，其中最著名的方法和方程如下 3．1 Ozawa法 通过对方程（23）变换，得Ozawa公式 （24） 方程（24）中的E，可用以下两种方法求得。 方法1 由于不同βi下各热谱峰顶温度Tpi处各α值近似相等，因此可用“”成线性关系来确定E值。令 这样由式（24）得线性方程组 解此方程组求出a，从而得E值。 Ozawa法避开了反应机理函数的选择而直接求出E值，与其它方法相比，它避免了因反应机理函数的假设不同而可能带来的误差。因此往往被其它学者用来检验由他们假设反应机理函数的方法求出的活化能值，这是Ozawa法的一个突出优点。 3．2 Phadnis法 式中 （25） 该方程由Phadnis等人提出。对于合适的机理函数，与成线性关系，由此求出E值，但无法求出A值。 3．3 Coats-Redfern近似式 取方程（23）右端括号内前二项，得一级近似的第一种表达式Coats-Redfern近似式 （26） 式中 并设，则有 积分方程（4-3），整理，两边取对数，得 当时， （27） 当时， （28） 上述两个方程都称为Coats-Redfern方程。 由于对一般的反应温区和大部分的E值而言，，所以方程（4-4）和（4-5）右端第一项几乎都是常数，当时，对作图，而时，对作图，都能得到一条直线，其斜率为（对正确的n值而言）。 3． 4 Mac Callum-Tanner近似式 该法无需对pu作近似处理，可以证明，对于一定的E值，-log pu与1/T为线性关系，并可表达为 而且，E对也是线性关系，可表达为 于是有 虽然u对E不是线性关系，但是logu对logE是线性关系，即 于是有 借助于附录A中列出的logpuu表计算出相应的常数后，代入上式，得 式中E 活化能，kcal/mol T 温度，K 上述方程称Mac Callum-Tanner近似式。 4． 计算结果判据 提出的选择合理动力学参数及最可几机理函数的五条判据是 1 用普适积分方程和微分方程求得的动力学参数E和A值应在材料热分解反应动力学参数值的正常范围内，即活化能E值在80250kJmol-1之间，指前因子的对数（lgA/s-1）值在730之间； 2 用微分法和积分法计算结果的线性相关系数要大于0.98； 3 用微分法和积分法计算结果的标准偏差应小于0.3； 4 根据上述原则选择的机理函数f（α）应与研究对象的状态相符； 5 与两点法、Kissinger法、Ozawa法和其它微积分法求得的动力学参数值应尽量一致。 5 常用的动力学机理函数 函数号 函数名称 机理 积分形式Gα 微分形式fα 1 抛物线法则 一维扩散, 1D, D1减速形α-t曲线 α2 2 Valensi方程 二维扩散, 园柱形对称, 2D, D2, 减速形α-t曲线 3 Jander方程 二维扩散, 2D, 4 Jander方程 二维扩散，2D，n2 5 Jander方程 三维扩散，3D， 6 Jander方程 三维扩散,球形对称，3D，D3，减速形α-t曲线，n2 7 G.-B方程* 三维扩散,球形对称，3D，D4，减速形α-t曲线 8 反Jander方程 三维扩散，3D 函数号 函数名称 机理 积分形式Gα 微分形式fα 9 Z.-L.-T.方程** 三维扩散，3D 10 Avrami-Erofeev方程 随机成核和随后生长，A4，S形α-t曲线，，m4 11 Avrami-Erofeev方程 随机成核和随后生长，A3，S形α-t曲线，，m3 12 Avrami-Erofeev方程 随机成核和随后生长， 13 Avrami-Erofeev方程 随机成核和随后生长，A2，S形α-t曲线，，m2 14 Avrami-Erofeev方程 随机成核和随后生长， 15 Avrami-Erofeev方程 随机成核和随后生长， 16 Mample单行法则，一级 随机成核和随后生长，假设每个颗粒上只有一个核心，A1，F1，S形α-t曲线，，m1 17 Avrami-Erofeev方程 随机成核和随后生长， 函数号 函数名称 机理 积分形式Gα 微分形式fα 18 Avrami-Erofeev方程 随机成核和随后生长， 19 Avrami-Erofeev方程 随机成核和随后生长， 20 Avrami-Erofeev方程 随机成核和随后生长， 21 P.-T.方程*** 自催化反应, 枝状成核, Au, B1 S形α-t曲线 22 Mampel Power法则幂函数法则 23 Mampel Power法则幂函数法则 24 Mampel Power法则幂函数法则 25 Mampel Power法则幂函数法则 相边界反应一维,R1, n1 1 26 Mampel Power法则幂函数法则 27 Mampel Power法则幂函数法则 n2 续表 函数号 函数名称 机理 积分形式Gα 微分形式fα 28 反应级数 29 30 收缩球状（体积） 相边界反应，球形对称，R3，减速形α-t曲线，，n3（三维） 31 32 收缩园柱体（面积） 相边界反应，园柱形对称，R2，减速形α-t曲线，，n2（二维） 33 反应级数 n2 34 反应级数 n3 35 反应级数 n4 36 二级 化学反应，F2，减速形α-t曲线 37 反应级数 化学反应 38 2/3级 化学反应 39 指数法则 E1，n1，加速形α-t曲线 40 指数法则 n2 续表 函数号 函数名称 机理 积分形式Gα 微分形式fα 41 三级 化学反应，F3，减速形α-t曲线 42 S-B方程**** 固相分解反应SB（m，n） 43 反应级数 化学反应， ROn， 44 J-M-A方程***** 随机成核和随后生长，An，JMA（n） 45 幂函数法则 P1，加速型α-t曲线 *, Ginstling-Brounstein方程 **, Zhuralev-Lesokin-Tempelman方程 ***,Prout-Tompkins方程 ****,šestk-Berggren方程 *****, Johnson-Mehl-Avrami方程 注函数No.1和27称谓不同，形式相同。