24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sailote

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 140.5
帖子: 28
在线: 18.6小时
虫号: 1615473
注册: 2012-02-13
专业: 传热传质学

[求助] 用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

回复此楼

» 猜你喜欢

导师想让我从独立一作变成了共一第一已经有8人回复
博士读完未来一定会好吗已经有23人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
申请2026年博士已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

初始化温度场的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人回复

1楼2014-05-14 15:02:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhaolei539

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 324.1
红花: 1
帖子: 76
在线: 21.1小时
虫号: 2870773
注册: 2013-12-14
性别: MM
专业: 制造系统与自动化

引用回帖:

4楼: Originally posted by woshihuolei at 2014-05-15 09:59:34
如果是前一时刻赋给当前时刻的话应该是用DEFINE_EXECUTE_AT_END的吧 DEFINE_ADJUST是每次迭代都计算一次

你好！想请教提取某个边界面的温度，赋值给另一个边界面
面上不同位置温度不一样不知道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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 15 个回答

sailote

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 140.5
帖子: 28
在线: 18.6小时
虫号: 1615473
注册: 2012-02-13
专业: 传热传质学

后来我怀疑是不是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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖