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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

enigmafwj

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 0.4
散金: 21
帖子: 106
在線: 35小時
蟲號: 1917147
注冊: 2012-07-30
專業(yè): 傳熱傳質(zhì)學(xué)

[求助] 【求助】FLUENT UDF 氣液兩相流汽蝕汽化潛熱(問題詳見內(nèi)容)

為液相、氣相質(zhì)量方程和能量方程添加源項，UDF程序如下，汽化、液化由壓力控制（壓力條件判斷），各物性參數(shù)用溫度擬合。問題是：編譯時，提示有語法錯誤(parse)，并且fv為未聲明的變量。請問如何解決？另：程序中有其他不妥之處，煩請指教，謝謝！
#include "udf.h"
#include "sg_mphase.h"

DEFINE_SOURCE(liq_m,c,l_t,dS,eqn) /*液相質(zhì)量輸運源項*/
{
Thread *m_t,*v_t;                /*計算區(qū)域指針*/
  m_t=THREAD_SUPER_THREAD(l_t);    /*混合區(qū)液相指針*/
  v_t=THREAD_SUB_THREAD(m_t,1);    /*單相區(qū)氣相指針*/
real m_dot_l;                      /*液相質(zhì)量輸運量*/
real stc;                         /*液氧表面張力系數(shù)*/
real fv;                         /*氣相質(zhì)量分?jǐn)?shù)*/
  fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧飽和蒸汽壓(考慮湍流影響)*/
  if(C_T(c,m_t)>154.58)
{
   stc=0;
   pres_sat=5.041e6;
}
  else if(C_T(c,m_t)>54.36)
{
   stc=0.038066*pow(1-C_T(c,l_t)/154.58,1.2136);
   pres_sat=133.322*pow(10,20.6695-526.97/C_T(c,l_t)-6.7062*log10(C_T(c,l_t))+0.012926*C_T(c,l_t)-C_T(c,l_t)*C_T(c,l_t)*(9.8832e-13));
}
  else
{
   stc=0.0225;
   pres_sat=146.041;
}
  pres_turb=0.39*C_R(c,l_t)*C_K(c,m_t);
  pres_vap=pres_sat+0.5*pres_turb;
if(C_P(c,m_t)<=pres_vap)          /*液相汽化*/
{
m_dot_l=-0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)))*(1-fv);
dS[eqn]=-0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)));
}
else                               /*氣相液化*/
{
m_dot_l=0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)))*fv;
dS[eqn]=0;
}
return m_dot_l;
}

DEFINE_SOURCE(vap_m,c,v_t,dS,eqn) /*氣相質(zhì)量輸運源項*/
{
Thread *m_t,*l_t;                /*計算區(qū)域指針*/
  m_t=THREAD_SUPER_THREAD(v_t);    /*混合區(qū)氣相指針*/
  l_t=THREAD_SUB_THREAD(m_t,0);    /*單相區(qū)液相指針*/
real m_dot_v;                      /*氣相質(zhì)量輸運量*/
real stc;                         /*液氧表面張力系數(shù)*/
real fv;                         /*氣相質(zhì)量分?jǐn)?shù)*/
  fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧飽和蒸汽壓(考慮湍流影響)*/
  if(C_T(c,m_t)>154.58)
{
   stc=0;
   pres_sat=5.041e6;
}
  else if(C_T(c,m_t)>54.36)
{
   stc=0.038066*pow(1-C_T(c,l_t)/154.58,1.2136);
   pres_sat=133.322*pow(10,20.6695-526.97/C_T(c,l_t)-6.7062*log10(C_T(c,l_t))+0.012926*C_T(c,l_t)-C_T(c,l_t)*C_T(c,l_t)*(9.8832e-13));
}
  else
{
   stc=0.0225;
   pres_sat=146.041;
}
  pres_turb=0.39*C_R(c,l_t)*C_K(c,m_t);
  pres_vap=pres_sat+0.5*pres_turb;
if(C_P(c,m_t)<=pres_vap)          /*液相汽化*/
{
m_dot_v=0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)))*(1-fv);
dS[eqn]=0;
}
else                               /*氣相液化*/
{
m_dot_l=-0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)))*fv;
dS[eqn]=-0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)));
}
return m_dot_v;
}

DEFINE_SOURCE(m_enrg,c,m_t,dS,eqn) /*混合相能量方程輸運源項*/
{
Thread *l_t,*v_t;                /*計算區(qū)域指針*/
  l_t=THREAD_SUB_THREAD(m_t,0);    /*混合區(qū)液相指針*/
  v_t=THREAD_SUB_THREAD(m_t,1);    /*混合區(qū)氣相指針*/
real m_dot;                      /*相變質(zhì)量*/
real hc;                         /*汽化潛熱*/
real source;                      /*能量方程源項*/
real stc;                         /*液氧表面張力系數(shù)*/
real fv;                         /*氣相質(zhì)量分?jǐn)?shù)*/
  fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧飽和蒸汽壓(考慮湍流影響)*/
  if(C_T(c,m_t)>154.58)
{
   stc=0;
   pres_sat=5.041e6;
   hc=0;
}
  else if(C_T(c,m_t)>54.36)
{
   stc=0.038066*pow(1-C_T(c,l_t)/154.58,1.2136);
   pres_sat=133.322*pow(10,20.6695-526.97/C_T(c,l_t)-6.7062*log10(C_T(c,l_t))+0.012926*C_T(c,l_t)-C_T(c,l_t)*C_T(c,l_t)*(9.8832e-13));
   hc=8040*pow(1-C_T(c,l_t)/154.58,0.201);
}
  else
{
   stc=0.0225;
   pres_sat=146.041;
   hc=2523;
}
  pres_turb=0.39*C_R(c,l_t)*C_K(c,m_t);
  pres_vap=pres_sat+0.5*pres_turb;
if(C_P(c,m_t)<=pres_vap)          /*液相汽化，溫度升高，源項為正*/
{
m_dot=0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)))*(1-fv);
source=m_dot*hc;
dS[eqn]=0;
}
else                               /*氣相液化，溫度降低，源項為負*/
{
m_dot=-0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)))*fv;
source=-m_dot*hc;
dS[eqn]=0;
}
return source;
}

回復(fù)此樓

1樓 2012-07-30 17:51:36

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

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

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17457
散金: 883
紅花: 18
帖子: 6206
在線: 1161.7小時
蟲號: 1061079
注冊: 2010-07-20
性別: GG
專業(yè): 流體力學(xué)

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

★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與！ 2012-07-30 21:50:07
enigmafwj: 金幣+2 2012-12-18 20:11:11

real fv;                         /*氣相質(zhì)量分?jǐn)?shù)*/
fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧飽和蒸汽壓(考慮湍流影響)*/

改為：
real pres_vap,pres_sat,pres_turb; /*液氧飽和蒸汽壓(考慮湍流影響)*/
real fv;                         /*氣相質(zhì)量分?jǐn)?shù)*/
fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
定義變量要在賦值語句開始之前全部定義完畢。

另外向兄臺請教一下，這兩句的thread具體指什么thread，感激不盡啊。
m_t=THREAD_SUPER_THREAD(v_t);    /*混合區(qū)氣相指針*/
l_t=THREAD_SUB_THREAD(m_t,0);    /*單相區(qū)液相指針*/
為啥第一句是“混合區(qū)氣相指針”，而第二句跑到“單相區(qū)液相指針”了呢

贊一下(1人)

回復(fù)此樓

2樓2012-07-30 21:19:25

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

愛問為什么

鐵蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 921
帖子: 62
在線: 54.8小時
蟲號: 2028631
注冊: 2012-09-25
專業(yè): 傳遞過程

引用回帖:

12樓: Originally posted by enigmafwj at 2012-08-11 10:37:58
朋友，這個UDF在單機多核下最后會出現(xiàn)這樣一個問題,能否麻煩你看看，懇請指教。
received suspend command for a pmi context that doesn't exist:unmatched id=2
unable to read the cmd header on the pmi con ...

之前遇到過這個問題，不過好久以前了忘了，好像是C_UDMI語句有錯，你試著改改這句

贊一下(1人)

回復(fù)此樓

16樓2014-01-07 15:37:26

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

普通回帖

enigmafwj

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 0.4
散金: 21
帖子: 106
在線: 35小時
蟲號: 1917147
注冊: 2012-07-30
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

2樓: Originally posted by kylafree at 2012-07-30 21:19:25
real fv; /*氣相質(zhì)量分?jǐn)?shù)*/
fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));
real pres_vap,pres_sat,pres_turb; /*液氧飽和蒸汽壓(考慮湍流影 ...

因為你所問的這兩句，是在氣相質(zhì)量方程源項里的，DEFINE_SOURCE(vap_m,c,v_t,dS,eqn)里的thread：v_t在該宏應(yīng)用于氣相質(zhì)量方程時指向了氣相區(qū)。因為之后要用到混合區(qū)氣相的參數(shù)，所以用THREAD_SUPER_THREAD(v_t)把指針恢復(fù)到混合區(qū)氣相并傳給m_t；
THREAD_SUB_THREAD(m_t,0)的意思是在m_t得到混合區(qū)氣相指針后，m_t作為最高區(qū)域(混合區(qū))的指針，可以用這個函數(shù)把單相區(qū)液相指針傳給l_t，才能用l_t調(diào)用到液相區(qū)的參數(shù)，0指的是主相(液相)，1指的是次相(氣相)。

贊一下

回復(fù)此樓

3樓2012-07-31 15:49:14

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

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17457
散金: 883
紅花: 18
帖子: 6206
在線: 1161.7小時
蟲號: 1061079
注冊: 2010-07-20
性別: GG
專業(yè): 流體力學(xué)

★ ★ ★
xiegangmai: 金幣+3, 鼓勵討論交流！ 2012-08-02 23:05:49

引用回帖:

3樓: Originally posted by enigmafwj at 2012-07-31 15:49:14
因為你所問的這兩句，是在氣相質(zhì)量方程源項里的，DEFINE_SOURCE(vap_m,c,v_t,dS,eqn)里的thread：v_t在該宏應(yīng)用于氣相質(zhì)量方程時指向了氣相區(qū)。因為之后要用到混合區(qū)氣相的參數(shù)，所以用THREAD_SUPER_THREAD(v_t)把 ...

感謝兄臺，受教了。我本來沒做多相流，但對這個有點興趣。
再問幾個小白的問題可以不。比如我按下面這樣理解對頭不：
氣相區(qū)----就是指這個區(qū)域的CELL只有氣相；
液相區(qū)----就是指這個區(qū)域的CELL只有液相；
混合區(qū)----就是指這個區(qū)域的CELL既有氣相也有液相；
混全區(qū)thread----是說混合區(qū)既然有混合thread，又有氣相thread，還有液相thead。氣相thread液相thread混合thread怎么聯(lián)系起來啊，糊涂了

贊一下

回復(fù)此樓

4樓2012-07-31 16:22:41

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

enigmafwj

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 0.4
散金: 21
帖子: 106
在線: 35小時
蟲號: 1917147
注冊: 2012-07-30
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

4樓: Originally posted by kylafree at 2012-07-31 16:22:41
感謝兄臺，受教了。我本來沒做多相流，但對這個有點興趣。
再問幾個小白的問題可以不。比如我按下面這樣理解對頭不：
氣相區(qū)----就是指這個區(qū)域的CELL只有氣相；
液相區(qū)----就是指這個區(qū)域的CELL只有液相；
混 ...

你看下UDF幫助的第3章第11節(jié)的區(qū)域結(jié)構(gòu)和第6章THREAD函數(shù)。

贊一下

回復(fù)此樓

5樓2012-08-02 09:54:37

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

enigmafwj

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 0.4
散金: 21
帖子: 106
在線: 35小時
蟲號: 1917147
注冊: 2012-07-30
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

照你的改完還是不行啊

贊一下

回復(fù)此樓

6樓2012-08-02 09:55:06

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

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17457
散金: 883
紅花: 18
帖子: 6206
在線: 1161.7小時
蟲號: 1061079
注冊: 2010-07-20
性別: GG
專業(yè): 流體力學(xué)

引用回帖:

6樓: Originally posted by enigmafwj at 2012-08-02 09:55:06
照你的改完還是不行啊...

你干脆把出錯代碼貼出來大家都可以看到

贊一下

回復(fù)此樓

7樓2012-08-02 16:09:10

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

enigmafwj

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 0.4
散金: 21
帖子: 106
在線: 35小時
蟲號: 1917147
注冊: 2012-07-30
專業(yè): 傳熱傳質(zhì)學(xué)

★
xiegangmai: 金幣+1, 鼓勵分享解決經(jīng)驗，金幣獎勵！ 2012-08-11 22:33:37

引用回帖:

7樓: Originally posted by kylafree at 2012-08-02 16:09:10
你干脆把出錯代碼貼出來大家都可以看到...

改好了。把變量聲明全部放前面就OK了

贊一下

回復(fù)此樓

8樓2012-08-02 21:54:26

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

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17457
散金: 883
紅花: 18
帖子: 6206
在線: 1161.7小時
蟲號: 1061079
注冊: 2010-07-20
性別: GG
專業(yè): 流體力學(xué)

引用回帖:

8樓: Originally posted by enigmafwj at 2012-08-02 21:54:26
改好了。把變量聲明全部放前面就OK了...

那恭喜啊，有時間給我講解下我的疑問��？

贊一下

回復(fù)此樓

9樓2012-08-02 22:00:06

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

enigmafwj

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 0.4
散金: 21
帖子: 106
在線: 35小時
蟲號: 1917147
注冊: 2012-07-30
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

9樓: Originally posted by kylafree at 2012-08-02 22:00:06
那恭喜啊，有時間給我講解下我的疑問��？

...

我剛學(xué)這些個東西，只是按自己的理解用它，你覺得這兩函數(shù)是否是這樣用的？

贊一下

回復(fù)此樓

10樓2012-08-03 16:01:53

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

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

返回列表

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

[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 271材料工程求調(diào)劑 +4	.6lL 2026-03-18	4/200	2026-03-18 20:05 by 楤哥
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 0703化學(xué)調(diào)劑 +4	pupcoco 2026-03-17	7/350	2026-03-18 12:14 by djl2006
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao

24小時熱門版塊排行榜

[求助] 【求助】FLUENT UDF 氣液兩相流 汽蝕 汽化潛熱(問題詳見內(nèi)容)

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

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

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

[求助] 【求助】FLUENT UDF 氣液兩相流汽蝕汽化潛熱(問題詳見內(nèi)容)

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