Question

FLuent中使用UDF求解边界面的压力

小弟编写了一个UDF来求解边界面的压力，但是得出的结果是错误的，求高手帮忙看看，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);
}

Answer 1

这个程序求的是面平均压强，重点检查id给的对不对

追问

id给的是对的，把F_P换成F_W时求面平均速度是对的，但是求平均压强时出现错误

追答

f_P好像取的是静压，所以你先要搞清楚你为什么判断udf算出的结果是错的，判断依据是否错了？

追问

我在边界条件上设置入口压力未200Pa，然后用上面的宏提取入口面的平均压力，在控制台显示出来的结构为-180左右，所以判断结果是错的

Answer 2

header file: #define "mem.h"
and so on.