| 查看: 2488 | 回复: 12 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
求尿素热解化学反应动力学参数估算方法 已有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之和。我现在想用这组试验数据拟合出上述九个化学反应的频率因子和活化能，请问要怎么做？我被这个问题困扰了好几个月了，一直没有找到很好的方法。  数据如下： 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  |
回复此楼» 收录本帖的淘帖专辑推荐
 ANSYS仿真计算 |
» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有244人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- SCR烟气脱硝各路气体流量问题
已经有7人回复
- SCR中尿素热解制备氨气
已经有5人回复
- 【求助】复旦考研药分
已经有4人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2014-02-28 11:16:56
2楼2014-02-28 10:25:26
4楼2014-02-28 18:09:25
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
感谢参与，应助指数 +1
|
能有高版本的1stopt软件做比较好。 MATLAB的话，参见： http://muchong.com/bbs/viewthread.php?tid=6425538&authorid=1122189 |
5楼2014-02-28 18:49:16
