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

[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多，而且没有下降的趋势

[李东岳]

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

[SSSSK]

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

[2019201300]

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

[李东岳]

@2019201300 厉害厉害，LES么？

[2019201300]

@李东岳 对的LES，用的WALE模型，搞出来有一段时间了，只不过最近在忙着开题，开完了看到还有人在问，然后就传上去吧，供大家参考。

[SSSSK]

@2019201300 超过上限访问不了了