燃烧学第三章作业.docx

第三章作业 使用编程语言求解常微分方程的方法，求解零维系统自 燃问题： k0 ――频率因子，初值取 100.0 E――活化能，初值取1.E5 R气体常数，初值取8.314 C——可燃混合物中反应物浓度，初值取 1.0 n――反应级数，初值取1.0 V ――容器体积，初值取1.0 Q可燃混合物的燃烧热，初值取 2.E7 T――容器内可燃混合物温度，初值取与 T0相等的数值。 h――散热系数，初值取5.0 S――容器壁散热面积，初值取6.0 T0――容器壁温度，初值取800.0 定解条件： t=0 时， T=T0。 以上方程描述了一个零维系统的温度从 t=0开始随时间变化的 过程。 a） 使用在第一章已经介绍过的Euler法，求解上述定解条件下的 常微分方程，获得系统温度 T随时间t的变化曲线。 b） 分别令系统初温 TO=3O0 400、500、600、700、800、900、 1000、1100（其他参数不变），获得不同初始温度下的系统升温曲线， 并讨论系统初温对热自燃过程的影响。 c） 分别令散热系数h= 1.0?10.0 （T0保持800K）,获得不同初 始温度下的系统升温曲线，并讨论散热系数对热自燃过程的影响。 问题求解过程如下： 一：程序如下： #include stdio.h #include math.h #include fstream #in clude iostream using namespacestd; int main（） { ofstream ofile; ofile.open（ f:\\myfile.txt ）; // 打开文件 double k0 = 100, R = 8.314, C = 1.0, n = 1.0, V = 1.0, yita = 5.0, S = 6.0; // 定义参 数 double rou = 1.0, cv = 4.02; double E = 1.0e5, Q = 2.0e7; double d, tao[101], T[101]; double h = 0.1; //设置时间间隔 int i; double f（ double x, double y[100]）; tao[0] = 0.0, T[0] = 800.0; // 设置初值 printf( %f,%f\n , T[0], tao[0]); for 仃[0] = 300; T[0]= 1100; T[0] = T[0] + 100) // 求解方程在 T0=300到 1100 时的温度变化值 { if仃[0]==800) //当初始温度为800度时改变 n使其在1.0至10变化 for (yita = 1.0; yita = 10; yita++) { ofile yita en dl; ofile \n endl; for (i = 0; i 100; i++) { d = (k0*exp(-E / R / T[i])*pow(C, n)*V*Q - yita*S*(T[i] - T[0])) / V / rou / cv; //求解温度随时间变化率 T[i + 1] = T[i] + h*d; // 计算下 一时刻温度值 tao[i + 1] = tao[0] + h*i; printf( %f,%f\n , T[i + 1], tao[i + 1]); ofile T[i] endl; // 输岀 T 至文件 myfile.txt | } } else for (i = 0; i 100; i++) // 若T0不 等于800，使yita=5.0，不改变yita值 { yita = 5.0; d = (k0*exp(-E / R / T[i])*pow(C, n)*V*Q - yita*S*(T[i] - T[0])) / V / rou / cv; T[i + 1] = T[i] + h*d; tao[i + 1] = tao[0] + h*i; printf( %f,%f\n , T[i + 1], tao[i + 1]); ofile T[i] en dl; } ofile \n en dl; } ofile.close(); // 关闭文件 二：程序求解结果 图一：反应温度与初始温度、散热系数关系图 6001(J图二：反应温度与散热系数关系图 TO=30O°C T0=^Q0n 600 1(J 图二：反应温度与散热系数关系图 TO=30O°C T0=^Q0nC —To=so(rc TO=600nC T0=?D0*C T0=80OnC T0=S(X)X TO=1DOOY TD-11D0T T0-800? yita=1.0 10=800, yita^2 0 U