|
发表于 2010-4-27 20:03:27
|
显示全部楼层
将r表示成圆盘上任一点与圆心的距离即可。
#include "udf.h"
DEFINE_PROFILE(heatflux, t, i)
{
face_t f;
real a,r;
real x[ND_ND],y[ND_ND],z[ND_ND];
begin_f_loop(f, t)
{
F_CENTROID(x,f,t);
y[0]=0.0; //假设圆心坐标为(0,0,0),若为其他坐标更改即可。
y[1]=0.0;
y[2]=0.0;
NV_VV(z,=,x,-,y);
r=NV_MAG(z);
F_PROFILE(f, t,i) = a*r;
}
end_f_loop(f, t)
} |
|