| 查看: 4664 | 回复: 14 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
用udf实现出口温度赋值给进口温度已有1人参与
|
||
|
最近在搞一个非稳态算例，前一时刻流体出口面的平均温度赋值给当前时刻的进口面。初始状态时，流体区域平均分成前后两个区域，前后的温度不一样。 我网上搞了一个udf(如下)，计算中，用monitor监控进口和出口的平均温度，每10个时间步长输出一次。发现两者温度曲线相差很大，如图所示（纵坐标为平均温度K，横坐标为时间步长）。请问是这个udf的问题还是怎啦？ real NV_VEC(A); real avg_temp; 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,14);/14是出口面的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_PROFLIE(inlet_T,t,i) {begin_f_loop(f,t) {F_PROFILE(f,t,i) = avg_temp; } end_f_loop(f,t) }  1212.JPG |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有9人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
13楼2015-12-14 13:06:15
|
后来我怀疑是不是DEFINE_PROFLIE这个有问题，所以把udf改成如下，发现结果像之前一样（上图所示） Thread *thread_in; Thread *thread_out; face_t f; Domain *domain; real avg_temp; real sum_T_A=0.0; real sum_A=0.0; real NV_VEC(A); DEFINE_ADJUST(adjust,d) { domain=Get_Domain(1); thread_in=Lookup_Thread(domain,15); /*15是进口面的id*/ thread_out=Lookup_Thread(domain,14); /*14是出口面的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,thread_in) { F_T(f,thread_in) = avg_temp; } end_f_loop(f,thread_in) } |
2楼2014-05-14 15:04:44
woshihuolei
银虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 270.2
- 红花: 4
- 帖子: 191
- 在线: 55.4小时
- 虫号: 2760275
- 注册: 2013-10-28
- 专业: 工程热物理相关交叉领域
4楼2014-05-15 09:59:34
5楼2014-05-15 13:11:22