udfÌí¼ÓÄÚÈÈÔ´ÎÊÌâÇóÖú£¬³ÌÐòÈçÏ£¬µ«µ¼ÈëfluentÖÐ×ÜÊÇÌáʾ³öÏÖ´íÎó
Error: received a fatal signal (Segmentation fault).
Error: received a fatal signal (Segmentation fault).
Error Object: #f
Ï£Íû´ó¼Ò°ïæ¿´¿´ÊDz»ÊdzÌÐòÓÐÎÊÌ⣬ллÁË
/***********************************************************************
UDF that adds energy source term
************************************************************************/
#include "udf.h"
#include "mem.h"
#include "metric.h"
#define Tair 298
#define Ts 500
DEFINE_SOURCE(e_energy,c,t,dS,eqn)
{
real a;
real b;
real Nu;
real qq,ave_q;
real source;
real temd,temh,volume;
real p_avg=0.0;
real vol_tot=0.0;
real v_avg=0.0 ;
begin_c_loop(c,t)
{
volume=C_VOLUME(c,t);
temd=C_P(c,t);
temh=C_U(c,t);
vol_tot+=volume;
p_avg+=temd*volume;
v_avg+=temh*volume;
C_UDMI(c,t,0)=C_U_M0(c,t);
}
end_c_loop(c,t)
p_avg/=vol_tot;
v_avg/=vol_tot;
a=0.8*(Ts-Tair)*pow(3.14*C_K(c,t)*C_R(c,t)*1.225*0.0006*v_avg/2.8,0.5);
b=3.14*0.025*0.025*340*1.732*p_avg;
qq=a*(pow(fabs(1.0/3.0+(temh-C_UDMI(c,t,0))),0.5)-0.577)/b;
Nu=0.664*pow(C_CP(c,t)*C_MU_L(c,t)/C_K_L(c,t),1.0/3)*pow(C_R(c,t)*C_U(c,t)*0.9/C_MU_L(c,t),0.5);
ave_q=Nu*C_K_L(c,t)*(Ts-Tair)*3.14*0.025*0.025/0.9; /*Å£¶ÙÀäÈ´¶¨ÂÉÇóµÃ*/
source=ave_q+qq;
dS[eqn]=0;
return source;
} |