更新一下
发现了问题的源头，是网格的问题。
我使用ICEM画网格，或者OpenFOAM的blockMesh工具画网格，都是先画背景网格，然后对近壁面进行多次加密。加密是八叉树类型的加密，所以在网格尺寸变化的地方，大网格的尺寸直接是小网格的二倍。这样的网格确实不好，可能一些拓扑结构上的不好导致了singleGraph程序进行插值时找不到插值的基点。
解决办法：
举例说距离壁面0~1的网格尺寸是1,1~2的网格尺寸是2，那么就使用singleGraph在0~1这个范围内提取一条线，在1~2再提取一条线，这时的type可以是uniform，也可以是midPoint或者midPointAndFace。一句话，就是分段提取信息。
我是重新画了网格， 老师说之前画的这网格根本不能用:mihu:

@希望先生 硫酸老师说他上传的版本坏掉了，缺东西，找时间再上传一个新的:135:

@adcpk 超算每个节点20核，是因为只有20核还是只用了20核？如果是只用了20核，其他核在跑别的作业吗？如果只有20核，建议换成15+15这样来测试。

一般情况下超算平台的CPU性能肯定没有你工作站的好，工作站肯定比超算新吧，快一些的话其实也可以理解。

这种边界条件不是特别稳定。很容易发散。所以总是报错。

@李东岳 感谢李老师！现在对这个求解器理解了一些，对于$\mathbf{U}_1$的控制方程, $K_d \mathbf{U}_2$一项作为显式项被直接纳入pEqn.H的界面通量的修正。

非常感谢！已更新

嘿，我还没关注这个历史。还真是。太折磨人了。

@韬智tz 在 OpenFOAM模拟空化气泡一些参数的问题 中说：

那最大气泡半径怎么得知呢

这个应该是随着模拟自动发展出来的吧？

你可以直接采用湍流模型来算试一下。用湍流模型计算层流，只要湍流粘度足够小，影响就很小。

@李东岳

sudo apt-get install libcgal-dev

这样就可以编译了

openFoam可以无量纲的，应该在etc/controlDict这个文件里里吧dimensionset 1改为0

@nanxuan 请问你解决这个问题了吗？能否分享一下？谢谢

@星星星星晴 谢谢星星前辈，这就去试试:146:

各位前辈好，我采用open foam并行计算，想用tecplot后处理。导入controlDict后只有模型显示，没有计算时间信息。另外在open foam中重组的时候出现以下报错Create time

--> FOAM Warning :
From function int main(int, char**)
in file reconstructPar.C at line 220
No times selected
不知道什么原因。请问有没有什么办法可以让我用tecplot后处理的。

@winsway_zero 千万级别笔记本开paraview比较难，还是得找台本地的服务器。

数据传回本地你可以写个脚本啥的自动去跑，不过这个得看集群的网络。我们是大量千万网格级别的计算随时回传的，感觉速度也还行。不过我们的计算是稳态不可压，算是最简单的了。

有效的后处理看你需求了，如果同类型的案例，只要一开始能够看到三维的展示就好，这样能明确知道哪些地方需要关注。后续写成各种function配到controldict里面，计算完成之后就直接跑一遍，这样生成的结果文件数据量小一些。

@李东岳 感谢东岳老师提供的思路，我在代码中把ShihQuadraticKE改了个名，目前在源文件中用Allwmake编译成功了。
不过孩子还是很不明白怎么单独拿出来编译这类非线性湍流模型，看来暂时只能通过这种方式代替了。

@bestucan 好的，谢谢老师

@wangfei9088 好的，谢谢，可以了

@chen_hao
解决了，改成了无量纲，就避免出问题：

volSymmTensorField Sij(twoSymm(tgradU())/omega_); volTensorField Wij(2.0*skew(tgradU())/omega_);

@cfdngu 收费也没办法，我之前找了一圈也没有找到比这个好用的