 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1598 | 回復(fù): 5 | |||
簡(jiǎn)單的快樂(lè)xf銅蟲(chóng) (初入文壇)
|
[求助]
關(guān)于水蒸發(fā)UDF問(wèn)題 已有2人參與
|
|
在網(wǎng)上看到的關(guān)于水蒸發(fā)的UDF，想問(wèn)如果改為其他種類(lèi)液體（如液化天然氣）需要改變程序哪些地方？不勝感激�。。� #include "udf.h" ‘包括常規(guī)宏 #include "sg_mphase.h" ‘包括體積分?jǐn)?shù)宏CVOF(C,T) #define T_SAT 373 ‘定義蒸發(fā)溫度100℃ #define LAT_HT 1.e3 ‘定義蒸發(fā)潛熱J/Kg DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) ‘液相質(zhì)量源項(xiàng)UDF { Thread *mix_th, *sec_th; ‘定義計(jì)算區(qū)域線指針 real m_dot_l; ‘定義液相質(zhì)量轉(zhuǎn)移 kg/(m2.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) ‘如果液相單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移 { m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dS[eqn] = -0.1*C_R(cell, pri_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, sec_th)*C_R(cell, sec_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(vap_src, cell, sec_th, dS, eqn) ‘氣相質(zhì)量源項(xiàng)UDF { Thread * mix_th, *pri_th; real m_dot_v; 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) ‘如果混合區(qū)單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量質(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; dS[eqn] = 0.; ‘由于是液相向氣相轉(zhuǎn)移，所以氣相的質(zhì)量源項(xiàng)對(duì)來(lái)自液相的質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零 } 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; ‘如果指向混合區(qū)的單元溫度小于蒸發(fā)溫度，氣相向液相的質(zhì)量轉(zhuǎn)移 dS[eqn] = -0.1*C_R(cell, sec_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_v; } 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, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT; } else { 相反，氣相向液相轉(zhuǎn)移則放熱 m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;} return LAT_HT*m_dot; 氣化潛熱與質(zhì)量轉(zhuǎn)移率相乘得能量源項(xiàng)W/m3 } |
回復(fù)此樓» 收錄本帖的淘帖專(zhuān)輯推薦
 博士科研新路歷程 |
» 猜你喜歡
283求調(diào)劑（080500）
已經(jīng)有3人回復(fù)
-
286求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有4人回復(fù)
-
308求調(diào)劑
已經(jīng)有6人回復(fù)
-
275求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
287求調(diào)劑
已經(jīng)有9人回復(fù)
-
274求調(diào)劑
已經(jīng)有23人回復(fù)
-
085602 化工專(zhuān)碩 338分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有16人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 請(qǐng)教一個(gè)關(guān)于閉式虹吸環(huán)路的二維帶相變的模擬
已經(jīng)有1人回復(fù)
- 關(guān)于用不同宏來(lái)描述揮發(fā)冷凝過(guò)程中的質(zhì)量傳輸問(wèn)題
已經(jīng)有5人回復(fù)
- 【求助】求助或討論一下關(guān)于能量傳遞源項(xiàng)中各個(gè)焓值的含義
已經(jīng)有3人回復(fù)
2樓2015-11-02 08:53:59
簡(jiǎn)單的快樂(lè)xf
銅蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 65.5
- 帖子: 36
- 在線: 31.7小時(shí)
- 蟲(chóng)號(hào): 3105188
- 注冊(cè): 2014-04-01
- 專(zhuān)業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
3樓2015-11-20 21:58:57
簡(jiǎn)單的快樂(lè)xf4樓2015-11-23 10:37:48
簡(jiǎn)單的快樂(lè)xf
銅蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 65.5
- 帖子: 36
- 在線: 31.7小時(shí)
- 蟲(chóng)號(hào): 3105188
- 注冊(cè): 2014-04-01
- 專(zhuān)業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
5樓2015-11-29 11:27:56
|
本帖內(nèi)容被屏蔽 |
6樓2017-04-10 17:42:32
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 簡(jiǎn)單的快樂(lè)xf 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書(shū)讀?/a> 2026-03-26 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
|
[考研] 277跪求調(diào)劑 +4 | 1915668 2026-03-27 | 5/250 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 安徽大學(xué)專(zhuān)碩生物與醫(yī)藥專(zhuān)業(yè)(086000)324分，英語(yǔ)已過(guò)四六級(jí)，六級(jí)521，求調(diào)劑 +3 | 美味可樂(lè)雞翅 2026-03-26 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 食品工程專(zhuān)碩求調(diào)劑 +3 | 小張zxy張 2026-03-26 | 3/150 |
|
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
|
[考研] 尋找調(diào)劑 +5 | 倔強(qiáng)芒? 2026-03-21 | 8/400 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
|
[考研] 一志愿國(guó)科過(guò)程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|