| 查看: 3352 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
Fluent，UDF中的UDMI使用问题求助！ 已有1人参与
|
||
|
请问各位，我在fluent得UDF编程中，用UDMI的过程中遇到了如下问题！发现UDMI中的值在DEFINE_PROFILE中调用错误，我用patch初始化了所有UDMI都为303.15，可是最后输出壁面平均温度却总是零点几K，应该是300K左右才对！ 我对UDMI的使用不熟悉，不知道能不能这样用。 还有我用DEFINE_EXECUTE_AT_END的原因是，我想让壁面温度在每个time step才更新一次，而不是每个迭代步更新！ 按照我下面写的程序，应该是UDMI的数每个时间步更新一次，然后赋值给壁面的free stream tem，可是实际运行总是走不通，求大家懂这方面的指点一下，搞了好几天这个程序，一直搞不定，脑袋都快蒙了。。 #include "udf.h" #include "mem.h" DEFINE_EXECUTE_AT_END(updatetimestepend) { Domain *d; Thread *t,*t_walltem_4,*t_walltem_3, face_t f; cell_t c; real NV_VEC(A4),NV_VEC(A3); d=Get_Domain(1); real sum_A4=0.0,sum_A3=0.0; real sum_T4=0.0,sum_T3=0.0; real ave_T4,ave_T3; t_walltem_4=Lookup_Thread(d,66);/*获取ID=66壁面指针*/ t_walltem_3=Lookup_Thread(d,67);/*获取ID=67壁面指针*/ begin_f_loop(f,t_walltem_4)/*求壁面66平均温度*/ { F_AREA(A4,f,t_walltem_4); sum_A4+=NV_MAG(A4); sum_T4+=NV_MAG(A4)*F_T(f,t_walltem_4); } end_f_loop(f,t_walltem_4) ave_T4=sum_T4/sum_A4; begin_f_loop(f,t_walltem_3)/*求壁面67平均温度*/ { F_AREA(A3,f,t_walltem_3); sum_A3+=NV_MAG(A3); sum_T3+=NV_MAG(A3)*F_T(f,t_walltem_3); } end_f_loop(f,t_walltem_3) ave_T3=sum_T3/sum_A3; thread_loop_c(t,d)/*将需要的值存在UDM中，在下面的PROFILE调用*/ { begin_c_loop(c,t) { C_UDMI(c,t,5)=C_UDMI(c,t,4)-0.2*(C_UDMI(c,t,4)-ave_T4); C_UDMI(c,t,4)=C_UDMI(c,t,3)-0.2*(C_UDMI(c,t,3)-ave_T3); C_UDMI(c,t,3)=303.15； } end_c_loop(c,t) } } DEFINE_PROFILE(tem_5,t,i)/*66壁面的Free stream temperature*/ { face_t f; cell_t c; real tem5; tem5=C_UDMI(c,t,5); begin_f_loop(f,t) { F_PROFILE(f,t,i)=tem5; } end_f_loop(f,t) } DEFINE_PROFILE(tem_4,t,i)/*65壁面的Free stream temperature*/ { face_t f; cell_t c; real tem4; tem4=C_UDMI(c,t,4); begin_f_loop(f,t) { F_PROFILE(f,t,i)=tem4; } end_f_loop(f,t) } DEFINE_PROFILE(tem_3,t,i)/*64壁面的Free stream temperature*/ { face_t f; cell_t c; real tem3; tem3=C_UDMI(c,t,3); begin_f_loop(f,t) { F_PROFILE(f,t,i)=tem3; } end_f_loop(f,t) } |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
国自然面上复盘~欢迎讨论
已经有15人回复
-
中青基了要发朋友圈吗？
已经有3人回复
3楼2018-12-18 18:55:39
2楼2018-12-18 00:12:25
4楼2018-12-18 20:35:42
10