UDF求解边界面编写求教

Firehole

编写了一个求解边界面的UDF，但是求解时发现结果是错的。求高手指点下，不胜感激。
#include "udf.h"
real NV_VEC(A);
real p;
real sum_P_A=0.0;
real sum_A=0.0;
Thread *thread_in;
face_t f;
Domain *domain;
DEFINE_ADJUST(inlet,d)
{
domain=Get_Domain(1);
thread_in=Lookup_Thread(domain,5);
begin_f_loop(f,thread_in)
{
  F_AREA(A,f,thread_in);
  sum_A+=NV_MAG(A);
  sum_P_A+=NV_MAG(A)*F_P(f,thread_in);
}
end_f_loop(f,thread_in);
p=sum_P_A/sum_A;
printf("P is:%g\n",p);
}

fty0083

把fluent的参量定义放到宏里面去，real型可以做全局变量

Thread *thread_in;
face_t f;
Domain *domain;

另外，先定义domain，再定义thread，调整下顺序，有时候这个说不定就出错了。

Firehole

我编的这个求一个边界面的平均压力的UDF没有问题是吧？只需要调整参量顺序就行了是不？