 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1601 | 回復: 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 } |
回復此樓» 收錄本帖的淘帖專輯推薦
 博士科研新路歷程 |
» 猜你喜歡
0856，材料與化工321分求調劑
已經(jīng)有6人回復
-
0703化學求調劑，各位老師看看我�。�！
已經(jīng)有5人回復
-
266分，求材料冶金能源化工等調劑
已經(jīng)有7人回復
-
張芳銘-中國農業(yè)大學-環(huán)境工程專碩-298
已經(jīng)有3人回復
-
070300化學求調劑
已經(jīng)有4人回復
-
086000調劑
已經(jīng)有3人回復
-
320分，材料與化工專業(yè)，求調劑
已經(jīng)有4人回復
-
265求調劑
已經(jīng)有7人回復
-
295材料工程專碩求調劑
已經(jīng)有4人回復
-
高溫高壓反應求助
已經(jīng)有4人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 請教一個關于閉式虹吸環(huán)路的二維帶相變的模擬
已經(jīng)有1人回復
- 關于用不同宏來描述揮發(fā)冷凝過程中的質量傳輸問題
已經(jīng)有5人回復
- 【求助】求助或討論一下關于能量傳遞源項中各個焓值的含義
已經(jīng)有3人回復
3樓2015-11-20 21:58:57
2樓2015-11-02 08:53:59
簡單的快樂xf4樓2015-11-23 10:37:48
5樓2015-11-29 11:27:56
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295求調劑 +4 | 1428151015 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京化工大學 070300 學碩 336分 求調劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 07化學280分求調劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 調劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿陜師大生物學071000，298分，求調劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研] 求調劑，一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
|
[考研] 284求調劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研]
 304材料求調劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研] 281求調劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
|
[考研] 303求調劑 +6 | 藍山月 2026-03-25 | 6/300 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
|
[考研] 290分調劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 生物技術與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 086003食品工程求調劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 318求調劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 269求調劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
|
[考研] 280分求調劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 求調劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 一志愿華中科技大學071000，求調劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|