固相反应动力学

第三节固相反应动力学一.一般动力学关系整个过程的速度将由其中速度最慢的一环控制。现以金属氧化反应M+1/2O2→MO为例（图5）说明之。图5由化学动力学和菲克第一定律，其反应速度V和扩散速度V分别为：V==KcV==D=DdtdQPdtdQDdxdccc0当过程达到平衡时，Vp=VD或K0=DC=c0V=Kc=cc0DK110011DcKc讨论：1、当扩散速度远大于化学反应速度时，即K＜＜D/δ，则V=kc=Vp最大（式中c0=c），说明化学反应速度控制此过程，称为化学动力学范围。2、当扩散速度远小于化学反应速度时，即K＞＞D/δ，即c=0，V=D说明扩散速度控制此过程，称为扩散动力学范围。。最大DVcDcc003、当扩散速度远和化学反应速度可相比拟时，则过程速度由上式确定，称为过渡范围，即V=最大最大DPVVDcKc1111100二化学动力学范围1.此过程的特点是：反应物通过产物层的扩散速度远大于接触面上的化学反应速度。过程总的速度由化学反应速度所控制。2.均相二元系统化学反应速度的表达式对于均相二元系统，化学反应速度的一般表达式是V=KncAmcBn对于反应过程中只有一个反应物的浓度是可变的，则上式可简化为：V=Kncn（n≠1，n＞0）令经过任意时间t，有X部分反应物消耗于反应，而剩下的反应物量为（c－X）。上式可写成：—积分并考虑到初始条件：t=0，x=0得：nnXcKdtXcd)()(dtKXcdXxtnn00)(tKcXcnXccnnnnn1111))(1()(或这里，n是反应级数。故给出除一级以外的任意级数的这类反应的动力学积分式。tKcXcnnnn111)(111讨论：零级反应：n=0，X=K0t二级反应：n=2，或tKcXc211tKXccX2）（一级反应：n=1,可求得：tkcxcdtkdxxcxcdtkxcxcdtx10101ln)()(3.非均相二元系统化学反应速度表达式对于均相二元系统，计算过程中未考虑接触面积的影响，实际上，固相反应层非均相反应，接触面积在反应过程是影响固相反应速度的一个重要参数。接触面积用F表示。非均相二元系统反应的一般速度方程：nBmAncFcKdtdxv当反应过程中只有一个反应物浓度可变时：下面，引入转化率G的概念转化率G=转化的反应物量（或消耗掉的反应物量）/原始反应物量消耗掉的反应物x即等于G反应物浓度C=1-G那么，二元系统非均相化学反应一般速度方程即可写成：dG/dt=knF（1-G）nnnFcKdtdxv4.接触面积F的计算转化程度：G=R0—x=R0（1—G）1/3或x=R0[1－（1－G）1/3]相应于每个颗粒的反应表面积F’与转化程度G的关系：F’=A’（1—G）2/3对于单位质量的物系，F=A（1—G）2/3303030(RxRR）其中对于球形颗粒A=，对于立方形A＝。03R018R5.化学反应控制范围的动力学方程将上式反应级数n代入数值，就可求出不同级数反应的微积分形式：1）零级反应n=0dG/dt=KA（1-G）2/3（1-G）0=K0（1-G）2/3nnnnGGAkGFkdtdG)1()1()1(3/2积分并考虑到初始条件：t=0，G=0，得：F0(G)=1-（1-G）2/3=K0t球形或立方体颗粒：F0(G)=1-（1-G）1/3=K0t园柱形颗粒：F0(G)=1-（1-G）1/2=K01t平板状颗粒：F0(G)=G=K02ttGdtKGdG0003/2)1(6.实验验证如何验证上述动力学是正确的？如果我们能够使扩散阻力很小，这时扩散很快，反应为化学反应所控制。实验上常采取降低反应物颗粒度，再加入助熔剂来使反应处于化学动力学范围。如NaCO3：SiO2=1：1，进行固相反应，其动力学是化学反应控制的一级反应。三扩散动力学范围1.过程特点扩散速度很慢，起控制作用，为整个固相反应中速度最慢的一步。在多数情况下,扩散速度往往起控制作用。图6平板扩散模型2.动力学方程（1）抛物线型速度方程——平板模型此方程可从平板扩散模型导出。如图6所示。若化学反应速度远大于扩散速度，则过程由扩散控制。经dt时间，通过AB层迁移的A物质量为dm，平板间接触面积为S；浓度梯度为dc/dx，则按菲克定律有：dxdcDSdtdm由于A物质迁移量dm是比例于sdx，故积分得：F4（G）=x2=2K4’Dt=K4t上式即为抛物线速度方程的积分式。说明反应产物层厚度与时间的平方根成比例。这是一个重要的基本关系，可以描述各种物理或化学的控制过程并有一定的精确度。xDKdtdx‘4xDSdtdm1实验验证：图7示出的金属镍氧化时的增重曲线就是一个例证。局限性：但是，由于采用的是平板模型，忽略了反应物间接触面积随时间变化的因素，使方程的准确度和适用性都受到局限。图7金属镍的氧化增重曲线（2）杨德方程——球体模型在材料生产中通常采用粉状物料作为原料，这时，在反应过程中，颗粒间接触界面积是不断变化的。为此，扬德在抛物线速度方程基础上采用了“球体模型”导出了扩散控制的动力学关系。A.扬德假设：a、反应物是半径为R的等径球粒；b、反应物A是扩散相，即A成分总是包围着B的颗粒，而且A、B同产物C是完全结束的，反应自球表面向中心进行；c）A在产物层中的浓度是线性的，而扩散层截面积一定，图8杨德模型B动力学方程的推导现令以B物质为基准的转化程度为G，则G=代入抛物线速度方程式得x2=R02[1-（1-G）1/3]2=K4tF5（G）=[1-（1-G）1/3]2=K4/R。2=K4t303030301)1(1)(RxRxRRVV23/1202])1(1[GRx微分得：其中k5（c是常数，q是活化能，R是气体常数），也称杨德速度常数。3/13/25)1(1)1(GGKdtdGC）实验验证对于反应BaCO3+SiO2→BaSiO3+CO2的实测结果示于图9。由图可见，随着反应温度的升高，反应强度规律也提高了，但都很好地符合杨德方程。0.0160.0120.0080.00401020304050607080[1-(1-G)1/3]2时间（分）890℃870℃830℃图9不同温度下BaCO3与SiO2的反应情况D、杨德方程的适用范围——反应初期、G较小时对碳酸盐和氧化物间的一系列反应进行实验研究，发现在反应初期都基本符合杨德方程式，而后偏差就愈来愈大。为什么会这样呢？原因是杨德方程虽然采用了球体模型，在计算产物厚度时考虑了接触界面的变化，即利用反应前后球体之体积差算出产物层厚度x。但在将x值代入抛物线方程时实质上又保留了扩散面积恒定的假设，这是导致其局限性的主要原因。（3）金斯特林格方程——三维球体模型金斯特林格采用了杨德的球状模形，但放弃了扩散截面不变的假设从而导出了更有普遍性的新动力学关系。A金斯特林格假设：a）假设反应A是扩散相，B是平均半径为R。的球形颗粒，反应沿B整个球表面同时进行，首先，A和B形成产物AB，厚度为x，x随反应进行而增厚b）A扩散到A-AB界面的阻力远小于通过AB层的扩散阻力，则A-AB界面上A的浓度不变为C。，因扩散控制则A在B-AB界面上的浓度为0图10金斯特林格模型B方程推导由于粒子是球形的，产物两侧界面A的浓度不变，故随产物层增厚，A在层内的浓度分布是r和时间t函数，即过程是一个不稳定扩散问题，可以用球面坐标情况下的菲克扩散方程描述：根据初始和边界条件r=Rt＞0，C(R0,t)=C0r=R0-xt＞0，C(R0-x,t)=0t=0，x=0]2[),(22）（rCrrCDttrxRrrCDdtdx0)(式中是比例常数，其中和分别是产物AB的比重和分子量，n是反应的化学计量常数，即和一个B分子化合所需的A的分子数，D是A在AB中的扩散系数。求解得：许多试验研究表明，金斯特林格方程具有更好的普遍性。n3/13/13/13/16)1(1)1()1(1)1(GGKGGKdtdGrF6(G)=1—2/3G—（1-G）2/3=K6t杨德方程与金斯特林格方程之比较,得：按上式令对G作图可得图11曲线。可见，当G值较小即转化程度较低时，说明两方程是基本一致的，反之，随G值增加，两式偏差越来越大。可见，杨德方程只是在转化程度较小时适用。当G值较大时，Kj将随G的增大而增大，而金斯特林格方程则在一定程度上克服了杨德方程的局限。3/13/23/1)1()1()1()()(GGKGKdtdGdtdGirirjrdtdGdtdG)/()图11反应产物层增厚速率与x/R0关系图12金斯特林格方程与杨德方程之比较表1部分重要的固相反应动力学方程控制范围反应类别动力学方翟的积分式A值对应于图13的曲线界面零级反应(对于球形颗粒)F0(G)=1一(1--G)1/3=K0t=0.2063(t/t0．5)0．20637化学零级反应(对于圆柱形颗粒)F1(G)=1一(1--G)1/2=K1t=0.2929(t/t0．5)0．29296反应零级反应(对于乎板试样)F2(G)=G)=K2t=0.5000(t/t0．5)0．50005控制一级反应(球形颗粒)F’3(G)=ln(1--G)=-K3t=0.6931(t/t0．5)0．69319扩散抛物线速度方程(干板试样)F4(G)=G2=K4’t=0.2500(t/t0．5)0．25001控制对圆柱形试样F7(G)=(1--G)ln（1-G）+G=K7t=0.1534(t/t0．5)0：15342杨德方程(球形试样)F5(G)=[1一(1--G)1/3]2=K5t=0.0426(t/t0．5)0．04263范围金斯特林格方程(球形试样)F6(G)=1一2/3G-(1--G)2/3=K6t=0.0367(t/t0．5)0．03674