Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam

[流散然]

@random_ran
感谢您这么长时间以来的耐心解答与帮助，问题我已解决，所有的症结都出在边界条件上，这一方法适用于矩形网格和圆形网格。具体的细节我不太方便透露，如果您对此有兴趣的话，我们可以加下联系方式（邮箱或微信等皆可）具体的聊下
最后再次对您一直以来的帮助致以诚挚的感谢，如果日后我因为算力不足转向DES，到时还请您多加指教！

[红豆沙]

@random_ran 你用fluent算的，，参考值是否设置正确，，这个对阻力系数有直接影响

[红豆沙]

@红豆沙 哈哈，看了老半天，讨论好详细，，

[hoversoar]

@random_ran 请问你圆形计算域的网格也是用ICEM画出来的吗？
我计算的Re大概在15w左右，边界层非常小，用meshing的结构化貌似不太成功。所以现在用了非结构化网格，但是结果也是偏大。还望请教一下！
谢谢！

[hoversoar]

您好！
我还有一些疑问想请教一下：
关于入口边界的K和nut值。

1. 我注意到一开始这个帖子你放的是在cylinder k=0.24，请问这个k您是怎么取的?
2. nut的值，这个帖子用了壁面函数。
   在另一个帖子里 你用的是fixedvalue 0. 没有使用壁面函数。
   请问在边界层足够小的情况下，是不是可以直接采取fixedvalue 0. 这样的设定？

[random_ran]

A video to show flow around circular cylinder in real life:
https://youtu.be/0ThQ_nD97hY?t=306

[好好科研]

@random_ran 您好！我最近也在做三维圆柱扰流的相关工作，我先开始做的是长方形流域的二维圆柱绕流，感觉还能差不多吻合升阻力系数，但是三维数值模拟很慢，看到您的帖子，开始转移到圆形流域，发现升阻力系数会随着流域大小的变化而变化，三维流域的 话，就算的很慢很慢，差不多要12个小时才要算1个时间单位。如果方便的话，可否把你的算例或者网格让我参考一下~x谢谢！！！这是我的邮箱~3310792254@qq.com~qq.com~祝好！

[炎冰灬灬]

@random_ran 你好。看了你的升阻力系数的图，我最近也出现了同样的问题。我算的是二维的刚性翼型的绕流问题，用的是SST k-w模型，pisoFoam求解器，算出来的升阻力系数的波动也是非常的大，不知道问题出在哪里，你能指点一二吗

[yfleung]

@random_ran 您好，请问您这个圆形流域的网格是设置的维的吗？我最近也在做一个雷诺数是3900的算例，也是采用的圆形计算网格，流域的直径设置为50D，圆的周长（paiD）方向是480个网格，半径方向是400个网格，计算的Y+为0.06，但是用ICEM划分的是二维的，在做雷诺数是200时候阻力系数Cd和文献里面差距不大，但是在雷诺数是3900时，Cd是0.5的样子，和文献里面差距很大。对了，我是用的直接数植模拟，没有用湍流模型，但是3900雷诺数情况下的阻力系数差特别大。不知道是什么原因

[yfleung]

@bingningmeng45 你好，请问您做圆柱绕流划分网格的话，要满足什么样的条件才是DNS啊

[落日酒店]

@yfleung 您好，我最近在做les模拟三维圆柱绕流，也出现了cd值只有0.6的问题，请问您解决这个问题了吗

[Jean]

@random_ran
老师您好，看了你的不少评论受益匪浅，我目前也在做圆柱绕流仿真，雷诺数为100000，但阻力系数升力系数一直不理想，自己不知道咋调试，能否帮忙看看问题在哪？

具体文件如下
0文件
k nut nuTilda p U

constant文件
transportProperties turbulenceProperties

system文件
controlDict createPatchDict decomposeParDict fvSchemes fvSolution

网格
中心区域为4D*4D

[李东岳]

@Jean 粗略一看，你需要多算一段时间看看

[Jean]

@李东岳 李老师您好！目前大概跑了20天左右，具体胜利阻力系数如下

但是升力系数似乎有些问题，不知道问题可能在哪，不知道要怎样处理，您能否帮忙看看，非常感谢🙏🙏🙏🙏🙏🙏🙏

[李东岳]

跑了20天？多少万网格？

[Jean]

@李东岳 李老师您好，网格4140024，y+=1.0

[李东岳]

400万网格为啥跑了20天，你多少核。

400万不行，debug不动，用的什么湍流模型，流场看起来稳么，有越界么

[Jean]

@李东岳 用了112核，目前case还可以跑，没有报debug，湍流模型用了DDES，是放到超算平台跑的，由于文件下载挺麻烦的，所以目前主要是看了这个力系数，关于云图啥的还没有查看

[李东岳]

圆柱扰流？你这个cl好像确实要起飞的节奏。。也降不下来。我不知道啥情况，你要不要看看云图有没有错误。你这跑20天太难debug了。你随便调整一个设置都要算20天。太难搞。按道理来说圆柱扰流这种特别简单不应该如此。

[李东岳]

不对，我把cd跟cl这两个线看混了。我刚才发的帖子作废。我觉得你可能需要继续多算一段时间