24小时热门版块排行榜    

查看: 1199  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

guoyhchina

金虫 (小有名气)

[交流] 请帮忙看下关于传热系数的UDF对不对,谢了 已有1人参与

我的问题是要写一个传热系数是壁温及壁面附近流场温度和速度函数的UDF:
      hc=(ts+ta)^0.3*U^0.7   
其中hc是传热系数,ts为壁面温度,ta为壁面附近流体温度,U为壁面附近流体速度

这个问题已困扰很久了,也在论坛里发贴问过,自己写了个UDF,请这方面的高手帮忙,看下我

写的这个UDF有什么问题?该怎样修改?不胜感谢。

/**********************************************************************/
/* UDF that implements a heat transfer coefficient hc*/
/* hc=(ts+ta)^0.3*U^0.7   其中ts为壁面温度,ta为壁面附近流体温度,U为壁面附近流体速度 */
/**********************************************************************/

#include "udf.h"
DEFINE_PROFILE(heat_transfer_coef,t,i)
{
   cell_t c;
   face_t f;
   Thread *t;
   real hc;

if (BOUNDARY_FACE_THREAD_P(t)) /*Most face values will be available*/

  real wt,ct,cv;
/*wt is wall temperature, ct and cv are cell temperature and velocity of the adjacent

wall,respectively*/

   wt=F_T(f,t);
    ct=C_T(c,t);
    cv=C_U(c,t);
    hc=power((wt+ct), 0.3)*power(cv, 0.7);

return hc;
}
回复此楼

» 猜你喜欢

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

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

guoyhchina

金虫 (小有名气)

引用回帖:
2楼: Originally posted by kylafree at 2012-08-22 14:39:02
幂函数为pow而非power,另外你自己测试了么?测试后有什么问题?

谢谢kylafree的提醒,我已把power改成pow了; 另外,“Thread *t;”  语句是多余的,因DEFINE_PROFILE(heat_transfer_coef,t,i)中已给定了指针,所以删除“Thread *t;” 。最后修改如下:

/**********************************************************************/
/* UDF that implements a heat transfer coefficient hc*/
/* hc=(wt+ct)^0.3*cv^0.7   其中wt为壁面温度,ct为壁面附近流体温度,cv为壁面附近流体速度 */
/**********************************************************************/
#include "udf.h"
DEFINE_PROFILE(heat_transfer_coef,t,i)
{
  cell_t c;
  face_t f;
  real hc;
  if(BOUNDARY_FACE_THREAD_P(t)) /*Most face values will be available*/
  real wt, ct, cv;
  wt = F_T(f,t);
  ct = C_T(c,t);
  cv = C_U(c,t);
  hc=pow((wt+ct), 0.3)*pow(cv, 0.7);
return hc;
}

但编译还是通不过去,有两处错误,提示如下:
if (BOUNDARY_FACE_THREAD_P(t))   有语法错误 (parse error);
wt=F_T(f,t);  未定义变量wt,但我明明已在前面定义了( real wt,ct,cv; )
请高手再次指点,谢谢
3楼2012-08-24 10:48:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

kylafree

至尊木虫 (知名作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
幂函数为pow而非power,另外你自己测试了么?测试后有什么问题?
2楼2012-08-22 14:39:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 7/350 2026-02-28 20:42 by 好好好1233
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 5/250 2026-02-28 20:11 by iwuli
[考博] 博士推荐 +5 花儿笑? 2026-02-21 6/300 2026-02-28 18:53 by nxgogo
[考研] 311求调剂 +7 南迦720 2026-02-28 7/350 2026-02-28 18:28 by leonnulll
[考研] 材料学调剂 +4 提神豆沙包 2026-02-28 4/200 2026-02-28 18:26 by houyaoxu
[考研] 285求调剂 +5 满头大汗的学生 2026-02-28 5/250 2026-02-28 18:10 by 材料专硕调剂;
[考研] 材料调剂 +3 爱擦汗的可乐冰 2026-02-28 3/150 2026-02-28 18:06 by houyaoxu
[考研] 化工专硕348,一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 290求调剂 +4 材料专硕调剂; 2026-02-28 5/250 2026-02-28 13:32 by houyaoxu
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[基金申请] 什么是人一生最重要的? +10 瞬息宇宙 2026-02-21 10/500 2026-02-27 08:46 by tfang
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
信息提示
请填处理意见