| 2 | 1/1 | 返回列表 |
| 查看: 475 | 回復(fù): 1 | |||
up2uggyy銅蟲 (初入文壇)
|
[交流]
【求助】水蒸發(fā) UDF 的應(yīng)用 已有1人參與
|
網(wǎng)上的水蒸發(fā)UDF如下,但是不知道如何使用.... ,這個(gè)UDF 應(yīng)該用在哪里?流體的屬性里面嗎?具體是哪項(xiàng)屬性?高手來指點(diǎn)下啊~~~~~#include "udf.h" /*包括常規(guī)宏*/ #include "sg_mphase.h" /*包括體積分?jǐn)?shù)宏*/ #define T_SAT 373 /*定義飽和溫度*/ #define LAT_HT 1.e3 /*定義水蒸汽潛熱*/ DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) /*定義液相源項(xiàng)*/ { Thread *mix_th, *sec_th; /*混合相、第二相 定義計(jì)算區(qū)域指針*/ real m_dot_l; /*蒸發(fā)冷凝速率 定義液相質(zhì)量轉(zhuǎn)移*/ 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){ m_dot_l = -0.02*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*如果液相單元的溫度高于蒸發(fā)溫度,液相向氣相的質(zhì)量轉(zhuǎn)移*/ dS[eqn] = -0.02*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*定義源項(xiàng)對(duì)質(zhì)量轉(zhuǎn)移偏導(dǎo)*/ } else { m_dot_l = 0.02*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th))/T_SAT; 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) /*定義氣相源項(xiàng)*/ { 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){ m_dot_v = 0.02*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; /*如果混合區(qū)單元的溫度高于蒸發(fā)溫度,液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移 */ dS[eqn] = 0.; /*由于是液相向氣相轉(zhuǎn)移,所以氣相的質(zhì)量源項(xiàng)對(duì)來自液相的質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零*/ } else { m_dot_v = -0.02*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = -0.02*C_R(cell, sec_th)* fabs(C_T(cell, sec_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)*/ { 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){ m_dot = -0.02*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*如果混合區(qū)的單元溫度高于蒸發(fā)溫度。質(zhì)量轉(zhuǎn)移由液相向氣相轉(zhuǎn) 移,吸熱,質(zhì)量轉(zhuǎn)移量前有負(fù)號(hào)*/ dS[eqn] = -0.02*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT; } else { m_dot = 0.02*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn] = -0.02*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;} /*相反,氣相向液相轉(zhuǎn)移則放熱*/ return LAT_HT*m_dot; /*能量大小氣化潛熱與質(zhì)量轉(zhuǎn)移率相乘得能量源項(xiàng)*/ } |
榮譽(yù)版主 (著名寫手)
小木蟲蛀蟲防疫站站長

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +7 | 自由煎餅果子 2026-03-16 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 17/850 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +4 | zxts12 2026-03-21 | 7/350 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|