24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

逆袭defy

新虫 (初入文坛)

应助: 2 (幼儿园)
帖子: 13
在线: 14.4小时
虫号: 3594110
注册: 2014-12-14
专业: 多相流热物理学

[交流] 利用udf统计出来的速度和直接report出来的速度差很多已有3人参与

我利用udf读取边界的速度，代码如下，
DEFINE_ADJUST(gasification, domain)
{

cell_t c;
      Thread *gas_thread;
      Thread *jinkou_thread;
      Thread *chukou_thread;
      Domain *mixure_domain,*solid_domain;
      real zhiliangfenshu1=0,tiji1=0;
      real zhiliangfenshu2=0,tiji2=0;
real chukousuduzonhe=0;
real chukousuduzonhe1=0;
real cellshuliang=0;

      solid_domain=Get_Domain(3);

gas_thread= Lookup_Thread(solid_domain, CHUKOUID);
chukou_thread = THREAD_T0(gas_thread);

begin_c_loop(c,chukou_thread)
   {
chukouwenduzonghe+=C_T(c,chukou_thread)*C_VOLUME(c,chukou_thread);
chukousuduzonhe+=fabs(C_U(c,chukou_thread)*C_VOLUME(c,chukou_thread));/*注意这里用的是x方向的速度作为出口速度，根据mesh的具体坐标而定*/
chukousuduzonhe1+=fabs(C_U(c,chukou_thread));
cellshuliang+=1;
tiji1+=C_VOLUME(c,chukou_thread);
}
end_c_loop(c,chukou_thread)
   chukousudu1=chukousuduzonhe1/cellshuliang;
chukousudu=chukousuduzonhe/tiji1;
Message ("出口固体速度: %g\n",chukousudu);
Message ("出口固体速度平均: %g\n",chukousudu1);
}

我用两种方法求速度，一个是把速度*体积加到一起然后除以总体积，另外一种是把所有的数量加一起，然后除以cell总数量，两种方法求出来的边界的速度很接近，但是和利用flunet直接report出来的速度差的很多，，但是这两种方法求出来的温度和质量分数和report出来的都是一样的，只有速度和report的不一样
求大侠指点下，已经想了好几天了还没有头绪

回复此楼