CFD中文网

请问有人能帮忙解答这个问题吗？

圆管DNS后处理中需要分析湍动能TKE和Reynolds shear stress（RSS），表达式分别如下：


其中Ui是圆柱坐标下三个方向上的速度。ui''表示Favre平均下的速度脉动。Favre平均就是：

那么我想请教下我这么推导是正确的吗？

或者说通过其他量就能计算得到TKE和RSS? 我试过Uprime2Meanxx来表示u''u''，画出来的趋势差不多。但不知道对不对。

另外我想了解圆柱坐标下的径向速度Ur如何得到呢，通过根号(UxUx+UyUy)吗？

以上问题还望解答，谢谢！

@cxzhmg 我之前跟同学讨论的是泄压可能是因为设置cyclic的效果就像是一道虚设的面，面上可以流进流出，而真正有阻挡作用的是wall壁面，因此channel就像是上下两个平板之前存在一个自由流动的空间一样，于是留着留着就泄压。但如果像圆管那样四周都是wall，就不会泄压。另外我采用了低马赫数假设，压力的变动不影响物性计算。因此无论压力是否泄压，我的物性计算还是正常的。

圆管DNS的后处理，有的论文里面图长这样：

纵坐标的u''是基于Favre平均下的脉动量。Favre平均就是：

所以我在想OpenFOAM在计算时会自动生成一个张量场UPrime2Mean，如果控制方程是可压缩的话，我看网上说可压的雷诺平均NS方程相当于Favre平均NS方程。

那么我不确定这个UPrime2Mean是基于雷诺平均下的量，还是基于Favre平均下的量？

如果是基于雷诺平均下的量，那么为了得到Favre平均下的速度脉动，是否需要在计算前先生成复合量ρu，然后计算出它的时均量，接着用ρu的时均量除以ρ的时均量，得到u的Favre平均量，然后再用u的瞬时量减去u的Favre平均量，得到Favre平均下的瞬时的脉动量。

但有个问题是作图需要时均的脉动量，而我只知道如何求得瞬时的脉动量 时均的脉动量要怎么办呢？ 自己手动将一个个瞬时量做时间平均，效率太低了。还望指点。

没有具体看流线，但我觉得应该是弯曲的，不过不清楚是否是规律弯曲的。

先前有研究竖直圆管的摩擦压降分析。我画出了圆管的轴向中心沿线分布，查看了压力场，发现整体趋势是笔直的，但局部却是成台阶状的分布

按理说不会出现这种情况吧。压力突然垂直向下，那不是压力梯度无穷大了吗？但是计算过程仍感觉没有问题，没报错。
不知道是什么引起的这种现象？

我想了解下OpenFOAM是否可以批量截取截面？比如轴向有1000个网格点，就截取1000个面，然后计算每个面上的某物理量平均值。我看cfdpost可以联合Python进行批处理。我尝试了foamDataToFluent，但导入fluent流场都错了，就做不了批处理。

我做的压力驱动流，也是两个平板间的channel flow，上下有温差，但不考虑浮力。瞬时温度场长这样：

@Exthan 一般有两种方法，查表法和多项式拟合

@东岳 谢谢老师。uTau的公式我知道，但不用把这个公式植入代码里面吗？如果需要的话，有点不太清楚怎么写入代码，还望赐教。
DNS圆管的算例目前还在算呢，所需要的资源太大了，算比较慢。

@东岳 老师 ，请问uTau不应该用公式植入吗？在沿管方向的壁面上，uTau应该是变化的吧。如果用初始值计算出来的uTau定值植入，我试了下，湍流模型没有达到自己预想的效果。

@东岳 好的谢谢！

如果要植入MK湍流模型的话，里面有一项
其中y+我按照了东岳老师提供的方法，那就是如下图：

编译成功，但按这个湍流模型跑的结果不对。

uTau为何设成1.003呢？如果要植入uTau的公式，那tauw这个如何得到呢？源代码里面相关的有tauw跟wallshearstress，但是不知道怎么引用。还望赐教。

@星星星星晴 然而也没其他办法了吧

目前速度场文件能用，温度场我批量地把低于301.15K的值统一改为301.15K（初始温度），然后温度场看起来正常了。压力场我是用之前跑过一小段的压力文件。四拼五凑跑起来了 不知道压力场用比较之前的文件，会不会对趋于充分发展的湍流造成影响，又给打回原形去了。

我目前用mapFields 以及 mapFieldsPar 两个工具，将750W的网格流场映射到7500W的网格上。
我先是用mapFieldsPar 映射，发现速度流场最大值从0.76变成1.1，然后温度场更是出现了150K这样的低温，本来都是301.15K以上的。

接着我用mapFields，默认的映射方法会报错。然后我改用-mapthod mapNearest，可以正常映射，但结果仍然跟之前的一样，流场变了。

然后自己改用了-mapthod cellPointInterpolation，结果速度场变得比较正常，最大值0.79，比0.76稍微大了一点，其他细节都跟源算例比较相像。但是温度场，压力场仍然出现很低的低值。

但是我点击了paraview的rescale to visible data range，流场就变得正常了，应该是忽略了一些细微的不怎么可见的温度/压力范围。映射前后对比如下：

接着我放大了细节，发现低温出现在了入口的边缘（图中蓝色），很细节的地方，映射不成功。

对于这种情况，可以忽略这种细节吗？或者还有其他的方法映射吗？试了好久，求助。谢谢！

@东岳 非常感谢耐心的解答和分析！颇有启发！
我有大致对比过advective跟inletoutlet的效果，如上面的图所示，不理解的是advective在出口处为何有个速度突增。但另一个边界条件的结果显得更为离谱，我简单plot了其壁面的速度分布，是疯狂地震荡。根据您以上的分析，advective似乎更合理。我如果有发表类似论文，乐意提及CFD-China对论文工作所做的贡献。

@东岳 好的谢谢！

前一段时间一直在算空气在加热圆管流动的DNS，目前算了两秒多，还在算，残差保持稳定，但数量级有点高，p跟U稍低于0.01.平均壁温还没有达到光滑分布。我对于速度出口设成 advictive 不太有把握，这个边界的原定目的是想消除反射波。另外对于DNS收敛的监测角度，有什么比较有效的判据吗？谢谢。

最近还在国家超算上用DNS算超临界CO2的圆管对流换热，热物性变化十分剧烈。后面看看情况。

@东岳 我看的DNS文章里面，都是用这个convective boundary作为速度出口，我的控制方程跟他们的一样。所以我就借鉴了。
另外我对比了advective跟inletoutlet两种边界对结果有何影响，对比如图示：


两张图的上面部分是inletlet，出现了很奇怪的现象。下面部分是advective，不太理解的是第一张图的下半部分，右端出口处有速度的突增，不知道怎么解释这种情况？谢谢。