这算什么现象？

classic1573

问题描述：图中模型为一方形填料塔，中间部分为填料层，入口气体仅看作是空气，想模拟填料层对气流分布的影响，之前把填料层作为多孔介质处理，模拟结果显示流场较为合理，但很多文献将填料层对压力降的影响作为动量方程的源项进行处理，于是按照公式编写了udf添加到填料层计算区域（图中的中间区域），在一个小模型里运行没问题，但是放在这个大模型中运行结果却是这样的，只有一边有气体通过，明显不对啊，什么原因？求解释
程序：
#include "udf.h"
#define C2 464
DEFINE_SOURCE(z_mom_source,c,t,ds,eqn)
{
real z[ND_ND];
real con,source;
C_CENTROID(z,c,t);
if(z[2]>=4.645&&z[2]<8.045)
{
con=C2*0.5*C_R(c,t)*(z[2]-4.645);
source=-con*fabs(C_W(c,t))*C_W(c,t)+C_P(c,t);
ds[eqn]=-2.*con*fabs(C_W(c,t));
return source;
}
}