关于power law udf的菜鸟问题

乌江挽歌

大家好！
我有做高分子溶液的理论研究，最近涉及到非牛顿流体的非经典模型，临时看教程用udf编了个非常简单的方程，但是在运算的时候怎么都不收敛。
程序如下：
#include "udf.h"
DEFINE_PROPERTY(cell_viscosity, cell, thread)
{
        real visco;
        real R;
        R = C_STRAIN_RATE_MAG(cell,thread);
        visco = 0.14*2.8284/pow(R,-0.56)-0.0018284;
        return visco;
        }

我试过所有的参数发现只有把“-0.56”改成正数才能收敛，但是我本来做的就是假塑性流体，并且用fluent自带的non-newtonian-power law模型（-0.2的flow index）就可以收敛.

所以也许是程序出了问题，希望有高手可以帮我看看
先谢谢啦

乌江挽歌

更正一下 中间一行打错了 是“visco = 0.14*2.8284/pow(R,0.56)-0.0018284;”