24小时热门版块排行榜    

查看: 2062  |  回复: 6

horstchen

银虫 (初入文坛)

[求助] fluent运用UDF三维弯管模拟, 如何给定进口圆周速度(初学者) 已有1人参与

本人对三维弯管进行数值模拟,模型如图1所示。
须给定进口速度,准备按轴向、径向和周向给定(velocity  specification method采用component   圆柱坐标系)。
径向速度取:0 m/s
轴向速度为常量
周向速度采用UDF 给定,具体速度公式图2所示。

之前曾采用如下UDF:
DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
  real x[ND_ND]; /* this will hold the position vector */
  real y;
  face_t f;
  begin_f_loop(f,thread)
{
  F_CENTROID(x, f, thread);
  y=x[1];
  F_PROFILE(f, thread, position) = 25*pow((y*y)/(1.44*1.44),0.25);
}
  end_f_loop(f, thread)
}

求虫友们帮帮忙,给点意见,看如何实现图2中进口圆周速度

fluent运用UDF三维弯管模拟, 如何给定进口圆周速度(初学者)
图1


fluent运用UDF三维弯管模拟, 如何给定进口圆周速度(初学者)-1
图2
回复此楼

» 猜你喜欢

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

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

mwm621

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
horstchen: 金币+8, ★★★很有帮助, pow((x*x+y*y)/(1.44*1.44),0.25) 应该是xc*XC+y*y 2014-06-04 09:29:17
DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
  real x[ND_ND]; /* this will hold the position vector */
  real xc,y;
  face_t f;
  begin_f_loop(f,thread)
{
  F_CENTROID(x, f, thread);
  xc=x[0];
  y=x[1];
  F_PROFILE(f, thread, position) = 25*pow((x*x+y*y)/(1.44*1.44),0.25);
}
  end_f_loop(f, thread)
}


楼主是说这个意思么?
2楼2014-06-02 11:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

horstchen

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by mwm621 at 2014-06-02 11:13:05
DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
  real x; /* this will hold the position vector */
  real xc,y;
  face_t f;
  begin_f_loop(f,thread)
{
  F_CENTROID(x, f, thread);
  xc=x; ...

fluent提示错误
invalid type for binary expression: array 3 of float * array 3 of float.

应该是
F_PROFILE(f, thread, position)=25*pow((x*x+y*y)/(1.44*1.44),0.25);
这句出错了
3楼2014-06-02 12:27:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mwm621

木虫 (小有名气)

引用回帖:
3楼: Originally posted by horstchen at 2014-06-02 12:27:12
fluent提示错误
invalid type for binary expression: array 3 of float * array 3 of float.

应该是
F_PROFILE(f, thread, position)=25*pow((x*x+y*y)/(1.44*1.44),0.25);
这句出错了...

real x[ND_ND]; /* this will hold the position vector */
...
{
  F_CENTROID(x, f, thread);
  xc=x[0];
  y=x[1];

是这个样子啊,x[ND_ND]后面必须有[]指明是哪个坐标上的点。。。
4楼2014-06-02 12:50:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

horstchen

银虫 (初入文坛)

引用回帖:
4楼: Originally posted by mwm621 at 2014-06-02 12:50:22
real x; /* this will hold the position vector */
...
{
  F_CENTROID(x, f, thread);
  xc=x;
  y=x;

是这个样子啊,x后面必须有[]指明是哪个坐标上的点。。。...

x[ND_ND]后面必须有[]指明是哪个坐标上的点

不是这个,fluent编译的UDF有x[ND_ND]

fluent提示 line 12 出问题,对应的是这句
F_PROFILE(f, thread, position)=25*pow((x*x+y*y)/(1.44*1.44),0.25)
5楼2014-06-02 13:57:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mwm621

木虫 (小有名气)

【答案】应助回帖

不知道这么做得到楼主需要的结果了么?
6楼2014-06-04 15:03:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

horstchen

银虫 (初入文坛)

引用回帖:
6楼: Originally posted by mwm621 at 2014-06-04 15:03:48
不知道这么做得到楼主需要的结果了么?

进口的周向速度与实际的不太符合,还需修改。
7楼2014-06-04 16:41:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 horstchen 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 7/350 2026-08-22 21:39 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见