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

[李东岳]

@falltime 是说$R{e}_{\tau }=u_{\tau }\delta /\nu$，这个公式只用从$u_{\tau }$算$R{e}_{\tau }$，而不能反过来用是么？

[falltime]

@李东岳 在 LES直流槽道边界层模拟，如何得到正则化速度u+以及正则化坐标y+？ 中说：

@falltime 是说$R{e}_{\tau }=u_{\tau }\delta /\nu$，这个公式只用从$u_{\tau }$算$R{e}_{\tau }$，而不能反过来用是么？

是的，算之前是不知道的。除非在给定压力梯度的情况下，通过压力梯度和壁面切应力的关系来算

[李东岳]

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

[coolhhh]

前期计算的${⟨u_{\tau }⟩}^{+}$会出现参差的现象，原来是之前在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 \mathrm{\partial }u/\mathrm{\partial }y{|}_{y=0}$ 啊

[SSSSK]

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

[coolhhh]

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

[SSSSK]

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

[学流体的小明]

@SSSSK
我$\mathrm{\Delta }{x}^{+}=\mathrm{\Delta }{z}^{+}=10$，壁面法向$\mathrm{\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]

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