openfoam添加湍流入口方法?

尚善若水

@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

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

疏影横斜水清浅

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

zzzzz7777

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

sxz0823

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

CFD中文网

openfoam添加湍流入口方法?