nwj0502

應(yīng)助: 5 (幼兒園)
金幣: 11135.3
散金: 6315
紅花: 38
帖子: 9032
在線: 1320.5小時(shí)
蟲號(hào): 771979
注冊(cè): 2009-05-16
性別: GG
專業(yè): 物理化學(xué)

[求助] fluent模擬堆積床內(nèi)焦炭顆粒燃燒，出現(xiàn)錯(cuò)誤提示“pressure correction”

大家好，我想用fluent模擬堆積床內(nèi)焦炭顆粒（直徑0.005m）燃燒過程傳熱傳質(zhì)，采用氣固雙能量方程，組分方程。
遇到的問題：不加氣固能量源項(xiàng)和氣體組分方程，給定固體初始溫度1000K，氣體初始溫度300K，模擬氣固換熱是沒有問題的；但是加上氣固能量源項(xiàng)后，點(diǎn)開始計(jì)算，給定氣體溫度1473K模擬點(diǎn)火，氣體想固體傳熱，當(dāng)固體溫度達(dá)到著火點(diǎn)923K時(shí)（這時(shí)開始發(fā)生焦炭燃燒反應(yīng)），fluent就會(huì)提示錯(cuò)誤“pressure correction”，是不是能量源項(xiàng)太大？需要對(duì)源項(xiàng)進(jìn)行線性化處理，具體怎么處理呢？還請(qǐng)fluent高手給予指點(diǎn)，小弟不勝感激。

模擬焦炭燃燒過程傳熱傳質(zhì)

焦炭燃燒反應(yīng)：C+O2→CO2
焦炭燃點(diǎn)：923K
反應(yīng)速率：Rc=3.14*（dc）^2*kc *CO2*Nc
               式中，dc是焦炭顆粒直徑，kc是綜合反應(yīng)速率常數(shù)，CO2是氧氣摩爾濃度，Nc是單位體積焦炭顆粒數(shù)；
初始條件和邊界條件：
Tg，0=1473K，Ts=340K，Co2,0= Co2，Co2,0= Cco2，Cn2,0= Cn2
點(diǎn)火溫度1473K(Tg=1473K)，時(shí)間60s；保溫溫度1043K(Tg=1043K)，時(shí)間90s；接下來抽風(fēng)溫度300K(Tg=300K)（壓力出口，-10KPa）

以下是我的部分udf代碼：
real kfc;                            //傳質(zhì)系數(shù),m/s；
real Sh,Re,Sc;                         //舍伍德數(shù)，顆粒雷諾數(shù)，施密特準(zhǔn)數(shù)；
real mu_lam=C_MU_L(c,t);             //氣體黏度；
real den_gas=C_R(c,t);                   //氣體密度，kg/m3；
real Nc=9.92*pow(10.0,5.0);             //單位體積焦炭顆粒數(shù)，1/m3;
real mole_total;                      //氣體總的摩爾數(shù)
real P_total;                         //氣體總壓,pa
real P_O2;                            //氧氣分壓,pa
real DO2;                            //氧氣擴(kuò)散系數(shù),m2/s;
real Co;                               //氧氣摩爾濃度，mol/m3;
real O2_mole_frac;                      //氧氣摩爾分?jǐn)?shù)；
real kc1;                            //綜合反應(yīng)速率常數(shù)，m/s;
real k1;                               //反應(yīng)速率常數(shù),m/s;
real R1;                               //焦炭完全燃燒的反應(yīng)速率，mol/m3s;

P_total=ABS_P(C_P(c,t),op_pres);          //氣體總壓?
mole_total=C_YI(c,t,0)/32.0+C_YI(c,t,1)/44.0+C_YI(c,t,2)/28.0; //
O2_mole_frac=C_YI(c,t,0)/(32.0*mole_total);
P_O2=P_total*O2_mole_frac;             //氧氣分壓
Co=P_O2/R*C_T(c,t);

DO2=2.00*pow(10.0,-4.0)*pow(0.001,1.78)/P_total;
Sc=mu_lam/(den_gas*DO2);             //施密特準(zhǔn)數(shù)舍伍德準(zhǔn)數(shù)
Re=ND_MAG(C_U(c,t),C_V(c,t),C_W(c,t))*DIAMETER*den_gas/mu_lam; //雷諾數(shù)
Sh=2.0/POROUSITY+0.75*pow(Re,0.5)*pow(Sc,0.333)/POROUSITY;    //舍伍德數(shù)
kfc=Sh*DO2/DIAMETER;                                  //傳質(zhì)系數(shù)
k1=6.53*pow(10.0,7.0)*pow(C_UDSI(c,t,0),0.5)*exp(-44000.0/(R*C_UDSI(c,t,0))); //反應(yīng)速率常數(shù)
kc1=(kfc+k1)/kfc*k1;

if (C_UDSI(c,t,0)<923.0)

      R1=0;

else

      R1=3.14*pow(DIAMETER,2.0)*kc1*Co*Nc;

return R1;    //
DEFINE_SOURCE(source_energy_gas,c,t,dS,eqn)
{
real source;
real delta_H=-2.79*pow(10.0,5.0);
real a=0.5;  //反應(yīng)熱分配系數(shù)
source=htc(c,t)*(C_UDSI(c,t,0)-C_T(c,t))-(1-a)*(delta_H*RC(c,t));
dS[eqn]=0;
return source;
}

DEFINE_SOURCE(source_energy_solid,c,t,dS,eqn)
{
real source;
real delta_H=-2.79*pow(10.0,5.0);
real a=0.5;  //反應(yīng)熱分配系數(shù)
source=htc(c,t)*(C_T(c,t)-C_UDSI(c,t,0))-a*(delta_H*RC(c,t));
dS[eqn]=0;
return source;
}

fluent模擬堆積床內(nèi)焦炭顆粒燃燒，出現(xiàn)錯(cuò)誤提示“pressure correction”

壓力修正.jpg

fluent模擬堆積床內(nèi)焦炭顆粒燃燒，出現(xiàn)錯(cuò)誤提示“pressure correction”-1

出現(xiàn)錯(cuò)誤.jpg

回復(fù)此樓

1樓 2014-06-16 00:00:16

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

nwj0502

至尊木蟲 (知名作家)

應(yīng)助: 5 (幼兒園)
金幣: 11135.3
散金: 6315
紅花: 38
帖子: 9032
在線: 1320.5小時(shí)
蟲號(hào): 771979
注冊(cè): 2009-05-16
性別: GG
專業(yè): 物理化學(xué)

這么晚大家都睡了吧，自己頂一下吧，希望好心人來幫助，在這個(gè)問題上糾結(jié)兩個(gè)禮拜了

贊一下

回復(fù)此樓

2樓2014-06-16 00:08:55

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

nwj0502

至尊木蟲 (知名作家)

應(yīng)助: 5 (幼兒園)
金幣: 11135.3
散金: 6315
紅花: 38
帖子: 9032
在線: 1320.5小時(shí)
蟲號(hào): 771979
注冊(cè): 2009-05-16
性別: GG
專業(yè): 物理化學(xué)

怎么沒人回復(fù)呢，熱心的蟲子幫幫唄

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

3樓2014-06-16 08:53:30

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

gearboy78

禁蟲 (小有名氣)

本帖內(nèi)容被屏蔽

4樓2014-06-16 09:20:43

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

nwj0502

至尊木蟲 (知名作家)

應(yīng)助: 5 (幼兒園)
金幣: 11135.3
散金: 6315
紅花: 38
帖子: 9032
在線: 1320.5小時(shí)
蟲號(hào): 771979
注冊(cè): 2009-05-16
性別: GG
專業(yè): 物理化學(xué)

引用回帖:

4樓: Originally posted by gearboy78 at 2014-06-16 09:20:43
具體調(diào)試是很費(fèi)時(shí)間和精力的事情，只能自己去做，不要期望別人能幫你。就好像寫一篇論文一樣，你覺得會(huì)有人幫你寫嗎？

我不是要?jiǎng)e人幫我調(diào)程序，因?yàn)閯倢W(xué)UDF不久，對(duì)UDF理解不深，我是想讓大家看看程序，源項(xiàng)部分處理有什么不對(duì)的地方

贊一下

回復(fù)此樓

5樓2014-06-16 09:47:07

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

gearboy78

禁蟲 (小有名氣)

本帖內(nèi)容被屏蔽

6樓2014-06-16 12:35:32

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

nwj0502

至尊木蟲 (知名作家)

應(yīng)助: 5 (幼兒園)
金幣: 11135.3
散金: 6315
紅花: 38
帖子: 9032
在線: 1320.5小時(shí)
蟲號(hào): 771979
注冊(cè): 2009-05-16
性別: GG
專業(yè): 物理化學(xué)

引用回帖:

6樓: Originally posted by gearboy78 at 2014-06-16 12:35:32
不調(diào)試直接看源碼的話只能看語法是否錯(cuò)誤，你編譯通過了語法應(yīng)該是沒問題的，其它就是系統(tǒng)問題了，不去加載調(diào)試很難看出問題。

嗯，謝謝你，最近自己也一直在調(diào)試，找問題

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

7樓2014-06-16 12:50:06

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

gearboy78

禁蟲 (小有名氣)

本帖內(nèi)容被屏蔽

8樓2014-06-16 13:23:42

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

longma135

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1186.6
散金: 15
帖子: 199
在線: 207小時(shí)
蟲號(hào): 2533271
注冊(cè): 2013-07-05
性別: GG
專業(yè): 能源化工

點(diǎn)火溫度1473K(Tg=1473K)，時(shí)間60s；保溫溫度1043K(Tg=1043K)，時(shí)間90s；接下來抽風(fēng)溫度300K(Tg=300K)（壓力出口，-10KPa）

我想問下上面這部分怎么做到的�。�

贊一下

回復(fù)此樓

9樓2014-06-18 11:01:14

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

nwj0502

至尊木蟲 (知名作家)

應(yīng)助: 5 (幼兒園)
金幣: 11135.3
散金: 6315
紅花: 38
帖子: 9032
在線: 1320.5小時(shí)
蟲號(hào): 771979
注冊(cè): 2009-05-16
性別: GG
專業(yè): 物理化學(xué)

引用回帖:

9樓: Originally posted by longma135 at 2014-06-18 11:01:14
點(diǎn)火溫度1473K(Tg=1473K)，時(shí)間60s；保溫溫度1043K(Tg=1043K)，時(shí)間90s；接下來抽風(fēng)溫度300K(Tg=300K)（壓力出口，-10KPa）

我想問下上面這部分怎么做到的�。�

點(diǎn)火，保溫，抽風(fēng)的邊界條件是通過UDF來完成的，具體幫助里有

贊一下

回復(fù)此樓

10樓2014-06-18 15:27:08

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

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

[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9	Ymlll 2026-03-18	12/600	2026-03-20 00:01 by 23Postgrad
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 本人考085602 化學(xué)工程專碩 +17	不知道叫什么！ 2026-03-15	19/950	2026-03-19 15:06 by 盡舜堯1
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +3	Ncdx123456 2026-03-19	3/150	2026-03-19 13:18 by houyaoxu
[考研] 材料080500調(diào)劑求收留 +4	一顆meteor 2026-03-13	4/200	2026-03-19 10:32 by 30660438
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425