关于UDF的问题，想求场协同角，请大家指点点

hitwuxing

我想用UDF求解场协同角（速度矢量与温度梯度的夹角），以下是我的程序，顺便交代一下，我是想用DEFINE_ADJUST和DEFINE_EXECUTE_AT_END两个宏进行比较一下，其中有一部分程序我是为了显示一些变量，然后自己检验一下计算结果合不合理。
DEFINE_EXECUTE_AT_END(cooperation_coff)
{
   real i=0,cos,coef;
   real T,T0,T1,T2,u,v,w,dT,dT_V,V;
   FILE *fp3=fopen(UDF_FILE_3, "w");
   int id=2;
   Domain *d=Get_Domain(1);
   Thread *t=Lookup_Thread(d,id);
   thread_loop_c(t,d)
     {
      cell_t c;
      if(THREAD_ID(t)==2)
         begin_c_loop(c,t)
          {
            T=C_T(c,t);
            T0=C_T_RG(c,t)[0];
            T1=C_T_RG(c,t)[1];
            T2=C_T_RG(c,t)[2];
            u=C_U(c,t);
            v=C_V(c,t);
            w=C_W(c,t);
            ND_SET(u,v,w,C_U(c,t),C_V(c,t),C_W(c,t));
            dT_V=T0*u+T1*v+T2*w;
            dT=sqrt(T0*T0+T1*T1+T2*T2);
            V=sqrt(u*u+v*v+w*w);
            cos=dT_V/(dT*V);
            coef=acos(cos);
            i++;
            fprintf(fp3,"%d %e %e %e %e %e %e %e %e %e %e %e\n",i,T0,T1,T2,u,v,w,dT,dT_V,V,cos,coef);
         }
       end_c_loop(c,t)
     }
  fclose(fp3);
}
DEFINE_ADJUST(cooperation_coff_1,d)
  {
    real i=0,cos,coef;
    real T,T0,T1,T2,u,v,w,dT,dT_V,V;
    FILE *fp4=fopen(UDF_FILE_4, "w");
    Thread *t;
    cell_t c;
    thread_loop_c(t,d)
       {
        if(FLUID_THREAD_P(t))
          begin_c_loop(c,t)
           {
             T=C_T(c,t);
             T0=C_T_G(c,t)[0];
             T1=C_T_G(c,t)[1];
             T2=C_T_G(c,t)[2];
             u=C_U(c,t);
             v=C_V(c,t);
             w=C_W(c,t);
             ND_SET(u,v,w,C_U(c,t),C_V(c,t),C_W(c,t));
             dT_V=T0*u+T1*v+T2*w;
             dT=sqrt(T0*T0+T1*T1+T2*T2);
             V=sqrt(u*u+v*v+w*w);
             cos=dT_V/(dT*V);
             coef=acos(cos);
             i++;
             fprintf(fp4,"%d %e %e %e %e %e %e %e %e %e %e %e\n",i,T0,T1,T2,T,u,v,w,dT,dT_V,V,cos);
             printf("%d %e %e %e %e %e %e %e %e %e %e %e\n",i,T0,T1,T2,u,v,w,dT,dT_V,V,cos,coef);
           }
         end_c_loop(c,t)
       }
  fclose(fp4);
}

飞翔鸟

买噶，还有人搞场协同这个东东，真没有意义啊。

hitwuxing

哥们不要这么说啊，中国的学术研究有多少不是从外国文献翻译过来的啊？好歹我还用个中国人搞出来的理论，虽然这个理论的确有点扯淡，对你来说没有意义，对我就不一样了，关系到我能不能毕业的问题！！！

sky321hit

没必要用UDF计算场协同角，用Custom Filed Functions就可以定义。
一、在Fluent空白窗口中输入
solve/set> expert
use conservative form of energy equation? [yes]
use alternate formulation for wall temperatures? [no]
Save cell residuals for post-processing? [no] yes
Keep temporary solver memory from being freed? [no] yes
Allow selection of all applicable discretization schemes? [no]
二、先迭代计算一两步以保证内存中储存了温度梯度
三、Define->Custom Filed Functions->
在对话框中找到Filed Functions下面的Temperature和Reconstructions dT/dx以及Velocity
和Reconstruction dX-Velocity/dx等等量，组成温度梯度矢量和速度梯度矢量的点积。
剩下的后处理我不说了，你应该都懂的吧。
Good Luck！

sky321hit

还有，你问的是场协同角，而不是（火积）耗散极值原理导出的场协同方程组的数值求解。你解决场协同方程源项和逆温度方程源项的处理了吗？解决了的话，麻烦通知一声哈，谢谢！

jsking

如何进行点积计算？

chenhuan1115

原帖由 sky321hit 于 2011-3-11 15:31 发表

                               
登录/注册后可看大图

没必要用UDF计算场协同角，用Custom Filed Functions就可以定义。
一、在Fluent空白窗口中输入
solve/set> expert
use conservative form of energy equation? [yes]
use alternate formulation for wall tempe ...

我试了下，但显示的好像有错误，肯定是哪设置错啦

zzu340

这个方法好，我喜欢。十分的感谢您！

chenhuan1115

原帖由 zzu340 于 2011-4-20 14:03 发表

                               
登录/注册后可看大图

这个方法好，我喜欢。十分的感谢您！

你试了吗可以做吗

tomato123

你好，请教下，我在Filed Functions没有找到Reconstructions dT/dx和Reconstruction dX-Velocity/dx，还有，协同角的公式中速度矢量与温度梯度的点积该怎么表示？这个问题困扰了我很久，劳烦阁下指点迷津，谢谢！

zersong

此贴甚好，必须收藏。

sky321hit

solve/set/expert回车

落月

您好，我也需要求取温度梯度和速度梯度场。按您的方法，收敛以后进行了如上设置。然后先迭代计算一两步以保证内存中储存了温度梯度
。在这里发现出现了问题。本来已经收敛的energy，残差曲线突然变大，不知道怎么回事。

落月

您好，我也需要求取温度梯度和速度梯度场。按您的方法，收敛以后进行了如上设置。然后先迭代计算一两步以保证内存中储存了温度梯度
。在这里发现出现了问题。本来已经收敛的energy，残差曲线突然变大，不知道怎么回事。

yaojie1986

我也为这个问题苦恼