| 4 | 1/1 | 返回列表 |
| 查看: 1367 | 回復(fù): 3 | ||
[求助]
udf編寫求換熱器管壁對流換熱系數(shù) 已有1人參與
|
上個(gè)帖子發(fā)完后,自己找了找錯(cuò)誤,新的程序如下,自己interpreted了一下沒錯(cuò)誤,每個(gè)程序后面有詳盡的解釋,請各位大神再給點(diǎn)意見,也造福小木蟲的其他蟲子們。![]() /************************************************************************* * 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)/*每次迭代都計(jì)算一次*/ { 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);/*定性溫度,按照常熱流壁面去計(jì)算*/ 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;/*水的動(dòng)力粘度*/ nu_wa=(0.1*(0.000000659-0.000000805)*(temp-30))+0.000000805;/*水的運(yùn)動(dòng)粘度*/ alpha_wa=(0.1*(0.000386-0.000305)*(temp-30))+0.000305;/*水的定壓比熱容*/ prandtl_wa=mu_wa*cp_wa/lamda_wa;/*銅管內(nèi)側(cè)的普朗特?cái)?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è)的努謝爾特?cái)?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 |
用戶注銷 (文壇精英)
新蟲 (初入文壇)

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +6 | leaflight 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
[考研] 321求調(diào)劑 +6 | 璞玉~~ 2026-03-25 | 7/350 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考博] 26申博 +3 | 加油沖! 2026-03-26 | 3/150 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|