内容发布更新时间 : 2024/5/31 21:21:51星期一 下面是文章的全部内容请认真阅读。

/* 使用时将此文档的扩展名改为.c 然后编译到fluent当中 每个参数都有说明 用户可以根据需要调整 调用UDM需要首先设置Define User-Defined Memory... 设置为1即可 此处只存储了一组数据*/

#include \

/* Author:wanghao SOUTH CHINA UNIVERSITY of TECHNOLOGY BEEL 2014*/

DEFINE_ON_DEMAND(on_demand_PMVcalc) {

/*定义real数据类型fluent会转译为float数据类型*/ real t_wh; /*空气温度*/

real speed_u_wh; /* u方向速度分量 */ real speed_v_wh; /* v方向速度分量*/ real speed_w_wh; /* w方向速度分量*/

real RH_wh=40; /* 输入相对湿度 水蒸气饱和时输入100 */

real pa_wh; /*物质质量分数C_YI(c,t,0) 后面通过相对湿度计算水蒸气分压力*/

real icl_wh=0.155; /*服装热阻 在计算时单位不是clo 此处将0.3clo换算成0.05㎡K/W*/

real fcl_wh; /*穿衣服人体外表面与裸体人体表面积之比*/ real tcl_wh; /*服装外表面温度*/ real tcl1_wh; real tcl2_wh; real temperary1_wh; real temperary2_wh;

real hc_wh; /*表面换热系数*/

real tr_wh; /*平均辐射温度 可以后期自己定义*/ real M_wh=69.78; /*新陈代谢量*/ real W_wh=0.0; /*机械做功*/

real a_wh; /*heat loss by radiation refer to */

real b_wh; /*heat loss by convection refer to */

real c_wh; /*heat loss diff.through skin refer to */

real d_wh; /*heat loss by sweating(comfort) refer to */

real e_wh; /*latent respiration heat loss refer to */

real f_wh; /*dry respiration heat loss refer to */

real L_wh;

real speed_wh; real PMV_wh;

Domain *d; /* declare domain pointer since it is not passed as an argument to the DEFINE macro */

Thread *t; /* 声明Thread指针 */ cell_t c;

d = Get_Domain(1); /*!!!!!这个要置于所有定义的变量之后 Get the domain using Fluent utility*/

thread_loop_c(t,d) /* Loop over all cell threads in the domain 因为loop是针对每个cell的所以针对cell的计算都要包含到loop里面*/ {

begin_c_loop(c,t) {

t_wh=C_T(c,t)-273.15; /*获取空气温度 ℃*/

speed_u_wh=C_U(c,t); /*获取u方向速度分量 m/s */ speed_v_wh=C_V(c,t); /*获取v方向速度分量 m/s */ speed_w_wh=C_W(c,t); /*获取w方向速度分量 m/s */

tr_wh=t_wh; /*此处把平均辐射温度近似处理为空气温度 ℃ 可以后期自己定义*/

pa_wh=RH_wh*10*exp(16.6536-4030.183/(t_wh+235)); /*计算水蒸气分压力*/

speed_wh=sqrt(pow(speed_u_wh,2.0)+pow(speed_v_wh,2.0)+pow(speed_w_wh,2.0)); /*计算绝对空气相对速度*/

if (icl_wh<0.078)

fcl_wh=1.00+1.290*icl_wh; else

fcl_wh=1.05+0.645*icl_wh; /*判断穿衣服人体外表面与裸体人体表面积之比*/

temperary1_wh=2.38*pow(fabs(tcl1_wh-t_wh-273),0.25); /*自然换热系数*/

temperary2_wh=12.1*pow(speed_wh,0.5); /*强迫换热系数*/ if (temperary1_wh

hc_wh=temperary1_wh; /*判断表面换热系数 并赋值*/

tcl1_wh=t_wh+273+(35.5-t_wh)/(3.5*icl_wh+0.1); /*计算服装表面温度的初始

值*/

tcl2_wh=(308.7-0.028*(M_wh-0)+icl_wh*fcl_wh*hc_wh*(t_wh+273)+icl_wh*0.0000000396*fcl_wh*pow((tr_wh+273),4)-icl_wh*0.0000000396*fcl_wh*pow(tcl1_wh,4))/(1+icl_wh*fcl_wh*hc_wh);/*计算中服装表面温度 K */

while (fabs(tcl1_wh-tcl2_wh)>0.001) /*通过判断差的绝对值进行循环逼近真实值 fabs表示取浮点数据类型的绝对值*/ {

tcl1_wh=tcl2_wh; /*循环过程中把上次的计算结果赋给tcl1_wh做下一次计算*/

temperary1_wh=2.38*pow(fabs(tcl1_wh-t_wh-273),0.25); /*自然换热系数*/

temperary2_wh=12.1*pow(speed_wh,0.5); /*强迫换热系数*/ if (temperary1_wh

hc_wh=temperary1_wh; /*判断表面换热系数 并赋值*/

tcl2_wh=(308.7-0.028*(M_wh-0)+icl_wh*fcl_wh*hc_wh*(t_wh+273)+icl_wh*0.0000000396*fcl_wh*pow((tr_wh+273),4)-icl_wh*0.0000000396*fcl_wh*pow(tcl1_wh,4))/(1+icl_wh*fcl_wh*hc_wh);/*计算中服装表面温度 K */ }

tcl_wh=tcl2_wh-273; /*最终的服装表面温度 ℃ */

a_wh=0.0000000396*fcl_wh*(pow((tcl_wh+273),4.0)-pow((tr_wh+273),4.0)); b_wh=fcl_wh*hc_wh*(tcl_wh-t_wh);

c_wh=0.00305*(5733-6.99*(M_wh-W_wh)-pa_wh);

if((M_wh-W_wh-58.15)>0) d_wh=0.42*(M_wh-W_wh-58.15); else d_wh=0; /*heat loss by sweating(comfort) refer to */ e_wh=0.000017*M_wh*(5867-pa_wh);

f_wh=0.0014*M_wh*(34-t_wh);

L_wh=M_wh-W_wh-(a_wh+b_wh+c_wh+d_wh+e_wh+f_wh);

PMV_wh=(0.303*exp(-0.036*M_wh)+0.028)*L_wh;

C_UDMI(c,t,0)=PMV_wh; /*调用UDM需要首先设置Define User-Defined Memory... 设置为1即可 此处只存储了一组数据 */ }

end_c_loop(c,t) } }