[求助] udf编写求换热器管壁对流换热系数 已有1人参与

上个帖子发完后，自己找了找错误，新的程序如下，自己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);/*数字为fluid_ID，需改*/ thread_out=Lookup_Thread(domain,20);/*数字为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);/*出口累加面积与对应温度的乘积*/ } 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;/*铜管当量直径*/   ruo_wa=(0.1*(992.2-995.7)*(temp-30))+995.7;/*水的密度*/   lamda_wa=(0.1*(0.635-0.618)*(temp-30))+0.618;/*水的导热系数*/   a_wa=(0.1*(0.000000153-0.000000149)*(temp-30))+0.000000149;/*水的热扩散率*/   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;/*铜管内侧的普朗特数*/   Re_wa=ruo_wa*v_wa*de/mu_wa;/*铜管内侧的雷诺数*/   inter_equ_wa=Re_wa*prandtl_wa*de/L;/*铜管内侧的中间数*/   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);/*铜管内侧的努谢尔特数,吕原丽*/    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

回复此楼