请大家帮我看看为啥我的这个动量源没作用

日照

#include "udf.h"
#define con 2.0
DEFINE_SOURCE(ymom_source,c,t,dS,eqn)
{
  real x[ND_ND];
  real r;
  real source;
begin_c_loop(c,t)
{
F_CENTROID(x,c,t);
r=sqrt(x[1]*x[1]+x[3]*x[3]);
if(x[2]>=-0.25&&x[2]<=0.25&&r<=2)
{ source=-con*r;
  dS[eqn]=con;
}
else
source=dS[eqn]=0;
end_c_loop(c,t)
}
return source;
}
程序能够编译，运行。而且在fluent里面已经选择了Y方向的动量源項使用UDF定义
但是计算后打开云图发现这个动量源根本没加上去
是程序中什么地方有毛病吗？请多多指教

binsce

源项宏不需要单元格的循环，它是自动循环的，我估计是这个原因，也不一定，拙见而已！

flyfish1234

另外，是否因为x=x[0],y=x[1],z=x[2],源项加错地方了呢？猜测而已，不妨一试