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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

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

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ì)量方程和能量方程添加源項(xiàng)，UDF程序如下，汽化、液化由壓力控制（壓力條件判斷），各物性參數(shù)用溫度擬合。問題是：編譯時，提示有語法錯誤(parse)，并且fv為未聲明的變量。請問如何解決？另：程序中有其他不妥之處，煩請指教，謝謝！
#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;                /*計算區(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;                /*計算區(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;                /*計算區(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;
}
這些東西直接在fluent里面設(shè)置行嗎，還是必須要編程？

回復(fù)此樓

» 猜你喜歡

南京大學(xué)化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
一志愿北京工業(yè)大學(xué)，324分求調(diào)劑已經(jīng)有6人回復(fù)
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑已經(jīng)有14人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑已經(jīng)有6人回復(fù)
26申博自薦已經(jīng)有6人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有10人回復(fù)
求調(diào)劑有機(jī)化學(xué)考研356分已經(jīng)有5人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有10人回復(fù)
材料與化工328分調(diào)劑已經(jīng)有7人回復(fù)

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

已閱回復(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的回帖

查看全部 6 個回答

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的回帖

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)助金幣可直接申請返還(點(diǎn)擊右下角“通知版主處理”)，重新開貼以提高求助成功率！為營造互助、和諧論壇共同努力，謝謝~

贊一下

回復(fù)此樓

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

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

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

查看全部 6 個回答

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

[考研] 281求調(diào)劑 +4	亞克西good 2026-03-26	6/300	2026-03-29 21:06 by 啊李999
[考研] 求調(diào)劑 +4	青春裁為三截 2026-03-29	4/200	2026-03-29 15:01 by 唐沐兒
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7	@taotao 2026-03-29	7/350	2026-03-29 12:03 by longlotian
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +7	小小甜面團(tuán) 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落睿可思
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9	一定上岸aaa 2026-03-27	13/650	2026-03-28 15:00 by 神馬都不懂
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito
[考研] 070300化學(xué)求調(diào)劑 +4	起個名咋這么難 2026-03-27	4/200	2026-03-27 21:39 by 83503孫老師
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6	kotoko_ik 2026-03-23	7/350	2026-03-27 12:29 by 惠州彭于晏
[考研] 求調(diào)劑 +3	劉柯@ 2026-03-24	4/200	2026-03-27 11:28 by shangxh
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte

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

24小時熱門版塊排行榜

renzhihao

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

» 猜你喜歡

renzhihao

jiachaolis

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

renzhihao

wolfghost

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

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