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

[SSSSK]

@coolhhh 感谢分享。dnsfoam是通过加了一个力得到的，这个力好像是按照wiener process加的，然后傅里叶变换到xy坐标。但这个加力的方式有个bug，他没法并行，所以像很多论文里都按照波数来加的那种方法的话就没法实现，也查了好多资料，就是查不到openfoam算forced各项同性湍流的。

[李东岳]

我之前一直想做dnsfoam解析，一直没倒出时间。dnsfoam这个跟icofoam啥的唯一的区别就是这个力，这个力的生成算法我还没细看，dnsfoam解析也一直没写出来。我记得是有什么傅里叶变化啥的东西在里面。

另外新版dnsfoam应该已经没了。变成了一个力源项。

[coolhhh]

@SSSSK 确实没看过对均匀各向同性湍流加force的，一般都是多是对类似channel395这种非均匀湍流加force。就是对均匀各向同性湍流加force，要想想是否理论上加force是否就能实现不衰减？或者加的force需要一个很复杂的形式才能实现，而不像channel395只加一个固定的压力梯度。

[coolhhh]

对均匀各向同性湍流加force的自保持，如果只是要实现到目标结构位置湍流特性是所想要的，可以对入口湍流乘以放大系数。均匀各向同性湍流衰减曲线的对数率形式接近一条直线，根据衰减曲线估计下在入口处对湍流增强洗漱，应该能调整出目标位置的所需要的湍流特性

[SSSSK]

@李东岳 新版的能像dnsfoam一样加白噪声力吗？原版的我稍微看了一下，是通过高斯分布得到白噪声然后再根据维纳过程施加和时间有关的力，最后傅里叶变换变成方程里的force，李老师你说的新版是of11和12吗？

[SSSSK]

@coolhhh 对这种湍流加的力和channel不一样，channel只需要一个固定压力梯度就行，但这个是从波数域加的，需要傅里叶变化，但傅里叶变换里面有个kmesh，就是这个不能并行

[LESKIN]

@2019201300文件过期呜呜呜，最近也在研究管道流动问题，想知道大佬怎么做的，可以重新传一下文件嘛

[2019201300]

@LESKIN 这个是个槽道，流向和展向都是周期性边界条件。通过网盘分享的文件：Channel180_Kong.rar
链接: https://pan.baidu.com/s/1TvKiFbKk-H0n-pm5LBmXyA?pwd=vnfi 提取码: vnfi

[LESKIN]

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

@LESKIN 这个是个槽道，流向和展向都是周期性边界条件。通过网盘分享的文件：Channel180_Kong.rar
链接: https://pan.baidu.com/s/1TvKiFbKk-H0n-pm5LBmXyA?pwd=vnfi 提取码: vnfi

太感谢了

[东方白杨]

您好，我也是用的上面说的论文里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 刚刚仔细找了一下，发现大佬后面又发了一下链接，能正常下载了，目前就第二个问题需要大佬解答