使用kOmegaSSTIDDES湍流模型计算大气边界层风场

simulationType      LES;

LES
{
    LESModel        kOmegaSSTIDDES;

    turbulence      on;

    printCoeffs     on;

    delta           IDDESDelta;

    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
    decayControl    true;//false;
    kInf            0;
    omegaInf        0;
    useSigma        false;
    kappa           0.41;
    CDESkom         0.82;
    CDESkeps        0.6;
    Ctrans          60;
    Cdt1            20;
    Cdt2            3;
    Cl              5;
    Ct              1.87;
    fe              true;
}

风速入口为对数廓线

            scalar kappa = 0.41;
            scalar Uref  = 2.0;
            scalar Zref  = 10.0;
            scalar z0    = 0.001;
            scalar zd    = 0.0;

            scalar Ustar = kappa*Uref/log((Zref+z0)/z0)；
            scalar z  = inletPatch.Cf()[i].z();
            Uf[i].x() = Ustar/kappa*log(((z-zd) + z0)/z0);
            Uf[i].y() = 0.0; 
            Uf[i].z() = 0.0;

结果为
4810b34c-ee97-480d-b96b-a90b3c278fbf-图片.png
其中，9个位置分别为10m、200m、400m、500m、600m、700m、800m、900m、990m（计算域流向1000m，垂向100m，展向100m），网格划分(250 25 50) z是垂向，膨胀率 (1 1 5)。
请问，上部分的速度随高度基本不变，是什么原因？怎么才可以使得风速随高度呈对数形态。

CFD中文网

疏影横斜水清浅