 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1368 | 回復(fù): 3 | ||
[求助]
udf編寫求換熱器管壁對流換熱系數(shù) 已有1人參與
|
上個帖子發(fā)完后，自己找了找錯誤，新的程序如下，自己interpreted了一下沒錯誤，每個程序后面有詳盡的解釋，請各位大神再給點意見，也造福小木蟲的其他蟲子們。 /************************************************************************* * THIS UDF PROGRAM IS PROGRAMMED BY LIU Wenjie ON 11.28.2016 * * IT IS USED FOR THE THERMAL BOUNDARY CONDITIONS OF HEAT EXCHANGER PIPE WALL'S CONVECTIVE HEAT TRANSFER COEFFICIENT（h）**************************************************************************/ #include "udf.h" #define v_wa 0.061 /*銅管（換熱管）入口流速*/ #define L 1.2 /*10銅管長度*/ #define W 0.006 /*銅管入口的寬度*/ #define H 0.008 /*銅管入口高度*/ #define T11 307.5 /*冷水入口處的水溫*/ #define cp_wa 4.174 /*水的定壓比熱容*/ real avg_temp; /************************************************************************* the UDF for average temperature ofoutlet求出口平均溫度 *************************************************************************/ DEFINE_ADJUST(lwj_adjust,d)/*每次迭代都計算一次*/ { real NV_VEC(A); real sum_T_A=0.0; real sum_A=0.0; Thread *thread_out;/*用于宏在編譯過程中的延伸*/ face_t f; cell_t cell; Domain *domain; domain=Get_Domain(11);/*數(shù)字為fluid_ID，需改*/ thread_out=Lookup_Thread(domain,20);/*數(shù)字為zone_ID，需改*/ begin_f_loop(f,thread_out) { F_AREA(A,f,thread_out); sum_A+=NV_MAG(A);/*求得出口的累加面積*/ sum_T_A+=NV_MAG(A)*F_T(f,thread_out);/*出口累加面積與對應(yīng)溫度的乘積*/ } end_f_loop(f,thread_out) avg_temp=sum_T_A/sum_A;/*求得出口的平均溫度*/ } /************************************************************************* the UDF for convective heat transfer coefficient of inner HEX wall *************************************************************************/ DEFINE_PROFILE(lwj_h,thread,index) { real temp; real prandtl_wa,Re_wa,NU_wa,inter_equ_wa; real U,S,de; real t; real h_wa;/*h*/ real ruo_wa,lamda_wa,a_wa,mu_wa,nu_wa,alpha_wa; t=RP_Get_Real("flow-time"  ;temp=0.5*(T11+avg_temp);/*定性溫度，按照常熱流壁面去計算*/ U=2*(W+H);/*銅管周長*/ S=W*H;/*銅管斷面積*/ de=4*S/U;/*銅管當(dāng)量直徑*/ ruo_wa=(0.1*(992.2-995.7)*(temp-30))+995.7;/*水的密度*/ lamda_wa=(0.1*(0.635-0.618)*(temp-30))+0.618;/*水的導(dǎo)熱系數(shù)*/ a_wa=(0.1*(0.000000153-0.000000149)*(temp-30))+0.000000149;/*水的熱擴(kuò)散率*/ mu_wa=(0.1*(0.0006533-0.0008015)*(temp-30))+0.0008015;/*水的動力粘度*/ nu_wa=(0.1*(0.000000659-0.000000805)*(temp-30))+0.000000805;/*水的運動粘度*/ alpha_wa=(0.1*(0.000386-0.000305)*(temp-30))+0.000305;/*水的定壓比熱容*/ prandtl_wa=mu_wa*cp_wa/lamda_wa;/*銅管內(nèi)側(cè)的普朗特數(shù)*/ Re_wa=ruo_wa*v_wa*de/mu_wa;/*銅管內(nèi)側(cè)的雷諾數(shù)*/ inter_equ_wa=Re_wa*prandtl_wa*de/L;/*銅管內(nèi)側(cè)的中間數(shù)*/ if (Re_wa<2300.0) { NU_wa=3.66; if(inter_equ_wa>10) NU_wa=1.86*pow(inter_equ_wa,1.0/3.0);/*銅管內(nèi)側(cè)的努謝爾特數(shù),呂原麗*/ else NU_wa=3.66+0.0668*inter_equ_wa/(1+0.04*pow(inter_equ_wa,2.0/3.0));/*NU,呂*/ } else if(Re_wa<10000.0) NU_wa=0.16*pow(Re_wa,2.0/3.0)*pow(prandtl_wa,1.0/3.0)*pow(1+de/L,2.0/3.0);/*NU*/ else NU_wa=0.023*pow(Re_wa,0.8)*pow(prandtl_wa,0.4);/*NU,呂*/ h_wa=NU_wa*lamda_wa/de; begin_f_loop(f,thread) { if(temp>C_T(cell,thread)) h_wa=0; F_PROFILE(f,thread,index)=h_wa; } end_f_loop(f,thread) } @wuming524 |
回復(fù)此樓» 猜你喜歡
352分-085602-一志愿985
已經(jīng)有3人回復(fù)
-
291求調(diào)劑
已經(jīng)有4人回復(fù)
-
330分求調(diào)劑
已經(jīng)有4人回復(fù)
-
309求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856求調(diào)劑
已經(jīng)有13人回復(fù)
-
0703化學(xué)
已經(jīng)有11人回復(fù)
-
329求調(diào)劑
已經(jīng)有7人回復(fù)
-
321求調(diào)劑
已經(jīng)有8人回復(fù)
-
0856材料化工調(diào)劑 總分330
已經(jīng)有12人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有4人回復(fù)
2樓2016-12-05 21:04:15
匿名
用戶注銷 (文壇精英)
- 應(yīng)助: 164 (高中生)
- 金幣: 9550.6
- 散金: 34344
- 紅花: 361
- 沙發(fā): 31
- 帖子: 20580
- 在線: 3633.3小時
- 蟲號: 0
- 注冊: 2011-12-15
- 專業(yè): 建筑環(huán)境與結(jié)構(gòu)工程
3樓2016-12-07 09:11:15
aaalancelot
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 85.5
- 帖子: 20
- 在線: 6.4小時
- 蟲號: 5295150
- 注冊: 2016-11-29
- 性別: GG
- 專業(yè): 流體力學(xué)
4樓2016-12-07 10:09:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 溪柳柳自搖 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+7
|
nnnnnnn5 2026-03-25 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +7 | 麻辣魷魚 2026-03-27 | 8/400 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個楠 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[有機(jī)交流]
有機(jī)合成求助
 20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|