@搬运工不好当

CFD界：相对于商业代码，OpenFOAM的特性是什么？

Henry：OpenFOAM是完全免费并且开源的（CFD界注：开源的意思即为使用者可以随意更改软件的内核）！这提供了HUGE的方便，不仅仅在费用上，也表现在灵活性、以及从软件本身来学习CFD上。OpenFOAM不仅仅让CFD工程师可以模拟算例，它也允许CFD工程师来学习CFD方法、开发CFD算法并完善，甚至某一天创造一个可以替代OpenFOAM的软件。
并且，许多CFD问题需要一些特别的功能、算法、方程、边界条件等。在商业软件中植入这些功能可能令人沮丧并且很可能失败。有些情况下可能有必要完全的重新写一个求解器、边界条件等。OpenFOAM提供了最大的灵活以及扩展性。这些对使用CFD做研究的人非常重要，因此我觉得OpenFOAM可能在某一天会替代闭源软件。

@wwzhao 在OpenFOAM-2.1的时候，程序里写的就是“CrankNicholson” ，所以controlDict里也应该写“CrankNicholson”，虽然这个拼写确实是不正确的。但是从OpenFOAM-2.3开始，这个名字改成了正确的拼写：“CrankNicolson”，相差一个字母。所以这就导致了2.1和2.3的算例会在这里不兼容。

@搬运工不好当

在对拉普拉斯项离散的时候，如果使用非正交网格会引入误差，因此需要非正交修正。所以nNonOrthCorr的作用就是消除非正交网格引入的误差。

你的arc定义错了，你觉得一个点能定义一条arc吗？

你得把错误内容贴出来，不然一句话看不出来。出错的原因太多了，而且不知道你修改了求解器，甚至case的哪部分...............

transportpProperties没有声明。

这几天放假没在学校，才看到回复。谢谢分享 。我研究一下。万分感谢

是的，可以忽略。

并没有用过pyfoam

@xpqiu 谢谢！

在liggghts/src下面找到这个fix_heat_gran_conduction.cpp文件，里面讲了颗粒和壁面之间的静态导热的计算，看看代码应该可以知道的。

用surfaceCheck检查STL文件的结果呢？

diagonal是对角矩阵求解器。

fvScalarMatrix rhoEqn ( fvm::ddt(rho) + fvc::div(phi) == fvOptions(rho) );

从这个方程里面可以看出，只有时间项是隐式离散，fvc和fvOptions都进入了矩阵右边。因此离散之后是一个对角阵。直接除过去就可以。并不需要迭代。因此迭代数为0.

这个原项是怎么弄得？fvOptions，之前我记得都是MRFSimpleFoam，SRFSimpleFoam，现在都整合到fvOptions里面了。原理是什么啊？

@mark 这种“提示文件找不到”的都是路径搞错了，我看你已经编译出来了，仔细查看查看路径。

把fvsolutuon里的tolerance由e-6改为e-8对结果影响大不。

像大家请教一个问题，算低速流体（10m/s以下）且求解域很小（1cm1cm1cm量级）当用可压求解器求解好像CFL不管设多小，算足够长时间都会发散，除开用不可压求解器以外还有什么其他解决方式么？

@李东岳 想明白了，多谢东岳大神！

是中国人，还没改国籍。

互联网的某个角落 or 提问。