| 查看: 2354 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
udf 根据监测点的压力值设置入口的质量流量
|
||
|
初学udf,请大家帮忙看下哪有问题，或者如果这样思路不对，应该如何做，已被折磨n天，跪求大神指教。 #include "udf.h" #include "mem.h" #define m_standard 8000. #define m_interval 500. DEFINE_ADJUST(mass_inlet,d) { Thread *t1,*t2; cell_t c; face_t f; real xc[ND_ND]; int i; real p=0.; thread_loop_c(t1,d) { /*坐标为(10,10)的点为监测点，循环找到一些和它接近的点， 把这些点的压力值赋给p,并认为最后一次赋的值即为监测点的值*/ begin_c_loop(c,t1) { C_CENTROID(xc,c,t1); if(sqrt(ND_SUM(pow(xc[0]-10.,2.), pow(xc[1]-10.,2.)))<0.001) p=C_P(c,t1); } end_c_loop(c,t1) } /*根据监测点的压力值，来改变入口的质量流量*/ begin_f_loop(f,t2) { if(p<=350000.) F_PROFILE(f,t2,i) = m_standard+m_interval; else if(p>350000.&&p<380000.) F_PROFILE(f,t2,i) = m_standard; else F_PROFILE(f,t2,i) = m_standard-m_interval; } end_f_loop(f,t2) } |
回复此楼» 猜你喜欢
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有4人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
什么时候开奖？
已经有10人回复
-
2026国自然函评费到账
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
让我中一个面上吧！
已经有14人回复
-
今天放榜吗？
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
5楼2013-05-16 20:36:33
shuai_1314
捐助贵宾 (小有名气)
- 应助: 33 (小学生)
- 金币: 1861.7
- 红花: 6
- 帖子: 167
- 在线: 163.1小时
- 虫号: 2415320
- 注册: 2013-04-12
- 性别: GG
- 专业: 内流流体力学
2楼2013-05-13 07:57:05
3楼2013-05-13 19:21:11
shuai_1314
捐助贵宾 (小有名气)
- 应助: 33 (小学生)
- 金币: 1861.7
- 红花: 6
- 帖子: 167
- 在线: 163.1小时
- 虫号: 2415320
- 注册: 2013-04-12
- 性别: GG
- 专业: 内流流体力学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
xiegangmai: 金币+2, 谢谢参与！ 2013-05-15 11:22:00
ltpreason: 金币+5, ★有帮助 2013-05-16 20:37:12
xiegangmai: 金币+2, 谢谢参与！ 2013-05-15 11:22:00
ltpreason: 金币+5, ★有帮助 2013-05-16 20:37:12
|
Thread *thread_name = Lookup_Thread(domain,zone_ID) 用这个命令，zone_ID就是你的各种边界和流动区域对应的那个编号，比如你的fluid编号是2，你的入口边界编号是5 Thread *t1 = Lookup_Thread(d,2)； Thread *t2 = Lookup_Thread(d,5)； 应该这样定义，多看看帮助文档搞清楚domain和thread |
4楼2013-05-14 23:09:57
