第三章作业。
使用编程语言求解常微分方程的方法,求解零维系统自燃问题:
k0——频率因子,初值取100.0
e——活化能,初值取
r——气体常数,初值取8.314
c——可燃混合物中反应物浓度,初值取1.0
n——反应级数,初值取1.0 v——容器体积,初值取1.0
q——可燃混合物的燃烧热,初值取
t——容器内可燃混合物温度,初值取与t0相等的数值。
h——散热系数,初值取5.0
s——容器壁散热面积,初值取6.0
t0——容器壁温度,初值取800.0
定解条件:t=0时,t=t0。
以上方程描述了一个零维系统的温度从t=0开始随时间变化的过程。
a)使用在第一章已经介绍过的euler法,求解上述定解条件下的常微分方程,获得系统温度t随时间t的变化曲线。
b)分别令系统初温t(其他参数不变),获得不同初始温度下的系统升温曲线,并讨论系统初温对热自燃过程的影响。
c)分别令散热系数h=1.0~10.0(t0保持800k),获得不同初始温度下的系统升温曲线,并讨论散热系数对热自燃过程的影响。
问题求解过程如下:
一:程序如下:
#include<>
#include<>
#include
#include
using namespace std;
int main()
ofstream ofile;
"f:\\打开文件。
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", t[0], tao[0]);
for (t[0] =300; t[0]<=1100; t[0] =t[0] +100求解方程在t0=300到1100时的温度变化值。
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", t[i + 1], tao[i + 1]);
ofile < ofile < 关闭文件。二:程序求解结果。 图一:反应温度与初始温度、散热系数关系图。 图二:反应温度与散热系数关系图。 图三:反应温度与初始温度关系图。 四:结果分析。 由图二易知,随着散热系数的增大,燃烧稳定性下降,从正常燃烧的状态变为不能燃烧的工况; 由图三易知,随着初始温度的增加,燃烧稳定性改善,从不能着火的状态变为能燃烧的工况。 v s 顺序执行下述两个动作 1.s值加1,即s s 1 2.如果s 0,则该进程继续运行 3.如果s 0,则唤醒等待信号量s阻塞队列中的头一个进程 把阻塞态改为就绪态 执行v操作的进程继续运行。procedure s var s semaphore begin s s 1 if s 0 then ... 1 顺序栈空 栈满条件2 链栈栈空 栈满条件。3 循环队列队空 队满条件,如何表示队列中数据元素的个数4 链队列队空 队满条件。5 以下运算实现在顺序栈上的进栈,请在 处用适当的语句予以填充。int push sqstacktp sq,datatype x if sp top sqstack max... 1.论述各类绿地的环境特点和树种的选择。一 高层建筑中的狭窄街巷绿地绿地内的环境特点 直射辐射量少,日照时间短 夏季气温偏低,冬季因受周围建筑物热辐射的影响,气温偏高 风速一般偏低,但有时会产生狭管效应,使风速增大。这些地方裸露土面极少,多为水泥铺装,严重阻碍了土壤与大气的水 气交换,且存在一定程度...第三章作业
第三章作业
第三章作业