圆柱绕流高雷诺数10e5

Afr1yne

把fvSolution中压力求解器改成PCG试试

lwf

谢谢，我试了一下，还是不行，我网格用ICEM画的，网格质量很高，就是不知道哪出的问题。还请大家指教一下
fvsolution

solvers
{ 
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0.05;
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    "(U|k|epsilon)"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0.1;
    }

    "(U|k|epsilon)Final"
    {
        $U;
        relTol          0;
    }
}

PIMPLE
{
    nNonOrthogonalCorrectors 0;
    nCorrectors         2;
}

fvscheme

ddtSchemes
{
    default         backWard;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;

    div(phi,U)      Gauss linearUpwind grad(U);

    turbulence      Gauss upwind;
    div(phi,k)      $turbulence;
    div(phi,epsilon) $turbulence;
    div(phi,nuTilda) $turbulence;
    div(phi,R)      $turbulence;
    div(R)          Gauss linear;

    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

李东岳

时间步长缩小

lwf

感谢李老师，我再试试

火丙三金

@lwf 怎么样兄弟

lwf

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

RANS的结果：

火丙三金

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

李东岳

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

lwf

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

火丙三金

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

xpqiu

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

估计你的网格是边界层首层厚度很小的那种吧，那么你用 standard k epsilon 这个湍流模型是不合适的，建议使用 kOmegaSST。
抛开 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

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

CFD中文网

圆柱绕流 高雷诺数10e5

圆柱绕流高雷诺数10e5