谢谢大神，我看到你原来的代码里有相关的。再次感谢。
http://www.cfd-china.com/topic/1228/owner-neighbour-网格寻址/2

在外挂硬盘上计算非常慢，不写log可以，写log的话，在机械硬盘上写log严重耽误了计算的时间。

另外，log文件有时候变的特别大，我这个算例时间步长比较小，写文件间隔比较大，导致结果没有多大，但是log上G。可以通过 > /dev/null来禁止写log，当然如果你确定不需要看log的话。

0_1525840011230_2018-05-09 12-25-51屏幕截图.png

约2G左右

https://www.sciencedirect.com/science/article/pii/S0010465517303375

See this dsmcFoam+ published on Computer Physics Communications.

5.5.1 fluentMeshToFoam
......
There is currently no support for embedded interfaces and refinement trees.
貌似对这个是不支持的。我也用过snappyHexMesh， 这种refinement的过渡处是会生成非结构网格的，也就是没有interface。现在snappyhexmesh是做了一些改进吗？

我的理解是fvc::div(phiHbyA)已经将原控制方程中的散度表示为高斯积分的求和形式了

是的。正确。

@yuan_neu 4如果从数学意义上来说没有错误。考虑到物理机理，如果模拟的是湍流，那么所谓的2D是错误的，因为湍流本身是三维的。

momentumPredictor为动量预测步骤，对应链接文本中的方程14。动量预测步骤可有可无。从经验来看，在粘度比较大的时候，方程比较刚性，若附加可能会引起发散。在多相流以及动网格算法中，通常此项关闭。具体原因，目前并没有见过相关文献讨论动量预测对收敛的严格证明。

目前尚不清楚动量预测的普适性方案，后续还有待系统研究。下面的算例是一个表示开启动量预测较好的算例。但不具有普适性意义。算例模拟一个静止的流体。纯测试算法，理论上模拟结果速度应该为0。模拟中考虑了重力。采用p方程计算。

首先，俩种结果的压力场（右图）均正确。均体现了由于重力导致的静水压力的渐进变化。但速度场都不是很正确。这是因为在附加重力场的时候，用的是p方程而不是p_rgh方程。类似讨论： http://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/37

然后，momentumPredictor开启的时候，假拟速度更好一些，如下图。

momentumPredictor off

0_1525236812110_pic2.png

momentumPredictor On

0_1525236812110_pic2.png

谢谢东岳老师考古！讲解非常清楚！

@东岳 谢谢您的回复。我的边界条件是进口速度fixvalue，进口压力zeroGradient。出口速度inletoutlet，出口压力0。因为我之前AMI几何形状不一致，后来才改的。是不是除以0就是因为AMI有一部分没有贴合上。

@东岳 谢谢东岳老师，这个据我所看文献现在of并没有植入，植入很麻烦。现在也不是我研究重点，等有时间在弄这个。

虚拟机肯定是没问题，跟Ubuntu真系统感觉没什么区别。
但是BSL下面的OpenFOAM能感觉出来有点慢，之前对比过，主要是来自于磁盘读写。如果你有兴趣的话，可以试试对比下，谢谢。

@浪迹天大 问题我已经解决了，就是说明一种现象

是的，cyclicAMI是CFD中的滑移网格，试试添加关键词lowWeightCorrection 0.2;。可以上传你的图片看看

参考interDyMFoam下面floatingObject算例，用topoSet生成中间的小方块

可以用输入边界法造波，比如IHFoam，waves2Foam等，现在的of也自带边界输入法造波了，只是功能还比较单一。

@东岳 在 DPMFoam求解器颗粒计算有颗粒跑出计算域外! 中说：

粒速度持续

假如设定反弹，颗粒轨迹与壁面有交叉，但还是跑出了计算域，正如之前碰到的情况，还是算法本身不够完善，openfoam5.0以经把颗粒追踪那一部分的算法完全更新了，颗粒追踪不会失效了

@东岳 对，后来查看了一下userguide 中的solver部分，里面PCG的全称是preconditioned conjugate gradient，就一下明白了

因为你肯定没回头去把手册或者视频看一遍:xinlei:

从哪算取决于starttime的设定和有没有那个时间节点的文件夹。

看是不是只能从0.8算，那就mapfields过来后，从0看看能不能算；算了之后0.8有没有更新。

至于看0用没用上，把0里的边界条件瞎胡改看看抱不报错不就知道了

高中学习和大学考试还是那种“我一恼”“我一生气”就做到了的模式。往后就是“临渊羡鱼不如退而结网”。重点是“退”，死磕的话。。。。不是有个成语叫“磕（客）死他乡”:wocao:

@dzw05 湍流模型中有一个运动黏度nu，正常情况下在后面赋值即可。我当时在网上查询如何提取升力系数和阻力系数时有篇帖子提示要在nu后面先再输入一个nu字符再赋值，由于之后的计算我都用的是这个错误的湍流模型文件，导致阻力系数求解一直出现问题。后来一个一个文件排查的时候发现是这里的问题。