关于simpleFoam中kEpsilon中k和epsilon的设定

[vivian]

kepsilon模型中k的选择和epsilon的选择根据公式算出，

但计算结果一直显示发散，如下图

我认为是由于k和epsilon的初始值选择不恰当造成的。不知道大家怎么选择这两个值？
（计算背景：来流0.1m/s的圆柱绕流，雷诺数设置为2000，各边界设置如下：）
压力设置：

internalField   uniform 0;

boundaryField
{
    INLET
    {
        type            freestreamPressure;
       
    }

    OUTLET
    {
        type            freestreamPressure;
      
    }
     TOP
    {
        type            zeroGradient;
        
    }
    BOTTOM
    {
       type            zeroGradient;
    }
    CYLI
    {
        type            slip;
    }
    frontAndBack
    {
        type            empty;
    }
}

速度设置：

internalField   uniform (0.1 0 0);

boundaryField
{
   INLET
    {        
           type         freestream;
        freestreamValue        uniform    (0.1 0 0);
    }

    OUTLET
    {
       type         freestream;
       freestreamValue        uniform    (0.1 0 0);
    }
     TOP
    {
       type             zeroGradient;
    }
    BOTTOM
    {
        type            zeroGradient;
    }
    CYLI
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    frontAndBack
    {
        type            empty;
    }
}

nut

internalField   uniform 0;

boundaryField
{
    INLET
    {
        type            freestream;
        freestreamValue uniform 0;
    }

    OUTLET
    {
        type            freestream;
        freestreamValue uniform 0;
    }
    
    TOP
    {
        type            zeroGradient;
        
    }
    BOTTOM
    {
       type            zeroGradient;
    }

    CYLI
    {
        type            kqRWallFunction;
        value           uniform 0;
    }

    frontAndBack
    {
        type            empty;
    }
}

[chpjz0391]

入口的特征长度多少，发散也可能由其他引起。松弛因子等等。

[vivian]

谢谢您的回复。我又查了一下，把特征长度L用我的进口长度乘以0.05可以运行了。之前我是看一个视频上选择乘以0.08.您说的松弛因子我一般取为0.8。请问这个的选取有什么规则吗？

[chpjz0391]

@vivian 特征长度有公式，圆形的特征长度就是其直径，长方形的好像是4倍的面积除以周长什么的我记不太清楚了。你自己搜一下。松弛因子是控制变量在每次迭代中的变化，可控制收敛速度，越接近1 收敛越快，但是收敛性相对较差。这个看你模拟了一般0.7左右就可以具体情况具体分析。举个例子如果p一直发散可以适当调低。等计算一段时间后再慢慢调高 加快收敛速度。

[vivian]

多谢您的回复