 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2001 | 回復(fù): 4 | ||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
bianboshen新蟲(chóng) (初入文壇)
|
[交流]
UDF源項(xiàng)編寫(xiě)中如何調(diào)用壁面溫度
|
|||
|
樓主新人，目前想寫(xiě)一個(gè)液滴的能量源項(xiàng)，但是在源項(xiàng)中不知道如何調(diào)用壁面溫度，代碼如下，可以編譯，但無(wú)法運(yù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=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; } } |
回復(fù)此樓» 收錄本帖的淘帖專(zhuān)輯推薦
 博士科研新路歷程 |
» 猜你喜歡
材料專(zhuān)碩英一數(shù)二306
已經(jīng)有6人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
265求調(diào)劑
已經(jīng)有7人回復(fù)
-
26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線
已經(jīng)有8人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有21人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有20人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有6人回復(fù)
-
085600材料與化工 求調(diào)劑
已經(jīng)有17人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有13人回復(fù)
-
317求調(diào)劑
已經(jīng)有6人回復(fù)
4樓2016-12-28 10:46:57
2樓2016-12-27 16:16:32
bianboshen
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 15
- 帖子: 13
- 在線: 3.6小時(shí)
- 蟲(chóng)號(hào): 2482995
- 注冊(cè): 2013-05-26
- 專(zhuān)業(yè): 流體力學(xué)
|
感謝回帖，我目前要調(diào)用確實(shí)是壁面溫度，將源箱只加載在與壁面相連的邊界層中的cell上，其他cell不不加。這樣是可以實(shí)現(xiàn)的么？ 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
3樓2016-12-27 23:44:14
bianboshen
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 15
- 帖子: 13
- 在線: 3.6小時(shí)
- 蟲(chóng)號(hào): 2482995
- 注冊(cè): 2013-05-26
- 專(zhuān)業(yè): 流體力學(xué)
|
這個(gè)是我重新寫(xiě)的，不知是否已經(jīng)實(shí)現(xiàn)我所說(shuō)的功能，另我這個(gè)源項(xiàng)加進(jìn)去是要模擬液滴與壁面之間的對(duì)流換熱，所以我的想法是給壁面液體設(shè)置一個(gè)正能量源項(xiàng)，給氣體設(shè)置一個(gè)負(fù)的能量源項(xiàng)。 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ā)表 | |
|---|---|---|---|---|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 085601專(zhuān)碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專(zhuān)家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽(tīng)風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 0856專(zhuān)碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|