[求助] 用fluent udf向一个圆盘区域添加动量源项，可是一直加不上，恳求大佬帮忙看一下？

如题， 拟用fluent udf向一个圆盘区域添加动量源项模拟一个叶片旋转，这里先在定常状态下测试往扇形区域加体积力，可是一直加不上，也不知道是不是代码有问题，烦请有会的看一下？ udf代码如下 使用的是cell zone condition里添加动量源项，compile没有问题 #include "udf.h" DEFINE_SOURCE(z_source, cell, thread, dS, eqn) {         real source;         real x0, y0, z0, z1, r, fai1,r1;            real x[3];                                           C_CENTROID(x, cell, thread);                 x0 = x[0];                 y0 = x[1];                 z0 = x[2];                 r = sqrt(x0 * x0 + y0 * y0);                 fai1 = atan(y0 / x0);                 z1 = z0;                 r1 = r / 1000;                 if (z1 > -200. && z1 < 200. && r>500. && r < 8500.  && fai1>0 &&fai1<1 )           //对应的扇形区域                 {                         source = -40* r1 * r1;                         dS[eqn] = 0;                                            }                 else                         {                                                     source = dS[eqn] = 0;                                                  }                         return source; }

回复此楼