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

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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

enigmafwj

捐助貴賓 (小有名氣)

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

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

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

DEFINE_SOURCE(liq_m,c,l_t,dS,eqn) /*液相質(zhì)量輸運(yùn)源項(xiàng)*/
{
Thread *m_t,*v_t;                /*計(jì)算區(qū)域指針*/
  m_t=THREAD_SUPER_THREAD(l_t);    /*混合區(qū)液相指針*/
  v_t=THREAD_SUB_THREAD(m_t,1);    /*單相區(qū)氣相指針*/
real m_dot_l;                      /*液相質(zhì)量輸運(yùn)量*/
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ì)量輸運(yùn)源項(xiàng)*/
{
Thread *m_t,*l_t;                /*計(jì)算區(qū)域指針*/
  m_t=THREAD_SUPER_THREAD(v_t);    /*混合區(qū)氣相指針*/
  l_t=THREAD_SUB_THREAD(m_t,0);    /*單相區(qū)液相指針*/
real m_dot_v;                      /*氣相質(zhì)量輸運(yùn)量*/
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) /*混合相能量方程輸運(yùn)源項(xiàng)*/
{
Thread *l_t,*v_t;                /*計(jì)算區(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;                      /*能量方程源項(xiàng)*/
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)          /*液相汽化，溫度升高，源項(xiàng)為正*/
{
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                               /*氣相液化，溫度降低，源項(xiàng)為負(fù)*/
{
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的回帖

kylafree

至尊木蟲 (知名作家)

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

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

引用回帖:

3樓: Originally posted by enigmafwj at 2012-07-31 15:49:14
因?yàn)槟闼鶈?wèn)的這兩句，是在氣相質(zhì)量方程源項(xiàng)里的，DEFINE_SOURCE(vap_m,c,v_t,dS,eqn)里的thread：v_t在該宏應(yīng)用于氣相質(zhì)量方程時(shí)指向了氣相區(qū)。因?yàn)橹笠玫交旌蠀^(qū)氣相的參數(shù)，所以用THREAD_SUPER_THREAD(v_t)把 ...

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

贊一下

回復(fù)此樓

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

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

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

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17507
散金: 883
紅花: 18
帖子: 6216
在線: 1161.7小時(shí)
蟲號(hào): 1061079
注冊(cè): 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));
定義變量要在賦值語(yǔ)句開始之前全部定義完畢。

另外向兄臺(tái)請(qǐng)教一下，這兩句的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的回帖

enigmafwj

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 0.4
散金: 21
帖子: 106
在線: 35小時(shí)
蟲號(hào): 1917147
注冊(cè): 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; /*液氧飽和蒸汽壓(考慮湍流影 ...

因?yàn)槟闼鶈?wèn)的這兩句，是在氣相質(zhì)量方程源項(xiàng)里的，DEFINE_SOURCE(vap_m,c,v_t,dS,eqn)里的thread：v_t在該宏應(yīng)用于氣相質(zhì)量方程時(shí)指向了氣相區(qū)。因?yàn)橹笠玫交旌蠀^(qū)氣相的參數(shù)，所以用THREAD_SUPER_THREAD(v_t)把指針恢復(fù)到混合區(qū)氣相并傳給m_t；
THREAD_SUB_THREAD(m_t,0)的意思是在m_t得到混合區(qū)氣相指針后，m_t作為最高區(qū)域(混合區(qū))的指針，可以用這個(gè)函數(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的回帖

enigmafwj

捐助貴賓 (小有名氣)

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

引用回帖:

4樓: Originally posted by kylafree at 2012-07-31 16:22:41
感謝兄臺(tái)，受教了。我本來(lái)沒(méi)做多相流，但對(duì)這個(gè)有點(diǎn)興趣。
再問(wèn)幾個(gè)小白的問(wèn)題可以不。比如我按下面這樣理解對(duì)頭不：
氣相區(qū)----就是指這個(gè)區(qū)域的CELL只有氣相；
液相區(qū)----就是指這個(gè)區(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的回帖

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

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

[考研] 279求調(diào)劑 +3	學(xué)而思兮知 2026-04-01	3/150	2026-04-01 23:28 by 欣喜777
[考研] 265求調(diào)劑 +4	林深溫瀾 2026-04-01	6/300	2026-04-01 22:30 by 林深溫瀾
[考研] 303分 0807學(xué)碩求調(diào)劑 +3	TYC3632 2026-04-01	3/150	2026-04-01 19:24 by lwk2004
[考研] 362求調(diào)劑 +13	西南交材料專碩3 2026-03-31	13/650	2026-04-01 17:38 by JYD2011
[考研] 08工科275分求調(diào)劑 +12	AaAa7420 2026-03-31	12/600	2026-04-01 15:45 by jp9609
[考研] 材料專碩306英一數(shù)二 +7	z1z2z3879 2026-03-31	7/350	2026-04-01 14:50 by ZXlzxl0425
[考研] 311求調(diào)劑 +10	李芷新1 2026-03-31	10/500	2026-04-01 14:38 by chenqifeng666
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語(yǔ)已過(guò)四六級(jí)，六級(jí)521，求調(diào)劑 +10	美味可樂(lè)雞翅 2026-03-26	12/600	2026-04-01 11:17 by syh9288
[考研] 環(huán)境工程 085701，267求調(diào)劑 +15	minht 2026-03-29	16/800	2026-04-01 10:13 by li_sujuan99
[考研] 土木304求調(diào)劑 +3	兔突突突， 2026-03-31	3/150	2026-04-01 09:42 by JourneyLucky
[考研] 土木304求調(diào)劑 +5	兔突突突， 2026-03-31	6/300	2026-04-01 09:37 by JourneyLucky
[考研] 085602 307分求調(diào)劑 +10	不知道叫什么！ 2026-03-26	10/500	2026-03-31 19:53 by Dyhoer
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +6	小小機(jī)靈蟲 2026-03-29	12/600	2026-03-31 16:58 by 記事本2026
[考研] 277跪求調(diào)劑 +8	1915668 2026-03-27	13/650	2026-03-31 14:58 by 王亮_大連醫(yī)科大
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3	080700調(diào)劑 2026-03-30	4/200	2026-03-31 12:00 by KLMY666
[考研] 262求調(diào)劑 +7	ZZ..000 2026-03-30	8/400	2026-03-31 10:05 by cal0306
[考研] 083000學(xué)碩274求調(diào)劑 +12	Li李魚 2026-03-26	12/600	2026-03-31 10:01 by cal0306
[考研] 281求調(diào)劑 +5	亞克西good 2026-03-26	7/350	2026-03-30 20:42 by dophin1985
[考研] 322求調(diào)劑 +10	宋明欣 2026-03-27	10/500	2026-03-30 18:47 by 544594351
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�