cavity基础上的幂律流体代码，但是有问题，大家看看问题出在哪里？

hongyan

cavity基础上，将牛顿流体改成幂律流体，代码在附件中。
迭代时出现error=-1.#IND00e+000,无法继续迭代。
在evolution中加入了以下表示幂律流体性质的代码
void evolution()
{
  for(i=1;i<NX;i++)
    for(j=1;j<NY;j++)
        {
                for (k=0;k<Q;k++)
                {
                two1[j]+=(f[j][k]-feq(k,rho[j],u[j]))*e[k][0]*e[k][0];
                two2[j]+=(f[j][k]-feq(k,rho[j],u[j]))*e[k][1]*e[k][1];
                two3[j]+=(f[j][k]-feq(k,rho[j],u[j]))*e[k][0]*e[k][1];
                two4[j]+=(f[j][k]-feq(k,rho[j],u[j]))*e[k][1]*e[k][0];  
                }
                D[j]=two1[j]*two1[j]+two2[j]*two2[j]+two3[j]*two3[j]+two4[j]*two4[j];
                W[j]=1.5*pow(D[j],0.5)/(rho[j]*tau_f[j]);
                niu[j]=m*pow(W[j],(t-1))/rho[j];
                tau_f[j]=3.0*niu[j]+0.5;
        }

[ 本帖最后由 hongyan 于 2014-1-10 12:26 编辑 ]

zhy888

你的松弛时间为什么是一维数组

uesoft

不懂这些。
但c/c++里面出现类似error=-1.#IND00e+000错误一般是因为某个变量没有初始化。

DelvonLiu

我也是这个问题 郁闷。。