ABL入口设置

sxz0823

是不是什么地方没有设置对呢，请教各位大佬！

sxz0823

这是实验值

allanZhong

似乎可以将initABL设置为false, 然后在value下面给出你的这一列数值，大概如下：

inlet
{
    type    atmBoundaryLayer;
    initABL false;
    value   nonuniform List<scalar> 
        5          // 数据点数量
        (
            1.0
            1.2
            1.1
            0.9
            1.3
        );
}

仅供参考

ZhangXF

是不是因为u取值的问题，能不能先拟合出符合实验数据的u，k，z0这三个关键参数，用这些拟合的参数来给openfoam算风廓线；
第二个问题也许是因为实验风廓线不是完美的对数廓线，导致不同参考高度和速度算出的风廓线不一样。

另外，我之前给ABL廓线的时候是用了codedFixedValue边界条件自己编了一段廓线生成代码给k、w和e用，你可以参考一下（我的环境是2.4.x）：

    xMin
    {
        type            codedFixedValue;
        value           $internalField; //default value
        redirectType    kAtm; //name of new BC type
        code
        #{
            const fvPatch& inletPatch = this->patch();
            scalarField& vf = *this; 

            forAll(vf, i)
            { 
                scalar z = inletPatch.Cf()[i].z();
                scalar z_Ground = 8;
                scalar Cmu = 0.11;
                scalar uStar = 0.618;
                scalar z0 = 0.03;
                scalar C1 = -0.1693041080673375; // C1和C2根据实验数据的k拟合得到
                scalar C2 = 3.668086709040146;
                vf[i] = sqr(uStar)/sqrt(Cmu) * sqrt(C1*log((z - z_Ground + z0)/z0) + C2);
            }
        #};
    }

w和e的边界条件代码同理

chenboyao

在OpenFOAM的算例中就有使用ABL的验证算例，OpenFOAM的Tutorial/verificationAndValidation/atmosphericModels/atmDownstreamDevelopment

CFD中文网

ABL入口设置