24小时热门版块排行榜    

查看: 386  |  回复: 1

greatwu1302

木虫 (初入文坛)

[交流] 【求助】请帮助解微分方程组 已有1人参与

% 物料平衡
dxdz = rhoB*r*Mm /(yA0 * MG);
% 热量衡算
dTdz = 4*U*(T-TJ)/(dt*MG*Cp)-(rhoB*r*H)/(MG*Cp);

其中r为反应速率,
% 反应速度, kmol/kg catalyst hr
r=yita*k*kH*PA*PH/(1+kH*PH+KVC*PVC);

PA=P0*yA0*(1-x);
PH=P0*yH0*(1-x);
PVC=P0*x;
% 速度常数, kmol/kg catalyst hr
if T<403
    k=9.21*(10^13)*exp(-27.3*1000/(8.314*T))+0.0476;
else
    k=1.714*(10^14)*exp(-27.55*1000/(8.314*T));
end
if T<403
    kH=19.32*exp(2.72*1000/(8.314*T));
else
    kH=1.09*10^20*exp(39.6*1000/(8.314*T));
end

KVC=1.783/10^6*exp(12900/(8.314*T))-2.53;
yita=0.295*exp(590/T)-120*exp(4.5*x);

对于微分方程组dzdx和dTdx怎么用ODE45来解呀,请高人帮助做这个程序
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
robert2020(金币+1):多谢解答,辛苦了! 2010-06-12 20:03:51
robert2020(金币+1):多谢应助,辛苦了! 2010-06-12 20:06:54
到下面的地方去下载《实用化工计算机模拟-MATLAB在化学工程中的应用》
你将受益匪浅

http://muchong.com/bbs/viewthread.php?tid=2083386&fpage=1
2楼2010-06-05 18:21:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 greatwu1302 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见