24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2124  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

19880520

新虫 (初入文坛)

[求助] 求尿素热解化学反应动力学参数估算方法 已有3人参与

小弟是搞内燃机尾气排放的,最近在研究SCR系统尿素结晶机理。现在想根据尿素的热重分析曲线(TGA曲线)来求化学反应动力学参数。
化学反应机理如下图:
求尿素热解化学反应动力学参数估算方法
其动力学方程如下:
反应速率
R(1)=k(1)*urea
R(2)=k(2)*NH4
R(3)=k(3)*NCO*H
R(4)=k(4)*urea*NCO*H
R(5)=k(5)*biuret
R(6)=k(6)*biuret*NCO*H
R(7)=k(7)*cya
R(8)=k(8)*cya*NCO*H
R(9)=k(9)*amme
其中k(i)可以用阿累尼乌斯方程计算,R=8.31,T为温度,e是自然对数
k(1)=A(1)*e^(Ea(1)/RT)
k(2)=A(1)*e^(Ea(2)/RT)
k(3)=A(3)*e^(Ea(3)/RT)
k(4)=A(4)*e^(Ea(4)/RT)
k(5)=A(5)*e^(Ea(5)/RT)
k(6)=A(6)*e^(Ea(6)/RT)
k(7)=A(7)*e^(Ea(7)/RT)
k(8)=A(8)*e^(Ea(8)/RT)
k(9)=A(9)*e^(Ea(9)/RT)
各组分浓度变化率
durea/dt=(-R(1)-R(4)+R(5))
   dNH3/dt = (R(2)+R(6))
    dHCNO/dt = R(3)
    dbiuret/dt= (R(4)-R(5)-R(6))
    dcya/dt= (R(6)-R(7)-R(8))
    damme/dt= (R(8)-R(9))
    dHCN/dt= R(9)
    dNH/dt= R(9)
    dCO2/dt= R(8)
    dNH4/dt= (R(1)-R(2))
    dNCO/dt= (R(1)-R(3)-R(4)+R(5)-R(6)+3*R(7)-R(8)+2*R(9))
    dH/dt=H_0+(R(2)-R(3)-R(4)+R(5)-R(6)+3*R(7)-R(8)+2*R(9))

我现在有这样一组热重分析TGA试验数据,温度升高率为10K/min,其中mass/initial mass为urea+biuret+cya+amme之和。我现在想用这组试验数据拟合出上述九个化学反应的频率因子和活化能,请问要怎么做?我被这个问题困扰了好几个月了,一直没有找到很好的方法。
求尿素热解化学反应动力学参数估算方法-1
数据如下:
temperature        mass/initial mass
349.422        0.999276
370.417        0.998046
389.63        0.992421
407.512        0.978029
424.514        0.950498
436.611        0.91416
446.498        0.855916
455.503        0.785992
463.618        0.71314
472.198        0.622788
480.786        0.523683
489.361        0.437708
502.854        0.347409
512.735        0.295
527.52        0.251398
547.622        0.250159
575.762        0.250467
594.969        0.250677
609.27        0.24354
624.93        0.217452
637.484        0.172367
650.497        0.114157
666.176        0.070564
686.315        0.032856
703.32        0.002407
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
9楼: Originally posted by 19880520 at 2014-03-01 20:10:07
谢谢1stopt真的是好方便啊。不过结果不太对,A和Ea应该都是大于零的。可惜没有这个软件啊,不知道有没有人会用matlab来编...

"都是大于零"这种限制条件开始就该说出来啊。
10楼2014-03-01 22:59:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

19880520

新虫 (初入文坛)

没啥金币,但是真心求教啊
2楼2014-02-28 10:25:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
各自的初始值都知道吗,比如urea,NH3,HCNO,biuret,cya...的起始值。
3楼2014-02-28 11:16:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

19880520

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by dingd at 2014-02-28 11:16:56
各自的初始值都知道吗,比如urea,NH3,HCNO,biuret,cya...的起始值。

一开始全部是尿素,即urea=1,其他都是0
4楼2014-02-28 18:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂 +6 luoziheng 2026-04-10 6/300 2026-04-10 23:22 by babysonlkd
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +27 努力奋斗112 2026-04-07 30/1500 2026-04-10 15:06 by Kilig0317
[考研] 277求调剂 +19 倪建设 2026-04-06 19/950 2026-04-10 09:24 by guosr9609
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 086000调剂 +4 十七sa 2026-04-07 4/200 2026-04-09 09:27 by 壹往無前
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 0854求调剂 +9 亨氏番茄沙司 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7 晚风不晚 2026-04-04 7/350 2026-04-06 14:06 by houyaoxu
[考研] 285求调剂 +8 AZMK 2026-04-04 11/550 2026-04-06 13:56 by BruceLiu320
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 341求调剂 +3 学无止境,冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
信息提示
请填处理意见