24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2123  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+3 chenny174 2026-04-10 17/850 2026-04-11 00:26 by 骑牛渡寒江
[考研] 还有化工二轮调剂的学校吗 5+13 化工人999 2026-04-09 47/2350 2026-04-11 00:20 by 骑牛渡寒江
[考研] 303求调剂 +7 SereinQ 2026-04-10 8/400 2026-04-10 21:08 by gong120082
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 085500求调剂材料 +8 易11122 2026-04-09 8/400 2026-04-09 23:15 by parmtree
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 求调剂 +7 chenxrlkx 2026-04-05 9/450 2026-04-09 09:04 by wj165256
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考博] 材料方向考博,求推荐 +3 言语aaa 2026-04-05 4/200 2026-04-08 22:22 by nxgogo
[考研] 土木水利专硕276分求调剂 +6 我想上学!!6 2026-04-05 9/450 2026-04-08 17:45 by 宋小宝HQ
[考研] 327求调剂 +12 Xxjc1107. 2026-04-06 12/600 2026-04-08 16:46 by luoyongfeng
[考研] 081700,311,求调剂 +17 冬十三 2026-04-04 18/900 2026-04-07 12:50 by Sammy2
[考研] 292求调剂 +4 lilllllxccc 2026-04-05 5/250 2026-04-07 09:29 by 纺大杨老师
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 第一志愿东南大学物理313,有科研竞赛获奖经历,希望物理复试调剂 +3 马内橙 2026-04-05 3/150 2026-04-06 10:32 by 蓝云思雨
[考研] 材料调剂 +14 壹贰贰亿 2026-04-04 14/700 2026-04-05 23:31 by 来看流星雨10
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
信息提示
请填处理意见