24小时热门版块排行榜    

查看: 3039  |  回复: 8

ytyyutianyun

木虫 (小有名气)

[求助] 介质为空气的圆柱绕流无法生成卡门涡街 已有1人参与

在做3d前,我先考虑2d没想到不像水,空气工况下,2d情况也无法获得卡门涡街,请虫友门传授一下经验,谢谢

我的工况是 300X300X1000mm, 圆柱直径1mm(没法改,我们要求就是这样小)。

画好网格后,用如下参数

我设置fluent参数为:2dpp  材料:空气。 来流速度为1.33m/s。

solver: unsteady--2nd-order-implicity
under-relaxation: pressure 0.1
residual convergence-criteria 1.0e-18 1.0e-18 1.0e-18
Solution controls 为: SIMPLEC, Pressure : second order, Momentum: first order upwind

时间步长一开始:0.001, 算了1100步,之后时间步长改为:5.0e-5, 算了1100步
得到如下图,






附件是我生成的cas
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 2dcross.7z
  • 2013-12-12 15:19:06, 165.71 K

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

onesupeng

金虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
ytyyutianyun: 金币+10, ★★★★★最佳答案, 谢谢你的热心帮助,非常感谢 2013-12-13 09:19:33
物理上说,雷诺数决定流动现象,你要有涡脱落,雷诺数需要Re>48~50的样子。由于商业软件一般耗散性都比较大,所以一般需要更大Re,才有涡脱落现象。

在圆柱绕流中,Re=rho*U*D/\mu,其中rho是流体密度,U是来流速度,\mu是流体的粘性系数。所以不要什么空气不空气,根据定义,Re一样,无论你用的什么介质,在低速流动中结果都是一致的。

另一方面,由于商业软件耗散比较大,所以如果你想看到更真实的流动现象,Solution controls里面的Momentum选项用 second order upwind 或者second order central scheme(如果有的话)。因为upwind本身是耗散性比较强的格式,它有助于数值稳定性,但是有时候却因为耗散性太大,使得本来是非定常的,都搞成定长流动了。
长期招收博士生,参见http://fsl-unsw.com
2楼2013-12-12 17:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytyyutianyun

木虫 (小有名气)

引用回帖:
2楼: Originally posted by onesupeng at 2013-12-12 17:00:26
物理上说,雷诺数决定流动现象,你要有涡脱落,雷诺数需要Re>48~50的样子。由于商业软件一般耗散性都比较大,所以一般需要更大Re,才有涡脱落现象。

在圆柱绕流中,Re=rho*U*D/\mu,其中rho是流体密度,U是来 ...

谢谢,也就是说,不是看材质,看粘性
比如我这个工况

Re=1.33*1e-3*1.225/1.7894e-05=91.06符号绕流条件

1.如果三维,可否选大涡模拟(Bounded Central Differencing)而不是second order upwind

2.在麻烦问一下,时间步长这样可行否
时间步长一开始:0.001, 算了1100步,之后时间步长改为:5.0e-5, 算了1100步
主要文献实验测到的200Hz, 因此 1/200=5e-3, 然后根在除以100.所以5.0e-5

谢谢
3楼2013-12-12 17:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytyyutianyun

木虫 (小有名气)

引用回帖:
2楼: Originally posted by onesupeng at 2013-12-12 17:00:26
物理上说,雷诺数决定流动现象,你要有涡脱落,雷诺数需要Re>48~50的样子。由于商业软件一般耗散性都比较大,所以一般需要更大Re,才有涡脱落现象。

在圆柱绕流中,Re=rho*U*D/\mu,其中rho是流体密度,U是来 ...

改了second order upwind还是不行,会不会我的网格问题

对长方体300X1000mm, 对短边分30,长边100,导致间隔为10mm
而圆柱直径1mm,分60,导致间隔约等于0.052mm。

会不会差别太大?
4楼2013-12-12 18:50:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

引用回帖:
3楼: Originally posted by ytyyutianyun at 2013-12-12 17:33:05
谢谢,也就是说,不是看材质,看粘性
比如我这个工况

Re=1.33*1e-3*1.225/1.7894e-05=91.06符号绕流条件

1.如果三维,可否选大涡模拟(Bounded Central Differencing)而不是second order upwind

2.在麻 ...

关于网格和dt。按照经验,壁面附近网格的高度(沿法向)大约是D/sqrt(Re)/n (n=3~5), sqrt是开方。

dt你也用无量纲的dt看看,即大约dt=0.001*D/U左右,这个你看时间离散几阶的了。时间大约模拟200*D/U的样子,才会得到比较稳定的结果
长期招收博士生,参见http://fsl-unsw.com
5楼2013-12-13 04:46:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytyyutianyun

木虫 (小有名气)

引用回帖:
5楼: Originally posted by onesupeng at 2013-12-13 04:46:40
关于网格和dt。按照经验,壁面附近网格的高度(沿法向)大约是D/sqrt(Re)/n (n=3~5), sqrt是开方。

dt你也用无量纲的dt看看,即大约dt=0.001*D/U左右,这个你看时间离散几阶的了。时间大约模拟200*D/U的样子,才 ...

1.壁面附近网格的高度是不是指边界层网格设置。

2.了解了,不过我的工况dt=1e-6, 总时间=0.2,要算20万步,计算机要算1年。因此可不可以分时间步长,比如时间步长一开始:0.001, 算了2000步,之后时间步长改为:5.0e-5, 算2100步,这样可行否?
6楼2013-12-13 09:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytyyutianyun

木虫 (小有名气)

引用回帖:
5楼: Originally posted by onesupeng at 2013-12-13 04:46:40
关于网格和dt。按照经验,壁面附近网格的高度(沿法向)大约是D/sqrt(Re)/n (n=3~5), sqrt是开方。

dt你也用无量纲的dt看看,即大约dt=0.001*D/U左右,这个你看时间离散几阶的了。时间大约模拟200*D/U的样子,才 ...

对了,有可能我误解了
引用回帖:
时间大约模拟200*D/U的样子

此时间是总时间呢还是我说的先用较大的时间步长算到的时间,谢谢
7楼2013-12-13 09:27:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

引用回帖:
7楼: Originally posted by ytyyutianyun at 2013-12-13 09:27:12
对了,有可能我误解了


此时间是总时间呢还是我说的先用较大的时间步长算到的时间,谢谢...

总时间
长期招收博士生,参见http://fsl-unsw.com
8楼2013-12-14 08:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ytyyutianyun

木虫 (小有名气)

引用回帖:
8楼: Originally posted by onesupeng at 2013-12-14 08:49:27
总时间...

谢谢
9楼2013-12-14 16:05:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ytyyutianyun 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +6 小熊joy 2026-03-14 6/300 2026-03-16 12:59 by Iveryant
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 311求调剂 +6 冬十三 2026-03-15 6/300 2026-03-16 08:00 by wang_dand
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
请填处理意见