版塊導(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í)熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者1173018735將贈(zèng)送您 5 個(gè)金幣

1173018735

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 33
帖子: 9
在線: 1.1小時(shí)
蟲號(hào): 35398980
注冊(cè): 2024-04-28
專業(yè): 工程熱物理與能源利用

[求助] fluent模擬流化床甲烷干重整反應(yīng)udf使用浮點(diǎn)溢出

我在復(fù)現(xiàn)一篇用歐拉雙流體方法復(fù)現(xiàn)流化床甲烷干重整反應(yīng)的論文，我找不出來有什么問題，udf中公式的表達(dá)是沒問題的，這個(gè)模型不使用UDF可以正常計(jì)算，但是一使用了UDF就會(huì)報(bào)錯(cuò)。而這個(gè)UDF在另外一個(gè)簡(jiǎn)單的模型里使用沒有問題。有大佬幫我看看什么問題嗎？
# include "udf.h"
# include "mem.h"

DEFINE_VR_RATE(rate, c, t, r, wk, yk, rate, rr_t)
{

#define UGC 8.314

real T = C_T(c, t);
real k1 = 1.29e6 * exp(-102065.0 / UGC / C_T(c, t));
real k2 = 0.35e6 * exp(-81030.0 / UGC / C_T(c, t));
real k3 = 6.95e3 * exp(-58893.0 / UGC / C_T(c, t));
real k4 = 5.55e9 * exp(-166397.0 / UGC / C_T(c, t));
real k5 = 1.34e15 * exp(-243835.0 / UGC / C_T(c, t));
real Kco2_1 = 2.61e-2 * exp(37641.0 / UGC / C_T(c, t));
real Kch4_1 = 2.60e-2 * exp(40684.0 / UGC / C_T(c, t));
real Kco2_2 = 0.5771 * exp(9262.0 / UGC / C_T(c, t));
real Kh2_2 = 1.494 * exp(6025.0 / UGC / C_T(c, t));
real Kch4_3 = 0.21 * exp(-567.0 / UGC / C_T(c, t));
real Kh2_3 = 5.18e7 * exp(-133210.0 / UGC / C_T(c, t));
real Kh2_4 = 1.83e13 * exp(-216145.0 / UGC / C_T(c, t));
real Kh2o_4 = 4.73e-6 * exp(97770.0 / UGC / C_T(c, t));
real Kch4_4 = 3.49 * exp(-0 / UGC / C_T(c, t));
real Kco_5 = 7.34e-6 * exp(100395.0 / UGC / C_T(c, t));
real Kco2_5 = 2.81e7 * exp(-104085.0 / UGC / C_T(c, t));
real Kp1 = 6.78e14 * exp(-259660.0 / UGC / C_T(c, t));
real Kp2 = 56.4971 * exp(-36580.0 / UGC / C_T(c, t));
real Kp3 = 2.95e5 * exp(-84400.0 / UGC / C_T(c, t));
real Kp4 = 1.3827e7 * exp(-125916.0 / UGC / C_T(c, t));
real Kp5 = 1.9393e9 * exp(-168527.0 / UGC / C_T(c, t));

real mCO = C_YI(c, t, 0);
real mH2 = C_YI(c, t, 1);
real mH2O = C_YI(c, t, 2);
real mCH4 = C_YI(c, t, 3);
real mCO2 = C_YI(c, t, 4);

real pOp = 101325;
real sum = (mCH4 / 16 + mCO2 / 44 + mCO / 28 + mH2 / 2 + mH2O / 18);
real PCH4 = mCH4 / 16 / sum * pOp;
real PCO2 = mCO2 / 44 / sum * pOp;
real PCO = mCO / 28 / sum * pOp;
real PH2 = mH2 / 2 / sum * pOp;
real PH2O = mH2O / 18 / sum * pOp;

if (!strcmp(r->name, "reaction-1"

)
{
*rate = k1 * Kco2_1 * Kch4_1 * (PCH4 / 101325) * (PCO2 / 101325) / pow((1 + Kco2_1 * (PCO2 / 101325) + Kch4_1 * (PCH4 / 101325)), 2) * (1 - pow(((PCO / 101325) * (PH2 / 101325)), 2) / (Kp1 * (PCH4 / 101325) * (PCO2 / 101325)));
C_UDMI(c, t, 0) = *rate;
}

else if (!strcmp(r->name, "reaction-2"

)
{
*rate = k2 * Kco2_2 * Kh2_2 * (PCO2 / 101325) * (PH2 / 101325) / pow((1 + Kco2_2 * (PCO2 / 101325) + Kh2_2 * (PH2 / 101325)), 2) * (1 - (PCO / 101325) * (PH2O / 101325) / (Kp2 * (PCO2 / 101325) * (PH2 / 101325)));
C_UDMI(c, t, 1) = *rate;
}

else if (!strcmp(r->name, "reaction-3"

)
{
*rate = k3 * Kch4_3 * ((PCH4 / 101325) - (pow((PH2 / 101325), 2) / Kp3)) / pow((1 + Kch4_3 * (PCH4 / 101325) + 1 / Kh2_3 * pow((PH2 / 101325), 1.5)), 2);
C_UDMI(c, t, 2) = *rate;
}

else if (!strcmp(r->name, "reaction-4"

)
{
*rate = k4 / Kh2o_4 * ((PH2O / 101325) / (PH2 / 101325) - (PCO / 101325) / Kp4) / pow((1 + Kch4_4 * (PCH4 / 101325) + ((PH2O / 101325) / (PH2 / 101325)) / Kh2o_4 + pow((PH2 / 101325), 1.5) / Kh2_4), 2);
C_UDMI(c, t, 3) = *rate;
}

else if (!strcmp(r->name, "reaction-5"

)
{
*rate = k5 / Kco_5 / Kco2_5 * ((PCO2 / 101325) / (PCO / 101325) - (PCO / 101325) / Kp5) / pow((1 + Kco_5 * (PCO / 101325) + ((PCO2 / 101325) / (PCO / 101325)) / (Kco_5 * Kco2_5)), 2);
C_UDMI(c, t, 4) = *rate;
}
*rr_t = *rate;
}@wuming524

回復(fù)此樓

» 猜你喜歡

085600材料與化工 292分求調(diào)劑已經(jīng)有3人回復(fù)
接受26屆調(diào)劑生已經(jīng)有8人回復(fù)
273求調(diào)劑已經(jīng)有7人回復(fù)
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向已經(jīng)有45人回復(fù)
334求調(diào)劑已經(jīng)有3人回復(fù)
求調(diào)劑推薦已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)
一志愿211，0860總分286食品工程專業(yè)求調(diào)劑已經(jīng)有6人回復(fù)
招調(diào)劑已經(jīng)有14人回復(fù)
誠(chéng)求調(diào)劑，323分有專利/科研/實(shí)習(xí)經(jīng)歷已經(jīng)有6人回復(fù)

1樓 2024-05-01 15:48:51

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 1173018735 的主題更新

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎(jiǎng)勵(lì)，但不允許灌水，必須填寫15個(gè)字符以上)

普通表情龍兔虎貓

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

[考研] 085600材料與化工 292分求調(diào)劑 +3	程晴之 2026-03-06	3/150	2026-03-06 20:25 by 吃骨頭的歡歡丫
[考研] 085701環(huán)境工程295求調(diào)劑，一志愿A區(qū)211 +6	等春來， 2026-03-06	7/350	2026-03-06 18:29 by wang_dand
[考研] 289求調(diào)劑 +3	馳光絕景 2026-03-04	4/200	2026-03-06 12:34 by 馳光絕景
[考研] 材料調(diào)劑 +7	ounce. 2026-03-04	8/400	2026-03-06 10:27 by pies112
[考研] 材料專碩調(diào)劑 +4	慕辰123 2026-03-05	7/350	2026-03-06 10:22 by 慕辰123
[考研] 274求調(diào)劑 +9	一個(gè)學(xué)習(xí)者 2026-03-04	9/450	2026-03-05 20:56 by 黑衣饅頭人
[考研] 一志愿天津大學(xué)085600 319分材料與化工金屬方向求調(diào)劑 +6	青科11 2026-03-02	6/300	2026-03-05 10:05 by oxidpl
[考研] 267調(diào)劑求助 +5	聰少OZ 2026-03-04	5/250	2026-03-05 09:38 by kakakapanpan
[考研] 求調(diào)劑院校 +6	云朵452 2026-03-02	12/600	2026-03-04 23:17 by 云朵452
[考研] 320材料與化工，求調(diào)劑 +6	鶴遨予卿 2026-03-04	8/400	2026-03-04 20:47 by wutongshun
[考研] 292求調(diào)劑 +9	yhk_819 2026-02-28	9/450	2026-03-04 16:06 by sslc1985
[考博] 26申博 +5	north， 2026-02-28	5/250	2026-03-04 10:01 by StarAura
[考研] 0703化學(xué)學(xué)碩 288分 +4	旅行中的紫葡萄 2026-03-03	4/200	2026-03-03 23:41 by ch1979
[考研] 化工專碩調(diào)劑 +4	利好利好. 2026-03-03	7/350	2026-03-03 21:30 by L135790
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5	Qwertyuop 2026-03-01	5/250	2026-03-03 17:50 by 梅子鴨鴨
[考研] 307求調(diào)劑 +6	wyyyqx 2026-03-01	6/300	2026-03-03 09:24 by 2235787770
[考研] 306分材料調(diào)劑 +5	chuanzhu川燭 2026-03-01	6/300	2026-03-02 19:51 by 張曉芳0105
[考研] 295求調(diào)劑 +8	19171856320 2026-02-28	8/400	2026-03-02 11:19 by yuchj
[考研] 一志愿鄭大材料學(xué)碩298分，求調(diào)劑 +6	wsl111 2026-03-01	6/300	2026-03-02 11:00 by ydudjddnd
[考研] 295復(fù)試調(diào)劑 +3	簡(jiǎn)木ChuFront 2026-03-01	3/150	2026-03-01 14:27 by zzxw520th

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

1173018735

[求助] fluent模擬流化床甲烷干重整反應(yīng)udf使用浮點(diǎn)溢出

» 猜你喜歡

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频