| 查看: 1398 | 回复: 6 | |||
leixiaobo713银虫 (正式写手)
|
[求助]
如何通过udf实现一段时间出口流体的平均温度，作为下个时间段入口流体的初始温度！ 已有1人参与
|
| 如题：我想在udf中实现 将一段时内出口流体的平均温度，作为下个时间段流体的初始温度，不知道如何实现，请大神指点, 谢谢!!! |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求热流体的出口温度
已经有10人回复
- fluent 6.3 13.0版本中多相流的比较
已经有20人回复
- fluent能量源项中温度T对时间的t的偏导数，即dT/dt如何用udf编写能量源项
已经有6人回复
- 怎样用UDF获得管道管内流体的中心温度？
已经有3人回复
死天
专家顾问 (著名写手)
-
专家经验: +60 - 应助: 399 (硕士)
- 金币: 3648.9
- 红花: 89
- 帖子: 1336
- 在线: 443.8小时
- 虫号: 2643057
- 注册: 2013-09-09
- 专业: 流体力学
- 管辖: 仿真模拟
2楼2013-09-27 10:25:43
leixiaobo713
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1228.5
- 散金: 148
- 红花: 4
- 帖子: 490
- 在线: 121.9小时
- 虫号: 1997539
- 注册: 2012-09-13
- 性别: GG
- 专业: 金属功能材料
3楼2013-09-27 11:20:41
leixiaobo713
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1228.5
- 散金: 148
- 红花: 4
- 帖子: 490
- 在线: 121.9小时
- 虫号: 1997539
- 注册: 2012-09-13
- 性别: GG
- 专业: 金属功能材料
4楼2013-09-27 11:21:49
leixiaobo713
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1228.5
- 散金: 148
- 红花: 4
- 帖子: 490
- 在线: 121.9小时
- 虫号: 1997539
- 注册: 2012-09-13
- 性别: GG
- 专业: 金属功能材料
|
#include "udf.h" real NV_VEC(A); real avg_temp; real sum_T_A=0.0; real sum_A=0.0; real time; Thread *thread_out; face_t f; Domain *domain; DEFINE_ADJUST(adjust,d) { domain=Get_Domain(1); thread_out=Lookup_Thread(domain,4);/*4是出口面的id*/ begin_f_loop(f,thread_out) { F_AREA(A,f,thread_out); sum_A+=NV_MAG(A); sum_T_A+=NV_MAG(A)*F_T(f, thread_out); } end_f_loop(f,thread_out) avg_temp=sum_T_A/sum_A; begin_f_loop(f,t) { time=RP_Get_Real("flow-time"  ; if(time<=0.5) F_PROFILE(f,t,i) =298; else if(time>0.5&&time<=1) F_PROFILE(f,t,i) = avg_temp; } end_f_loop(f,t) } 加载时老师出现 begin_f_loop(f,t)中的 t: undeclared variable，这个是为什么，请您予以指点！ |
5楼2013-09-27 11:24:31
leixiaobo713
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1228.5
- 散金: 148
- 红花: 4
- 帖子: 490
- 在线: 121.9小时
- 虫号: 1997539
- 注册: 2012-09-13
- 性别: GG
- 专业: 金属功能材料
6楼2013-11-22 11:05:39
malasong9471
木虫 (正式写手)
- 应助: 58 (初中生)
- 金币: 7327
- 红花: 29
- 帖子: 639
- 在线: 117.1小时
- 虫号: 4100890
- 注册: 2015-09-25
- 性别: GG
- 专业: 工程热力学
【答案】应助回帖
|
少了DEFINE_PROFILE(inlet_temperature,t,i) { ——————————————————————————————————————————————————————— #include"udf.h" real avg_temp=303; real NV_VEC(A); int ID1=40; real sum_T_A=0.0; real sum_A=0.0; Thread *thread_out; face_t f; Domain *domain; DEFINE_ADJUST(adjust,d) { domain=Get_Domain(1); thread_out=Lookup_Thread(domain,ID1);/*40是流体ID*/ begin_f_loop(f,thread_out) { F_AREA(A,f,thread_out); sum_A+=NV_MAG(A); sum_T_A+=NV_MAG(A)*F_T(f,thread_out); } end_f_loop(f,thread_out) avg_temp=sum_T_A/sum_A; } DEFINE_PROFILE(inlet_temperature,t,i) { begin_f_loop(f,t) { F_PROFILE(f,t,i)=avg_temp+3; } end_f_loop(f,t) } |
7楼2016-09-27 19:52:28
