24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1507  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

snecu

新虫 (小有名气)

[交流] UDF 求面上压力积分,与直接Fluent得到的不同。已有2人参与

写了一个UDF求面上的静压积分,但是与fluent自带的积分结果不同,不知道程序哪里出错了。求指教。
#include "udf.h"
#define rvol1_ID 22
DEFINE_EXECUTE_AT_END(pressureudf)
{
FILE *fp1;
real time1=CURRENT_TIME;
Domain *d;
Thread *thread;
face_t f;
double A[ND_ND];
double pressure=0.;
d= Get_Domain(1);
thread=Lookup_Thread(d, rvol1_ID);
begin_f_loop(f,thread)
{
F_AREA(A,f,thread);
pressure+=F_P(f,thread)*NV_MAG(A);
}
end_f_loop(f,thread)
fp1=fopen("pressure223.txt","a";
fprintf(fp1,"%g %g\n",time1,pressure);
fclose(fp1);

}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 snecu 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见