[原创]

vvvms

[watermark]Tin；为固体内节点温度
T ；为固体表面温度
T；（t）为环境温度
DEFINE_PROFILE(T_profile, t, i)
{
real pi=3.14;
face_t tf;
cell_t c0;
Thread *t0;
real time=RP_Get_Real("flow-time");
real T_t=273.15+0.33+4.55*(0.96*sin(2*pi*(time-9)/24)+0.14*sin(4*pi*(time-9)/24));/*insert the function of time-dependent outdoor air temperature*/
real U_t=1.4; /*insert the function of time-dependent outdoor air moisture potential*/
real A1_ext=lambda_1/Dgrid+hq_1, A1_int=lambda_4/Dgrid+hq_4; /*T boundary equation coefficient*/
real C1_ext, C1_int, C2_ext, C2_int;
    begin_f_loop(tf, t)
{
   real xw[ND_ND];
               F_CENTROID(xw,tf,t);
   c0=F_C0(tf, t);
               t0=F_C0_THREAD(tf, t);
   C1_ext=lambda_1*C_UDSI(c0,t0,Temp)/Dgrid+hq_1*T_t; /*T boundary equation constant*/
   C1_int=lambda_4*C_UDSI(c0,t0,Temp)/Dgrid+hq_1*T_indoor; /*T boundary equation constant*/
               
   if(xw[0]<=x1+0.25*Dgrid) F_PROFILE(tf, t, i)=C1_ext/A1_ext;
               else if(xw[0]>=x2-0.25*Dgrid) F_PROFILE(tf, t, i)=C1_int/A1_int;
   
}
    end_f_loop(tf, t)
  
}

vvvms

没人回复啊 ？

vvvms

高手帮帮忙啊?