燃燒學(xué)第三章作業(yè)

1、第三章作業(yè)使用編程語(yǔ)言求解常微分方程的方法，求解零維系統(tǒng)自 燃問(wèn)題：EdT 5(-麗)0-小(卩-7；)drV pc vk()頻率因子，初值取1()().()E活化能，初值取1.E5R氣體常數(shù)，初值取8.314C可燃混合物中反應(yīng)物濃度，初值取1.()n反應(yīng)級(jí)數(shù)，初值取1()V容器體積，初值取1.()Q可燃混合物的燃燒熱，初值取2.E7T容器內(nèi)可燃混合物溫度，初值取與TO相等的數(shù)值。h散熱系數(shù)，初值取5.()S容器壁散熱面積，初值取6.()TO容器壁溫度，初值取800.()定解條件：t=0時(shí)，T=T0o以上方程描述了一個(gè)零維系統(tǒng)的溫度從t=()開(kāi)始隨時(shí)間變化的 過(guò)程。勾使用在第一章已經(jīng)介紹過(guò)的E

2、ukr法，求解上述定解條件下的 常微分方程，獲得系統(tǒng)溫度T隨時(shí)間t的變化曲線。b) 分別令系統(tǒng)初溫 T0=300 400、500、600、700、800、900、1()0()、 1100 (其他參數(shù)不變)，獲得不同初始溫度下的系統(tǒng)升溫曲線，并討 論系統(tǒng)初溫對(duì)熱自燃過(guò)程的影響。c) 分別令散熱系數(shù)h = l.()1().() (TO保持800K),獲得不同初始 溫度下的系統(tǒng)升溫曲線，并討論散熱系數(shù)對(duì)熱自燃過(guò)程的影響。問(wèn)題求解過(guò)程如下:一：程序如下：#includc#includc#includc#includc using iiamcspacc std;int main()ofsircnm of

3、ile；ofilc.openf,f:myfile.txtM); /打開(kāi)夭件double kO = 100,8.314, C 二 1.0, n= 1.0, V= 1.0, yi恰=50, S = 60;定義參數(shù)double rou = 1Q cv = 4.02;double E = 1.0c5 Q = 2.0c7;double d, iao101, T101;/設(shè)應(yīng)時(shí)間何隔誅置初值/求解方程龍T()=300到1100廿的溫皮變/當(dāng)初始溫度為800度時(shí)改變耳使double h = 0；int i;double f(doublc x, double y100); too0 = 0.0, T0 = 8

4、00.0;T0, taoOJ);for (T0 = 300; T0= 1100; T0 = T0 + 100) 化值if(T0=800)其在l0至10變化fc)r (vita = 1.0; yita = 10; vita+)ofile yita cndl;精品ofilc ,n* cndl;for (i = 0; i Ti + 1, taoi + 1);ofile Ti endl;else計(jì)尊下一時(shí)/輸出T至文for (i = 0; i w(C, n)V*Q - yita*S*(Ti T0) / V / rou / cv; Ti + 1 = Ti + h*d；tao|i + 1 = tao0 +

5、 h*i; printf(H%f%fnH, Ti + 1, taop + 1);ofilc Ti cndl;ofile vv endl;ofilc.closeQ;關(guān)閉文伸二：程序求解結(jié)果圖一：反應(yīng)溫度與初始溫度、散熱系數(shù)關(guān)系圖圖二：反應(yīng)溫應(yīng)與散熱系數(shù)關(guān)系圖TO二800,yita=1.0 T0=800,yita=2.0TO二800,yita=3.0TO二800,yita=4.0 T0=800,yita=5.0T0=800,yita=6.0TO二800,yita=7.0TO二800, yita=8.0 TO二800,yita=9.0T0=800,yita=10.0200018001600140012001000800600400200T0=300CT0400CTO二500CTO二600CT0=700oCT0=80OCTO二900CT0=1000CT0=1100C0r-0246 8 10圖三：反應(yīng)遍度與初始遍廈關(guān)系圖四：結(jié)果分析由圖二易知，隨著散熱系數(shù)的增大，燃