| 1 | 1/1 | 返回列表 |
| 查看: 1930 | 回復(fù): 0 | |||
llongtan木蟲 (小有名氣)
老木蟲
|
[求助]
求助:fluent的沸騰UDF是否可以應(yīng)用的三維中
|
|
在網(wǎng)上看的沸騰的UDF一般都是在二維的模型上應(yīng)用,請(qǐng)問(wèn)一下他可不可以應(yīng)用到三維模型上。如果不能需要怎么修改。 謝謝了。 下面是網(wǎng)上找到的沸騰UDF: #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 } [ Last edited by llongtan on 2012-6-29 at 09:49 ] |
fluent模擬及問(wèn)題 |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北化求調(diào)劑 +4 | Jsman 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +4 | Shern—- 2026-03-28 | 4/200 |
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|