24小时热门版块排行榜

返回列表

溪柳柳自摇

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 2748
红花: 2
帖子: 223
在线: 40.2小时
虫号: 1269257
注册: 2011-04-18
专业: 建筑环境与结构工程

[求助] 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

回复此楼