《燃烧学》全套课件.pptx

下载文档

14
0
约1.5千字
约 477页
2024-01-08 发布于湖北
举报
版权申诉
保障服务

《燃烧学》全套课件.pptx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;示例：使用FORTRAN语言求解常微分方程的方法，求解零维系统自燃问题;k0——频率因子，初值取100.0

E——活化能，初值取1.E5

R——气体常数，初值取8.314

C——可燃混合物中反应物浓度，初值取1.0

n——反应级数，初值取1.0

V——容器体积，初值取1.0

Q——可燃混合物的燃烧热，初值取2.E7

T——容器内可燃混合物温度，初值取与T0相等的数值。

?——散热系数，初值取5.0

S——容器壁散热面积，初值取6.0

T0——容器壁温度，初值取800.0

定解条件：?=0时，T=T0。以上方程描述了一个零维系统的温度从?=0开始随时间变化的过程。

;Euler法求解常微分方程;programself_ignite

implicitnone

real::tao0,T0,tao,T,tao1,T1,dtao,taomax,Tmax

integer::i

i=0

tao0=0.

T0=800.

dtao=0.1

tao1=tao0

T1=T0

taomax=1000.

Tmax=1.e30

; open(01,file=t.dat)

dowhile(abs(tao)=taomax.and.abs(T)=Tmax)

i=i+1

tao=tao+dtao

T=T1+dtao*f(tao1,T1)

print*,i,tao,f(tao,T),T

write(01,*)tao,f(tao,T),T

tao1=tao

T1=T

enddo

close(01)

; contains

functionf(tao,T)

implicitnone

realf,tao,T,q1,q2,V,rho,cv

realk0,c,n,E,R,Q,yita,S

k0=100.

c=1.

n=1.

E=100.e3

R=8.314

V=1.

Q=2E7

yita=5

S=6.

rho=1.

cv=1000.

q1=k0*c**n*exp(-E/(R*T))*V*Q

q2=yita*S*(T-T0)

f=(q1-q2)/(V*rho*cv)

endfunction

endprogramself_ignite;求解结果：变化初始温度T0的情况;求解结果：变化散热系数?的情况;;;;;自选作业：;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;