|
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
#include "udf.h"
#define WALL1_ID 5
real avg_temp=285.0;
DEFINE_ADJUST(wall1_avg_temp, domain)
{
real NV_VEC(A);
real sum_T_A=0.0,sum_A=0.0;
face_t f;
Thread*f_thread;
f_thread = Lookup_Thread(domain, WALL1_ID);
begin_f_loop(f, f_thread)
{
F_AREA(A,f,f_thread);
sum_A+=NV_MAG(A);
sum_T_A+=NV_MAG(A)*F_T(f, f_thread);
}
end_f_loop(f, f_thread)
avg_temp=sum_T_A/sum_A;
}
DEFINE_PROFILE(Wall2_T,t,position)
{
face_t f;
FILE *fp;
int i;
float temp1;
fp=fopen("t.txt","r");
begin_f_loop(f,t)
{
fscanf(fp,"%f",&temp1);
F_PROFILE(f,t,position) =avg_temp-temp1;
}
end_f_loop(f,t)
fclose(fp);
} |
|