24小时热门版块排行榜    

查看: 1275  |  回复: 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的回帖

kylafree

至尊木虫 (知名作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
幂函数为pow而非power,另外你自己测试了么?测试后有什么问题?
2楼2012-08-22 14:39:02
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 guoyhchina 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 7/350 2026-08-25 03:24 by BZKMTicpDhFj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见