udf比热宏：DEFINE_SPECIFIC_HEAT的用法

yc9316

Ansys 12的udf有个自定义比热的改进，DEFINE_SPECIFIC_HEAT（name，T，Tref，h，yi），可是用这个宏自定义的过程中，我发现没有办法按不同区域的方法定义cp，比如：
DEFINE_SPECIFIC_HEAT(proof,T,Tref,h,yi)
{cell_t c;
Thread *t;
Domain *d;
real cp,x[ND_ND];
d=Get_Domain(1);

thread_loop_c(t,d)
        {
          begin_c_loop(c,t)
          {C_CENTROID(x,c,t);
                  if(x[1]<0.3)
                  cp=1000;
                  else
                  cp=2000;
          }
          end_c_loop(c,t)
        }
return cp;
}
返回的值只有1000或2000两种答案。
当网格总高度y<0.3时，整个区域的cp=1000;
当网格总高度y>0.3时，整个区域的cp=2000;
那么C_CENTROID(x,c,t);
          if(x[1]<0.3)
没有按c，t循环。不知道是怎么回事？望高人指点。

yc9316

比热的udf难道不能按区域划分么？求助各位啊。

gearboy

原帖由 yc9316 于 2010-10-25 02:01 发表

                               
登录/注册后可看大图

比热的udf难道不能按区域划分么？求助各位啊。

明显不可以。传入的参数是温度T，Cp应该写成T的函数，而不是位置的函数。除非先对整个c循环，然后通过温度相等来判断以获取当前的c号码，但是这样存在风险，如果区域中有两个cell温度相等，那就结果不对了。还是直接将Cp写成温度的函数比较好。

begin_c_loop(c,t)
          {
if(C_T(c,t)==T) then C_CENTROID(x,c,t);


          }
          end_c_loop(c,t)

yc9316

真是感谢gearboy，帮我答疑解惑好多次了。今天才明白原来DEFINE宏后面（）里的参数类型决定了它对于数据的传递方式。
但是对于一个计算区域内有不同的比热表达式，对于比热宏的定义还是很困难表达，只能在找关于T，h或者Yi的表达式了，不知道这样理解对不对？

swimming02

没试过，如果不行的话可以分成两个fluid zone，分别将比热设定为1000和2000

yc9316

Define material界面中没有按不同zone的物性参数进行定义，不知道您说的是哪个定义方式？

swimming02

你可以人为的将计算域分为两个区

yc9316

哦，没接触过这个功能，不知道是哪个操作命令

lovewb

你讲的画区是指，先算一个区在一定比热下， 然后打开别一个区算别一个比热下？？？

swimming02

我是指在gambit中事先定义两个不同的流体区域，然后在fluent中在这两个区的流体参数分别进行设置

lovewb

哦，明白

yc9316

明白了，谢谢大家

jntan

原帖由 yc9316 于 2010-10-24 10:13 发表

                               
登录/注册后可看大图

Ansys 12的udf有个自定义比热的改进，DEFINE_SPECIFIC_HEAT（name，T，Tref，h，yi），可是用这个宏自定义的过程中，我发现没有办法按不同区域的方法定义cp，比如：
DEFINE_SPECIFIC_HEAT(proof,T,Tref,h,yi)
{cel ...

你这个UDF整体思路是对的，只是你还没有定义t 和 d的关系。在你开始loop之前，加上一个
t = Lookup_Thread(d, Zone_ID[x])，具体指出一个cell thread，你接下来的loop才能用。