| 查看: 5257 | 回复: 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 |
回复此楼» 猜你喜欢
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
40+中年人失业后，还能找到好工作吗？
已经有5人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有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人回复
woshihuolei
银虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 270.2
- 红花: 4
- 帖子: 191
- 在线: 55.4小时
- 虫号: 2760275
- 注册: 2013-10-28
- 专业: 工程热物理相关交叉领域
4楼2014-05-15 09:59:34
|
后来我怀疑是不是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
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
200
