用UDF传递流量对气体浓度场造成影响？

gjmfreeboy

DEFINE_ADJUST(avg_liuliang_Q, domain)
{

real NV_VEC(P);
face_t f;
Thread *t;
real sum_mass=0.0,sum_P=0.0,f_vec=0.0;

t=Lookup_Thread(domain,WALL1_ID );
begin_f_loop(f,t)
{
f_vec=sqrt(F_U(f,t)*F_U(f,t)+F_V(f,t)*F_V(f,t)+F_W(f,t)*F_W(f,t));

       F_AREA(P,f,t);
        sum_P+=NV_MAG(P);  
   sum_mass+=NV_MAG(P)*f_vec*F_R(f,t);
}
end_f_loop(f,t)
avg_temp5=sum_mass;
}


DEFINE_PROFILE(liuliang_Q,t,position)
{
        face_t f;
        begin_f_loop(f,t)
        {
                F_PROFILE(f,t,position) = avg_temp5;
  }
        end_f_loop(f,t)
}
这是我的传递函数，目的是把出口的流量赋予进口。编译的时候没有出现问题，但是迭代求解的时候发现对已经初始化的气体浓度场影响很大，观测点的气体浓度场不断地上下波动。请问这是什么问题呢？

mingtongseo

做等高手回复

mayo_0

看ADJUST里流量算对了么？
计算方法是错误的。

gjmfreeboy

恩，谢谢。那应该怎么算？请高手指点