| 查看: 1624 | 回复: 2 | |||
[交流]
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); } |
回复此楼» 猜你喜欢
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
星星在我家门
木虫 (著名写手)
热爱互联网
- 应助: 21 (小学生)
- 金币: 2024.1
- 散金: 1232
- 红花: 8
- 沙发: 1
- 帖子: 1070
- 在线: 173.8小时
- 虫号: 3229413
- 注册: 2014-05-23
- 专业: 传递过程
2楼2018-09-03 10:39:00
MoonLee777
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2274.6
- 红花: 17
- 帖子: 839
- 在线: 78小时
- 虫号: 5947368
- 注册: 2017-03-12
- 性别: GG
- 专业: 流体力学
3楼2018-09-03 11:08:47
20