热分析动力学基础知识资料下载.pdf
《热分析动力学基础知识资料下载.pdf》由会员分享,可在线阅读,更多相关《热分析动力学基础知识资料下载.pdf(16页珍藏版)》请在冰点文库上搜索。
由于f()和G()分别为机理函数的微分形式和积分形式,它们之间的关系为:
d/)(d1)
(1)(GGf(4)k与反应温度T(绝对温度)之间的关系可用著名的Arrhenius方程表示:
)/exp(RTEAk(5)式中:
A表观指前因子;
E表观活化能;
R通用气体常数。
方程
(2)(5)是在等温条件下出来的,将这些方程应用于非等温条件时,有如下关系式:
tTT0(6)即:
/tddT式中:
T0DSC曲线偏离基线的始点温度(K);
加热速率(Kmin-1)。
于是可以分别得到:
非均相体系在等温与非等温条件下的两个常用动力学方程式:
)E/RT)f(Atddexp(/(等温)(7))/exp()(ddRTEfAT(非等温)(8)动力学研究的目的就在于求解出能描述某反应的上述方程中的“动力学三因子”动力学研究的目的就在于求解出能描述某反应的上述方程中的“动力学三因子”E、A和f(E、A和f()对于反应过程的DSC曲线如图所示。
在DSC分析中,值等于Ht/H0,这里Ht为物质A在某时刻的反应热,相当于DSC曲线下的部分面积,H0为反应完成后物质A的总放热量,相当于DSC曲线下的总面积。
二、微分法21Achar、Brindley和Sharp法:
对方程)/exp()(ddRTEfAT进行变换得方程:
)/exp(dd)(RTEATf(9)对该两边直接取对数有:
RTEATflndd)(ln(10)由式(11)可以看出,方程两边成线性关系。
通过试探不同的反应机理函数、不同温度T时的分解百分数,进行线性回归分析,就可以试解出相应的反应活化能E、指前因子A和机理函数f().22Kissinger法22Kissinger法Kissinger在动力学方程时,假设反应机理函数为nf)1()(,相应的动力学方程表示为:
nRTEAet)1(dd/(11)该方程描绘了一条相应的热分析曲线,对方程(12)两边微分,得tAeteAttnRTERTEnd)1(ddd)1(dddd/tnAetTRTEeAnRTERTEndd)1(dd)1()()1(1/2/tnAetTRTEtnRTEdd)1(dddd1/2RTEneAnRTtTEt/12)1(dddd(12)在热分析曲线的峰顶处,其一阶导数为零,即边界条件为:
T=Tp(13)0ddddtt(14)将上述边界条件代入(13)式有:
RTEnpeAnRTtTE/1p2)1(dd(15)Kissinger研究后认为:
与无关,其值近似等于1,因此,从方程(16)可变换为:
1p)1(nnp/2pRTEAeRTE(16)对方程(15)两边取对数,得方程(18),也即Kissinger方程:
pikkk2pi1lnlnTREERATi,i=1,2,4(17)方程(18)表明,2pilnTi与pi1T成线性关系,将二者作图可以得到一条直线,从直线斜率求Ek,从截距求Ak,其线性相关性一般在0.9以上。
23两点法Kissinger法是在有假定条件下得到的简化方程。
如果我们不作任何假设,只是利用数学的方法进行,可以得到两点法。
由方程
(2)、(5)知)(ddfAetRTE(18)方程(19)两边对T微分,得2/)()(ddRTEefAeAfTdtdRTERTE(19)当T=Tp时,反应速率达到最大,=p,从边界条件有:
0,ddddppTTTt我们得到第一个方程:
0)(2p/ppRTEefARTE(20)方程(20)两边对T微分,得RTERTERTEefRTAEefAeAfTt/22222/22)(3)()(dddd4222222)()(TRERTEeffARTE(21)这相当于对DSC曲线求二阶导,为的是求DSC曲线的拐点。
在DSC曲线的拐点处,我们有边界条件:
0,ddddpi22TTTt将该条件代入方程(22),从而得到第二个方程iiRTEiRTEiefRTAEefA/22222)(3)(+4222222)()(iiRTEiTRERTEeffAi=0(22)联立方程(21)和(22),即得到只与反应温度T、机理函数f()有关的方程如下:
021)()f(YE,422iiEUTRERTeDCBmmRTfTAm2EREe式中:
22RmmiTffB222R3immiTTffC422R1mmiiTfffDmimiTTTTUR通过解方程就可求出非等温反应动力学参数E和A的值。
在该方法中,只需要知道升温速率,拐点的温度Ti、分解百分数i,峰顶的温度Tm、分解百分数m,就可以试算不同的f(),以求解出对应于该f()时的活化能E值、指前因子A值。
三积分法对于积分法,tkG)(则由方程(8)求积分得TTTTRTEATRTEAfG00d)/exp(d)/exp()(d)(0)()(d2uueRAEupRAEuueRAEuuu(23)式中:
RTEuuuuup);
()exp()(对P(u)的不同处理,构成了一系列的积分法方程,其中最著名的方法和方程如下:
31Ozawa法31Ozawa法通过对方程(23)变换,得Ozawa公式:
RTERGAE4567.0315.2)(loglog(24)方程(24)中的E,可用以下两种方法求得。
方法1:
由于不同i下各热谱峰顶温度Tpi处各值近似相等,因此可用“T1log”成线性关系来确定E值。
令:
REaLiTyZii4567.0),2,1(/1logpii315.2)(logRGAEb这样由式(24)得线性方程组),2,1(LibayZii解此方程组求出a,从而得E值。
Ozawa法避开了反应机理函数的选择而直接求出E值,与其它方法相比,它避免了因反应机理函数的假设不同而可能带来的误差。
因此往往被其它学者用来检验由他们假设反应机理函数的方法求出的活化能值,这是Ozawa法的一个突出优点。
3232Phadnis法Phadnis法RTEuRTEeERTueREupRETeFK/22T0/)(d式中2)(FKueupuTERTfGdd)()(2(25)该方程由Phadnis等人提出。
对于合适的机理函数,)()(fG与TTdd2成线性关系,由此求出E值,但无法求出A值。
33Coats-Redfern近似式取方程(23)右端括号内前二项,得一级近似的第一种表达式Coats-Redfern近似式:
RTEuuRTEeERTERTuueREuueREupRETe/232T0/21221)(d(26)式中:
uueuueuPuuCR2112)(23并设nf)1()(,则有RTEneERTERTA/2021)1(d积分方程(4-3),整理,两边取对数,得当时,1nRTEERTEARnTn21ln)1()1(1ln21(27)当时,1nRTEERTEART21ln)1ln(ln2(28)上述两个方程都称为Coats-Redfern方程。
由于对一般的反应温区和大部分的E值而言,121,1ERTRTE,所以方程(4-4)和(4-5)右端第一项几乎都是常数,当1n时,)1()1(1ln21nTn对T1作图,而时,1n2)1ln(lnT对T1作图,都能得到一条直线,其斜率为RE(对正确的n值而言)。
334MacCallum-Tanner近似式4MacCallum-Tanner近似式该法无需对p(u)作近似处理,可以证明,对于一定的E值,-logp(u)与1/T为线性关系,并可表达为:
Tauup)(log而且,E对也是线性关系,可表达为:
abEya于是有TbEyuup)(log虽然u对E不是线性关系,但是logu对logE是线性关系,即:
EcAulogloglog于是有TbEyAEupc)(log借助于附录A中列出的logp(u)u表计算出相应的常数后,代入上式,得:
TEEupMT001.0217.0449.04828.0)(log4357.0TEEup001.0217.0449.04828.0MT4357.010)(式中:
E活化能,kcal/molT温度,K上述方程称MacCallum-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曲线21212Valensi方程二维扩散,园柱形对称,2D,D2,减速形-t曲线)1ln()1
(1)1ln(3Jander方程二维扩散,2D,21n2121)1(1212121)1
(1)1(44Jander方程二维扩散,2D,n=2221)1(112121)1
(1)1(5Jander方程三维扩散,3D,21n2131)1(1213132)1
(1)1(66Jander方程三维扩散,球形对称,3D,D3,减速形-t曲线,n=2231)1(113132)1
(1)1(237G.-B方程(*)三维扩散,球形对称,3D,D4,减速形-t曲线32)1(3211311)1(238反Jander方程三维扩散,3D2311)1(131321)1()1(23函数号函数名称机理积分形式G()微分形式f()9Z.-L.-T.方程(*)三维扩散,3D2311)1(131341)1()1(2310Avrami-Erofeev方程随机成核和随后生长,A4,S形-t曲线,41n,m=441)1ln(43)1ln()1(411Avrami-Erofeev方程随机成核和随后生长,A3,S形-t曲线,31n,m=331)1ln(32)1ln()1(312Avrami-Erofeev方程随机成核和随后生长,52n52)1ln(53)1ln()1(2513Avrami-Erofeev方程随机成核和随后生长,A2,21)1ln(21)1ln()1(2S形-t曲线,21n,m=214Avrami-Erofeev方程随机成核和随后生长,32n32)1ln(31)1ln()1(2315Avrami-Erofeev方程随机成核和随后生长,43n43)1ln(41)1ln()1(3416Mample单行法则,一级随机成核和随后生长,假设每个颗粒上只有一个核心,A1,F1,S形-t曲线,1n,m=1)1ln(117Avrami-Erofeev方程随机成核和随后生长,23n23)1ln(21)1ln()1(32函数号函数名称机理积分形式G()微分形式f()18Avrami-Erofeev方程随机成核和随后生长,2n2)1ln
(1)1ln()1(2119Avrami-Erofeev方程随机成核和随后生长,3n3)1ln
(2)1ln()1(3120Avrami-Erofeev方程随机成核和随后生长,4n4)1ln(3)1ln()1(4121P.-T.方程(*)自催化反应,枝状成核,Au,B1(S形-t曲线)1ln)1(22MampelPower法则(幂函数法则)41n4143423MampelPower法则(幂函数法则)31n3132324MampelPower法则(幂函数法则)21n2121225MampelPower法则(幂函数法则)相边界反应(一维),R1,n=111)1(1126MampelPower法则(幂函数法则)23n23213227MampelPower法则(幂函数法则)n=22121续表函数号函数名称机理积分形式G()微分形式f()28反应级数41n41)1(143)1(42930收缩球状(体积)相边界反应,球形对称,R3,减速形-t曲线,31n,n=3(三维)31)1(131)1(1332)1(332)1(3132收缩园柱体(面积)相边界反应,园柱形对称,R2,减速形-t曲线,21n,n=2(二维)21)1(121)1(1221)1(221)1(33反应级数n=22)1(11)1(2134反应级数n=33)1(12)1(3135反应级数n=44)1(13)1(4136二级化学反应,F2,减速形-t曲线1)1
(2)1(37反应级数化学反应1)1(12)1(382/3级化学反应21)1(23)1(239指数法则E1,n=1,加速形-t曲线ln40指数法则n=22ln21续表函数号函数名称机理积分形式G()微分形式f()41三级化学反应,F3,减速形-t曲线2)1(3)1(2142S-B方程(*)固相分解反应SB(m,n)nm)1(43反应级数化学反应,RO(n),nR11nn1)1(11n)1(44J-M-A方程(*)随机成核和随后生长,An,JMA(n)n/1)1ln(nn11)1ln()1(45幂函数法则P1,加速型-t曲线n/1nnn/)1()(*,Ginstling-Brounstein方程*,Zhuralev-Lesokin-Tempelman方程*,Prout-Tompkins方程*,estk-Berggren方程*,Johnson-Mehl-Avrami方程注:
函数No.1和27称谓不同,形式相同。