版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

enigmafwj

捐助貴賓 (小有名氣)

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

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

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

DEFINE_SOURCE(liq_m,c,l_t,dS,eqn) /*液相質量輸運源項*/
{
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;                      /*液相質量輸運量*/
real stc;                         /*液氧表面張力系數(shù)*/
real fv;                         /*氣相質量分數(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) /*氣相質量輸運源項*/
{
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;                      /*氣相質量輸運量*/
real stc;                         /*液氧表面張力系數(shù)*/
real fv;                         /*氣相質量分數(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;                      /*相變質量*/
real hc;                         /*汽化潛熱*/
real source;                      /*能量方程源項*/
real stc;                         /*液氧表面張力系數(shù)*/
real fv;                         /*氣相質量分數(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;
}

回復此樓

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

多相流UDF

傳熱傳質

» 猜你喜歡

0856求調(diào)劑已經(jīng)有13人回復
0703化學已經(jīng)有11人回復
329求調(diào)劑已經(jīng)有7人回復
321求調(diào)劑已經(jīng)有8人回復
0856材料化工調(diào)劑總分330 已經(jīng)有12人回復
309求調(diào)劑已經(jīng)有3人回復
070300化學354求調(diào)劑已經(jīng)有4人回復
330分求調(diào)劑已經(jīng)有3人回復
299求調(diào)劑已經(jīng)有8人回復
一志愿北京理工大學本科211材料工程294求調(diào)劑已經(jīng)有6人回復

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

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

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

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

kylafree

至尊木蟲 (知名作家)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

查看全部 22 個回答

kylafree

至尊木蟲 (知名作家)

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

【答案】應助回帖

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

real fv;                         /*氣相質量分數(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;                         /*氣相質量分數(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人)

回復此樓

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

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

enigmafwj

捐助貴賓 (小有名氣)

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

引用回帖:

2樓: Originally posted by kylafree at 2012-07-30 21:19:25
real fv; /*氣相質量分數(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; /*液氧飽和蒸汽壓(考慮湍流影 ...

因為你所問的這兩句，是在氣相質量方程源項里的，DEFINE_SOURCE(vap_m,c,v_t,dS,eqn)里的thread：v_t在該宏應用于氣相質量方程時指向了氣相區(qū)。因為之后要用到混合區(qū)氣相的參數(shù)，所以用THREAD_SUPER_THREAD(v_t)把指針恢復到混合區(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指的是次相(氣相)。

贊一下

回復此樓

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

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

kylafree

至尊木蟲 (知名作家)

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

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

引用回帖:

3樓: Originally posted by enigmafwj at 2012-07-31 15:49:14
因為你所問的這兩句，是在氣相質量方程源項里的，DEFINE_SOURCE(vap_m,c,v_t,dS,eqn)里的thread：v_t在該宏應用于氣相質量方程時指向了氣相區(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)系起來啊，糊涂了

贊一下

回復此樓

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

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

查看全部 22 個回答

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

[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +9	辰沐5211314 2026-03-26	9/450	2026-03-29 01:12 by 我是小康
[考研] 調(diào)劑求院校招收 +5	鶴鯨鴿 2026-03-28	5/250	2026-03-29 01:10 by 我是小康
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +4	來財aa 2026-03-25	4/200	2026-03-28 07:41 by 棒棒球手
[考研] 求調(diào)劑推薦材料 304 +15	荷包蛋hyj 2026-03-26	15/750	2026-03-28 04:13 by fmesaito
[考研] 331環(huán)境科學與工程求調(diào)劑 +3	熠然好運氣 2026-03-27	3/150	2026-03-28 04:11 by fmesaito
[考研] 086000調(diào)劑 +3	7901117076 2026-03-26	3/150	2026-03-27 21:34 by Jianing_Mi
[考研] 化學調(diào)劑 +4	愛吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 305求調(diào)劑 +5	哇盧卡庫 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 333求調(diào)劑 +3	question挽風 2026-03-23	3/150	2026-03-27 11:29 by 不吃魚的貓
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 0703化學338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 329求調(diào)劑 +7	鈕恩雪 2026-03-25	7/350	2026-03-27 04:28 by wxiongid
[考研] 085602化學工程求調(diào)劑。 +4	平樂樂樂 2026-03-26	4/200	2026-03-26 17:57 by fmesaito
[考研] 297求調(diào)劑 +6	田洪有 2026-03-26	6/300	2026-03-26 15:55 by 不吃魚的貓
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 309求調(diào)劑 +4	gajsj 2026-03-25	5/250	2026-03-26 00:27 by Dyhoer
[考研] 考研一志愿蘇州大學初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5	超贊的 2026-03-24	5/250	2026-03-24 21:07 by 星空星月
[考研] 一志愿北化315 求調(diào)劑 +3	akrrain 2026-03-24	3/150	2026-03-24 19:35 by 了了了了。。

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

24小時熱門版塊排行榜

enigmafwj

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

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

» 猜你喜歡

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

kylafree

kylafree

【答案】應助回帖

enigmafwj

kylafree

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

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