求教fluent 中PR方程计算气体混合物的udf编写

huaer_03

以下是我的求密度程序，请高手指点，本人不胜感激！另外，请问fluent中是否可以计算双层for循环？
#define N 7
DEFINE_PROPERTY(density,c,t)
{
  int i,j;
/*中间省略部分程序，程序中其它数据计算正确，但是以下for循环总是迭代错误*/
for(i=0;i<N;i++)
   {
     for(j=0;j<N;j++)
     {
      k[j]=sqrt(a*a[j]);
      r[j]=y*y[j];
      m=m+k[j]*r[j];
      }
    }
     do
     {
       v=V;
       V=R*C_T(c,t)/p+n-m*(v-n)/(p*v*(n+v)+p*n*(v-n));
     }while(fabs(V-v)<0.001);
    rho_mix=Mmix/V*1000;
    return rho_mix;
}