@浪迹天大 非常感谢您的帮助，我已经完美解决这个问题。

http://www.cfd-china.com/topic/1548/

好像跟这个帖子有关。气动噪声算法上我还不太理解。目前只知道是计算压力小波的传递？

你用的是非结构网格？很有意思 哈哈

在你画出来流线图之后，在你原始的.foam基础上，再做contour plot，然后同时显示流线和contour。

你试试运行这个命令：

sudo chmod u+x /home/lnq/OpenFOAM/

有时候相间的物质转换会造成局部高压，然后碰到wall就会反弹。。。
一般non reflective boundary都是用于出口入口的，那么它
可以用在wall上吗。

他是不是想截这个图 :mihu:

替代文字

而且用的是压力基 两相流

谢谢！经过测试，完美地解决了我的问题！

谢谢大神，我看到你原来的代码里有相关的。再次感谢。
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并没有植入，植入很麻烦。现在也不是我研究重点，等有时间在弄这个。