24小时热门版块排行榜    

查看: 1177  |  回复: 1
【奖励】 本帖被评价1次,作者nba24增加金币 0.6

nba24

新虫 (小有名气)


[资源] 两相流----传质UDF -----

#include "udf.h"
#include "mem.h"                                                    #define q 0.40
#define c 0.22
#define m 0.20
#define INT_YI 0.03
#define UDF_FILENAME "udf_loc_velo"
real AVG_YI;


DEFINE_ADJUST(outlet_avg_yi, d)
{

real sum_yi_A=0.0;
real sum_yi=0.0;
int j=0;
face_t f;
Thread *thread_out=Lookup_Thread(d,4);
  begin_f_loop(f,thread_out)
{
  sum_yi_A+=F_YI(f,thread_out,0);
  j++;      
}
end_f_loop(f,thread_out)

  AVG_YI=sum_yi_A/j;


printf("%f\n",sum_yi_A);/* write out kount to data file */
printf("%f\n",AVG_YI);
printf("%d\n",j);
}


DEFINE_PROFILE(inlet,thread_in,index)
{
    real inlet_yi=0.0;
    face_t inlet_face;
  
    real time=RP_Get_Real("flow-time";   
    begin_f_loop(inlet_face,thread_in)
        {
          if(time==0)
        {
         inlet_yi=INT_YI;
          F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
        }
        
          else
        {
          inlet_yi=AVG_YI;
          F_PROFILE(inlet_face,thread_in,index)=inlet_yi+0.2*0.5;
        }
      }
        end_f_loop(inlet_face,thread_in)   

     printf("%f\n",inlet_yi);                     
     }



大家看看    多提问题
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2012-11-03 10:41   回复  
三星好评  顶一下,感谢分享!
相关版块跳转 我要订阅楼主 nba24 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见