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

圆管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圆管的算例目前还在算呢，所需要的资源太大了，算比较慢。