24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1503  |  回复: 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的回帖

星星在我家门

木虫 (著名写手)

热爱互联网


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问楼主毕业了吗?这问题最后怎么解决的啊
于无声听惊雷
2楼2018-09-03 10:39:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MoonLee777

新虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2楼: Originally posted by 星星在我家门 at 2018-09-03 10:39:00
请问楼主毕业了吗?这问题最后怎么解决的啊

不要这样做,求面上的值一定是要由网格中心点的值计算得来,这是有限体积法的基础,

发自小木虫IOS客户端
3楼2018-09-03 11:08:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 snecu 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见