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

[qjh888]

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

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

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

[qjh888]

为了方便，我把我的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的求解器...