24小时热门版块排行榜

返回列表

liyujie1245

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 611.3
帖子: 17
在线: 41.5小时
虫号: 2395140
注册: 2013-04-02
专业: 电机与电器

[求助] 在fluent里面的udf中profile宏考虑传热模块怎么写udf程序已有1人参与

我刚刚开始使用udf还不是很熟练，我做的是一个灭弧室的模型，考虑到栅片的吸热过程。
http:// $\"在fluent里面的udf中profile宏考虑传热模块怎么写udf程序\"$
我描述一下：就是灭弧室右端是出气口，中间有5个栅片（铁），其余部分是空气。在在此模型中考虑栅片从空气中吸热，使栅片的温度提高，但是我还需要考虑栅片的热导率和比热随温度的变化（这些我在udf中已经完成），但是导师还希望我考虑栅片表面的蒸发潜热。我查看了蒸发潜热的工具书，它和栅片表面的温度有关，
   在T温度下，则单位时间单位面积上带走的热量为：Q=[4.0117*(10^25)*(10^((13.27-19710/T)*T^(-1.27)]/(T^0.5)
但是我在以往看udf里面编写profile宏的时候，都是直接定义边界压强或者温度等于常数，而我觉得我的问题比较特殊，因为它是一个根据边界温度来改变能量的一个过程。
而且因为热量的空气栅片的传热过程默认的是couple设置,我也觉得不好改能量传递的边界，所以我很困惑。
我现在不了解的有这么几个问题：
问题1:我看以前一些人发的profile宏的程序中肯定有：
begin_f_loop(f,thread_solid)
{
tem_solid=F_UDMI(f,t,0);
F_PROFILE(f,thread_solid)=.....
}
  end_f_loop(f,thread_solid)
}
我想知道这个F_PROFILE究竟是如何判断它是设置的压力值还是温度值。（因为按照我的想法我得设置成能量值)
问题2：上面的程序那些帖子里面解释说：设置边界温度为多少多少。但是在我的模型里面有很多的边界：栅片与空气的接触边界、右端出气口边界以及其他三个边的密闭边界。我不知道怎么写udf的profile宏程序才能找到我需要的那个“栅片与空气的接触边界”，再对它进行设置。
下面是我参考了以前几个帖子里面的程序进行的修改：但是加载进去总是出错，我不知道是c语言编译的问题，还是我设置的思想方法不对，希望小木虫里面做过有关传热仿真的朋友能够指导一下，谢谢啦！
下面的我修改的程序：
程序1：
#include "udf.h"
DEFINE_PROFILE(hnt_temperature,thread)
{
face_t f;
real tem_solid;
Thread *thread_gas=Lookup_Thread(domain,8); /*8是交界面流体侧面的id*/;
Thread *thread_solid=Lookup_Thread(domain,7); /*7是交界面固体侧面的id*/;
pvl=133.33*pow(10,(-19710/F_T(f,thread)))*pow(10,13.27)*pow(F_T(f,thread),-1.27);    /* the unit is Pa */
begin_f_loop(f,thread_solid)
{
tem_solid=F_UDMI(f,t,0);
F_PROFILE(f,thread_solid)=-He*(3.00882e23)*pvl/pow(F_T(f,thread),0.5);    /*液态条件下蒸发带走的热量 */;
}
  end_f_loop(f,thread_solid)
}
//程序1我觉得可以定位到我想要的栅片与气体接触的边界，但是还是显示函数有错误，我不知道是方法问题还是C语言上面的语言错误//
程序2 ：
#include "udf.h"
DEFINE_PROFILE(Plate_vapour,thread,position)
{
real He=5.1e-19; /* evaporation heat of an atom for copper material */
real x[ND_ND],xx,yy,pvs,pvl,pv,Fev;
  face_t f;
  begin_f_loop(f,thread)
{
  F_CENTROID(x,f,thread);
  yy=x[1];
  xx=x[0];
  pvl=133.33*pow(10,(-19710/F_T(f,thread)))*pow(10,13.27)*pow(F_T(f,thread),-1.27);    /* the unit is Pa */
  if(((yy=0.008)||(yy=0.010)||(yy=0.012)||(yy=0.014)||(yy=0.016)||yy=0.018)||(yy=0.020)||(yy=0.022)||(yy=0.024)||(yy=0.026))&&(0.34<=xx<=0.46))
{
  F_PROFILE(f,thread,position)=-He*(3.00882e23)*pvl/pow(F_T(f,thread),0.5);  /*蒸发带走的热量 */
  }
else if((xx=0.0034)||(xx=0.0046))&&((0.008<=yy<=0.010)||(0.012<=yy<=0.014)||(0.016<=yy<=0.018)||(0.020<=yy<=0.022)||(0.024<=yy<=0.026)))
  {
  F_PROFILE(f,thread,position)=-He*(3.00882e23)*pvl/pow(F_T(f,thread),0.5);    /*蒸发带走的热量 */
  }
else
  {
  F_PROFILE(f,thread,position)=0;
  }
}end_f_loop(f,thread)
}
//程序2 是我听师兄说可以通过坐标来确定边界，然后再写能量函数。//
再次感谢小木虫里面有做过有关传热仿真的朋友能够指导一下，谢谢啦！

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有119人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有6人回复

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

fluent用udf输入进口速度，求指点！！！已经有7人回复
急，在FLUENT中如何同时调用两个UDF？已经有11人回复
FLUENT三维建模后入口速度为抛物线怎么写UDF？已经有8人回复
fluent udf加载后无法运行的问题已经有6人回复
入口速度和温度都随时间变化，怎么写UDF？已经有10人回复
fluent UDF 修改物性已经有13人回复
关于材料属性定义的UDF 已经有6人回复
fluent中UDF编译的问题！已经有11人回复
fluent udf 编译问题 DEFINE_TRANSIENT_PROFILE(rotation_rate_ramp,time) 已经有8人回复
fluent中UDF编写自定义矢量——UDF大牛请进已经有7人回复
FLUENT中编写UDF时指针变量是否要初始化已经有4人回复
fluent 中二维模型中，材料属性随着温度变化的udf源代码？已经有14人回复
fluent udf请教已经有16人回复
fluent UDF 语句 pow函数值导出已经有5人回复
Fluent中UDF自定义了一个热流密度已经有12人回复
fluent中通过UDF获得系统时间已经有15人回复
【求助】FLUENT UDF 气液两相流汽蚀汽化潜热(问题详见内容) 已经有21人回复
求助！！！我想要在fluent中针对一个面进行积分，那么UDF应该用什么宏函数？怎么用？已经有5人回复
好心人帮帮忙吧，为什么这个UDF导入fluent以后，在边界条件处没有显示~~~ 已经有17人回复

1楼 2014-08-24 00:06:37

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

liyujie1245

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 611.3
帖子: 17
在线: 41.5小时
虫号: 2395140
注册: 2013-04-02
专业: 电机与电器

http:// $\"在fluent里面的udf中profile宏考虑传热模块怎么写udf程序-1\"$

赞一下

回复此楼

2楼2014-08-24 00:10:42

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

tank_fire

铜虫 (小有名气)

应助: 11 (小学生)
金币: 476.7
红花: 3
帖子: 107
在线: 60.9小时
虫号: 470732
注册: 2007-11-29
专业: 多相流热物理学

【答案】应助回帖

看不见图，虽然我不搞传热，但是我觉得你应该分别写2个UDF，一个用在固体界面（假设wall a），一个用在液体界面（wall b）。在Boundary conditions里面在固体界面wall a里thermal项heat flux（或者其他导热的计算方式）加载第一个udf，wall b的加载第二个，这样就不用考虑界面位置的问题，只要loop

赞一下

回复此楼

3楼2014-08-27 23:51:57

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