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

[东方白杨]

您好，我也是用的上面说的论文里preturb生成初始湍流场方法，可是目前有点问题。一开始初始场的速度型线确实是抛物线，但是算50s就开始变形，算到1000秒速度型线就非常离谱了
下图是50s的

下图是1000s的

[东方白杨]

是我搞错了。。。湍流就应该这样，抛物线型速度型线是层流吧

[李东岳]

是的，你得湍流的没问题，就是时间平均多搞搞，目前不对称

[东方白杨]

@李东岳 嗯嗯，老师您知道Eugene de Villiers论文里的LES数据在哪里能下载吗？我想对比一下

[东方白杨]

一个星期，终于把这个算例搞明白了，纪念一下，感谢上面前辈们的分享！

[李东岳]

@东方白杨 我去，厉害啊，老铁，你这个厉害。

上面帖子当时我没回，后来我就忘了

[东方白杨]

@李东岳 方法都是一样的，只不过算的更多时间更久，疯狂造超算资源（滑稽滑稽）

[风云5091]

@学流体的小明 大佬，请问你们y+、u+沿某一条轴上的分布曲线是怎么画出来的啊，openfoam自带的后处理函数yPlus貌似只能得到壁面第一层网格上的y+，其余部分都设置的是0

[学流体的小明]

@风云5091 我们这些都是槽道流里面的实际坐标，除以特征长度，得到的无量纲化坐标y+。
你如果只做验证的话，OF输出的yPlus足够了，应该可以在paraview显示出来。
OF输出的yPlus确实只存在于那个边界patch上。

[风云5091]

@学流体的小明 谢谢大佬解答。我是做可压缩湍流（平板）边界层内的热化学非平衡效应的，最近要验证openfoam里的求解器，希望能算的准边界层内的流动，应该要用到整个流场里面的u+、y+，我看您上面发了一个网盘链接有直接算摩擦速度的matlab程序，但是链接显示过期了，麻烦能再发下吗，我想学习参考一下。另外大佬之前有用过rhocentralFoam吗，我用的hy2Foam非平衡程序是基于这个求解器开发的，不知道这个求解器算边界层怎么样？

[风云5091]

@风云5091 刚刚仔细找了一下，发现大佬后面又发了一下链接，能正常下载了，目前就第二个问题需要大佬解答