CFD中文网

@学流体的小明 线比较多，分成了两张图。第一个数字代表起始时刻，第二个代表结束时刻。20s开始采集的在25s时仅采集了5s就远大于150s开始采集的值，并且采集到80和120的值基本不再变化了。而从150s开始采集的，会从一个特别小的值开始，也是逐渐增大，但即使采集了100s也不如另一个采集到25s时的值

@学流体的小明 20s以后Umean统计的Umean和100s统计的都是一样的，只有urms不一样而且差很大

@学流体的小明 uMEAN区别是不大，主要的区别再Uprime2mean

@学流体的小明 我感觉是不是从100s开始统计的时候，其实是把前面的时间步也考虑了，在平均下来就会很小。从20s统计的时候，前面没统计的时间步数比100s开始的要少得多，所以urms的值能更快的趋于定值。另外统计的时候我不是重新算的，是直接设置的starttime 20；应该没有你说的那种波动吧。

用openfoam计算完全发展channel时，以前都是跑到100s左右开始统计，发现urms会从一个非常小的值开始增加，最后统计大概200s后区域一个定值变化很小。最近发现20s左右就已经发展完全了，就从20s开始统计，这个时候发现即使只统计5s，得到的urms也远大于从100开始统计到200的值。这个是为什么？感觉很奇怪

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

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

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

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

@coolhhh 大佬你用of算过forced的这种湍流没？

@李东岳 李老师，跑了以下您的算例，有一点很奇怪，为什么网格加密之后反而对不上了呢，您的算例里网格是50个，这个时候我用v2对的挺好的，但变成500之后反而结果很差

@2019201300 感谢感谢

@2019201300 超过上限访问不了了

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

@李东岳 有一个，差的特别大，感觉是我设置的问题，但设置我是参照李老师发的那个设的 ？感觉网格问题不大，图里面两个方法都是一套网格，280w左右，LES就对的很好，RAS差好多。对的是摩擦雷诺数180的，实际算出来的LES有179，在只改湍流模型不改动别的参数的情况下，RAS的有260多，而且没有下降的趋势

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

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

@学流体的小明 好吧我看错了那个u是rms值，所以如果用这个的话得自己修改程序才能求剪切力得平均啊？

@coolhhh 那剪切力还是得自己手动求平均么23333

@学流体的小明 postchannel命令里面出来的u就是x方向的剪切力吗，之前一直不确定，所以都用的是postprocess命令直接看的wallshearstress