基于A(yíng)rrhenius定理的化學(xué)動(dòng)力學(xué)數值計算法

第22卷第3期中國科技核心期刊2011年5月CHEMICAL RESEARCHhxy@henu. edu. en基于 Arrhenius定理的化學(xué)動(dòng)力學(xué)數值計算法鐘巍(西北核技術(shù)研充所陜西西安710024)摘要:基于 Arrhenius定理建立了一種新的化學(xué)動(dòng)力學(xué)數值模擬方法.將其與化學(xué)動(dòng)力學(xué)過(guò)程契合,能較好地體現化學(xué)動(dòng)力學(xué)過(guò)程.將該方法對簡(jiǎn)單一級反應、平行反應和復雜的綜合反應進(jìn)行模擬計算,模擬結果與準確解相對誤差小于0.5%關(guān)鍵詞: Arrhenius定理;化學(xué)動(dòng)力學(xué);數值算法中圖分類(lèi)號:0643.12文獻標志碼:A文章編號:1008-1011(2011103-0056-05Arrhenius theory-based numerical algorithm forchemical kineticsZHONG WeiNorthwest Institute of Nuclear Technology, Xian 710024, Shanri, China)Abstract: A numerical simulation method, based on Arrhenius theory, was proposed in relationto calculation of chemical kinetic equations. The simulation method has been found to be con-sistent of chemical kinetics and can well reflect the process of chemical kinetics. When it wasadopted to solve the first-order reaction, parallel reaction and complicated comprehensive reac-tion the relative errors between the numerical solutions and exact solutions were less thanKeywords: Arrhenius theory; chemical kinetics; numerical algorithm化學(xué)動(dòng)力學(xué)是物理化學(xué)的一個(gè)重要組成部分,是一門(mén)研究化學(xué)反應速率和反應機理的科學(xué).研究化學(xué)反應機理通常需要先求解相關(guān)的動(dòng)力學(xué)微分方程組.由于化學(xué)動(dòng)力學(xué)學(xué)科的發(fā)展以及數學(xué)學(xué)科的發(fā)展,求解微分方程組的解析方法已經(jīng)有很多,常用的有本征值和本征向量法,分離變量法和消元法2等,但大部分化學(xué)動(dòng)力學(xué)方程組都很復雜,很難求得解析解,于是數值求解方法常被用于化學(xué)動(dòng)力學(xué)方面的研究.常見(jiàn)的數值求解方法有 Euler方法,四階 Runge-Kutta方法,數值積分法,數值微分法和遺傳算法, MonteCarlo方法等. Euler方法等純數學(xué)方法能很好地解出動(dòng)力學(xué)方程,但數學(xué)性很強,難以體現解法的物理意義.遺傳算法等是比較新型的算法,具體的應用也還在探索當中, Monte carlo方法已經(jīng)在分子動(dòng)力學(xué)中有詳細的應用,但其在化學(xué)動(dòng)力學(xué)中的應用遠遠不及分子動(dòng)力學(xué).本文作者受到 Arrhenius定理的啟發(fā),提出了一種處理化學(xué)反應動(dòng)力學(xué)問(wèn)題的數值方法,可以計算典型的化學(xué)反應動(dòng)力學(xué)問(wèn)題.該方法不僅有數學(xué)算法的運算特點(diǎn),模擬過(guò)程與化學(xué)動(dòng)力學(xué)過(guò)程相似,具有明確的物理意義,而且還能夠解決各種復雜化學(xué)動(dòng)力學(xué)問(wèn)題模擬方法介紹中國煤化工Arrhenius定理指出:在恒定濃度下,基元反應的速率與CN以朵系可用下面的數收稿日期:2010-11-30作者簡(jiǎn)介:鐘巍(1986—),男,研究方向為化學(xué)反應動(dòng)力學(xué)過(guò)程的數值模擬及其應用,E-mail:lengshui222@163.com.鐘巍:基于 Arrhenius定理的化學(xué)動(dòng)力學(xué)數值計算法學(xué)形式表示,即:k= Ae -E,rT上式中k為反應速率常數,R為理想氣體通用常數,T為反應體系的熱力學(xué)溫度,A和E是兩個(gè)與反應溫度及濃度無(wú)關(guān)的常數.觀(guān)察該表達式可以發(fā)現,方程式右邊實(shí)質(zhì)上是考慮了溫度等影響的一個(gè)指標,用來(lái)表征反應速率常數k的大小.于是,通過(guò)上面的啟發(fā)得到了下面的化學(xué)動(dòng)力學(xué)數值模擬方法某化學(xué)反應的速率計算公式即反應動(dòng)力學(xué)微分方程表示如下:出f i (t,ci, L, cn)i=1, 2,L,N(1)中t表示反應時(shí)間,c;為組分的濃度,f為組分i的化學(xué)動(dòng)力學(xué)表達式,稍作變形為,f (t,C1, L,CN)i=1,2,L,N(2)本文提供的模擬方法和上述 Arrhenius定理是相符合的,事實(shí)上(2)式左邊是某個(gè)時(shí)間步反應掉的物質(zhì)與原物質(zhì)的比值,這可以看作是數值模擬時(shí)化學(xué)反應速率常數的相對值,右邊則可以看作是該速率常數相對值的迭代計算方法.于是,通過(guò)可以將組分進(jìn)行等分,給每等分設計出一個(gè)合理的考慮了溫度等影響的一個(gè)指標,使其分別與相對速率常數值相比較以判斷該組分是否發(fā)生化學(xué)反應.得到具體的計算機模擬解法求解各組分的濃度過(guò)程如下:1)將組分分成N等分(N越大,計算結果越精確),各等分依次賦值i/Ni=1,2,L,N表征該等分發(fā)生反應的能力,即數值模擬中考慮溫度等影響的一個(gè)指標,下面簡(jiǎn)稱(chēng)為數值指標;2)求組分的濃度c(t).令k=f(t,c),dt,若=0,則組分濃度c;(t)不變,即c(t)=c(0);否則將k,|與所有等分的數值指標進(jìn)行比較,記數值指標小于|k|的等分數為S,c(t)=c1(0)(3)用最新迭代所得的值代替組分的初始濃度,重復步驟2),直到反應終點(diǎn);4)輸出計算結果在選擇迭代初始值時(shí),如果(2)式中c為0,則先用其他方法對第一個(gè)時(shí)間步估算,得出第一時(shí)間步的估計值作為初始值,迭代也從第一時(shí)間步以后開(kāi)始,若第一時(shí)間步仍為0,則估算到第二時(shí)間步,依此類(lèi)推.2模擬結果與分析下面使用上述模擬方法對簡(jiǎn)單一級反應,平行反應和復雜的綜合反應進(jìn)行數值模擬,驗證數值算法的可行性,對于下面反應中的準確解的得到辦法,本文不作贊述2.1簡(jiǎn)單一級反應考慮簡(jiǎn)單一級反應如下AP其動(dòng)力學(xué)微分方程式如下:取反應速率系數k=0.5L·moll·s-1,物質(zhì)A的初始濃度為cA(0)=1.5mol/L按照本文提供的方法編程模擬該問(wèn)題,時(shí)間步長(cháng)取為M=1/2000s,組分等分數N=10000,模擬解與解析解以及相對誤差見(jiàn)表1衰1物質(zhì)A與解析解的Table 1 The relative error of中國煤化工t/(s)0.10CNMHG 0.30解析解1.46301.42681.39161.35731.32371.2911模擬解1.46311.42691.39181.3240相對誤差(%)0.0060.0070.0140.0150.0220.0152011年模擬結果如圖1所示,由表1和圖1可見(jiàn),該方法模擬結果和解析解的擬合程度很好2.2平行反應乙烯酮可由高溫下醋酸裂解制備,副產(chǎn)甲烷.反應如下(CH=CO+H OCH+ C已知溫度為916℃,k1=4.65s-1,醋酸的初始濃度為1mol/L,產(chǎn)物初始濃度均為0mol/L.計算反應0.5s后醋酸、乙烯酮和甲烷的濃度,假設A,B,C分別代表CH3COOH,CH2=CO,CH4,化學(xué)反應速率方程如下:kICA(t)+kzCA(t)kcA(t)(4dce (t)kpcA(t)按照本文提供的方法模擬該問(wèn)題,時(shí)間步長(cháng)取為1/2000s,組分等分數N=100000,模擬結果與解析解及相對誤差見(jiàn)表2(只給出產(chǎn)物C).由表2可見(jiàn),該方法精度很高,且當時(shí)間增大時(shí),收斂性良好衰2產(chǎn)物C與解析解的相對誤差Table 2 The relative error of product Ct/(s)0.41.0解析解0.18530.29320.35610.39280.41420.4266模擬解0.18560.29380.35680.41420.4261相對誤差/(%)0.160.1012模擬結果如圖2所示,由圖2可見(jiàn),本文方法模擬結果和解析解能很好的吻合棱擬孵08摩解0.30.0.20.3040.5060.7080.900圖1簡(jiǎn)單一級反應物質(zhì)濃度隨時(shí)間的變化圖圖2平行反應物質(zhì)濃度隨時(shí)間變化Fig. I Plot of concentrations against time forig. 2 Plot of concentrations against time forfirst-order kinetics2.3綜合反應以制備膽甾醇異硫脲甲苯磺酸滃鹽( Cholesteryl Tosylate)的反應為例分析綜合反應的情形.針對該反應R.G. Pearson提出的反應機理如下:B+C中國煤化工反應條件為,[A]=0.00425mol/L,[B]=[C]=[D]=0CNMHGK, =2.96 L-min-1,k3=0.0021Lmol-min-1,計算到350min時(shí)各物質(zhì)的濃度反應速率方程為:第3期鐘巍:基于 Arrhenius定理的化學(xué)動(dòng)力學(xué)數值計算法d[Al=-(k1+k3)[A].dt dt= k,.,=kLAJ-k2LBJRd[d]k2[B][C]+k[A],=k2[B]2+k[A]這是一個(gè)復雜的綜合反應,求解析解過(guò)程很復雜.為了體現本文提供方法的優(yōu)點(diǎn),作者將本文方法得到的數值解與陸兆仁用 Monte Carlo方法得到的數值解進(jìn)行了對比.模擬結果對比如表3.用該方法模擬時(shí),考慮到驗證算法在長(cháng)時(shí)間下的穩定性,且該反應速率常數值偏小,于是將組分等分數取為N=1×106,時(shí)間步長(cháng)取為2s.由表3可見(jiàn),雖然各組分濃度本身值比較小,三個(gè)反應速率系數差距較大,計算的時(shí)間也更長(cháng),但本文方法仍與準確解(作者使用傳統數值方法求得)相對誤差很小,而 Monte carlo方法針對這類(lèi)復雜的反應誤差較大衰3產(chǎn)物D與解析解的相對誤差Table 3 The relative error of product D80260MC方法相對誤差/(%)6.06245.70554.69824.56424.8904.7701本文方法相對誤差/(%)0.09170.0650.0504模擬結果如圖3所示,左圖為 Monte carlo方法數值模擬結果,右圖為本文方法模擬結果.從圖3也可以看到本文的方法得到的數值解精度更高物C解值2.5Monte carto方法本文方法圖3綜合反應物質(zhì)濃度隨時(shí)間變化Fig 3 Plot of concentrations against time for the comprehensive reaction3總結1)受 Arrhenius定理的啟發(fā),提出了一種新的化學(xué)反應動(dòng)力學(xué)的數值模擬解法,用該方法可以得到精度較高的數值解2)該方法之所以能得到精度很高的模擬結果,原因在于可以任意提高組分等分數N的值,從而可以使得模擬隨著(zhù)N值的增大而接近化學(xué)反應的真實(shí)情況,具有明確中國煤化工3)現代化學(xué)反應動(dòng)力學(xué)理論研究在微觀(guān)角度已經(jīng)深入到CNMHG學(xué)反應涉及成千上萬(wàn)的基元反應,用穩定高效的數值方法研究這些反應很有必要,因此,該方法具有一定的應用價(jià)值2011年參考文獻[1]許越,化學(xué)反應動(dòng)力學(xué)[M].北京:化學(xué)工業(yè)出版社,2008:2-16[2]黃雪征.化學(xué)反應動(dòng)力學(xué)的計算與計算機模擬[D].2004,北京化工大學(xué)[3]關(guān)治,陸金甫,數值分析基礎[M].高等教育出版社,2010[4]王建國,李永旺,陳涌英,一些簡(jiǎn)單和復雜反應的 Monte Carlo模擬[]化學(xué)通報,1993(12);55-58[5]韓德剛,高盤(pán)良.化學(xué)動(dòng)力學(xué)基礎[M].北京:北京大學(xué)出版社,2001:3-8.[6] ESPENSON J H. Chemical kinetics and reaction mechanism[M]. New York: McGraw-Hill Inc, 1995: 15-18.[7]英徐根,應用計算物理化學(xué)[M.上海:上?？茖W(xué)技術(shù)文獻出版社,1989:283-321[8] PEARSON R G. Mechanism of formation of cholesterylisothiuromium tosylate form thiourea and cholesteryl to slate[J]J Am Chem Soc,1951,72:4149-4153.[9]陸兆仁.化學(xué)動(dòng)力學(xué)方程的計算機模擬求解法[]紡織高?；A科學(xué)學(xué)報,1998,11(3):221-225.(上接第47頁(yè))[6] WIPF P, CUNNINGHAM A. A solid phase protocol of the Biginelli dihydropyrimidine synthesis suitable for combinatorialchemistry[J]. Tetrahedron Lett, 1995, 36(43):7819-7822[7] SUN Qi, WANG Yi Qiang, GE Ze Mei, et al. Highly efficient solvent-free synthesis of dihydropyrimidinones catalyzed byzinc chloride[j]. Synthesis, 2004,(7): 1047-1051[8] SU Wei Ke, LI Jian Jun, ZHENG Zhi Guo, et al. One- pot synthesis of dihydropyrimidiones catalyzed by strontium(II)tri-flate under solvent-free conditions[J]. Tetrahedron Lett, 2005, 46(36):6037-6040.[9] GHOSH R, MAITI S, CHAKRABORTY A. In(OTf), catalyzed one-pot synthesis of 3, 4-dihydropyrimidin-2(IH)-ones[ J mol Catal A: Chem, 2004, 217(1):47-50.[10] RANU B C, HAJRA A, JANA U. Indium(IlIchloride-catalyzed one- pot synthesis of dihydropyrimidinones by a three-component coupling of 1, 3-dicarbonyl compounds, aldehydes, and urea: an improved procedure for the Biginelli reaction]. J Org Chem,2000,65(19):6270—6272[11] RAMALINGAN C, KWAKY W. Tetrachlorosilane catalyzed multicomponent one-step fusion of biopertinent pyrimidineheterocycles[J]. Tetrahedron, 2008, 64(22):5023-5031.[12] DILIP C B, JADAB C S. Three component condensations catalyzed by iodine alumina for the synthesis of substituted 3dihydropyrimidin-2(1H)-ones under microwave irradiation and solvent-free conditions [J]. Tetrahedron Lett,2005,46(7):1159—1160[13] HASHEN S, MAHBOUBEH J. Al, O,/MeSO, H: a novel and recyclable catalyst for one-pot synthesis of 3, 4-dihydro-pyrimidinones or their sulfur derivatives in Biginelli condensation[JJ. Synth Commun, 2009, 39(6):958-979[14] MANDAL P K, MISRA K. HCIO. -SiO, catalyzed multicomponent reactions for the synthesis of privileged heterocyclicstructures[J]. Lett Org Chem, 2006, 3(11):848-853.中國煤化工CNMHG