| 5 | 1/1 | 返回列表 |
| 查看: 1874 | 回復: 9 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
FLUENT中UDF編譯始終報錯,氣液兩相流問題,求大神指點 已有3人參與
|
|||
|
#include "udf.h" //常規(guī)宏 #include "sg_mphase.h" //系統(tǒng)頭文件,體積分數(shù)宏 #define T_SAT 373 //沸騰溫度373K #define LAT_HT 2257600 //氣化潛熱2257.6kJ/kg DEFINE_SOURCE(pri_src,cell,pri_th,dS,eqn) //液相質(zhì)量源項UDF { Thread *mix_th; Thread *sec_th; //定義指針 real m_dot_l; //液相質(zhì)量轉(zhuǎn)移kg/(㎡·s) mix_th=THREAD_SUPER_THREAD(pri_th); //混合區(qū)主相為液相指針 sec_th=THREAD_SUB_THREAD(mix_th,1); //單相控制區(qū)的氣相指針,氣相為第二相 if(C_T(cell,mix_th)>=T_SAT) //混合區(qū)液相溫度高于蒸發(fā)溫度,液相向氣相質(zhì)量轉(zhuǎn)移 { m_dot_l=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_R(cell,mix_th)-T_SAT)/T_SAT; //液相向氣相質(zhì)量轉(zhuǎn)移 dS[eqn]=-0.1*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; //定義源項為對質(zhì)量轉(zhuǎn)移偏導 } else //混合區(qū)液相溫度低于蒸發(fā)溫度 { m_dot_l=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; //氣相向液相轉(zhuǎn)移,液相增加 dS[eqn]=0;//液相質(zhì)量源項對質(zhì)量的偏導為零 } return m_dot_l; } DEFINE_SOURCE(sec_src,cell,sec_th,dS,eqn) //氣相質(zhì)量源項UDF { Thread *mix_th; Thread *pri_th; real m_dot_v; mix_th=THREAD_SUPER_THREAD(sec_th); //混合區(qū)第二項為氣相 pri_th=THREAD_SUB_THREAD(mix_th,0); //液相為第一相 if(C_T(cell,mix_th)>=T_SAT) //假設混合區(qū)溫度高于蒸發(fā)溫度,液相向氣相的質(zhì)量轉(zhuǎn)移 { m_dot_v=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; //液相向氣相轉(zhuǎn)移 dS[eqn]=0;//氣相的質(zhì)量源項對來自液相的質(zhì)量轉(zhuǎn)移的偏導為零 } else { m_dot_v=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; //氣相向液相質(zhì)量轉(zhuǎn)移 dS[eqn]=-0.1*C_R(cell,sec_th)*fabs(C_T(cell,sec_th)-T_SAT)/T_SAT; //因為氣相向液相轉(zhuǎn)移,所以氣相的質(zhì)量源項對自身的質(zhì)量轉(zhuǎn)移的偏導不為零 } return m_dot_v; } DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn) //混合模型源項UDF { Thread *pri_th; Thread *sec_th; real m_dot; pri_th=THREAD_SUB_THREAD(mix_th,0); //指向混合區(qū)液相的指針 sec_th=THREAD_SUB_THREAD(mix_th,1); //指向混合區(qū)氣相的指針 if(C_T(cell,mix_th)>=T_SAT) //假設混合區(qū)溫度高于蒸發(fā)溫度,液相向氣相質(zhì)量轉(zhuǎn)移,吸熱(質(zhì)量轉(zhuǎn)移為負) { m_dot=-0.1*C_VOF(cell,pri_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; //吸熱,液相向氣相轉(zhuǎn)移 dS[eqn]=-0.1*C_VOF(cell,pri_th)*fabs(C_R(cell,pri_th)/T_SAT; } else { m_dot=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; //放熱,氣相向液相轉(zhuǎn)移 dS[eqn]=0.1*C_VOF(cell,pri_th)*fabs(C_R(cell,pri_th)/T_SAT; return LAT_HT*m_dot; //氣化潛熱*質(zhì)量轉(zhuǎn)移率=能量源項W/m3 } cpp -I"D:\fluent\Fluent.Inc\fluent6.3.26/src" -I"D:\fluent\Fluent.Inc\fluent6.3.26/cortex/src" -I"D:\fluent\Fluent.Inc\fluent6.3.26/client/src" -I"D:\fluent\Fluent.Inc\fluent6.3.26/multiport/src" -I. -DUDFCONFIG_H="<udfconfig.h>" "H:\UDF\UDF .c" Error: H:\UDF\UDF .c: line 21: parse error. 一直報錯,求大神指點 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +181 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學求調(diào)劑,各位老師看看我。! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 化學調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 299求調(diào)劑 +6 | 嗯嗯嗯嗯2 2026-03-27 | 6/300 |
|
|
[考研]
|
孅華 2026-03-22 | 16/800 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 機械學碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[有機交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 一志愿吉大化學322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考研] 一志愿山東大學藥學學碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學學碩304分,有大創(chuàng)項目,四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|