|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
我想计算interior面(二维问题)的截面平均总压,编写的udf见下,一初始化就出现进去错误,已尝试重启fluent、初始化后重新加载udf、分母(sum_A)设置成常数,但都不对。。还有我看有大侠说可能是网格问题--有些面未进行投影?这个不是太懂。。我认为最可能的问题就是udf问题,但是实在看不出来了,烦请各位大侠帮帮忙了。。。小女子在此拜谢了啊。。
#include"udf.h"
real con=0.019895;
real rho=998.2;
real in_total_p;
real sum_A=0.0;
real sum_p_A=0.0;
real NV_VEC(A),NV_VEC(vel_vec);
face_t f_in;
Thread *thread_in;
Domain *domain;
DEFINE_ADJUST(adjust,d)
{
real dym_p;
real vel_in;
domain=Get_Domain(1);
thread_in=Lookup_Thread(domain,11);
begin_f_loop(f_in,thread_in)
{
NV_DS(vel_vec,=,F_U(f_in,thread_in),F_V(f_in,thread_in),F_W(f_in,thread_in),*,1.0);
vel_in=NV_MAG2(vel_vec);
F_AREA(A,f_in,thread_in);
sum_A+=NV_MAG(A);
dym_p=0.5*rho*vel_in;
sum_p_A+=NV_MAG(A)*(F_P(f_in,thread_in)+dym_p);
}
end_f_loop(f_in,thread_in)
in_total_p=sum_p_A/sum_A;
} |
|