圆柱绕流 高雷诺数10e5

[lwf]

@火丙三金 我改完后跑出来的结果感觉不对，没有涡脱落，fluent倒是能跑出来涡脱落，改到崩溃，看着不合理，水平不行...
LES的结果：

RANS的结果：

[火丙三金]

@lwf fluent可以分享一下嘛 我也来改一改！

[李东岳]

@lwf openfoam那个你要算充分长的时间，前期需要发展

[lwf]

@火丙三金 fluent的简单，把网格导进去设置一下就行，我把网格贴上来。ICM-11_10wan.zip
foam里面低雷诺数倒是正常

[火丙三金]

@lwf 我来试着调一调 兄弟没试下圆形区域嘛

[xpqiu]

你这个算例有几处不太合理的地方：

1. 估计你的网格是边界层首层厚度很小的那种吧，那么你用 standard k epsilon 这个湍流模型是不合适的，建议使用 kOmegaSST。
2. 抛开 k-epsilon 模型的合理性不谈，你这个 k epsilon nut 的 Inlet 边界条件设置也是不合理的。首先，你的入口速度是 9m/s，k 的 Inlet 值为 0.375，意味着入口湍流度为约 5.5%，这个值是偏大的。而 epsilon 的如何值为 0.07，意味着入口的 nut 值为 0.09*0.375^2/0.07=0.1808，超过流体粘度的 10000倍，这个是严重偏大的值。

[lwf]

@李东岳 李老师，LES和RANS都需要前期充分发展时间吗，怎么fluent里面很快就能出现涡脱落？

[lwf]

@xpqiu 非常感谢提醒，请问有可以计算模型边界初始值的参考资料吗

[lwf]

@火丙三金 计算域圆形矩形不大吧，圆形也试过，现在就是我设置的问题，模型和边界条件不合理

[火丙三金]

@lwf OK 我拿fluent试试

[李东岳]

@lwf 都需要。你可以把你的RANS的nut截图发上来看一下。是URANS么

[lwf]

@李东岳 李老师，因为不知道怎么设初始值合适，所以我当时是参照foam里面其他案例设置的nut，用的piso，这应该属于URANS吧，我看论坛以前讨论过RANS是用稳态模拟，URANS是瞬态模拟。

nut：

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    INLET
    {
        type            calculated;
        value           uniform 0;
    }

    OUTLET
    {
        type            calculated;
        value           uniform 0;
    }

    WALL
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    CYLINDER
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    frontAndBackPlanes
    {
        type            empty;
    }
}

[李东岳]

nut的瞬态云图截图看一下

[lwf]

@李东岳 好的
nut的截图：

[lwf]

这是之前跑的，我后面把湍流模型和初始边界条件改掉，然后跑个时间长点的试试

[李东岳]

形态倒是问题不大，数值0.3是太大了，湍流模型预测的nut太大，不对