这里的errors是误差的意思，而不是错误。可参考cfd-online上的这个帖子。

谢谢回答

fvm::xxx(u)，其中的u是需要求出来的，fvc::xxx(u)，其中的u为当前时间步的值，其返回一个场。所以，需要你要求某个场u，就用fvm。

@李东岳 MIFX采用cut-Cell的方法处理几何模型，一般都是用来计算几何结构比较简单的几何体，但是网站上也挂出来了比如旋风分离器等的例子。其处理方法和基于MP-PIC方法的Barracuda软件一致，原则上复杂网格都可以处理的。

Liu et al., 2016 "A Comprehensive Benchmark Suite for Simulation of Particle Laden Flows Using the Discrete Element Method with Performance Profiles from the Multiphase Flow with Interface eXchanges (MFiX) Code" 一文中对MFIX-DEM的并行能力进行了验证，采用了Weak Scaling和Strong Scaling两种评估方法，结果证明在核数达到1000左右，MFIX-DEM code的并行能力还很不错。但是注意到，在做自由堆积等的时候，文章回避了负载均衡的问题，也就是如果流化床下半部分非常密集，而上半部分非常稀疏，这是对就需要一种新的分区方式，使得每个核都能计算差不多的颗粒数。现在的MFIX好像还达不到这个功能，但我知道有人开发出来了。文章还认为，当颗粒数达到10(8)，即使用1000核，每个核也有10(5)个颗粒，还是非常有挑战性的。气固两相流的计算时间主要花费在了DEM颗粒彭碰撞计算，颗粒和流体之间的数据交换，曳力的求解。
看到文章的一点心得，和大家分享一下。

不好意思啊，我不是做颗粒流的，所以这个问题无法帮你。

我建议你阅读相应的源代码，可以用find命令查找具体模型的源码：

find $FOAM_SRC -iname 'coneInjuection'

@搬运工不好当 说：

OpenFOAM2.3.0中有buoyantPressure这个边界条件吗?没有的话要用哪个边界条件来代替 buoyantPressure呢？

OpenFOAM 2.3.0 中没有buoyantPressure，替代的边界条件为fixedFluxPressure，详细可参考cfd-online的这个帖子。

@搬运工不好当

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.