| 查看: 1612 | 回复: 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); } |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
星星在我家门
木虫 (著名写手)
热爱互联网
- 应助: 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
50