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

[SSSSK]

@coolhhh 最近刚接触这种湍流，有些地方有点疑惑，decay的话是不给力，最后让他一直衰减到0吗。然后forced就是给他一个力维持，最后让他平衡，给的力不同的话对结果影响会不会很大，如果我算forced的话，直接用dnsfoam里面的那种force可以吗

[coolhhh]

@SSSSK

1. 我算的是均匀衰减湍流，没有额外的force，最后湍动能衰减到0。以下是书上关于均匀湍流的阐述：

引自《张兆顺, 崔桂香, 许春晓. 湍流理论与模拟[M]. 北京: 清华大学出版社, 2005》31页

2. dnsfoam我没用过，没了解过是怎么加的，是否类似of自带算例channel395，施加了一个梯度源项？我理解湍流能够自保持，NS方程右边是合力，如果没有施加外力，耗散项不等于0，使得湍流不断衰减。channel395算例加了梯度源项，这个力刚好抵消耗散项，等到充分发展状态时候，NS方程右边合力等于0，这个时候就能保持不变发展。根据充分发展时候合力等于0，施加的力大小不同，那么耗散项的大小也不同，耗散项与速度梯度有关，也就对应了那几个平均速度剖面，雷诺应力剖面的不同。

[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上。