24小时热门版块排行榜    

查看: 1778  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +7 明远求学 2026-03-01 7/350 2026-03-02 19:22 by zhukairuo
[考博] 博士自荐 +4 kkluvs 2026-02-28 5/250 2026-03-02 19:19 by 轻松不少随
[考研] 290分材料工程085601求调剂 数二英一 +4 llx0610 2026-03-02 4/200 2026-03-02 19:19 by caszguilin
[考研] 一志愿中石油(华东)本科齐鲁工业大学 +3 石能伟 2026-03-02 3/150 2026-03-02 18:54 by caszguilin
[考研] 290求调剂 +6 ErMiao1020 2026-03-02 6/300 2026-03-02 18:14 by lature00
[考研] 299求调剂 +4 kkcoco25 2026-03-02 4/200 2026-03-02 18:04 by barlinike
[考研] 清华大学 材料与化工 353分求调剂 +4 awaystay 2026-03-02 5/250 2026-03-02 17:33 by yeahyou
[考研] 一志愿华中科技大学,化学专业344分,求调剂 +3 邢xing1 2026-03-02 3/150 2026-03-02 17:32 by houyaoxu
[考研] 考研复试调剂,过国家线的同学都可报名 +4 黑!在干嘛 2026-02-28 5/250 2026-03-02 16:48 by chock1337
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 14/700 2026-03-02 16:09 by 今夏不夏
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 338求调剂 +3 18162027187 2026-03-02 3/150 2026-03-02 13:12 by houyaoxu
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师!
[考研] 0857调剂 +4 一ll半 2026-02-28 5/250 2026-03-02 02:33 by 908055542
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
信息提示
请填处理意见