deadhunter
2011-05-19, 21:31
clc;clear;
function dxdl=Equations(l,x)
global dxdl, l,Fa0,s,d,l,x(1),x(2),cp,cpCH4,cpCO,cpCO2,cpH2,cpO2,cpN2,cpH2O,FCO,FCO2,FCH4,FH2,FO2,FN2,FH2O,Fi,A,B,C,D,k,Z,Kp,detH,r;
d=3.6;s=d^2*3.14/4;Fa0=44.348;
cpCH4=25.360+168.678*e-4*x(2)+713.121*e-7*x(2)^2-408.371e-10*x(2)^3;
cpCO=26.5366+7.6831e-3*x(2)-1.1719e-6*x(2)^2;
cpCO2=26.748+42.258e-3*x(2)-14.547e-6*x(2)^2;
cpN2=27.016+5.812e-3*x(2)-0.289e-6*x(2)^2;
cpH2O=30.204+9.933e-3*x(2)+1.117e-6*x(2)^2;
cpH2=29.062-0.820e-3*x(2)+1.9903e-6*x(2)^2;
cpO2=25.594+13.257e-3*x(2)-4.205e-6*x(2)^2;
Fi*cp=FCH4*cpCH4+FCO*cpCO+FCO2*cpCO2+FH2*cpH2+FO2*cpO2+FN2*cpN2+FH2O*cpH2O;
A=FCO/(FCH4+FCO2+FCO+FO2+FN2+FH2+FH2O);
B=FH2O/(FCH4+FCO2+FCO+FO2+FN2+FH2+FH2O);
C=FCO2/(FCH4+FCO2+FCO+FO2+FN2+FH2+FH2O);
D=FH2/(FCH4+FCO2+FCO+FO2+FN2+FH2+FH2O);
k=1.4869e4*exp(-15873/1.987*x(2));
Kp=exp(5025.163/x(2)-0.0936*ln(x(2))+1.4555e-3*x(2)-2.4887e-7*x(2)^2-5.2894);
detH=4.187*(10000+0.219*x(2)-2.845e-3*x(2)^2+0.9703e-6*x(2)^3);
r=k*((A-A*x(1))*(B-A*x(1))-(C+A*x(1))*(D+A*x(1))/Kp);
dxdl1=r*s/Fa0;
dxdl2=r*(detH)*s/Fi*cp;
(l,x)=ode45(@Equations,[0:0.02:10],[0;593.5]);
这是计算反应转化率、温度随反应器长度的关系的程序，现在不知道问题出在哪里，运行不了，还望各位大神帮助
function dxdl=Equations(l,x)
global dxdl, l,Fa0,s,d,l,x(1),x(2),cp,cpCH4,cpCO,cpCO2,cpH2,cpO2,cpN2,cpH2O,FCO,FCO2,FCH4,FH2,FO2,FN2,FH2O,Fi,A,B,C,D,k,Z,Kp,detH,r;
d=3.6;s=d^2*3.14/4;Fa0=44.348;
cpCH4=25.360+168.678*e-4*x(2)+713.121*e-7*x(2)^2-408.371e-10*x(2)^3;
cpCO=26.5366+7.6831e-3*x(2)-1.1719e-6*x(2)^2;
cpCO2=26.748+42.258e-3*x(2)-14.547e-6*x(2)^2;
cpN2=27.016+5.812e-3*x(2)-0.289e-6*x(2)^2;
cpH2O=30.204+9.933e-3*x(2)+1.117e-6*x(2)^2;
cpH2=29.062-0.820e-3*x(2)+1.9903e-6*x(2)^2;
cpO2=25.594+13.257e-3*x(2)-4.205e-6*x(2)^2;
Fi*cp=FCH4*cpCH4+FCO*cpCO+FCO2*cpCO2+FH2*cpH2+FO2*cpO2+FN2*cpN2+FH2O*cpH2O;
A=FCO/(FCH4+FCO2+FCO+FO2+FN2+FH2+FH2O);
B=FH2O/(FCH4+FCO2+FCO+FO2+FN2+FH2+FH2O);
C=FCO2/(FCH4+FCO2+FCO+FO2+FN2+FH2+FH2O);
D=FH2/(FCH4+FCO2+FCO+FO2+FN2+FH2+FH2O);
k=1.4869e4*exp(-15873/1.987*x(2));
Kp=exp(5025.163/x(2)-0.0936*ln(x(2))+1.4555e-3*x(2)-2.4887e-7*x(2)^2-5.2894);
detH=4.187*(10000+0.219*x(2)-2.845e-3*x(2)^2+0.9703e-6*x(2)^3);
r=k*((A-A*x(1))*(B-A*x(1))-(C+A*x(1))*(D+A*x(1))/Kp);
dxdl1=r*s/Fa0;
dxdl2=r*(detH)*s/Fi*cp;
(l,x)=ode45(@Equations,[0:0.02:10],[0;593.5]);
这是计算反应转化率、温度随反应器长度的关系的程序，现在不知道问题出在哪里，运行不了，还望各位大神帮助