版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

19880520

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 25.3
帖子: 30
在線: 32.2小時(shí)
蟲(chóng)號(hào): 2550410
注冊(cè): 2013-07-17
專業(yè): 傳熱傳質(zhì)學(xué)

[求助] 求尿素?zé)峤饣瘜W(xué)反應(yīng)動(dòng)力學(xué)參數(shù)估算方法已有3人參與

小弟是搞內(nèi)燃機(jī)尾氣排放的，最近在研究SCR系統(tǒng)尿素結(jié)晶機(jī)理�，F(xiàn)在想根據(jù)尿素的熱重分析曲線（TGA曲線）來(lái)求化學(xué)反應(yīng)動(dòng)力學(xué)參數(shù)。
化學(xué)反應(yīng)機(jī)理如下圖：

其動(dòng)力學(xué)方程如下：
反應(yīng)速率
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）可以用阿累尼烏斯方程計(jì)算，R=8.31,T為溫度，e是自然對(duì)數(shù)
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))

我現(xiàn)在有這樣一組熱重分析TGA試驗(yàn)數(shù)據(jù)，溫度升高率為10K/min，其中mass/initial mass為urea+biuret+cya+amme之和。我現(xiàn)在想用這組試驗(yàn)數(shù)據(jù)擬合出上述九個(gè)化學(xué)反應(yīng)的頻率因子和活化能，請(qǐng)問(wèn)要怎么做？我被這個(gè)問(wèn)題困擾了好幾個(gè)月了，一直沒(méi)有找到很好的方法。

求尿素?zé)峤饣瘜W(xué)反應(yīng)動(dòng)力學(xué)參數(shù)估算方法-1

數(shù)據(jù)如下：
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

回復(fù)此樓

1樓 2014-02-27 22:55:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

19880520

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 25.3
帖子: 30
在線: 32.2小時(shí)
蟲(chóng)號(hào): 2550410
注冊(cè): 2013-07-17
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

8樓: Originally posted by dingd at 2014-03-01 14:55:37
1stOpt試試，太化時(shí)間，大概算了下，不知理解對(duì)否，供參考。有時(shí)間自己慢慢算。

Constant R=8.31;
ConstStr
k1=a1*exp(ea1/(r*t)),k2=a1*exp(ea2/(r*t)),k3=a3*exp(ea3/(r*t)),
k4=a4*exp(ea ...

謝謝1stopt真的是好方便啊。不過(guò)結(jié)果不太對(duì)，A和Ea應(yīng)該都是大于零的�？上](méi)有這個(gè)軟件啊，不知道有沒(méi)有人會(huì)用matlab來(lái)編

贊一下

回復(fù)此樓

9樓2014-03-01 20:10:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 13 個(gè)回答

19880520

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 25.3
帖子: 30
在線: 32.2小時(shí)
蟲(chóng)號(hào): 2550410
注冊(cè): 2013-07-17
專業(yè): 傳熱傳質(zhì)學(xué)

沒(méi)啥金幣，但是真心求教啊

贊一下

回復(fù)此樓

2樓2014-02-28 10:25:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

各自的初始值都知道嗎，比如urea，NH3，HCNO，biuret，cya...的起始值。

贊一下

回復(fù)此樓

3樓2014-02-28 11:16:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

19880520

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 25.3
帖子: 30
在線: 32.2小時(shí)
蟲(chóng)號(hào): 2550410
注冊(cè): 2013-07-17
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

3樓: Originally posted by dingd at 2014-02-28 11:16:56
各自的初始值都知道嗎，比如urea，NH3，HCNO，biuret，cya...的起始值。

一開(kāi)始全部是尿素，即urea=1，其他都是0

贊一下

回復(fù)此樓

4樓2014-02-28 18:09:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 13 個(gè)回答

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9	Ymlll 2026-03-18	12/600	2026-03-20 00:01 by 23Postgrad
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5	Doleres 2026-03-19	5/250	2026-03-19 20:14 by 制度的
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢(mèng)人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu