| 查看: 3661 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
出口压力的UDF编写问题 已有4人参与
|
|||
|
小弟刚接触CFD不久，在fluent计算二相流中想设置outlet为压力出口。模型二维矩形网格，出口处y轴0.5m以上为空气，压力为零，y轴0.5以下为水，压力随水深变化。这个UDF是参照一个三维类似模型压力出口改的，可是怎么都没法compile成功，求大神帮忙看看有没有什么新人小白的错误什么的…… #include "udf.h" #define p 998.2 #define g 9.81 DEFINE_PROFILE(pressure_profile,thread,postion) { real x[ND_ND]; face_t f; real y; begin_f_loop(f,thread)/*loops over all face in the thread passed in the defing macro argument */ { F_CENTROID(x,f,thread); y=x[1] if (y>=0.5) {F-PROFILE(f,thread,position)=0;} else {F_PROFILE(f,thread,psotion)= p*g*(0.5-y);} end_f_loop(f,thread) } } |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有13人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一个FLUENT里面关于边界条件设置的问题
已经有5人回复
- UDF中如何编写体积分数的梯度
已经有9人回复
- 多孔介质燃烧编写UDF气固两相控制方程的问题
已经有5人回复
- 编UDF将一个边界上的温度作为另一个边界条件的温度值
已经有13人回复
- fluent求压降，设置边界条件如何设置。
已经有8人回复
- fluent中UDF编译的问题！
已经有6人回复
- 用UDF编写化学反应热作为源项
已经有12人回复
- FLUENT中编写UDF时指针变量是否要初始化
已经有4人回复
- UDF编译时出现语法错误
已经有13人回复
- 全部身家求大神帮忙解决fluent压力脉冲怎么编写UDF的问题
已经有8人回复
- 如何用UDF编写比热容随温度变化的程序，我试了一天了，总也不成功！
已经有4人回复
- fluent出口边界条件和回流设置问题
已经有9人回复
- 关于udf编译问题
已经有6人回复
- 如何编写UDF来读取txt格式的一组数据？
已经有11人回复
- UDF 编写变上限积分
已经有8人回复
- fluent中用udf编写三维速度边界条件出现错误，大家来指导一下啊
已经有9人回复
- udf编译出错！急
已经有10人回复
- 求助：UDF模型的编译问题
已经有15人回复
- fluent中如何编写风剖面的UDF
已经有6人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
我还是第一次见到 写的这么粗心的UDF，变量名就重复了几次，错误都不一样，LZ可真得认真点才行，还好这个程序简单。 我修改了一下，这个应该是能运行的 #include "udf.h" #define p 998.2 #define g 9.81 int position; DEFINE_PROFILE(pressure_profile,thread,position) { real x[ND_ND]; face_t f; real y; begin_f_loop(f,thread)/*loops over all face in the thread passed in the defing macro argument */ { F_CENTROID(x,f,thread); y=x[1]; if (y>=0.5) {F_PROFILE(f,thread,position)=0.0;} else {F_PROFILE(f,thread,position)= p*g*(0.5-y);} } end_f_loop(f,thread) } [ Last edited by 奔跑鱼 on 2014-3-11 at 23:21 ] |
8楼2014-03-11 23:01:40
yegavin
至尊木虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 11776.7
- 散金: 32
- 红花: 7
- 帖子: 1587
- 在线: 744.9小时
- 虫号: 1408483
- 注册: 2011-09-19
- 性别: GG
- 专业: 多相流热物理学
2楼2014-03-10 20:49:32
3楼2014-03-11 09:18:07
4楼2014-03-11 12:20:27
