24小时热门版块排行榜    

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

dante1991

新虫 (初入文坛)

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

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

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


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

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 6/300 2026-10-26 02:42 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:38 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:05 by cLOojpOxkX0G
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:03 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 01:50 by cLOojpOxkX0G
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:46 by cLOojpOxkX0G
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-26 01:24 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 7xHjyB3bAiwg 2026-10-24 7/350 2026-10-26 01:11 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:41 by cLOojpOxkX0G
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:32 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 7xHjyB3bAiwg 2026-10-23 9/450 2026-10-26 00:28 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:19 by cLOojpOxkX0G
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 GVGRDgNfv7P4 2026-10-23 10/500 2026-10-26 00:15 by cLOojpOxkX0G
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 GVGRDgNfv7P4 2026-10-23 5/250 2026-10-26 00:08 by cLOojpOxkX0G
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 WqtG0Aem8Ufl 2026-10-24 4/200 2026-10-25 20:14 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:11 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-24 23:35 by XcOYG53d3x4p
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 15:41 by Equinoxhua
[考博] 27秋申博自荐 +4 且听风吟时 2026-10-21 4/200 2026-10-23 13:26 by 大大大肥虫
信息提示
请填处理意见