| 24小時(shí)熱門(mén)版塊排行榜 |
| 2 | 1/1 | 返回列表 |
| 查看: 854 | 回復(fù): 1 | |||
kobe8875新蟲(chóng) (初入文壇)
|
[交流]
【求助】用species transport模型模擬蒸發(fā)UDF的疑問(wèn)
|
|
我用某教材里的下面這個(gè)UDF模擬130攝氏度(或更高)的熱空氣加熱80攝氏度的水,熱空氣采用組分運(yùn)輸模型,多項(xiàng)流采用VOF模型(因?yàn)橛凶杂身?yè)面),含有氧氣、氮?dú)狻⑺魵狻?br />
但是模擬的結(jié)果空氣中產(chǎn)生不出水蒸氣,而用FLUENT自帶的沸騰蒸發(fā)模型可以產(chǎn)生。 現(xiàn)在初步判斷可能是UDF中幾個(gè)指針有問(wèn)題,但是不知道怎么改,或者是源項(xiàng)加載的地方不對(duì)。 其中,air_src加給空氣中的H2O組分,liq_src加給水,enrg_src加給混合物 麻煩各位幫助,謝謝。。 #include "udf.h" //包括常規(guī)宏 #include "sg_mphase.h" #define T_SAT 373 //定義蒸發(fā)溫度100℃ #define LAT_HT 2.3e6 //定義蒸發(fā)潛熱J/Kg DEFINE_SOURCE(air_src,cell,pri_th,dS,eqn) //氣相質(zhì)量源項(xiàng)UDF { Thread * mix_th, *sec_th; real m_dot_a; 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ì)量質(zhì)量轉(zhuǎn)移 { m_dot_a = 0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn] = 0; //由于是液相向氣相轉(zhuǎn)移,所以氣相的質(zhì)量源項(xiàng)對(duì)來(lái)自液相的質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零 } else { m_dot_a = -0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; //如果指向混合區(qū)的單元溫度小于蒸發(fā)溫度,氣相向液相的質(zhì)量轉(zhuǎn)移,氣相失 dS[eqn] = -0.1*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; //由于是氣相向液相轉(zhuǎn)移,所以氣相的質(zhì)量源項(xiàng)對(duì)自身的質(zhì)量轉(zhuǎn)移的偏導(dǎo)不為零 } return m_dot_a; } DEFINE_SOURCE(liq_src,cell,sec_th,dS,eqn) //液相質(zhì)量源項(xiàng)UDF { Thread *mix_th, *pri_th; //定義計(jì)算區(qū)線指針 real m_dot_l; //定義液相質(zhì)量轉(zhuǎn)移 kg/(m2.s) mix_th = THREAD_SUPER_THREAD(sec_th); //指向混合區(qū)的主相即液相的指針 pri_th = THREAD_SUB_THREAD(mix_th,0); //指向單相控制區(qū)的氣相的指針,氣相為主相 if(C_T(cell, mix_th) >= T_SAT) //如果液相單元的溫度高于蒸發(fā)溫度,液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移 { m_dot_l = -0.1*C_VOF(cell,sec_th)*C_R(cell, sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn] = -0.1*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; //定義源項(xiàng)對(duì)質(zhì)量轉(zhuǎn)移偏導(dǎo) } else { m_dot_l = 0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; //如果指向混合區(qū)液相的單元溫度小于蒸發(fā)溫度,氣相向液相的質(zhì)量轉(zhuǎn)移,液相得 dS[eqn] = 0; //由于是氣相向液相轉(zhuǎn)移,所以液相的質(zhì)量源項(xiàng)對(duì)質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零 } return m_dot_l; } DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn) //混合模型能量源項(xiàng)UDF { Thread *pri_th, *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)移由液相向氣相轉(zhuǎn)移,吸熱,質(zhì)量轉(zhuǎn)移量前有負(fù)號(hào) { m_dot = -0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn] = -0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)/T_SAT; } else { //相反,氣相向液相轉(zhuǎn)移則放熱 m_dot = 0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = 0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT; } return LAT_HT*m_dot; //氣化潛熱與質(zhì)量轉(zhuǎn)移率相乘得能量源項(xiàng)W/m3 } [ Last edited by kobe8875 on 2013-10-18 at 14:31 ] |
新蟲(chóng) (初入文壇)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專(zhuān)碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 294求調(diào)劑材料與化工專(zhuān)碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 材料專(zhuān)碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|