LES直流槽道边界层模拟，如何得到正则化速度u+以及正则化坐标y+？

李东岳

@falltime nice！一针见血。厉害厉害！

coolhhh

前期计算的$\left \langle u_\tau \right \rangle ^+$会出现参差的现象，原来是之前在of8计算复制了一个算例的wallShearStress文件，文件中设置executeControl writeTime; 这导致了最终计算的平均只有writeTime次数，而不是每个timeStep的结果。参考of2206一个算例的设置就没问题了

of2206中算例设置

wallShearStress
{
    type                wallShearStress;
    libs                (fieldFunctionObjects);
    patches             ( bottomWall topWall );
    writePrecision      10;
    writeToFile         yes;
    log                 yes;
    executeControl      timeStep;
    executeInterval     1;
    writeControl        writeTime;
}

SSSSK

@学流体的小明想问一下，你这个里面计算摩擦速度得方法3，我输出的剪切力有个最大值有个最小值，我是该用流动方向得最大值吗？

SSSSK

@coolhhh 想问一下你统计使用的哪个命令啊

coolhhh

@SSSSK 我用最笨的方法，sample沿着横向设置所有的纵向线，然后后处理求平均。of.com版本有自带函数columnAverage求展向平均，自带算例有使用示例，但不能对边界面的物理量（wallShearStress）平均。

学流体的小明

@SSSSK 应该用空间平均后的数据，最后就是一个只随y坐标变化的函数，比如u(y)。建议把这个帖子好好翻一翻，好多东西应该都讨论过。

学流体的小明在 LES直流槽道边界层模拟，如何得到正则化速度u+以及正则化坐标y+？中说：

是底面平均之后的，postChannel工具会把槽道按照xz线（x和z坐标固定，y坐标变化）求平均

SSSSK

@学流体的小明 postchannel命令里面出来的u就是x方向的剪切力吗，之前一直不确定，所以都用的是postprocess命令直接看的wallshearstress

SSSSK

@coolhhh 那剪切力还是得自己手动求平均么23333

SSSSK

@学流体的小明好吧我看错了那个u是rms值，所以如果用这个的话得自己修改程序才能求剪切力得平均啊？

学流体的小明

@SSSSK Uf是平均速度啊，都知道Uf(y)了，直接算剪切力 $ \mu\partial u/\partial y \big |_{y=0} $ 啊

SSSSK

@学流体的小明好吧哈哈，直接看剪切力习惯了

coolhhh

@SSSSK 对，没找到有计算壁面展向平均的函数

SSSSK

@学流体的小明大佬，你这个最后是多少万网格对上的？我算了好多组也都是你这种情况，而且如果换成ras得话结果更差，但网格多了又算得太慢了感觉

学流体的小明

@SSSSK
我$ \Delta x^+=\Delta z^+=10 $，壁面法向$ \Delta y^+=2 - 10 $

你可以看看这篇文章
https://www.sciencedirect.com/science/article/pii/S0021999117304059?via%3Dihub

李东岳

@SSSSK

而且如果换成ras得话结果更差

结果很差么？我看没啥问题啊，还算的非常快。你有RANS的结果么？

SSSSK

@李东岳有一个，差的特别大，感觉是我设置的问题，但设置我是参照李老师发的那个设的？感觉网格问题不大，图里面两个方法都是一套网格，280w左右，LES就对的很好，RAS差好多。对的是摩擦雷诺数180的，实际算出来的LES有179，在只改湍流模型不改动别的参数的情况下，RAS的有260多，而且没有下降的趋势 LESRASDNS对比.jpg

李东岳

这种流动算RANS没必要整280万啊，28个网格就行了。你这个280w不好debug，跑都得跑1小时嗨。你看看研究研究吧。应该是湍流模型没用好。

SSSSK

@李东岳谢谢李老师，小网格我也试过，我在调试调试，肯定是哪没设对

2019201300

空槽道的网格，U+Y+对应结果放下边了。非常感谢李老师和@学流体的小明的帮助，再就是多看看前边提到的那篇文献，很重要的，认真看完就都没问题了。
网格文件附下，需要自取链接：https://pan.baidu.com/s/1DfaqXnRQMle7JvyHv0WI1A?pwd=kexh
提取码：kexh

李东岳

@2019201300 厉害厉害，LES么？

CFD中文网

LES直流槽道边界层模拟，如何得到正则化速度u+以及正则化坐标y+？