openfoam添加湍流入口方法?

cccrrryyy

@dyb 把两个结合起来肯定是有办法的，就相当于自己写一个边界条件了。但是我感觉这样需要的代码工作量可能很大，因为两部分代码里面都包含一些实际上你不想要管的东西，比如什么构造函数啊之类的，还要去好好的想每个类里面哪些数据和函数应该做成protected，哪些做成public之类的，种种问题。除非你是一个编程高手很熟悉这些东西，不然真的可能工作量太大了。实际上对inflowGenerator来说你就是需要给一个速度的profile，我觉得这个profile你按照之前那个post里面的方法去给就好了，根据位置去确定每个cell的速度是多少的时候你用atmBoundaryLayerInletVelocity里面提到的那个公式即可。

dyb

@cccrrryyy 谢谢老师

dyb

@星星星星晴谢谢老师，我按照老师的方法已经弄出来了

zxx

@cccrrryyy https:Eugene de Villiers的扰动在这里//github.com/wyldckat/perturbU

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
e595f49f-347a-4180-82e2-b6295cd45a15-图片.png

zzzzz7777

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

疏影横斜水清浅

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

R      uniform (0.8 0 0 0.4 0 0.4);

zzzzz7777

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

CFD中文网

openfoam添加湍流入口方法?