openfoam添加湍流入口方法?

[cccrrryyy]

@zxx 感谢！

[浪迹天大]

我把代码移植到OF-上了，并附上了我的算例以供参考：
https://github.com/ZhangYanTJU/boundaryConditions

[浪迹天大]

@浪迹天大 在 openfoam添加湍流入口方法? 中说：

我把代码移植到OF-上了，并附上了我的算例以供参考：
https://github.com/ZhangYanTJU/boundaryConditions

少打了一个字。。。OF-7

[dyb]

@浪迹天大 您好老师，有没有关于LeMOS的参数设置的计算方法，或者文献，感谢老师

[yhdthu]

我看到DFSEM入口有bug？

讨论地址：https://www.cfd-online.com/Forums/openfoam-solving/177711-turbulentdfseminlet.html

有问题的是11式，本来是个无量纲参数，根号出来其实有量纲，不知有人测试过没？

[尚善若水]

@cccrrryyy 有个问题想请教下您。我也是做射流的，我是平板声速射流，来流空气为超声速（1800m/s）。这种情况下，应该如何确定积分尺度呢？我的射流孔直径为2mm，试验壁面边界层为1.2mm。按照cfd online的介绍，0.22thicknessBL，和取射流口直径的的10%基本一致，不知道这样取对不对呢？另外我网格最小尺寸0.3mm0.3mm*0.3mm，取0.2mm积分尺度会不会使计算速度很慢呢。ps：最近测过Jiawei Wan博士的TinF边界，2千万网格在720processor上计算，计算速度很慢，所以有此担心。

[cccrrryyy]

@尚善若水 不知道你说的太慢，是慢到什么程度。积分尺度这个东西我自己也拿捏不准，关键是这个东西不一定对结果影响特别大。毕竟，它只是入口方法中的一个参数而已。

它和计算速度的关系我没有关注过，为什么它会对速度有影响呢。

[尚善若水]

@cccrrryyy 按道理来讲不应该有很大影响，但是我的计算确实变慢了很多，我在想是不是vorton太多了。我打算直接用turbulentInlet了。

[cccrrryyy]

@尚善若水 边界条件里面有inject vorton的时间我记得，应该没有很久吧。我觉得对计算速度应该没啥影响。

[sgggggg]

@dyb 你好我可以请问一下这个边界的几个参数的设置方法么，就是LField、RField和refField，我拜读过你的硕士论文但是里面我看没有做细致介绍，所以可以请教一下么，我也是搞风电现在也需要这个和les做耦合

[dyb]

这个参数设置，有两种，一种是均匀，一种是自定义风切入口风速，

 type            decayingTurbulenceInflowGenerator;
        direction       1;
        LField          nonuniform List<scalar>  
(
0.025
0.025


...

9.75
9.75
);
        RField          uniform (1 0 0 1 0 1   );
        refField        nonuniform List<vector> 

(
( 2.4495403099953 0 0 )
( 2.4495403099953 0 0 )

....

( 5.48123361668986 0 0 )
( 5.48123361668986 0 0 )

        );

[sgggggg]

@dyb 谢谢你的回答，我还有个问题你这边L一般就是取特征长度的10%来做是么，这个是按边界网格为基准计算么

[疏影横斜水清浅]

请问在使用turbulentDFSEMInlet入口类型计算过程中，随着计算的进行，后半部分的风场慢慢的变得不是对数廓线了，请问是什么原因，有没有大佬知道。
我的初始边界信息：
计算域：xyz为1011（m）
U：

inlet
    {
        type            turbulentDFSEMInlet;
        delta           1;
        U
        {
            type        mappedFile;
            mapMethod   nearest;
        }
        R               uniform (0.8 0 0 0.4 0 0.4);
        L               uniform 0.004;
        d               1;
        nCellPerEddy    1;
        scale           1;
        value           $internalField;
    }
    outlet
    {
        type            advective;
        value           $internalField;
    }
    "(topWall|front|back)"
    {
        type            zeroGradient;
    }
    "(bottomWall|Model)"
    {
        type            noSlip;
    }

输入文件：
point：

(
( 0 0 0.01 )
( 0 0 0.05 )
( 0 0 0.1 )
( 0 0 0.15 )
( 0 0 0.2 )
( 0 0 0.25 )
( 0 0 0.3 )
( 0 0 0.35 )
( 0 0 0.4 )
( 0 0 0.45 )
( 0 0 0.5 )
( 0 0 0.55 )
( 0 0 0.6 )
( 0 0 0.65 )
( 0 0 0.7 )
( 0 0 0.75 )
( 0 0 0.8 )
( 0 0 0.85 )
( 0 0 0.9 )
( 0 0 0.95 )
( 0 0 0.99 )
)

U：

(
( 0 0 0 )
( 1.048097372 0 0 )
( 1.499488341 0 0 )
( 1.763535131 0 0 )
( 1.95087931 0 0 )
( 2.096194744 0 0 )
( 2.2149261 0 0 )
( 2.31531203 0 0 )
( 2.402270279 0 0 )
( 2.47897289 0 0 )
( 2.547585713 0 0 )
( 2.609653562 0 0 )
( 2.666317069 0 0 )
( 2.718442442 0 0 )
( 2.766702999 0 0 )
( 2.811632503 0 0 )
( 2.853661247 0 0 )
( 2.893141184 0 0 )
( 2.930363858 0 0 )
( 2.965573488 0 0 )
( 2.992431707 0 0 )
)

计算3600s后，平均风速廓线，U1U2U3U4U5U6U7U8U9的位置为0.1，2，4，5，6，7，8，9，9.9

[zzzzz7777]

@疏影横斜水清浅 为啥你的输入文件没有雷诺应力的信息呀 可以讨论下

[疏影横斜水清浅]

@zzzzz7777 我给了一个相同的雷诺应力，不是廓线。

R      uniform (0.8 0 0 0.4 0 0.4);

[zzzzz7777]

@疏影横斜水清浅 这样给可以就不符合湍流边界层的规律了啊

[疏影横斜水清浅]

@zzzzz7777 大气边界层的雷诺应力应该怎么给，也有人说这个初始条件不重要，给的差不多后面也可以发展出正常的

[zzzzz7777]

我是用其他人算的数据差值给的 RANS算不准雷诺应力

[sxz0823]

@sgggggg 硕士论文我也想拜读一下，可以加个联系方式吗2132306535