 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2245 | 回復(fù): 5 | ||
renzhihao新蟲 (小有名氣)
|
[求助]
想模擬單個(gè)氣泡的潰滅，想用fluent編程把壓力導(dǎo)入進(jìn)去，用的vof模型。
|
|
我在做空化，想模擬單個(gè)氣泡的潰滅，但是現(xiàn)在要求內(nèi)部流體的壓力是個(gè)隨著位置變化的函數(shù)，想用fluent編程把壓力導(dǎo)入進(jìn)去，用的vof模型。 下面的為我從別的帖子里面看到的。想問下這些東西需要在fluent里面設(shè)置還是編程？ 為液相、氣相質(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; } 這些東西直接在fluent里面設(shè)置行嗎，還是必須要編程？ |
回復(fù)此樓» 猜你喜歡
材料專業(yè)求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有8人回復(fù)
-
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑
已經(jīng)有6人回復(fù)
-
請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎？
已經(jīng)有4人回復(fù)
-
一志愿南航材料專317分求調(diào)劑
已經(jīng)有3人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271
已經(jīng)有4人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有8人回復(fù)
-
化學(xué)308分求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有11人回復(fù)
jiachaolis
銅蟲 (小有名氣)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 184.1
- 紅花: 2
- 帖子: 134
- 在線: 51.5小時(shí)
- 蟲號(hào): 3307685
- 注冊(cè): 2014-07-04
- 性別: GG
- 專業(yè): 動(dòng)力學(xué)與控制
2樓2016-04-24 20:01:44
renzhihao
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1682.5
- 帖子: 121
- 在線: 15.8小時(shí)
- 蟲號(hào): 3808757
- 注冊(cè): 2015-04-14
- 專業(yè): 海洋工程
3樓2016-04-25 16:41:57
renzhihao
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1682.5
- 帖子: 121
- 在線: 15.8小時(shí)
- 蟲號(hào): 3808757
- 注冊(cè): 2015-04-14
- 專業(yè): 海洋工程
4樓2016-04-26 16:53:41
wolfghost
版主 (知名作家)
- 博學(xué)EPI: 432
- 應(yīng)助: 101 (高中生)
- 貴賓: 7.987
- 金幣: 36647.2
- 散金: 11283
- 紅花: 198
- 沙發(fā): 8
- 帖子: 7791
- 在線: 1298小時(shí)
- 蟲號(hào): 573180
- 注冊(cè): 2008-06-13
- 性別: GG
- 專業(yè): 電化學(xué)
- 管轄: 有獎(jiǎng)問答
5樓2016-05-25 13:50:12
聆聽風(fēng)聲love
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 821
- 散金: 100
- 紅花: 1
- 帖子: 273
- 在線: 170.9小時(shí)
- 蟲號(hào): 3694947
- 注冊(cè): 2015-02-24
- 性別: GG
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
6樓2016-07-03 11:41:17
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 renzhihao 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +4 | jhybd 2026-03-23 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料292調(diào)劑 +3 | 橘頌思美人 2026-03-23 | 3/150 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +7 | 用戶573181 2026-03-20 | 11/550 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 9/450 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +5 | vv迷 2026-03-22 | 5/250 |
|
|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|