有遇见过simpleFoam达不到稳态的情况吗?



  • 最近再用snappyHexMesh画网格,网格画完了checkMesh也通过了,但是simpleFoam达不到稳态。用了几个probe来检查,发现速度一直在震荡。 同样的system下面的设置,之前跑其他的算例都没问题的。
    我贴一下我的fvSchemes

    
    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2.3.1                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       dictionary;
        location    "system";
        object      fvSchemes;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    ddtSchemes
    {
        default         steadyState;
    }
    
    gradSchemes
    {
        default         cellLimited leastSquares 1; //Gauss linear;
    }
    
    divSchemes
    {
        default         none;
        div(phi,U)      bounded Gauss linearUpwindV cellLimited leastSquares 1;
        div(phi,k)      bounded Gauss linearUpwind celllimited leastSquares 1;
        div(phi,epsilon) bounded Gauss linearUpwind celllimited leastSquares 1;
        div(phi,R)      bounded Gauss linear corrected;
        div(R)          Gauss linear;
        div(phi,nuTilda) bounded Gauss upwind;
        div((nuEff*dev2(T(grad(U))))) Gauss linear;
    }
    
    laplacianSchemes
    {
        default         Gauss linear corrected;
    }
    
    interpolationSchemes
    {
        default         linear;
    }
    
    snGradSchemes
    {
        default         corrected;
    }
    
    fluxRequired
    {
        default         no;
        p               ;
    }
    
    
    // ************************************************************************* //
    

  • 网格教授 OpenFOAM教授 管理员

    simpleFoam达不到稳态的情况吗?

    发现速度一直在震荡

    这是可能的。有些物理过程,用稳态模拟就是会一直震荡。不过,用snappyHexMesh做的网格,也有可能是网格的问题,snappyHexMesh做出来的网格某些部位质量很差。但如果几何不是特别复杂,一般不会影响到全局。

    另外,如果你要降低震荡,应该增加粘性,用低阶格式,比如迎风试试。



  • @东岳 好的谢谢东岳,我先用一阶迎风格式试试