版塊導(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 氣液兩相流汽蝕汽化潛熱(問題詳見內(nèi)容)

為液相、氣相質(zhì)量方程和能量方程添加源項(xiàng)，UDF程序如下，汽化、液化由壓力控制（壓力條件判斷），各物性參數(shù)用溫度擬合。問題是：編譯時(shí)，提示有語法錯(cuò)誤(parse)，并且fv為未聲明的變量。請(qǐng)問如何解決？另：程序中有其他不妥之處，煩請(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ù)此樓

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

多相流UDF

傳熱傳質(zhì)

» 猜你喜歡

材料調(diào)劑已經(jīng)有9人回復(fù)
311（085601）求調(diào)劑已經(jīng)有5人回復(fù)
307分求調(diào)劑已經(jīng)有6人回復(fù)
085601英二數(shù)二求調(diào)劑總分325 已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有10人回復(fù)
生醫(yī)工0831調(diào)劑求推薦已經(jīng)有6人回復(fù)
生物學(xué)296求調(diào)劑已經(jīng)有12人回復(fù)
一志愿武理材料工程302調(diào)劑環(huán)化或化工已經(jīng)有4人回復(fù)
286分調(diào)劑已經(jīng)有12人回復(fù)
復(fù)試調(diào)劑已經(jīng)有3人回復(fù)

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

求助�。�！我想要在fluent中針對(duì)一個(gè)面進(jìn)行積分，那么UDF應(yīng)該用什么宏函數(shù)？怎么用？已經(jīng)有5人回復(fù)
FLUENT求幫助啊啊啊啊啊已經(jīng)有8人回復(fù)
好心人幫幫忙吧，為什么這個(gè)UDF導(dǎo)入fluent以后，在邊界條件處沒有顯示~~~ 已經(jīng)有17人回復(fù)
fluent中用udf編寫三維速度邊界條件出現(xiàn)錯(cuò)誤，大家來指導(dǎo)一下啊已經(jīng)有9人回復(fù)
FLUENT源項(xiàng)UDF非穩(wěn)態(tài)迭代的一個(gè)疑問請(qǐng)教已經(jīng)有8人回復(fù)
有關(guān)fluent的udf中的hook 已經(jīng)有6人回復(fù)
【求助】氣液分離兩相流動(dòng)模擬fluent gambit出入口的設(shè)置已經(jīng)有11人回復(fù)
【求助】單機(jī)多核fluent并行udf 已經(jīng)有13人回復(fù)
【求助】fluent計(jì)算中，udf編譯通過，初始化出錯(cuò) 已經(jīng)有11人回復(fù)
【求助】如何改變Fluent中的曳力模型（UDF）已經(jīng)有12人回復(fù)

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

已閱回復(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)槟闼鶈柕倪@兩句，是在氣相質(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的回帖

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

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17502
散金: 883
紅花: 18
帖子: 6215
在線: 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));
定義變量要在賦值語句開始之前全部定義完畢。

另外向兄臺(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的回帖

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17502
散金: 883
紅花: 18
帖子: 6215
在線: 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)槟闼鶈柕倪@兩句，是在氣相質(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)，受教了。我本來沒做多相流，但對(duì)這個(gè)有點(diǎn)興趣。
再問幾個(gè)小白的問題可以不。比如我按下面這樣理解對(duì)頭不：
氣相區(qū)----就是指這個(gè)區(qū)域的CELL只有氣相；
液相區(qū)----就是指這個(gè)區(qū)域的CELL只有液相；
混合區(qū)----就是指這個(gè)區(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小時(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)，受教了。我本來沒做多相流，但對(duì)這個(gè)有點(diǎn)興趣。
再問幾個(gè)小白的問題可以不。比如我按下面這樣理解對(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ā)表

[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-30	11/550	2026-03-31 16:34 by 邱gl
[考研] 材料科學(xué)與工程調(diào)劑 +13	深V宿舍吧 2026-03-30	14/700	2026-03-31 16:28 by hypershenger
[考研] 物理學(xué)調(diào)劑 +4	小羊36 2026-03-30	4/200	2026-03-31 16:16 by lishahe
[考研] 考研材料工程351分調(diào)劑 +3	整個(gè)好的 2026-03-31	3/150	2026-03-31 16:14 by 今天星期八LQ
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +5	GENJIOW 2026-03-31	5/250	2026-03-31 13:58 by limeifeng
[考研] 286求調(diào)劑 +5	丟掉懶惰 2026-03-27	8/400	2026-03-31 11:27 by Delta2012
[考研] 289求調(diào)劑 +6	BrightLL 2026-03-29	6/300	2026-03-31 11:22 by oooqiao
[考研] 0817化工學(xué)碩調(diào)劑 +7	努力上岸中！ 2026-03-31	7/350	2026-03-31 09:58 by nalakaiqi
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4	橙子（勝意） 2026-03-30	4/200	2026-03-31 00:44 by Linzejun
[考研] 285求調(diào)劑 +6	AZMK 2026-03-29	9/450	2026-03-30 21:02 by dophin1985
[考研] 342求調(diào)劑 +4	加油a李zs 2026-03-26	4/200	2026-03-30 16:39 by 晶體之美
[考研] 284求調(diào)劑 +14	junqihahaha 2026-03-26	15/750	2026-03-30 14:12 by 探123
[考研] 生物技術(shù)與工程 +7	1294608413 2026-03-25	8/400	2026-03-30 11:36 by 唐沐兒
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4	IceGuo 2026-03-26	5/250	2026-03-29 01:25 by griffith2014
[考研] 0856，材料與化工321分求調(diào)劑 +12	大饞小子 2026-03-27	13/650	2026-03-28 10:56 by self2008
[考研] 275求調(diào)劑 +10	jjjjjjjjjjl 2026-03-27	10/500	2026-03-27 23:47 by barnett0632
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 085600，材料與化工321分，求調(diào)劑 +9	大饞小子 2026-03-27	9/450	2026-03-27 14:30 by mmm just
[考研] 085600，材料與化工321分調(diào)劑 +4	大饞小子 2026-03-27	6/300	2026-03-27 14:11 by 松花缸1201

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

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

enigmafwj

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

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

» 猜你喜歡

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

enigmafwj

kylafree

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

kylafree

enigmafwj

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

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

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