This questions is not clear. What is 液力变矩器?

@搬运工不好当 说：

求速度方程之后为什么没有求压力方程，而

Hi，

我觉得能量方程放在压力方程之后和之前结果是相同的。就像湍流求解，标量传输放在压力方程之前和压力方程之后也是相同的。

你可以自己测试一下并反馈？

@李东岳 谢谢东岳大神。祝论坛越来越好。

@李东岳 东岳兄，这个是一个基于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的这个帖子。