fluent-计算pmv、ppd、dr值自定义函数代码(udf)

/*使用时将此文档的扩展名改为.c然后编译或解释到fluent当中每个参数都有说明用户可以根据需要调整调用UDM需要首先设置NumberofUser-DefinedMemoryLocations(注意不是存储在Node节点)设置为3,此处存储了三组数据，分别为PMV,PPD,DR*//*userdefinememory0变量为PMV*//*userdefinememory1变量为PPD*//*userdefinememory2变量为DR*/#includeudf.h/*Author:wanghaoSOUTHCHINAUNIVERSITYofTECHNOLOGYBEEL2014*//*Author:zhanghongliangDALIANUNIVERSITYofTECHNOLOGYFIE2019*//*增加计算DR值(吹风感系数)和PPD值代码2019*/DEFINE_ON_DEMAND(PMV_PPD_DR_calc){/*定义real数据类型fluent会转译为real数据类型*/#if!RP_HOSTrealt_wh;/*空气温度,℃*/realspeed_u_wh;/*u方向速度分量,m/s*/realspeed_v_wh;/*v方向速度分量,m/s*/realspeed_w_wh;/*w方向速度分量,m/s*/realspeed_wh;/*三个方向的和速度,m/s*/realre_zhhl;/*雷诺数*/realrho_zhhl;/*空气密度*/realmu_zhhl;/*动力粘度系数*/realdiam_zhhl=3.0;/*房间特征长度,置换通风取房间高度,根据实际情况修改*/realtu_zhhl;/*湍动强度,10%~60%之间,混合通风默认40%,置换通风默认20%*/realRH_wh=40.0;/*输入相对湿度水蒸气饱和时输入100*/realpa_wh;/*水蒸气分压力后面通过相对湿度计算水蒸气分压力,pa*/realicl_wh=0.155;/*服装热阻在计算时单位不是clo此处将0.3clo换算成0.05㎡K/W*/realfcl_wh;/*穿衣服人体外表面与裸体人体表面积之比*/realtcl_wh;/*服装外表面温度,℃*/realtcl1_wh;realtcl2_wh;realtemperary1_wh;realtemperary2_wh;realhc_wh;/*表面换热系数,w/m2K*/realtr_wh;/*平均辐射温度可以后期自己定义,℃*/realM_wh=69.78;/*新陈代谢量,W*/realW_wh=0.0;/*机械做功,W*/reala_wh;/*heatlossbyradiationrefertoENISO7730:2005page24*/realb_wh;/*heatlossbyconvectionrefertoENISO7730:2005page24*/realc_wh;/*heatlossdiff.throughskinrefertoENISO7730:2005page24*/reald_wh;/*heatlossbysweating(comfort)refertoENISO7730:2005page24*/reale_wh;/*latentrespirationheatlossrefertoENISO7730:2005page24*/realf_wh;/*dryrespirationheatlossrefertoENISO7730:2005page24*/realL_wh;realPMV_wh;/*计算所得PMV值*/realPPD_zhhl;/*计算所得PPD值*/realDR_zhhl;/*计算所得DR(吹风感系数）值*/Domain*d;/*declaredomainpointersinceitisnotpassedasanargumenttotheDEFINEmacro*/Thread*t;/*声明Thread指针*/cell_tc;/*声明cell_t变量，用于循环宏*/d=Get_Domain(1);/*!!!!!这个要置于所有定义的变量之后GetthedomainusingFluentutility,单相流取值为1*/thread_loop_c(t,d)/*Loopoverallcellthreadsinthedomain因为loop是针对每个cell的所以针对cell的计算都要包含到loop里面*/{begin_c_loop(c,t){/*计算PMV值--------------------------------------------------------------------------------------------*/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.0*exp(16.6536-4030.183/(t_wh+235.0));/*计算水蒸气分压力*/speed_wh=sqrt(pow(speed_u_wh,2.0)+pow(speed_v_wh,2.0)+pow(speed_w_wh,2.0));/*计算绝对空气相对速度*/if(icl_wh0.078)fcl_wh=1.0+1.290*icl_wh;elsefcl_wh=1.05+0.645*icl_wh;/*判断穿衣服人体外表面与裸体人体表面积之比*/tcl1_wh=t_wh+273.0+(35.5-t_wh)/(3.5*icl_wh+0.1);/*计算服装表面温度的初始值*/tcl2_wh=(308.7-0.028*(M_wh-W_wh)+icl_wh*fcl_wh*hc_wh*(t_wh+273.0)+icl_wh*0.0000000396*fcl_wh*pow((tr_wh+273.0),4.0)-icl_wh*0.0000000396*fcl_wh*pow(tcl1_wh,4.0))/(1.0+icl_wh*fcl_wh*hc_wh);/*计算中服装表面温度K*/temperary1_wh=2.38*pow(fabs(tcl2_wh-t_wh-273.0),0.25);/*自然换热系数*/temperary2_wh=12.1*pow(speed_wh,0.5);/*强迫换热系数*/if(temperary1_whtemperary2_wh)hc_wh=temperary2_wh;elsehc_wh=temperary1_wh;/*判断表面换热系数并赋值*/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),0.25);/*自然换热系数*/temperary2_wh=12.1*pow(speed_wh,0.5);/*强迫换热系数*/if(temperary1_whtemperary2_wh)hc_wh=temperary2_wh;elsehc_wh=temperary1_wh;/*判断表面换热系数并赋值*/tcl2_wh=(308.7-0.028*(M_wh-0.0)+icl_wh*fcl_wh*hc_wh*(t_wh+273.0)+icl_wh*0.0000000396*fcl_wh*pow((tr_wh+273.0),4.0)-icl_wh*0.0000000396*fcl_wh*pow(tcl1_wh,4.0))/(1.0+icl_wh*fcl_wh*hc_wh);/*计算中服装表面温度K*/}tcl_wh=tcl2_wh-273.0;/*最终的服装表面温度℃*/a_wh=0.0000000396*fcl_wh*(pow((tcl_wh+273.0),4.0)-pow((tr_wh+273),4.0));b_wh=fcl_wh*hc_wh*(tcl_wh-t_wh);c_wh=0.00305*(5733.0-6.99*(M_wh-W_wh)-pa_wh);if((M_wh-W_wh-58.15)0.0)d_wh=0.42*(M_wh-W_wh-58.15);elsed_wh=0.0;/*heatlossbysweating(comfort)refertoENISO7730:2005page24*/e_wh=0.000017*M_wh*(5867.0-pa_wh);f_wh=0.0014*M_wh*(34.0-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;/*计算PPD值--------------------------------------------------------------------------------------------*/PPD_zhhl=100.0-95.0*exp(-0.03353*pow(PMV_wh,4.0)-0.2179*pow(PMV_wh,2.0));/*计算DR值---------------------------------------------------------------------------------------------*/if(speed_wh0.05)speed_wh=0.05;rho_zhhl=C_R(c,t);mu_zhhl=C_MU_L(c,t);re_zhhl=rho_zhhl*speed_wh*diam_zhhl/mu_zhhl;/*计算雷诺数*///tu_zhhl=0.16*pow(re_zhhl,-0.125);tu_zhhl=20.0;DR_zhhl=(34.0-t_wh)*pow((speed_wh-0.05),0.62)*(0.37*speed_wh*tu_zhhl+3.14);if(DR_zhhl100.0)DR_zhhl=100.0;/*UDM设置---------------------------------------------------------------------------------------------*//*调用UDM需要首先设置NumberofUser-DefinedMemoryLocations，该值设置为3*/C_UDMI(c,t,0)=PMV_wh;C_UDMI(c,t,1)=PPD_zhhl;C_UDMI(c,t,2)=DR_zhhl;}end_c_loop(c,t)}#endif}