看了这个，有点懂了，问题应该在于离散格式和梯度限制器啥的。
Finite Volume Discretisation in OpenFOAM Best Practice Guidelines

Hi,

并没有用过这个，不过cfd-online有关waves2Foam的讨论已经多达2000贴。从常识来看，附加湍流波高减小是可以理解的。

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的这个帖子。

@搬运工不好当

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 提问。

@搬运工不好当 每个边界都要赋值，格式大概是
nonuniform List<scalar>
xxx
(
aaa
bbb
ccc
.....
);
xxx 要等于边界面元的数目，括号里的数字的个数要等于xxx。

另一个用法是vtk可以存储部分数据，而不是整个的场。大型算丽非常有效。

另外有些paravirw的后处理需要vtk做源source.

@dyj19901127 hi.

我觉得一个好的习惯就是禁止复制构造函数。除非你每个类都写，但一般也用不到。henry应该是采取不写的做法吧。

这个在c++ primer plus 和effictive c加加里面提到过。书并不在手边你可以看看那个类里面存在指针的时候，复制构造函数怎么泄露的。我更觉得这是个和编程有关的问题和cfd关系不大？

@搬运工不好当 你先wclean一下试试:）

@搬运工不好当 shm生成的网格有的时候非常的tricky，不错如果你采用shm成功的生成了mesh okay的网格，openfoam大部分求解器都能算收敛。of里面单相求解器已经非常稳健了。

如果对于现有几何，一定要生成一个非常好的网格，我建议你

简化几何 换专业的网格生成程序 慢慢的，充分的调节shm里面的网格控制参数，主要是网格密度，迭代次数，甚至删除边界层，这个步骤需要很长时间。

话说回来，shm对于复杂几何生成的是混合网格。虽然单相求解器对混合网格非常稳健，但某些多相求解器还算不了。

@hongfu2233 不，湍流粘度并不是一个直接求解的量。并没有湍流粘度的方程，因此不需要初始值。

@liliang2009223 我都那个丸不是很了解。这俩个互溶么？液体

@mark dgdt分析已经上传