利用udf统计出来的速度和直接report出来的速度差很多

terrorblade1992

我利用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的不一样
求大侠指点下，已经想了好几天了还没有头绪


我的联系方式qq 1032993302，望不吝赐教