为了方便，我把我的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;
}
CFD中文网

qjh888
