中性大气环境湍流动能的自保持 | 附有算例下载

coolhhh

@李东岳李老师，我尝试按照自己的理解，用fluent尝试算了Yang2009算例，参照OpenFOAM的设置，结果能对上，但和原论文结果有点区别：

按照Yang2009设置，顶面设置为slip wall，结果基本能自保持，但接近顶部的k无法自保持，猜测是顶面边界条件设置问题。

fluent中slip wall设置如下，把壁面剪切应力设置为0。

把顶面设置为symmetry，结果就基本一样。实在不知道为何OpenFOAM算出来k衰减那么多

coolhhh

似乎找到原因了，是要把算例0文件夹中的U、k、epsilon中用codedFixedValue编写的边界name命名为不同名称，基于李老师算例文件修改后如下，用OpenFOAM-9直接Allrun
Yang2009OFCase.tar

    inlet
    {
        type            codedFixedValue;
        value           uniform 0.1; //default value
        name            linearTBC1; //name of new BC type

修改后的算例结果如下，近壁面有点偏差，但总体是保持住了

李东岳

完美解决！

太感谢了！这个小bug让我卡了一个月。太无语了

。。。。。

李东岳

@coolhhh 在中性大气环境湍流动能的自保持 | 附有算例下载中说：

把顶面设置为symmetry，结果就基本一样

另外，symmetry与slip的差异极小。在openfoam这面，symmetry是一种几何性的限定，所有的场各种场都会严格的对称。但是slip只是速度会这样，其他场不会，比如HbyA。但是这个是fluent的结果，那fluent为什么这样，就不知道了。

coolhhh

@李东岳或许是我这边fluent有某个设置错了，杨老师论文中顶面用的slip结果不会出现k衰减。我再试了fluent中顶面设置no-slip，结果是比较符合预期，接近顶面壁面平均速度减小，湍动能增大

SHUKK

@李东岳李老师，我最近尝试复现Yang2017的kOmegaSST，按照文献的设置和参数设置，k都无法保持稳定。不知道老师您是否复现出来了。

李东岳

@SHUKK kOmegaSST没有尝试。你可以把你的算例打包发上来我给你看看

SHUKK

@李东岳好的，谢谢李老师！y2017kOmegaSST.zip

chszkc

@Vortex 前辈您好，请问可以向您请教一下LESGO温度分层算例的相关设置吗？

疏影横斜水清浅

请问有没有人做大涡模拟的中性大气环境湍动能自保持，大涡模拟的有没有参考文献？

疏影横斜水清浅

@李东岳老师您好，请问Yang2009中的垂向网格设置，第一层网格必须很小吗，和y+差不多。如果将垂向网格的尺寸变大，算例会发散，有人遇到过这种情况吗？

李东岳

@疏影横斜水清浅

如果将垂向网格的尺寸变大，算例会发散，有人遇到过这种情况吗？

我提供的这个算例么

疏影横斜水清浅

@李东岳是的

李东岳

@疏影横斜水清浅明天我试一下我记得我放大过很大的网格也能算

疏影横斜水清浅

@李东岳测试的网格如下，网格尺寸比较大

convertToMeters 1;

vertices
(
    (0 0 0)//0
    (5000 0 0)//1
    (5000 500 0)//2
    (0 500 0)//3
    (0 0 1)//4
    (5000 0 1)//5
    (5000 500 1)//6
    (0 500 1)//7
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (500 50 1) simpleGrading (1 10 1)
);

edges
(
);

boundary
(
    inlet
    {
        type       patch;
        faces
        (
             (3 0 4 7)
        );
    }
    outlet
    {
        type       patch;
        faces
        (
             (1 2 6 5)
        );
    }
    topWall
    {
        type       patch;
        faces
        (
             (2 3 7 6)
        );
    }
    flatWall1
    {
        type       wall;
        faces
        (
             (0 1 5 4)
        );
    }
    front
    {
        type empty;
        faces
        (
             (4 5 6 7)
        );
    }
    back
    {
        type empty;
        faces
        (
             (3 2 1 0)
        );
    }
);

mergePatchPairs
(
);

李东岳

你这个不是网格尺寸变大，是网格尺度变大。网格尺寸放大，是计算域不变，网格变粗。

yang本身是小尺度，你变成了大气环境大尺度，类似工业放大。或许会有问题。

疏影横斜水清浅

@李东岳这样呀，那我明白了，HW07那个可以放大计算域，但是那个的缺点是湍动能不随高度变化，所以我想着Yang09的也可以放大的，试了试，每次都发散

李东岳

@疏影横斜水清浅是，HW07本身计算雨很大，可以放大计算域，Yang2009这个我记得我尝试也不行，现在有点记不清了。

疏影横斜水清浅

@李东岳我明白了，HW07的湍动能和耗散率都是定值，不随高度变化，而Yang09的乘了一个系数，这个系数是个开平方数，而里面的数随着高度的增加会减小，当高度超过一定高度（与粗糙度系数有关）时，这个数就变成了负数，那么开方运算就不可以了，程序报错。

李东岳

@疏影横斜水清浅厉害厉害！非常厉害！从算法上解决了问题。有没有什么解决办法呢。

CFD中文网

中性大气环境湍流动能的自保持 | 附有算例下载

完美解决！