 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1607 | 回復: 5 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
關于水蒸發(fā)UDF問題 已有2人參與
|
|||
|
在網(wǎng)上看到的關于水蒸發(fā)的UDF，想問如果改為其他種類液體（如液化天然氣）需要改變程序哪些地方？不勝感激！�。� #include "udf.h" ‘包括常規(guī)宏 #include "sg_mphase.h" ‘包括體積分數(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) ‘液相質量源項UDF { Thread *mix_th, *sec_th; ‘定義計算區(qū)域線指針 real m_dot_l; ‘定義液相質量轉移 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ā)溫度，液相向氣相的質量質量轉移 { 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; ‘ 定義源項對質量轉移偏導 ， } 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ā)溫度，氣相向液相的質量轉移，液相得 dS[eqn] = 0.; ‘由于是氣相向液相轉移，所以液相的質量源項對質量轉移的偏導為零 } return m_dot_l; } DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) ‘氣相質量源項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ā)溫度，液相向氣相的質量質量轉移 { 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.; ‘由于是液相向氣相轉移，所以氣相的質量源項對來自液相的質量轉移的偏導為零 } 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ā)溫度，氣相向液相的質量轉移 dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; 由于是氣相向液相轉移，所以氣相的質量源項對自身的質量轉移的偏導不為零 } return m_dot_v; } DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) ‘混合模型能量源項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ā)溫度。質量轉移由液相向氣相轉移，吸 熱，質量轉移量前有負號 { 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 { 相反，氣相向液相轉移則放熱 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; 氣化潛熱與質量轉移率相乘得能量源項W/m3 } |
回復此樓» 收錄本帖的淘帖專輯推薦
 博士科研新路歷程 |
» 猜你喜歡
307求調劑
已經(jīng)有20人回復
-
285求調劑
已經(jīng)有8人回復
-
材料與化工（0856）304求B區(qū)調劑
已經(jīng)有8人回復
-
370求調劑
已經(jīng)有3人回復
-
070300化學279求調劑
已經(jīng)有5人回復
-
福建理工大學材料學院先進合金團隊招收考研調劑學生
已經(jīng)有4人回復
-
英一數(shù)一408，總分284，二戰(zhàn)真誠求調劑
已經(jīng)有5人回復
-
083000環(huán)境科學與工程調劑，總分281
已經(jīng)有4人回復
-
求調劑
已經(jīng)有3人回復
-
食品工程專碩一志愿中海洋309求調劑
已經(jīng)有10人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 請教一個關于閉式虹吸環(huán)路的二維帶相變的模擬
已經(jīng)有1人回復
- 關于用不同宏來描述揮發(fā)冷凝過程中的質量傳輸問題
已經(jīng)有5人回復
- 【求助】求助或討論一下關于能量傳遞源項中各個焓值的含義
已經(jīng)有3人回復
|
本帖內容被屏蔽 |
6樓2017-04-10 17:42:32
2樓2015-11-02 08:53:59
3樓2015-11-20 21:58:57
簡單的快樂xf4樓2015-11-23 10:37:48
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學與工程調劑 +9 | 深V宿舍吧 2026-03-30 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 22408 359分調劑 +4 | Qshers 2026-03-27 | 7/350 |
|
|
|
[考研] 288資源與環(huán)境專碩求調劑，不限專業(yè)，有學上就行 +11 | lllllos 2026-03-30 | 11/550 |
|
|
|
[考研] 317求調劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
|
[考研] 生物學學碩，一志愿湖南大學，初試成績338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
|
[考研] 一志愿：西北大學，英一數(shù)一408-284分求調劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
|
[考研] 324求調劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
|
[考研] 求調劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
|
[考研] 0703化學/290求調劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
|
[考研] 環(huán)境科學與工程334分求調劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
|
[考研] 337求調劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
|
[考研] 求調劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
|
[考研] 289求調劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
|
[考研] 311（085601）求調劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
|
[考研] 347求調劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 086502化學工程342求調劑 +6 | 阿姨復古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 求調劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 266求調劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 305求調劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[碩博家園] 北京林業(yè)大學碩導招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|