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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

837729744

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 455.7
散金: 181
紅花: 1
帖子: 152
在線: 117.6小時
蟲號: 1806845
注冊: 2012-05-09
專業(yè): 可再生與替代能源利用中的

[求助] 幫忙看下這個udf有問題嗎？初始化沒問題，點計算出錯

幫忙看下這個udf有問題嗎？初始化沒問題，點計算出錯

#include "udf.h"                                              /*包括常規(guī)宏*/
#include "sg_mphase.h"                                     /*包括體積分數(shù)宏*/
#define  LAT_HT 2260000                                        /*定義水蒸汽潛熱*/
#define  pi 3.1415                                           /*圓周率常數(shù)*/
#define  R 8314                                           /*gas constant 8.314j/mol.K*/
#define  moleweight 0.018                                  /*SI unit water mole weight 0.018kg/mol */
#define A 7.19621
#define B 1730.63
#define C 233.426
DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)                /*定義液相源項*/
{
  Thread *mix_th, *sec_th;                                  /*混合相、第二相定義計算區(qū)域指針*/
  real m_dot_l;                                              /*蒸發(fā)冷凝速率  定義液相質(zhì)量轉(zhuǎn)移*/
  real P_SAT ;
  mix_th = THREAD_SUPER_THREAD(pri_th);                      /*指向混合區(qū)的主相即液相的指針*/
  sec_th = THREAD_SUB_THREAD(mix_th, 1);                   /* 指向單相控制區(qū)的氣相的指針*/
  P_SAT=exp(A-B/(C_T(cell,pri_th)+C-273.15))*1e3;    /*定義飽和蒸汽壓與液相溫度的關(guān)系式，即antoine公式*/
if(C_P(cell, sec_th)< P_SAT)
{m_dot_l = -(P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); /*定義液相向氣相轉(zhuǎn)移的速率，如果氣相的壓力小于一定溫度下的飽和蒸汽壓，液相向氣相發(fā)生質(zhì)量轉(zhuǎn)移*/
dS[eqn]= -2*pi*R*sqrt(moleweight)*(P_SAT-P_SAT*C_T(cell,pri_th)*B/pow((C_T(cell,pri_th)+C-273.15),2)-C_P(cell, sec_th));                /*定義源項對質(zhì)量轉(zhuǎn)移偏導(dǎo)，fabs為取絕對值*/
}
else
{m_dot_l = (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); /*定義液相向氣相轉(zhuǎn)移的速率，如果氣相壓力大于一定溫度下的飽和蒸汽壓，氣相向液相發(fā)生質(zhì)量轉(zhuǎn)移*/
dS[eqn] = 0;                      /*由于氣相向液相轉(zhuǎn)移，所以液相的質(zhì)量源項對質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零*/
}
return m_dot_l;
}
DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn)
{
  Thread *mix_th, *pri_th;
  real m_dot_v;
  real P_SAT ;
  mix_th = THREAD_SUPER_THREAD(sec_th);
  pri_th = THREAD_SUB_THREAD(mix_th, 0);
  P_SAT=exp(A-B/(C_T(cell,pri_th)+C-273.15))*1e3;
if(C_P(cell, sec_th)< P_SAT)
{
  m_dot_v = (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th));
  dS[eqn] =0;
}
else
  {m_dot_v = -(P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th));
  dS[eqn]=-2*pi*R*sqrt(moleweight)*(P_SAT-P_SAT*C_T(cell,pri_th)*B/pow((C_T(cell,pri_th)+C-273.15),2)-C_P(cell, sec_th));
  }
return m_dot_v;
}
DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
{
Thread *pri_th,*sec_th;
real m_dot;
real P_SAT;
pri_th=THREAD_SUB_THREAD(mix_th,0);
sec_th=THREAD_SUB_THREAD(mix_th,1);
P_SAT=exp(A-B/(C_T(cell,pri_th)+C-273.15))*1e3;
if(C_P(cell, sec_th)< P_SAT)
{
  m_dot =- (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th));
  dS[eqn]=-2*pi*R*sqrt(moleweight)*(P_SAT-P_SAT*C_T(cell,pri_th)*B/pow((C_T(cell,pri_th)+C-273.15),2)-C_P(cell, sec_th));
}
else
{
m_dot = (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th));
dS[eqn]=2*pi*R*sqrt(moleweight)*(P_SAT-P_SAT*C_T(cell,pri_th)*B/pow((C_T(cell,pri_th)+C-273.15),2)-C_P(cell, sec_th));
}

return LAT_HT*m_dot;

}

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有7人回復(fù)
275求調(diào)劑已經(jīng)有6人回復(fù)
招08考數(shù)學已經(jīng)有7人回復(fù)
303求調(diào)劑已經(jīng)有6人回復(fù)
306求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
考研調(diào)劑已經(jīng)有4人回復(fù)
285求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學070300 學碩336求調(diào)劑已經(jīng)有7人回復(fù)
318求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

FLUENT中編寫UDF時指針變量是否要初始化已經(jīng)有4人回復(fù)
fluent14編譯UDF時總是出錯，求高手啊�。�！難道沒人遇到過嗎？？？已經(jīng)有11人回復(fù)
求助：模型初始化時出錯已經(jīng)有5人回復(fù)
fluent udf請教已經(jīng)有16人回復(fù)
一個UDF并行報錯已經(jīng)有17人回復(fù)
導(dǎo)入UDF時出這種錯誤！！已經(jīng)有4人回復(fù)
UDF 改為并行計算已經(jīng)有9人回復(fù)
FLUENT UDF 編譯通過，初始化出錯的原因？？？求助已經(jīng)有4人回復(fù)
求助�。�！我想要在fluent中針對一個面進行積分，那么UDF應(yīng)該用什么宏函數(shù)？怎么用？已經(jīng)有5人回復(fù)
好心人幫幫忙吧，為什么這個UDF導(dǎo)入fluent以后，在邊界條件處沒有顯示~~~ 已經(jīng)有17人回復(fù)
怎樣解決FLUENT received fatal signal (ACCESS_VIOLATION)的錯誤？已經(jīng)有20人回復(fù)
udf編譯出錯！急已經(jīng)有10人回復(fù)
udf編程通過，但是在FLUENT中運行出錯，麻煩udf高手幫忙看看已經(jīng)有4人回復(fù)
使用UDF后，初始化錯誤已經(jīng)有11人回復(fù)
【求助】fluent初始化的時候出現(xiàn)問題，請大家指點，在線等，多謝各位了已經(jīng)有4人回復(fù)
【求助】fluent計算中，udf編譯通過，初始化出錯已經(jīng)有11人回復(fù)
【求助】udf出錯已經(jīng)有13人回復(fù)

1樓 2013-09-05 10:38:31

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

回帖支持 ( 顯示支持度最高的前 50 名 )

轉(zhuǎn)石千仞

金蟲 (正式寫手)

擊破長空

應(yīng)助: 24 (小學生)
金幣: 941
散金: 164
紅花: 13
帖子: 673
在線: 161.6小時
蟲號: 1163587
注冊: 2010-12-06
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

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

★ ★ ★ ★ ★ ★ ★ ★
837729744: 金幣+8, ★★★很有幫助 2013-09-06 08:43:19

引用回帖:

5樓: Originally posted by 837729744 at 2013-09-05 11:03:37
操作壓力設(shè)為一個大氣壓，那么初始化時Gauge pressure的壓力應(yīng)該設(shè)為101325pa還是0pa（我想讓我的模型處于一個大氣壓下）。我感覺是不是壓力這塊出問題了...

不是的。
設(shè)置參考壓力后，你初始為0，即絕對壓力為一個大氣壓。

贊一下(1人)

回復(fù)此樓

故善戰(zhàn)人之勢，如轉(zhuǎn)圓石于千仞之山者，勢也

6樓2013-09-06 08:34:10

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

轉(zhuǎn)石千仞

金蟲 (正式寫手)

擊破長空

應(yīng)助: 24 (小學生)
金幣: 941
散金: 164
紅花: 13
帖子: 673
在線: 161.6小時
蟲號: 1163587
注冊: 2010-12-06
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

11樓: Originally posted by 837729744 at 2013-09-07 08:37:34
還有就是都是在什么情況下，才用求導(dǎo)啊（dS=），明白為什么要對上式求導(dǎo)，求導(dǎo)的作用是什么？...

是為了加快收斂。
你應(yīng)該去看相關(guān)理論些的書，你要明白，并不是所有答案，別人都可以告訴你的。你要自己去尋找。

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

837729744

故善戰(zhàn)人之勢，如轉(zhuǎn)圓石于千仞之山者，勢也

12樓2013-09-07 10:16:23

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

普通回帖

837729744

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 455.7
散金: 181
紅花: 1
帖子: 152
在線: 117.6小時
蟲號: 1806845
注冊: 2012-05-09
專業(yè): 可再生與替代能源利用中的

錯誤代碼：Error: Floating point error: invalid number

Error Object: ()
把所有質(zhì)量轉(zhuǎn)移速率和求導(dǎo)的式子改成為=0，（m_dot=0;dS[eqn]=0),可以進行計算，是不是式子的問題��？求指導(dǎo)

贊一下

回復(fù)此樓

2樓2013-09-05 10:39:39

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

轉(zhuǎn)石千仞

金蟲 (正式寫手)

擊破長空

應(yīng)助: 24 (小學生)
金幣: 941
散金: 164
紅花: 13
帖子: 673
在線: 161.6小時
蟲號: 1163587
注冊: 2010-12-06
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
837729744: 金幣+8 2013-09-05 11:04:03
1592203609: 金幣+3, 謝謝回帖 2013-09-05 16:01:15

引用回帖:

2樓: Originally posted by 837729744 at 2013-09-05 10:39:39
錯誤代碼：Error: Floating point error: invalid number

Error Object: ()
把所有質(zhì)量轉(zhuǎn)移速率和求導(dǎo)的式子改成為=0，（m_dot=0;dS=0),可以進行計算，是不是式子的問題��？求指導(dǎo)

這種情況下，建議你對計算的數(shù)據(jù)進行檢查。
你可以采用程序截斷的形式，用Message插入可能出問題的位置前后，看是否能輸出語言與數(shù)據(jù)，這樣可以檢測你的程序運行到哪里出現(xiàn)了浮點錯誤。
你還可以用Error語句來控制程序的停止，你可以讓程序運行到任意位置停止，以便不讓程序在每一個控制容積里運行。
出現(xiàn)浮點錯誤的原因，是你計算出現(xiàn)極大或極小值，比0等等。
調(diào)試程序需要耐心與經(jīng)驗，祝你好運！

贊一下

回復(fù)此樓

故善戰(zhàn)人之勢，如轉(zhuǎn)圓石于千仞之山者，勢也

3樓2013-09-05 10:51:03

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

轉(zhuǎn)石千仞

金蟲 (正式寫手)

擊破長空

應(yīng)助: 24 (小學生)
金幣: 941
散金: 164
紅花: 13
帖子: 673
在線: 161.6小時
蟲號: 1163587
注冊: 2010-12-06
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

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

★ ★
1592203609: 金幣+2, 謝謝回帖 2013-09-05 16:01:29

引用回帖:

你的求導(dǎo)的式子，可能產(chǎn)生了極大值，“變化“太大，你的方程立馬發(fā)散。

贊一下

回復(fù)此樓

故善戰(zhàn)人之勢，如轉(zhuǎn)圓石于千仞之山者，勢也

4樓2013-09-05 10:52:09

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

837729744

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 455.7
散金: 181
紅花: 1
帖子: 152
在線: 117.6小時
蟲號: 1806845
注冊: 2012-05-09
專業(yè): 可再生與替代能源利用中的

引用回帖:

4樓: Originally posted by 轉(zhuǎn)石千仞 at 2013-09-05 10:52:09
你的求導(dǎo)的式子，可能產(chǎn)生了極大值，“變化“太大，你的方程立馬發(fā)散。...

操作壓力設(shè)為一個大氣壓，那么初始化時Gauge pressure的壓力應(yīng)該設(shè)為101325pa還是0pa（我想讓我的模型處于一個大氣壓下）。我感覺是不是壓力這塊出問題了

贊一下

回復(fù)此樓

5樓2013-09-05 11:03:37

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

837729744

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 455.7
散金: 181
紅花: 1
帖子: 152
在線: 117.6小時
蟲號: 1806845
注冊: 2012-05-09
專業(yè): 可再生與替代能源利用中的

引用回帖:

3樓: Originally posted by 轉(zhuǎn)石千仞 at 2013-09-05 10:51:03
這種情況下，建議你對計算的數(shù)據(jù)進行檢查。
你可以采用程序截斷的形式，用Message插入可能出問題的位置前后，看是否能輸出語言與數(shù)據(jù)，這樣可以檢測你的程序運行到哪里出現(xiàn)了浮點錯誤。
你還可以用Error語句來控 ...

Message插入和用Error語句來控制程序的停止，我不會啊，難嗎？怎么弄啊？

贊一下

回復(fù)此樓

7樓2013-09-06 08:56:51

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

轉(zhuǎn)石千仞

金蟲 (正式寫手)

擊破長空

應(yīng)助: 24 (小學生)
金幣: 941
散金: 164
紅花: 13
帖子: 673
在線: 161.6小時
蟲號: 1163587
注冊: 2010-12-06
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

7樓: Originally posted by 837729744 at 2013-09-06 08:56:51
Message插入和用Error語句來控制程序的停止，我不會啊，難嗎？怎么弄��？...

建議你去看幫助文件。里面有，兩條語句而已，不難的。

贊一下

回復(fù)此樓

故善戰(zhàn)人之勢，如轉(zhuǎn)圓石于千仞之山者，勢也

8樓2013-09-06 18:32:36

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

837729744

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 455.7
散金: 181
紅花: 1
帖子: 152
在線: 117.6小時
蟲號: 1806845
注冊: 2012-05-09
專業(yè): 可再生與替代能源利用中的

現(xiàn)在可以計算了，如果轉(zhuǎn)移質(zhì)量速率公式是這個：m_dot = -(P_SAT-C_P(cell,pri_th))*sqrt(moleweight/(2*pi*R*C_T(cell,sec_th)))。不知道對這個公式的質(zhì)量項求導(dǎo)是什么（dS[eqn]=。。。）。原先dS[eqn]公式估計也不對。

贊一下

回復(fù)此樓

9樓2013-09-06 19:10:55

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

轉(zhuǎn)石千仞

金蟲 (正式寫手)

擊破長空

應(yīng)助: 24 (小學生)
金幣: 941
散金: 164
紅花: 13
帖子: 673
在線: 161.6小時
蟲號: 1163587
注冊: 2010-12-06
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

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

★
1592203609: 金幣+1, 謝謝回帖 2013-09-07 21:16:40

引用回帖:

9樓: Originally posted by 837729744 at 2013-09-06 19:10:55
現(xiàn)在可以計算了，如果轉(zhuǎn)移質(zhì)量速率公式是這個：m_dot = -(P_SAT-C_P(cell,pri_th))*sqrt(moleweight/(2*pi*R*C_T(cell,sec_th)))。不知道對這個公式的質(zhì)量項求導(dǎo)是什么（dS=。。。）。原先dS公式估計也不對。...

對變量求導(dǎo)，你求解的是什么方程，方程中將哪個量作為變量，你就對該求量求導(dǎo)。

贊一下

回復(fù)此樓

故善戰(zhàn)人之勢，如轉(zhuǎn)圓石于千仞之山者，勢也

10樓2013-09-06 20:10:00

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

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

返回列表

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

[考研] 306求調(diào)劑 +4	來好運來來來 2026-03-22	4/200	2026-03-22 12:43 by 素顏傾城1988
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	4/200	2026-03-22 10:25 by QGZDSYS
[基金申請] 山東省面上項目限額評審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 【考研調(diào)劑】化學專業(yè) 281分，一志愿四川大學，誠心求調(diào)劑 +11	吃吃吃才有意義 2026-03-19	11/550	2026-03-21 18:23 by 學員8dgXkO
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 070300化學319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 324求調(diào)劑 +5	lucky呀呀呀鴨 2026-03-20	5/250	2026-03-20 22:30 by 促天成
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 081700化工學碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 一志愿福大288有機化學，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy

24小時熱門版塊排行榜

[求助] 幫忙看下這個udf有問題嗎？初始化沒問題，點計算出錯

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

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

» 本帖已獲得的紅花（最新10朵）

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

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

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

[求助] 幫忙看下這個udf有問題嗎？初始化沒問題，點計算出錯