| 查看: 6027 | 回復(fù): 21 | ||||
enigmafwj捐助貴賓 (小有名氣)
|
[求助]
【求助】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; } |
多相流UDF | 傳熱傳質(zhì) |
|
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ū)液相指針”了呢 |
鐵蟲 (小有名氣)
捐助貴賓 (小有名氣)
|
因?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指的是次相(氣相)。 |
|
感謝兄臺(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ā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過,有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|