kOmegaSSTIDDES湍流模型风速廓线

[疏影横斜水清浅]

使用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;  

结果为

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

[coolhhh]

@疏影横斜水清浅 kOmegaSSTIDDES没用过，但用过WALE计算，一般入口只是平均风的话，都会保持的很好。看起来网格非常粗糙，不确定是不是全尺寸网格问题导致的，可以做个缩尺的模型验证一下。还有如果算全尺寸的话，计算域底面的yPlus应该量级很大，这也是我一直的疑惑点，目前的壁面函数一般要求使用范围是1000以内，壁面函数能否直接用于全尺寸模型

[李东岳]

正如楼上大佬所说，网格粗糙，看起来像是网格的问题。边界层没有捕获到精确的速度梯度。

[疏影横斜水清浅]

@coolhhh 是的，yplus很大，三四千，这个是我在测试的，我是想做全尺寸的，所以想用的是IDDES，IDDES是否还要考虑yplus，我不是很清楚。测试算例，所以使用的网格很粗。

[coolhhh]

@疏影横斜水清浅 算全尺寸，如果只是关注风速，对结果影响不大。现在很多算城市风环境的，基本yPlus都是不满足的。但如果算全尺寸结构表面风压，这个我没算过不确定影响多大。另外算全尺寸入口加的脉动风，因为网格太大脉动基本会被过滤很多，结果也会有点差

[疏影横斜水清浅]

@李东岳 我先将网格细化做一下测试，排除或者确认一下网格的原因。