关于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;”
页:
[1]