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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

renzhihao

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1682.5
帖子: 121
在線: 15.8小時
蟲號: 3808757
注冊: 2015-04-14
專業(yè): 海洋工程

[求助] 想模擬單個氣泡的潰滅，想用fluent編程把壓力導(dǎo)入進(jìn)去，用的vof模型。

我在做空化，想模擬單個氣泡的潰滅，但是現(xiàn)在要求內(nèi)部流體的壓力是個隨著位置變化的函數(shù)，想用fluent編程把壓力導(dǎo)入進(jìn)去，用的vof模型。

下面的為我從別的帖子里面看到的。想問下這些東西需要在fluent里面設(shè)置還是編程？
為液相、氣相質(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                               /*氣相液化，溫度降低，源項為負(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;
}
這些東西直接在fluent里面設(shè)置行嗎，還是必須要編程？

回復(fù)此樓

» 猜你喜歡

086003食品工程求調(diào)劑已經(jīng)有3人回復(fù)
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑已經(jīng)有12人回復(fù)
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑已經(jīng)有12人回復(fù)
材料調(diào)劑已經(jīng)有5人回復(fù)
081700 調(diào)劑 267分已經(jīng)有8人回復(fù)
請教下大家 2026年國家基金申請是雙盲審嗎？已經(jīng)有5人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷已經(jīng)有10人回復(fù)
一志愿南航材料專317分求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)學(xué)碩求調(diào)劑已經(jīng)有5人回復(fù)

1樓 2016-04-21 17:27:53

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

jiachaolis

銅蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 184.1
紅花: 2
帖子: 134
在線: 51.5小時
蟲號: 3307685
注冊: 2014-07-04
性別: GG
專業(yè): 動力學(xué)與控制

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

樓主，請把你相變的模型公式貼出來，這樣就好查錯誤了

贊一下(1人)

回復(fù)此樓

付出就會有收獲

2樓2016-04-24 20:01:44

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

renzhihao

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1682.5
帖子: 121
在線: 15.8小時
蟲號: 3808757
注冊: 2015-04-14
專業(yè): 海洋工程

引用回帖:

2樓: Originally posted by jiachaolis at 2016-04-24 20:01:44
樓主，請把你相變的模型公式貼出來，這樣就好查錯誤了

怎么查錯誤

回復(fù)此樓

3樓2016-04-25 16:41:57

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

renzhihao

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1682.5
帖子: 121
在線: 15.8小時
蟲號: 3808757
注冊: 2015-04-14
專業(yè): 海洋工程

引用回帖:

2樓: Originally posted by jiachaolis at 2016-04-24 20:01:44
樓主，請把你相變的模型公式貼出來，這樣就好查錯誤了

怎么查錯誤

回復(fù)此樓

4樓2016-04-26 16:53:41

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

wolfghost

版主 (知名作家)

博學(xué)EPI: 432
應(yīng)助: 101 (高中生)
貴賓: 7.987
金幣: 36647.2
散金: 11283
紅花: 198
沙發(fā): 8
帖子: 7791
在線: 1298小時
蟲號: 573180
注冊: 2008-06-13
性別: GG
專業(yè): 電化學(xué)
管轄: 有獎問答

樓主，長期無人應(yīng)助金幣可直接申請返還(點擊右下角“通知版主處理”)，重新開貼以提高求助成功率！為營造互助、和諧論壇共同努力，謝謝~

贊一下

回復(fù)此樓

大刀向鬼子的頭上砍去！ChopJapstoDeath!

5樓2016-05-25 13:50:12

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

聆聽風(fēng)聲love

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 821
散金: 100
紅花: 1
帖子: 273
在線: 170.9小時
蟲號: 3694947
注冊: 2015-02-24
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

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

質(zhì)量源項的方程出自哪里?

贊一下

回復(fù)此樓

努力！奮斗！

6樓2016-07-03 11:41:17

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

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

返回列表

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

[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
[考研] 361求調(diào)劑 +3	Glack 2026-03-22	3/150	2026-03-23 22:03 by fuyu_
[考研] 材料專碩英一數(shù)二306 +8	z1z2z3879 2026-03-18	8/400	2026-03-23 20:49 by baobaoye
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5	葵梓衛(wèi)隊 2026-03-18	7/350	2026-03-23 16:26 by lingjue
[考研] 材料與化工考研調(diào)劑 +4	孅華 2026-03-22	4/200	2026-03-23 16:13 by 一休哥FU
[考研] 招08考數(shù)學(xué) +6	laoshidan 2026-03-20	14/700	2026-03-23 14:37 by 15614359529
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	9/450	2026-03-23 12:57 by yqdszhdap－
[考研] 求調(diào)劑材料學(xué)碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 287求調(diào)劑 +8	晨昏線與星海 2026-03-19	9/450	2026-03-22 17:01 by i_cooler
[考研] 生物學(xué)071000 329分求調(diào)劑 +5	我愛生物生物愛?/a> 2026-03-17	5/250	2026-03-22 16:42 by tcx007
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne

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

24小時熱門版塊排行榜

renzhihao

[求助] 想模擬單個氣泡的潰滅，想用fluent編程把壓力導(dǎo)入進(jìn)去，用的vof模型。

» 猜你喜歡

jiachaolis

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

renzhihao

renzhihao

wolfghost

聆聽風(fēng)聲love

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

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

[求助] 想模擬單個氣泡的潰滅，想用fluent編程把壓力導(dǎo)入進(jìn)去，用的vof模型。