|
发表于 2014-5-21 15:00:02
|
显示全部楼层
如果是轴对称,需要乘以2PI。轴对称情况下,如果有线段AB,fluent汇报的是线段AB绕轴旋转一周后形成的三维旋转体表面积。如果是用UDF求,那么fluent中汇报的值应该等于2*PI*sum(Ai),其中sum表示求和,Ai为线段AB上每一小段的长度。
real sum=0;
face_t f;
real Ai;
Thread*tf;
real NV_VEC(A);
tf=Lookup_Thread(domain,Line_AB_ID);
begin_f_loop(f,tf)
{
F_AREA(f,tf,A);
Ai=NV_MAG(A);
sum+=2*M_PI*Ai;
}
end_f_loop(f,tf)
此时,sum计算结果就是fluent后处理中汇报的值。 |
|