UDF错误-----UDF编湍流模型，有重谢！！！！！111

shangdong

/* Set the turbulent viscosity */
  thread_loop_c (t, domain)
  {   
    if (FLUID_THREAD_P(t))
    {
      begin_c_loop(c,t)
      {
      real RO=C_R(c,t);
      real K=C_UDSI(c,t,TKE);
      real W=C_UDSI(c,t,TDR);
   real DUDX=C_DUDX(c,t);
   real DUDY=C_DUDY(c,t);
            real DVDX=C_DVDX(c,t);
            real DVDY=C_DVDY(c,t);
            real vort_mag=fabs(DVDX-DUDY);
      real MU_T;
      real MU_T1=alphx_(c,t)*RO*K/W;
      real MU_T2=0.31*RO*K/vort_mag/BLF_2(c,t);
      if(MU_T1<MU_T2)
        MU_T=MU_T1;
      else
        MU_T=MU_T2;
    错误处：  C_MU_T(c,t)=MU_T;
      }
   end_c_loop(c,t)
    }
  }
}
这是本人用UDF编的湍流模型中的一段，编译时出现以下错误：
Deleted old libudf\ntx86\2ddp\libudf.dll
已复制 1 个文件。
已复制 1 个文件。
(system "copy E:\Fluent.Inc\fluent6.2.16\src\makefile_nt.udf libudf\ntx86\2ddp\makefile")已复制 1 个文件。
0
(chdir "libudf")()
(chdir "ntx86\2ddp")()
SST.c
..\..\src\SST.c(291) : error C2106: ';='; : left operand must be l-value
Done.
程序中以标出错误处：错误处：C_MU_T(c,t)=MU_T;
此段程序是用来定义湍流模型中的湍流粘性的。
不懂请大侠指教，如成功，有重谢--------愿与他分享我便得这个湍流模型UDF程序。

Joywe

好像MU_T没有初始化啊，呵呵..........
是不是可以重谢啦.....

shangdong

楼上那位兄弟能不能说清楚一点，我还是不明白？什么初始化MU_T???

Joywe

     real MU_T;
     real MU_T1=alphx_(c,t)*RO*K/W;
     real MU_T2=0.31*RO*K/vort_mag/BLF_2(c,t);
     if(MU_T1<MU_T2)
       MU_T=MU_T1;
     else
       MU_T=MU_T2;
   错误处：  C_MU_T(c,t)=MU_T;
仔细看看，MU_T是多少啊？

shangdong

还是不明白？？？？？？？？？？？？？？？
真的！！！！！！！！！！

Joywe

完了，一看就没有诚意贴出来

Joywe

俺看错了，呵呵。。。。。。贴出来讨论一下也好啊
我的邮箱fluent_udf@126.com

hp270388026

[这个贴子最后由hp270388026在 2005/10/26 09:04pm 第 2 次编辑]

简单：
打个比方：
int a=1;//正确
int a=f(x);//正确
f(x)=1;  //wrong

看你怎么谢我。
给我宣传一下软件
http://219.232.54.3/cgi-bin/LB5000/topic.cgi?forum=51&topic=78&show=0

shangdong

换一个，这是Stanford大学UDF编译ke湍流模型的错误：
DEFINE_ADJUST(turb_adjust, domain)
{
  Thread *t;
  cell_t c;
/* Set the turbulent viscosity */
  thread_loop_c (t, domain)
    if (FLUID_THREAD_P(t))
    {
      begin_c_loop(c,t)
      {
出错处：C_MU_T(c,t) = C_R(c,t)*C_MU*f_mu(c,t)*sqrt(C_UDSI(c,t,TKE))/C_UDSI(c,t,TDR);
      }
  end_c_loop(c,t)
    }
错误显示：
已复制 1 个文件。
已复制 1 个文件。
(system "copy E:\Fluent.Inc\fluent6.2.16\src\makefile_nt.udf libudf\ntx86\2ddp\makefile")已复制 1 个文件。
0
(chdir "libudf")()
(chdir "ntx86\2ddp")()
kemodel.c
..\..\src\kemodel.c(81) : error C2106: ';='; : left operand must be l-value
Done.

附上Stanford的UDF的PPT，以及我复制的程序。
这个错误还是这样子的，好像不是你说的那样，请众位高手下下来，到自己的机器上调一下，看一下是不是这个错误，

shangdong

传上去这是资料

hp270388026

你去查查这个宏的定义，我查过了。
别管他什么狗屁美国大学。

hp270388026

DEFINE_PRANDTL_K(user_pr_k, c, t)
{
real pr_k, alpha;
real mu = C_MU_L(c,t);
real mu_t = C_MU_T(c,t);
alpha = rng_alpha(1., mu + mu_t, mu);
pr_k = mu_t/((mu+mu_t)*alpha-mu);
return pr_k;
}
DEFINE_PRANDTL_D(user_pr_d, c, t)
{
real pr_d, alpha;
real mu = C_MU_L(c,t);
real mu_t = C_MU_T(c,t);
alpha = rng_alpha(1., mu + mu_t, mu);
pr_d = mu_t/((mu+mu_t)*alpha-mu);
return pr_d;
}
看看比如人的

hp270388026

DEFINE_PRANDTL_K(user_pr_k, c, t)
{
real pr_k, alpha;
real mu = C_MU_L(c,t);
real mu_t = C_MU_T(c,t);
alpha = rng_alpha(1., mu + mu_t, mu);
pr_k = mu_t/((mu+mu_t)*alpha-mu);
return pr_k;
}
DEFINE_PRANDTL_D(user_pr_d, c, t)
{
real pr_d, alpha;
real mu = C_MU_L(c,t);
real mu_t = C_MU_T(c,t);
alpha = rng_alpha(1., mu + mu_t, mu);
pr_d = mu_t/((mu+mu_t)*alpha-mu);
return pr_d;
}

shangdong

高手先受累调一下我的那个程序，有什么改正再告诉我吧，谢谢！！！！！！！！

shangdong

这位大哥说的好像与我讲的不是一回事
利用ADJUST宏命令来在每一步迭代的开始重新计算C_MU_T
难道不能直接修改？？？