柴油机颗粒氧化动力学反应分析方法比较
2024-03-16
来源:步旅网
第32卷第1期 西华大学学报(自然科学版) 2013年1月 Vo1.32,No.1 Journal of Xihua University・Natural Science Jan.2013 ・新能源汽车与低碳运输・ 柴油机颗粒氧化动力学反应分析方法比较 . 孟忠伟,杨冬,闫妍 (西华大学交通与汽车工程学院,汽车测控与安全四川省重点实验室,四川成都610039) 摘要:采用热重法(TGA)分析柴油机颗粒的氧化动力学性能,其动力学参数的处理方法较多,具有典型代表 的是积分法(Coats—Redfern)、微分法(ABSW)和Arrhenius法。文章采用上述3种方法,研究了改变柴油机颗粒初 始质量m 对活化能E及指前因子A 的影响。实验结果表明:当取反应级数n=1时,利用试样总质量和干碳黑质 量分别作为初始质量进行比较,当使用积分法(Coats—Redfern)处理数据时,活化能及指前因子都将增大;而使用 微分法(ABSW)时,活化能不会发生改变,指前因子会增大;使用Arrhenius方法时,活化能只有微小的改变,指前因 子有一定增大。 关键词:柴油机颗粒;热重分析法;活化能;指前因子 中图分类号:TK421 文献标志码:A 文章编号:1673—159X(2013)01—0051—05 doi:10.3969/j.issn.1673—159x.2013.01.009 The Comparison of Analysis Method of Diesel Particulate Oxidation Kinetics MENG Zhong—wei,YANG Dong,YAN Yan (Vehicle Measurement,Control and Safety Key Laboratory of Sichuan Province,School fo Transportation and Automotive Engineering,Xihua University,Chengdu 610039 China) Abstract:There are several thermo—gravimetric analysis methods used to analysis the oxidative combustion of diesel particulate for TGA data processing.The most typical ones are the integration(Coats—Redfern),the differentila(ABSW)and the Arrhenius method. Based on these three methods,the effects of changing the initial sample mass on the activation energy and the pre-exponential factor are investigated and compared.The results show that if the total soot mass instead of the dry soot mass is used in calculation,the value of activation energy and pre—exponential factor will be increased obviously when using the Coats—Redfern method.The value of activation energy will be unchanged and the value of the pre—exponential factor will be increased when using the ABSW method.The value of acti— ration energy will be slightly changed and the pre—exponential factor will be increased when using the Arrhenius method. Key words:diesel particulate;thermo—gravimetric analysis;activation energy;pre-exponential factor 柴油机排放颗粒物中,70%的粒径小于 物进行氧化燃烧再生可以有效实现碳烟的净化_3 J。 0.3 m,主要是由碳烟颗粒以及吸附在其表面的可 为了使碳烟颗粒能够在排气温度下更易燃烧,减小 溶性有机物(SOF)、碳氢液滴和硫酸液滴组成 。 其活化能是一种非常行之有效的方法。 碳烟颗粒容易被人体吸入肺泡且不易排除,其上面 采用热重法对柴油机颗粒进行升温氧化,并根 吸附的SOF物质能致癌 J,为此,排放法规对柴油 据实验数据利用热分析动力学可以计算出相应的 车的颗粒物排放进行了严格的控制。通过将颗粒 活化能和指前因子等动力学参数。柴油机颗粒的 收稿日期:2012-l1-25 资金项目:国家自然科学基金资助项目(51106130);四川省教育厅青年基金项目(1OZB078,12ZB138);四川省教育厅科技创新团队项目 (KYTD201003);西华大学重点科研基金项目(Zl120320,Z1120319);汽车测控与安全四川I省重点实验室开放基金项目 (szJJ2011—001);西华大学流体及动力机械省部共建教育部重点实验室学术成果培育资助项目(SBZDPY—l1—18);四川省重 点科技项目(2011JYZ014)。 作者简介:孟忠伟(1980一),男,副教授,博士,主要研究方向为柴油机颗粒物捕集器。E mail:mengzw@mail.xhu.edu.cn. 52 西华大学学报・自然科学版 热重氧化分析可以分为2个阶段,第1阶段在150 ~400℃温度区间,主要是SOF等成分的燃烧,第2 阶段在600~700℃温度区间内,主要是干碳烟的燃 烧 J。目前已开展的研究工作有:常仕英等利用 Arrhenius法分析了碳烟的不同氧化阶段和燃烧速 率 ;何百磊等分别利用Coats—Redfern、ABSW和 Arrhenius法这3种方法对碳烟氧化的活化能和指 前因子进行了计算,并探究了传质对热重分析实验 结果的影响 ;张广龙等利用Coats—Redfern法计 算了当使用不同类型的催化剂时碳烟的活化能和 指前因子 。然而,上述研究都将碳烟的总质量作 为初始质量进行计算,而没有选取干碳烟质量来研 究其在600~700℃温度区间内氧化动力学参数。 其实在这一温度区间,剩余的样品质量主要是干碳 烟的质量。因此,有必要深入比较初始质量的选取 对碳烟氧化动力学参数的影响。本文利用碳黑和 柴油的混合物模拟柴油机碳烟颗粒进行热重实验, 分别利用试样总质量和干碳黑质量计算碳黑的氧 化动力学参数,并进行比较,是运用典型的动力学 公式 进行分析,这为柴油机碳烟颗粒氧化动力学 分析提供一定的科学依据。 1 实验方法 1.1热重实验 柴油机颗粒中大部分物质主要是干碳烟颗粒, 大多数情况下SOF的质量占颗粒总质量的40%左 右 j。本实验用碳黑(德固赛Printex—U碳黑,其 比表面积S =100 m /g)与柴油的混合物来模拟柴 油机颗粒物,配制了柴油质量分数为40%(碳黑质 量分数为60%)的实验试样。利用TGAQ600热重 分析仪,在空气氛围中,以70 K/min的升温速率将 试样从室温(25 oC左右)升温到800℃,TGA记录 下不同反应时刻的样品质量及样品温度。对实验 得出的数据进行热重动力学分析,可以得到O/一T (转化率一温度)曲线、失重率等曲线。 根据O/一T曲线可以发现,混合试样中柴油最 先反应燃烧,当温度到达520℃以后碳黑才开始氧 化燃烧。本文分别选取室温时(试样未开始反应 时)混合试样的总质量 。 以及升温到520℃时碳 黑剩余质量m。 。作为初始质量m。,比较分别采用 上述3种氧化动力学方法对分析碳黑氧化活化能和 指前因子的影响。 1.2动力学参数 首先定义任意时刻t或者温度 时试样的转化 率(也称失重率) OL =—: I二— (1)J m0一m 式中:m 为试样起始质量;m为任意时刻时的质量; m 为终止质量。 相对反应速率表达式: r= =A …p(_ ).(1 S ・P : (2) 本征反应速率常数: k=A'exp(一 ) (3) 式中:碳黑比表面积S :S.,o/(m。一m ),(m ・ g ),其中S 是初始表面积,m ;E为表观活化能, J・tool~;A 为总指前因子,s・m~;P。 为氧气 分压l6 J。 2 实验结果与分析 试样以70 K/min的升温速率进行的TGA实验 得出的质量变化曲线和相对反应速率曲线如图1和 图2所示。根据曲线可以很明显地看出实验试样的 2个反应阶段,当升温到300 oC左右柴油几乎燃烧 完全。相对于柴油机颗粒物可以认为是SOF的氧 化燃烧,相对反应速率r在200℃左右达到最大。 在经过一段时间的升温后,在350~500℃试样的相 对反应速率趋近于0,试样在这段时间内并没有明 显的质量变化。当温度达到520℃左右时,碳黑开 始反应燃烧,在600~700℃区间内开始快速反应, 试样质量迅速减少,相对反应速率r在升温到 650 oC左右时达到最大。 料 K 图1转化率、失重率与温度曲线 通过图2可以获得柴油和碳黑燃烧的3个特征 温度:起燃温度、峰值温度和燃尽温度。起燃温度 通常是指试样反应掉的质量达到试样总质量10%时 第1期 孟忠伟,等:柴油机颗粒氧化动力学反应分析方法比较 53 ∞ 如 加 m 0 始质量m 对拟合直线的斜率和截距依然有比较大 的影响,当初始质量m 减小,拟合出的直线的斜率 和截距值会明显增大。利用碳黑的质量m。 。来拟 喜 合曲线时,其斜率和截距都明显大于采用试样总质 量m。一 得出的拟合曲线,说明利用纯碳黑质量 m。 。得出的活化能和指前因子将大于采用试样总 200 400 600 800 l000 m( 图2相对反应速率曲线 的反应温度,用 。表示;峰值温度是反应速率达到 最大时的反应温度,用 来表示;燃尽温度是指颗 粒物燃烧结束的温度,通常是反应掉的质量达到总 质量90%时的反应温度,用 表示。表1得出了 柴油和碳黑在反应中的这3个温度值。从柴油的特 征温度能看出,碳黑的快速氧化燃烧的温度区间是 110—240℃。从碳黑的特征温度也能看出,碳黑的 快速氧化燃烧的温度区间是610—710℃。 表1柴油、碳黑的特征温度值 ℃ 2.1 积分法(Coats—Redfern)的数据处理结果 将反应温度按照一定的升温速率 进行升温, 其表达式为: = dT (4) 联合公式(2)并进行移项积分可得: A'SPoⅧ(1 a, o:2f r。e p(一E/RT)dr(5)一Ot)“ .“ 公式两边分别取对数,取反应级数n=1得: g 【 】一 E 【o … lg(一In(1一 )/ )与1/T呈线性关系,然后对 曲线进行线性拟合,通过拟合直线的斜率和截距就 可以求得活化能E和指前因子A 。 分别用总质量m。 和碳黑的质量m 。作为 起始质量m。对实验数据进行曲线拟合,结果如图3 所示。根据公式(6)可以看出,公式两边经过2次 取对数后,公式(6)不能进行进一步的简化,初始质 量m。不能被通过数学方法舍去或者简化,因此,初 质量m。 得出的数据。 图3积分法的拟合曲线 2.2微分法(ABSW法)的数据处理结果 将公式(2)分离变量取对数有: ・n【 一 ㈩ 取 =1时有: n[ 一 (8) 根据高数知识,在 内取任一微元得到:d1 d = (9) — 带入公式(8)有: 一 (10) 联合公式(1)化简可得: ・ 一面E (11) 指前因子的计算公式: = (12) 式中,口。是所拟合出曲线的截距值。 n[ ]与・/ 呈线性关系, 拟合出相应直线,根据直线的斜率求活化能E,根据 直线截距求指前因子A 。 分别用总质量m 和碳黑的质量m 作为 54 西华大学学报・自然科学版 起始质量m 对实验数据进行曲线拟合,从公式 (11)可以看出,在利用In[ 亏 】对 ln(卜mmz )=In EA S Po2 ) 一面E(18) 利用lnk一1/T图可以拟合出相应的直线,由其 1/T作图时,公式经过变化化简后完全消去初始质 量m ,表明用该种方法拟合出的曲线,其斜率和截 距与反应试样的初始质量m 并没有直接的关系。 因此,公式不管是用总质量还是利用碳黑的质量作 为反应的初始质量m 进行曲线的拟合来计算活化 能和指前因子,碳黑活化能的计算结果都不会产生 变化。根据指前因子的计算公式(12),改变初始质 量m 将略微改变指前因子的数值。 2.3 Arrhenius方法的数据处理结果 对公式(2)和(3)整理并取对数有: n— 一丽/7,(13S a oP 0,(一 )“1 R?’ ),根据高数知识,在d ̄/dt内取任一微元,有: 垫: 二 ,14) dt t,一t 、 带人公式(9)有: ln =ln[ : _==_ I__)雨】= l 一 RT f 15) 联立公式(1)、(11)化简可得: n = n[ 二¨_ l_ 瓣]: lnA 一 /5 (16) 根据TGA实验中反应物的质量变化可知,在最 后反应阶段剩余质量很少,即m 趋近于0,因此为 了简化公式这里可以暂时忽略不计,从而公式变为: n = n[L S :P— o,(2—1)m0J-t= t] = n4 一 c 7 斜率和截距分另U可以求得活化能E和指前因子 A 值。 分别用总质量m。 和碳黑的质量m。 。作为 起始质量m 进行曲线拟合,如图4所示。分析公式 (18)可以得出,利用Arrhenius法作图时,反应试样 的初始质量m。对活化能E的计算无影响,对指前 因子A 有一定的影响,增大初始质量m。,指前因子 的值会增大。由于在进行公式的推导简化时忽略 了反应结束时的质量m ,实际上在反应结束后,还 会有一些很少的试样中的杂质等没有燃烧完全,对 结果会造成一定的影响。 1/丁TK。。1 图5 Arrhenius方法的拟合曲线 3 活化能的计算结果及分析 在进行TGA实验时,扩散的影响对碳黑本征活 化能的影响非常大¨ 。但是通过研究发现,可以选 择合适的坩埚尺寸和试样初始质量等因素来减小 甚至避免扩散的影响 卜 J。可以使用高度较低、直 径相对较大的坩埚,同时试样的初始质量应该偏 小¨ ,一般在2.5~5 mg范围内比较适宜。本文用 的TGAQ600的坩埚尺寸为西6.5 mm×4 mm,H/D =0.62,试样的初始质量为2.5 mg左右。因此,本 文实验可忽略扩散对活化能计算的影响,认为采用 碳黑质量m 加计算得到的活化能为本征活化能; 而采用总质量m。 计算得到的活化能为考虑柴油 和碳黑氧化相互影响的一种综合活化能。 表2列出了利用上述3种方法所计算出的活化 能和指前因子的值。利用ABSW和Arrhenius方法 无论是利用试样总质量m。 ,还是式样中碳黑质量 //Z。 。作为计算活化能的初始质量m。,所计算出的 活化能值几乎没变化,指前因子的值会略有增大; 而利用Coats—Redfern方法不仅活化能会增大,指 第1期 孟忠伟,等:柴油机颗粒氧化动力学反应分析方法比较 55 前因子也会增大。这是因为根据上文对各个公式 的变化推导结果可以得出:利用ABSW和A ̄henius 方法计算出的活化能是碳黑的本征活化能,不管试 碳烟的本征活化能。 参监测,2006,22(3):23—26. 考文献 [1]邹建国,钟秦.柴油机排放颗粒物组成分析[J].中国环境 样初始质量m 如何变化,活化能的计算结果都不 会发生改变;而使用Coats—Redfern法将总质量 [2]郝吉明,段雷.燃烧源可吸人颗粒物的物理化学特征[M]. 北京:科学出版社,2008:19—25. m。 作为初始质量m。时,所计算出的活化能是碳 黑与柴油混合物的综合活化能。由于柴油的影响, 试样的表观活化能会低于碳黑的本征活化能。3种 [3]王逢瑚,郭秀荣,马岩,等.柴油车尾气微粒捕集器技术研究 现状及发展趋势[J].小型内燃机与摩托车,2010,39(1):93—96. [4]Stanmore B R,Brilhac J F,Gilot P.The Oxidation of Soot:a 方法所计算出的碳黑本征活化能数值比较接近,且 都在合理误差范围内,从而验证上述分析的合理性。 表2计算出的动力学参数 不同初始质量动力学参数Coats-Redfem ABSW Arrhenius 4 结论 1)根据Ot—T曲线分析可知,柴油+碳黑模拟 的柴油机颗粒物的氧化主要分为2个阶段:第1阶 段为l10—240℃,主要为柴油(模拟SOF组分)的 氧化阶段;第2阶段为610~710℃,主要为碳黑(模 拟干碳烟组分)的氧化燃烧阶段。 2)推导并建立了3种典型的计算活化能的热 分析动力学模型,通过计算结果对比验证了它们的 实用性和正确性。 3)在取反应级数 =1时,当改变计算公式初 始质量m。后,使用积分法(Coats—Redfem)所计算 的活化能值和指前因子则都会增加,使用微分法 (ABSW)和A ̄henius方法所计算出的活化能值几 乎不会改变,但是计算出的指前因子则略有增加。 4)积分法(Coats—Redfern法)不仅能计算本征 活化能,还能根据需求计算试样的综合活化能,而 微分法(ABSW)和Arrhenius方法只能计算试样中 Review of Experiments,Mechanisms and Models[J].Carbon,2001,39 (15):2248—2265. [5]常仕英,吴庆伟,杨冬霞,等.柴油车碳烟的燃烧特性及动力 学研究[J].内燃机学报,2009,27(3):256—258. [6]何百磊.柴油机颗粒物捕集器再生用CaMnO3系列催化剂 的研究[D].北京:清华大学,2006. [7]张广龙,荆军航,宁智,等.碳烟微粒催化氧化过程动力学分 析[J].内燃机学报,1995,13(1):47—51. [8]胡荣祖,史启祯.热分析动力学[M].北京:科学出版社, 2001:10—25. [9]张延峰,宋崇林,成存玉,等.车用柴油机排气颗粒物中有机 组分和无机组分的分析[J].燃烧科学与技术,2004,10(3):198 —201. [10]Neeft J P A,Makkee M,Moulijn J A,et a1.The Effects of Heat and Mass Transfer in Thermogravimetrieal Analysis:A Case Study Towards the Catalytic Oxidation of Soot[J].Thermochimica Acta,1996, 287(2):261—278. [11]Zouaoui N,Brilhac J F,Mechati F,et a1.Study of Experi— mental and Theoretical Procedures when Using Thermogravimetric Analy— sis to Determine Kinetic Parametersof Carbon Black Oxidation[J].Jour- hal of Thermal Analysis and Calorimetry,2010,lO2(3):837—849. [12]Sharma H N,Pahalagedara L,Joshi A,et a1.Expeirmental Study of Carbon Black and Diesel En ̄ne Soot Oxidation Kinetics Using hTermogravimetric Analysis[J].Energy&Fuels,2012,26(9):5613 —5625. [13]Gilot P,Bonnefoy F,Marcuceilli F,et a1.Determination of Kinetic Data for Soot Oxidation:Modeling of Competition between Oxygen Difusion and Reaction during Thennogravimetirc Analysis[J].Combus— tion and Flame,1993,95(1—2):87—100. (编校:夏书林)