openfoam添加湍流入口方法?

尚善若水

@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添加湍流入口方法?