24小时热门版块排行榜    

查看: 623  |  回复: 1
本帖产生 1 个 ICEPI ,点击这里进行查看

dante1991

新虫 (初入文坛)

[求助] 欧拉法求解微分方程求解反应速率系数已有1人参与

问题不复杂,H2S+CO2热解反应,根据实验结果推算出反应的速率系数。
小弟选择了四个可逆反应描述该过程,并且写出主要物质的质量守恒方程。
然后根据不同停留时间下的组分浓度去计算。
根据文献,由欧拉法可以求解。
小弟不太懂这个,手动积分求解,利用相邻两点组分浓度带入,可以求出各反应的正向速率系数,但是发现随着时间增长,反应趋于平衡,组分浓度变化不大时,求出的k值也逐渐减小。按道理同一温度下k值是不变的啊。
那么如何用欧拉法求解呢?

欧拉法求解微分方程求解反应速率系数
H2S-CO2热解结果图


欧拉法求解微分方程求解反应速率系数-1
质量守恒方程
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mbchen

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dante1991: 金币+20, 有帮助, 说明很详细 2017-01-09 08:22:53
冰点降温: 金币+3, ICEPI+1, 谢谢专家解答,辛苦了。解答很详细,具有普遍参考意义。 2017-01-09 09:40:47
无论怎么样的化学反应,总可以写出一阶微分方程dx(t)/dt=f(x(t)),其中x(t)为某产物的浓度。
然后用Euler法求一阶微分方程dx(t)/dt=f(x(t))的解,并且已知起始条件x(t=0)。

Euler法的原理如下:
将时间看成离散的,即t=0,h,2h,3h,...
原微分方程dx(t)/dt=f(x(t))近似于如下差分方程:
[x(t+h)-x(t)]/h = f(x(t)),即
x(t+h) = x(t) +h* f(x(t))
该差分方程在t=0,h,2h,...时分别为如下几个方程
x(h) = x(0) +h* f(x(0))
x(2h) = x(h) +h* f(x(h))
x(3h) = x(2h) +h* f(x(2h))
……
这样就依次得到上述微分方程(和起始条件)的近似解x(h),x(2h),x(3h),...。
2楼2017-01-08 18:29:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dante1991 的主题更新
信息提示
请填处理意见