foam-extend求解稳态问题,如何才能使最终解残差下降



  • 各位朋友,最近我在计算一个稳态的Laval喷管的稳态问题。可是最后的马赫盘的位置一直在非常微小幅度的振荡,导致最后的残差一直居高不下。我尝试过如下的几种方法:

    1。 waveTransmissive改为稳态适用版本;
    2。 减少co数和松弛因子。
    3。 加密网格。

    我用的是kOmegaSST, 求解器比较小众,是transonicMRFDyMFoam. 请问各位在稳定稳态求解方面有什么好的经验?
    谢谢各位!



  • 为了方便,我把我的fvScheme贴在这儿:

    ddtSchemes
    {
        default          none;//EulerLocal physDeltaT CoDeltaT;
        ddt(rho)       EulerLocal physDeltaT CoDeltaT;
        ddt(rhoU)      EulerLocal physDeltaT CoDeltaT;
        ddt(rhoE)      EulerLocal physDeltaT CoDeltaT;
        ddt(rho,omega) EulerLocal physDeltaT CoDeltaT;
        ddt(rho,k)     EulerLocal physDeltaT CoDeltaT;
        p              steadyState;   // Needed for waveTransmissive;
        phiMin           steadyState;   // Needed for waveTransmissive;
        phiMax          steadyState;    // Needed for waveTransmissive;
        
    }
    
    gradSchemes
    {
        default         Gauss linear;
        grad(p)         Gauss linear;
        grad(rhoSlope)  Gauss linear;
        grad(pSlope)    Gauss linear;
        grad(USlope)    Gauss linear;
        grad(TSlope)    Gauss linear;
        snGradCorr(U)   Gauss linear;
        snGradCorr(h)   Gauss linear;
        snGradCorr(omega) Gauss upwind phi;
        snGradCorr(k)   Gauss upwind phi;
        grad(omega)     Gauss upwind phi;
        grad(k)         Gauss upwind phi;
    }
    
    divSchemes
    {
        default                         none;
        div(phi,U)                      Gauss linearUpwind grad(U);
        div(phi,k)                      Gauss QUICK;//linearUpwind;
        div(tau)                        Gauss linear;
        div((tau&U))                    Gauss linear;
        div(phi,omega)                  Gauss QUICK;//linearUpwind;
    
    }
    
    laplacianSchemes
    {
        default         Gauss linear limited  0.33;
    }
    
    interpolationSchemes
    {
        default               none;
        interpolate(rho) linear corrected;
    }
    
    snGradSchemes
    {
        default         orthogonal;
    }
    
    fluxRequired
    {
        default         no;
    }
    
    


  • 最近看了几个extend的问题,暂时没用过extend的求解器...


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2