24小时热门版块排行榜    

查看: 3564  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金啊基金 +3 longfie172 2026-08-20 3/150 2026-08-21 00:04 by wuchongjun
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 言之凿凿,明天官网会通知放榜? +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见