udf动量源求助

cici0725

在管道中设置一段类似于泵的功能的边界条件，有文献说用动量源项，不用进出口边界条件。
本人编的动量源项只在源项附近有流体，整个计算域都没有流体，如图所示。
还有用动量源后，残差曲线的连续性方程精度到10-2次方就不下降了，流场不收敛，请求帮助！！
这是我自己编的动量源程序，不知是否能实现流体的循环。多谢帮助！
#include "udf.h"
#define CON 20.0
DEFINE_SOURCE(mou_source,cell,thread,ds,equ)
{
real vel=0.45;
real x[ND_ND];
real y;
real z;
real source=0;
C_CENTROID(x,cell,thread);
y=x[1];
z=x[2];
if((((x[0]+0.0305621)*(x[0]+0.0305621)+(y-0.1733262)*(y-0.1733262))<=0.000049) && (z<=0.24 && z>=0))   
source=CON*vel;
ds[equ]=0;
return source;
}