| 查看: 2009 | 回复: 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) } |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
6楼2013-05-16 20:38:09
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