求助——如何利用udf添加能量源相

fpw0412

现在遇到这个的一个问题，我想在基台上方加一个圆形的温度源相，要求圆内的温度在3000K左右，圆外的温度fluent默认或300K也可。需要自定义源相，其编程如下：
/**************************************************************/
/* UDF for specifying an heat source term */
/**************************************************************/
#include "udf.h"
#define C1 0.028
#define C2 0.030
#define PI 3.14159
#define T_SAT 300.

DEFINE_SOURCE(energy_source, cell, thread, dS, eqn)
{
real r[ND_ND];
real d, temp,source,x,y,co;
temp=C_T(cell,thread);

x=r[0];
y=r[1];
d= sqrt(x*x+(y-C1)*(y-C1));
co=1.+cos(PI*d/C2);
if(d<=(C2/2))
  {
           source =C_CP(cell,thread)*C_R(cell,thread)*(temp-T_SAT)*co;
   dS[eqn]=C_CP(cell,thread)*C_R(cell,thread)*co;
  }
   else
    source=0;
return source;
}
导入到fluent里面，显示正确，却没有计算结果， 我很少苦恼，不知道有哪位大侠帮我指导下啊，某将感激不尽！

[ 本帖最后由 fpw0412 于 2010-8-20 10:56 编辑 ]

fpw0412

另外，能量源相对谁求导啊？望大侠们多多指导啊，急！！

fpw0412

上面编的不是很清楚
热源方程如下
H=hs/2*[1+cos(pi*d/c0)]   条件d<=c0
H=0   条件  d>c0
其中c0=0.030，hs是常数为30000，pi=3.14
d= sqrt(x*x+(y-C1)*(y-C1));c1=0.028

fpw0412

#include "udf.h"
#define C1 0.028
#define C2 0.030
#define PI 3.14159
#define hs 300000000
DEFINE_SOURCE(heat_source, c, t, dS, eqn)
{
real r[ND_ND];
real d, source, x, y, co;
x= r[0];
y= r[1];
d= sqrt(x*x+(y-C1)*(y-C1));
co=1.+cos(PI*d/C2);
if(d<=(C2/2))
  {
           source =hs/2*co;
   dS[eqn]=0;
  }
   else
    source=0;
return source;
}

swimming02

你只是导入到fluent了，还应该把源项加到相应的流体区域中。

caicainide

原帖由 swimming02 于 2010-8-20 23:43 发表

                               
登录/注册后可看大图

你只是导入到fluent了，还应该把源项加到相应的流体区域中。

怎么添加就是source terms 里面添加吗？

fangdian

能量源好像对温度求偏导。

simwon

C_CENTROID(r,c,t);
x=r[0];
y=r[1];
不然xy无法付值