24小时热门版块排行榜    

查看: 2534  |  回复: 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的回帖

huab1984666

新虫 (著名写手)

【答案】应助回帖

用MATLAB 很简单的,只要有初始值
春风又绿江南岸,明月何时照我还。
12楼2014-11-26 19:43:05
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 范进中举一文的中心思想 +4 炎黄贵胄 2026-08-22 5/250 2026-08-24 02:09 by 流流伤
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见