 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2002 | 回復: 4 | ||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
[交流]
UDF源項編寫中如何調(diào)用壁面溫度
|
||||
|
樓主新人，目前想寫一個液滴的能量源項，但是在源項中不知道如何調(diào)用壁面溫度，代碼如下，可以編譯，但無法運行。 DEFINE_SOURCE(liq_enrg_src, c, sec_th, dS, eqn) { Thread *mix_th, *pri_th, *t; Domain *domain=Get_Domain(1); int zone_ID=18;/*wall ID*/ real m_dot_d=0; real kd=0; real Cnt=0; real a; real b; real d=0; real vg=0; real Reg=0; real tr=0; real Enrg_Source=0; face_t f; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th,0); t=Lookup_Thread(domain,zone_ID); begin_f_loop(f, t) { Cnt=(1.0-C_VOF(c,pri_th))*C_R(c,sec_th); a=Cnt/C_R(c,pri_th); if(a<0.3) { kd=0.18*sqrt(0.026/(C_R(c,pri_th)*0.008)); } else { kd=pow(a,-0.85)*sqrt(0.026/(C_R(c,pri_th)*0.008)); } m_dot_d=kd*Cnt; dS[eqn]=0; vg=sqrt(pow(C_U(c,pri_th),2)+pow(C_V(c,pri_th),2)+pow(C_W(c,pri_th),2)); Reg=C_R(c,pri_th)*vg*0.008/C_MU_EFF(c,pri_th); d=0.00796*0.026/C_R(c,pri_th)/pow(C_VOF(c,pri_th)*vg,2)*pow(Reg,2/3)*pow(C_R(c,pri_th)/C_R(c,sec_th),-1/3)*pow(C_MU_EFF(c,pri_th)/C_MU_EFF(c,sec_th),2/3); tr=3.14*sqrt(C_R(c,sec_th)*d*d*d/16/0.026); b=pow(18*pow(C_K_L(c,pri_th)*tr,3)*C_R(c,pri_th)*LAT_HT*pow(m_dot_d,5)/pow(C_R(c,sec_th),4)/C_MU_EFF(c,pri_th)/(1.0-C_VOF(c,pri_th))/(F_T(f,t)-T_SAT),0.25); Enrg_Source=(F_T(f,t)-T_SAT)*b*4/0.008; if(C_T(c,mix_th)<T_SAT) { return 0; } else { return Enrg_Source; } } |
回復此樓» 收錄本帖的淘帖專輯推薦
 博士科研新路歷程 |
» 猜你喜歡
一志愿重慶大學085700資源與環(huán)境總分308求調(diào)劑
已經(jīng)有6人回復
-
268求調(diào)劑
已經(jīng)有7人回復
-
302求調(diào)劑
已經(jīng)有11人回復
-
306求調(diào)劑
已經(jīng)有4人回復
-
材料學學碩080502 337求調(diào)劑-一志愿華中科技大學
已經(jīng)有4人回復
-
求調(diào)劑
已經(jīng)有6人回復
-
313求調(diào)劑
已經(jīng)有3人回復
-
南昌大學材料專碩311分求調(diào)劑
已經(jīng)有6人回復
-
316求調(diào)劑
已經(jīng)有6人回復
-
346求調(diào)劑[0856]
已經(jīng)有7人回復
|
感謝回帖，我目前要調(diào)用確實是壁面溫度，將源箱只加載在與壁面相連的邊界層中的cell上，其他cell不不加。這樣是可以實現(xiàn)的么？ 發(fā)自小木蟲Android客戶端 |
3樓2016-12-27 23:44:14
2樓2016-12-27 16:16:32
4樓2016-12-28 10:46:57
|
這個是我重新寫的，不知是否已經(jīng)實現(xiàn)我所說的功能，另我這個源項加進去是要模擬液滴與壁面之間的對流換熱，所以我的想法是給壁面液體設置一個正能量源項，給氣體設置一個負的能量源項。 DEFINE_SOURCE(liq_enrg_src, c, sec_th, dS, eqn) { Thread *mix_th, *pri_th, *t; Domain *domain=Get_Domain(1); int zone_ID=18;/*wall ID*/ real m_dot_d=0; real kd; real Cnt; real a; real Area[ND_ND]; real b; real d; real vg; real Reg; real tr; real Enrg_Source=0; face_t f; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th,0); t=Lookup_Thread(domain,zone_ID); begin_f_loop(f, t) { if(c == F_C0(f,t)) { Cnt=(1.0-C_VOF(c,pri_th))*C_R(c,sec_th); a=Cnt/C_R(c,pri_th); if(a<0.3) { kd=0.18*sqrt(0.026/(C_R(c,pri_th)*0.008)); } else { kd=pow(a,-0.85)*sqrt(0.026/(C_R(c,pri_th)*0.008)); } m_dot_d=kd*Cnt; vg=sqrt(pow(C_U(c,pri_th),2)+pow(C_V(c,pri_th),2)+pow(C_W(c,pri_th),2)); Reg=C_R(c,pri_th)*vg*0.008/C_MU_L(c,pri_th); d=0.00796*0.026/C_R(c,pri_th)/pow(C_VOF(c,pri_th)*vg,2)*pow(Reg,2/3)*pow(C_R(c,pri_th)/C_R(c,sec_th),-1/3)*pow(C_MU_L(c,pri_th)/C_MU_L(c,sec_th),2/3); tr=3.14*sqrt(C_R(c,sec_th)*d*d*d/16/0.026); b=pow(18*pow(C_K_L(c,pri_th)*tr,3)*C_R(c,pri_th)*LAT_HT*pow(m_dot_d,5)/pow(C_R(c,sec_th),4)/C_MU_EFF(c,pri_th)/(1.0-C_VOF(c,pri_th))/(F_T(f,t)-T_SAT),0.25); F_AREA(Area,f,t); Enrg_Source=(F_T(f,t)-T_SAT)*b/C_VOLUME(c,mix_th)/(1-C_VOF(c, pri_th))*NV_MAG(Area); dS[eqn]=0; //Enrg_Source=2; return Enrg_Source; } } end_f_loop(f, t) } |
5樓2016-12-28 13:28:51
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
|
[考研] 354求調(diào)劑 +5 | Tyoumou 2026-03-18 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研]
 295復試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 生物學調(diào)劑招人�。。� +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 0703化學336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|