| 查看: 4821 | 回复: 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 |
回复此楼» 猜你喜欢
26申博
已经有3人回复
-
求调剂，总分315，考的生物医药，一志愿湖南师范大学。调剂到任何专业都可以
已经有3人回复
-
材料专硕306英一数二
已经有7人回复
-
材料工程专硕274一志愿211求调剂
已经有6人回复
-
318求调剂
已经有5人回复
-
一志愿南京大学，080500材料科学与工程，调剂
已经有4人回复
-
332求调剂
已经有5人回复
-
271求调剂
已经有14人回复
-
275求调剂
已经有4人回复
-
302求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 初始化温度场的udf
已经有14人回复
- 解释型UDF如何调用温度梯度
已经有9人回复
- 入口速度和温度都随时间变化，怎么写UDF？
已经有10人回复
- 编UDF将一个边界上的温度作为另一个边界条件的温度值
已经有13人回复
- 用udf实现coupled边界固体区域中温度最值提取
已经有8人回复
- udf如何实现温度热源的自然变化？
已经有12人回复
- fluent中UDF编译的问题！
已经有6人回复
- 如何用UDF获取某个温度值下的坐标
已经有7人回复
- fluent 中二维模型中，材料属性随着温度变化的udf源代码？
已经有14人回复
- UDF编译时出现语法错误
已经有13人回复
- 一个fluent UDF问题
已经有4人回复
- udf更新壁面温度求助
已经有15人回复
- DEFINE_SOURCE(energy,c,t,ds,eqn)条件赋值的相关问题
已经有8人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- UDF 热源随温度变化的处理
已经有4人回复
- UDF添加动能源项，发现温度没什么变化
已经有20人回复
- 激光焊接udf，加载温度场不对
已经有7人回复
- fluent能量源项中温度T对时间的t的偏导数，即dT/dt如何用udf编写能量源项
已经有6人回复
- 求助气固两相流的边界条件问题
已经有8人回复
- 使用UDF后，初始化错误
已经有11人回复
- 【求助】fluent计算中，udf编译通过，初始化出错
已经有11人回复
|
后来我怀疑是不是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
3楼2014-05-14 20:05:52
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
woshihuolei
银虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 270.2
- 红花: 4
- 帖子: 191
- 在线: 55.4小时
- 虫号: 2760275
- 注册: 2013-10-28
- 专业: 工程热物理相关交叉领域
6楼2014-05-15 18:44:35
zhaolei539
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 324.1
- 红花: 1
- 帖子: 76
- 在线: 21.1小时
- 虫号: 2870773
- 注册: 2013-12-14
- 性别: MM
- 专业: 制造系统与自动化
7楼2014-06-12 20:37:53
8楼2014-06-13 16:17:58
zhaolei539
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 324.1
- 红花: 1
- 帖子: 76
- 在线: 21.1小时
- 虫号: 2870773
- 注册: 2013-12-14
- 性别: MM
- 专业: 制造系统与自动化
|
你好！想请教提取某个边界面的温度，赋值给另一个边界面 面上不同位置 温度不一样 不知道UDM怎么存储的？我这样写行不行？ #include "udf.h" DEFINE_PROFILE(hnt_temperature,thread) { face_t f; real tem_liquid; real tem_solid; Thread *thread_liquid=Lookup_Thread(domain,14); /*14是交界面液体侧面的id*/; Thread *thread_solid=Lookup_Thread(domain,15); /*15是交界面固体侧面的id*/; begin_f_loop(f,thread_liquid) { temp=F_T(f,thread_liquid); F_UDMI(f,t,0)=temp_liquid; } end_f_loop(f,thread_liquid) begin_f_loop(f,thread_solid) { tem_solid=F_UDMI(f,t,0); F_PROFILE(f,thread_solid)=tem_solid; } end_f_loop(f,thread_solid) } 麻烦您帮忙看看 谢谢！ |
9楼2014-06-15 09:39:39
zhaolei539
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 324.1
- 红花: 1
- 帖子: 76
- 在线: 21.1小时
- 虫号: 2870773
- 注册: 2013-12-14
- 性别: MM
- 专业: 制造系统与自动化
|
楼主能不能把你的代码发给我 我的问题跟你类似，但是代码导进去之后有错误，您帮我看看我的这个？给我参考一下？卡这里一个月了，不胜感激~~ #include "udf.h" DEFINE_PROFILE(hnt_temperature,thread) { face_t f; real tem_liquid; real tem_solid; Thread *thread_liquid=Lookup_Thread(domain,14); /*14是交界面液体侧面的id*/; Thread *thread_solid=Lookup_Thread(domain,15); /*15是交界面固体侧面的id*/; begin_f_loop(f,thread_liquid) { temp=F_T(f,thread_liquid); F_UDMI(f,t,0)=temp_liquid; } end_f_loop(f,thread_liquid) begin_f_loop(f,thread_solid) { tem_solid=F_UDMI(f,t,0); F_PROFILE(f,thread_solid)=tem_solid; } end_f_loop(f,thread_solid) } 麻烦您帮忙看看 谢谢！ |
10楼2014-06-15 18:30:21
50