@D-Benjamin
谢谢你的回复！
一开始我表述的问题，我是想取网格节点上的值，不是网格中心的值。
这个表是东岳前辈翻译的openfoam3.01版本的用户指南。

找到一个这个，用face应该是网格点上的值？

@Mania 还有另外一种方法是安装两个版本的openfoam，然后后处理的时候用低版本的操作

@Mania 抱歉昨天忙了一天没有看到消息，没回复。
他的版本是比较低的版本，大概是2.3.0左右的，所以高版本下会出问题，你可以找下那个连接下有没有高版本的
或者另外一个帖子下有其他的链接，你可以试一下

@Mania 我的意思是直接在第一次计算的时候就用大括号，不知道会不会有影响。
有钝体的话，matlab作图确实会出现自动插值的图，可能会有问题。当然也是可以解决的，就是在有钝体的地方你要放上网格信息，并把值赋值，做成白色就可以。
你是想做展向平均吗，做展向平均的话，cfdonline上有别人写的程序

@Mania 不知道版本有没有影响，在2017版本里面公式内的变量需要用大括号{}，你可以看下tecplot help，下面这张图是17版本的help的提示

所以你换这个写的话，应该是{V34}=...，你可以试一下。
另外，其实你这个做三个面的平均的话，完全可以用matlab，非常简单，几行代码的事

一年多了还是没有解决这个算例，不知道有没有人做过相类似的，重新描述下这个算例，如图所示一个后台阶流动，但是是在图中左下（0,0,0）位置为原点绕Z轴旋转，主要流动参数，
雷诺数$Re_b=\frac{U_b(H-h)}{\nu}=5600$，其中$U_b$是体积平均速度；
旋转数$U_b=\frac{\Omega(H-h)}{U_b}$，$\Omega$是角速度；
台阶前段槽道流部分壁面剪切雷诺数$Re_\tau=\frac{u_\tau(H-h)}{\nu}=360$，这个与槽道流DNS结果是完全一致的，也就是说台阶前是完全发展的槽道流。目前存在以下问题：
1、求解器选择，SRF类还是需要自己修改求解器？
2、合理的进口条件，不旋转的相对好给，但是旋转的case，速度进口边界条件和初始值不知道如何给合理。
3、无量纲化的流向速度为什么会小于1，如下图所示。以我的经验，槽道流里面最大速度比体积平均速度肯定是大的，也就是$\frac{U_{max}}{U_b}$是大于1的，但是为什么在旋转流动里面就会于1呢？

参考文献
1、DNS原始文献：Turbulent flow over a backward-facing step. Part 1. Effects of anti-cyclonic system rotation
2、图片来自：Development and Application of Rotation and Curvature Correction，是篇博士论文

@liujunCFD 湍流层流粘度比这种说法来自哪，没有听过。湍流粘度比是一个估计值，指的是湍流粘度比流体粘度 $\mu_t/\mu $，一般取10就可以，这里有个解释可以看下

@liujunCFD CFDonline计算工具，可以自己对照一下
https://www.cfd-online.com/Tools/turbulence.php

@chszkc 在我的理解里是不行的，设定了cyclic边界后就是周期性边界条件了，在ployMesh里面的boundary里面，inlet/outlet就必须设定为cyclic，而zerogradient适用于wall和patch，不能设定为cyclic，两者矛盾，应该计算不了。

@Stan同学 是的，这个跟paraview是不一样的，paraview有internalMesh这个选项，tecplot是没有的
你做振动的话，使用FSI吗？

@Stan同学 不太清楚你算的什么算例，tecplot显示出来的应该是只有边界，要处理内部场的话可以尝试做个slice，或者处理一个Q准则等内部流场的涡结构，显示一下看看有没有结果