 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 853 | 回復(fù): 1 | |||
kobe8875新蟲 (初入文壇)
|
[交流]
【求助】用species transport模型模擬蒸發(fā)UDF的疑問
|
|
我用某教材里的下面這個(gè)UDF模擬130攝氏度（或更高）的熱空氣加熱80攝氏度的水，熱空氣采用組分運(yùn)輸模型，多項(xiàng)流采用VOF模型（因?yàn)橛凶杂身撁妫醒鯕�、氮�(dú)�、水蒸氣�?br />
但是模擬的結(jié)果空氣中產(chǎn)生不出水蒸氣，而用FLUENT自帶的沸騰蒸發(fā)模型可以產(chǎn)生。 現(xiàn)在初步判斷可能是UDF中幾個(gè)指針有問題，但是不知道怎么改，或者是源項(xiàng)加載的地方不對。 其中，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)對來自液相的質(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)對自身的質(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)對質(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)對質(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ù)號 { 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 ] |
回復(fù)此樓» 猜你喜歡
307求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
-
301求調(diào)劑
已經(jīng)有10人回復(fù)
-
303求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料工程（專）一志愿985 初試335求調(diào)劑
已經(jīng)有4人回復(fù)
-
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑
已經(jīng)有3人回復(fù)
-
初始318分求調(diào)劑（有工作經(jīng)驗(yàn)）
已經(jīng)有3人回復(fù)
-
二本跨考鄭大材料306英一數(shù)二
已經(jīng)有3人回復(fù)
-
化學(xué)求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- VOF模型靜止壁面上速度不為0
已經(jīng)有10人回復(fù)
- 哪位大俠會(huì)fluent定義化學(xué)反應(yīng)
已經(jīng)有3人回復(fù)
- FLUENT加入VOF模型后的一些問題
已經(jīng)有15人回復(fù)
- 用UDF定義了volumetric reaction rate 還怎么在fluent里設(shè)置阿累尼烏斯定律參數(shù)
已經(jīng)有11人回復(fù)
- 采用VOF模型對磨削加工中的流場進(jìn)行仿真，仿真結(jié)果不對，摸索中，望高手相助��！
已經(jīng)有28人回復(fù)
- fluent中的species transport中如何正確輸入反應(yīng)參數(shù)
已經(jīng)有5人回復(fù)
- ansys13.0啟動(dòng)fluent出錯(cuò)
已經(jīng)有30人回復(fù)
- udf獲得多相流中某一項(xiàng)中某一組分的含量
已經(jīng)有14人回復(fù)
- Fluent，VOF模型，怎么在追蹤的界面上施加反沖壓
已經(jīng)有7人回復(fù)
- UDF的編譯通不過
已經(jīng)有8人回復(fù)
- vof模型和species transport模型到底沖突不
已經(jīng)有18人回復(fù)
- fluent-udf自定義標(biāo)量方程邊界條件
已經(jīng)有5人回復(fù)
- 在氣相混合中如何選擇正確的物理模型
已經(jīng)有7人回復(fù)
- fluent 關(guān)于prePDF模型求助
已經(jīng)有13人回復(fù)
- 組分輸運(yùn)模型的出口組分設(shè)置？
已經(jīng)有26人回復(fù)
- 【求助】通氣空泡多相流（氣、水蒸汽、水）空化模型的UDF如何入手啊？？？？
已經(jīng)有4人回復(fù)
- 【求助】求UDF高手幫忙看看我的問題
已經(jīng)有18人回復(fù)
- 【求助】用fluent模擬超重力氣液流體動(dòng)力模型和傳質(zhì)過程，中間有化學(xué)吸收反應(yīng)，fluent
已經(jīng)有23人回復(fù)
kobe8875
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 188
- 帖子: 19
- 在線: 8.7小時(shí)
- 蟲號: 2677845
- 注冊: 2013-09-24
- 專業(yè): 多相流熱物理學(xué)
2樓2013-10-17 12:20:08
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 kobe8875 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +9 | 戲精丹丹丹 2026-03-17 | 9/450 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +5 | 77chaselx 2026-03-20 | 5/250 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|