@李东岳 东岳兄，这个是一个基于OpenFOAM和LAMMPS的求解器，有一篇文章叫：Sun_2016_SediFoam A general-purpose, open-source CFD-DEM solver for particle-laden，可以一看。这个求解器其实功能是和CFDEM差不多，也就是求解密相流动的。其不属于CFDEM一部分。

Hi，

这和你的求解方程有关，例如：

fvm::ddt(rho) + fvc::div(U)
和
fvm::laplacian(nu, p) = fvc::div(phi)

是对称矩阵。更进一步的，第一个是对角阵。

fvm::ddt(U) + fvm::div(U)

则是非对称矩阵。

of用的是自己的容器

@xiaohsu2006 是根据温度和压力算出来的，　constant 参数后面的值是根据什么给的呢？

“\”好像是Linux的语法。 我的理解，是换行的意思，表示后面还有代码需要执行。你看最后一行就没有""了

musgs= rho*nusgs alphasgs = musgs/prt

对于mvCoeff 和　mvCoeff 两个变量，当他是可压缩的时候，成为密度的函数，该怎样修改呢？？

很感谢！！！

omegaWallFunction 这个边界条件继承自 fixedValueFvPatchField类，其构造函数如下
0_1460115657665_upload-47e7e05b-ec74-4419-8ff9-d55730af0dfa

这个构造函数里，需要注意两点，第一是传给基类fixedValueFvPatchField的参数，也就是 fixedValueFvPatchField<scalar>(p, iF, dict)，这一句里，基类会通过参数iF 去读取 类似 value uniform 2; 这样的语句，来讲当前边界的值初始化为 2。不过，看下面最后一句，this->operator==(patchInterField());，这里将当前边界的值重新赋值为临近边界的第一层网格中心的值。

所以，综合起来，value uniform 2; 这一句不起实际作用。

@搬运工不好当

这种情况很常见，如果使用第三方网格生成程序比如ICEM做的网格需要使用periodic严格定义才可以。如果你确定你的面试匹配的话，可以加大tolerance来忽略这个提醒，也不会有错误。这就是网格的原因。

这里的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的作用就是消除非正交网格引入的误差。