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

@李东岳

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里面，计算完成之后就直接跑一遍，这样生成的结果文件数据量小一些。

@rachel0096 我没看过他这个 但基本都大同小异 其实 写一个ibm在openfoam就用不了多少代码 我也就看了几篇文章 找准一篇你觉得好实现的对着做 你可以先用最简单的方法 比如在pimplefoam里 对于2d圆柱绕流 写一个函数检测小于半径r范围内的cell 这样就检测出一个圆柱范围内的流体cell 标记为1 其他流体的cell标记为0 然后在压力方程那里 根据固体cell和desired velocity的速度差算出体积力 加给phi那里 进入压力迭代 这就是一个最简单的ibm 这个弄对了再加上二阶插值分布函数那些来考虑更高精度的界面识别

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

@bestucan 好的，谢谢老师

@wangfei9088 好的，谢谢，可以了

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

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

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

@tens 感谢！这正是我想要的

如果已知流场内某一点的坐标（a，b，c），想通过插值获得这一点的流场物理量，可以使用interpolate（）函数，这一函数的定义如下： template<class Type> Type Foam::AveragingMethods::Basic<Type>::interpolate ( const barycentric& coordinates, const tetIndices& tetIs ) const { return data_[tetIs.cell()]; }

代码中coordinates、tetIs分别表示什么呢？
如果已知点（a，b，c）如何获取coordinates、tetIs呢？
谢谢！

@李东岳 十分感谢

@李东岳 正在学习，谢谢李老师！

@李东岳 好的李老师，里边的数据已经拷出来了。多谢。

@cccrrryyy 感谢感谢！给我指了一个好的修改方向！我做过RANS模型的对比，确实影响不大，但我看文献里有的是Smagorinsky算的很好，比如文献1，有的就是算的误差比较大，比如文献2，一头雾水:136: 还是得自己测试呀！

@sibo 我靠，前辈，我也是做结冰的