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



  • kepsilon模型中k的选择和epsilon的选择根据公式算出,
    0_1476111011982_upload-aaccde47-b559-4023-bc98-ee7f8100e9bf
    但计算结果一直显示发散,如下图
    0_1476111104742_upload-6e7b27bb-d32b-492b-97c2-b6c212a2815a
    我认为是由于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;
        }
    }
    


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



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



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



  • 多谢您的回复


Log in to reply