而且用的是压力基 两相流

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

谢谢大神，我看到你原来的代码里有相关的。再次感谢。
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字符再赋值，由于之后的计算我都用的是这个错误的湍流模型文件，导致阻力系数求解一直出现问题。后来一个一个文件排查的时候发现是这里的问题。

paraivew据我所知不能。但在C++里一切皆有可能，但是可能你需要自己写点代码。

@dzw05 明白了，也就是说，water并不是一个真正原始的量，这个量的值实际上是由alpha1来指定的，那么我想输出water，直接就输出alpha1就行了。我刚刚试了一下，确实是这样的。多谢啦~~~~

@bestucan 那也是学长，厉害了~

@东岳 这个模型是根据实际的物体建的，所以模型重心就直接取了实际物体的重心坐标，没有数学公式。

@dzw05 谢谢您的解答，讲解的很明白，我还要看一看，如果有新想法，会回复您的。

明白了，谢谢！

@东岳 谢谢东岳的耐心:baobao: :baobao: 感动！

@肖恩曼 看了下paraview手册，实现的逻辑一样，就是用python实现。画图可以直接调用paraview，比用其他绘图软件方便。

但是我设置数组变量的时候也是您这样设置的，然后这样怎么能够分别设置边界条件呢？

我知道了，你要在创建数组的时候，自动给outlet，inlet之类的创建边界条件？这没法实现。只能设定为MUST_READ，然后手动在0文件夹里面设定。或者，就是在求解器内部强制赋值某边界条件。

而且我要的变量场是因变量，就是需要在0文件中设置初始时刻的边界条件。

这和上文不是相符的么？

@东岳 好的，李老师。我仔细看看ICEM。

嗯嗯，好的，谢谢岳哥！

@bestucan 谢谢你~我去仔细研读一下~

只能看出来是发散了，无法更深的debug，对结果影响很严重。

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

学生在用simpleFoam 求解NACA0012翼型的升力系数和阻力系数时，所得结果Cd100多，Cl 20多。学生采用速度1.6m/s，攻角8度，弦长1m，SA模型。U文件中入口和出口类型均为freestream，翼型为fixedValue;P文件中入口和出口类型均为 freestreamPressure;nut和nuTlida文件中参数为0.14。controlDict文件密度输入1.225,lref1,aref0.08。学生是刚学习openfoam，很多东西还不是很清楚。请各位老师帮忙看一下问题出在哪里，不胜感谢！

@东岳 大致明白了，这里相当于是用新的算法代替了传统的散度操作，从而在能达到原有求散度目的的情况下避免了直接求散度造成的麻烦。多谢前辈！！！

@东岳
谢谢，我也做了类似的仿真，相比sonicfoam， rhocentralfoam似乎结果得到的时间更短。

对于简单的问题，例如无黏流流经wedge一类，两个求解器的结果其实是相似的，而且较为满意。

但是，不知道您有没有发现，在这个喷管的模拟中，在打洞的地方，喷管的收缩段似乎出现了速度为0的区域，很像边界层的感觉。

然而我做的是inviscid flow我设置mu=0。我看了rhocentralfoam的代码，感觉rhocentralfoam是能做无粘的，但相似的现象仍然有。不知道是哪里的问题。

@yhdthu 现在很多实验研究都证实了这种压力波动，并且影响到空泡的周期行为

@yuan_neu https://github.com/OpenFOAM/OpenFOAM-dev/tree/master/tutorials/multiphase/compressibleInterFoam/laminar/climbingRod 最新版本代码已经把climbingRod算例换位置了。所以曾经的连接失效了。5版本就有着个求解器compressibleInterFoam，不过没这个算例climbingRod

@东岳 谢谢。这个好像是针对2.2版本的。改一下还是可以用在2.3版本。

感觉是个非常简单的小问题，目前你应该已经解决了。。

@散漫守望2016 openfoam3.0.1下，自行设置了一个算例，入口速度2m/s， 利用simpleFoam求解器求解流场的话，得到的流场从数值大小角度来看是合理的；而同样的初始及边界条件（网格也采用同一套网格），换个求解器DPMFoam单纯求解流场（关闭颗粒发射），得到的区域内最大流场速度却达到了惊人的百米量级 ，也用其它版本的of算过，DPMFoam求解出的流场数值都在百米量级左右。结果确实怪异！

DPMFoam求解器求解流场也有什么限制吗？

你可以看下这篇文章
Direct numerical simulation of near nozzle diesel jet evolution with full temporal-spatial turbulence inlet profile

我先做一个从水中剥离一氧化二氮的模拟，准备修改bubbleColumnEvaporatingDissolving。 但是我的收敛有问题，然后我有些不太理解的地方，在初始值的文件夹里water.liquid, water.gas, air.gas, air.liquid到底代表什么？water.liquid + air.liquid = 1?

还有我的T一直不收敛。谢谢大家帮忙

罗宾边界条件是对流扩散方程的一般形式的绝缘边界条件。这里，边界处的对流和扩散通量总和为零：

0_1522589151720_a1f0a086-969a-48dd-ac44-d23c64e49ea4-image.png

其中D是扩散常数，u是边界处的对流速度，c是浓度。第二项是菲克扩散定律的结果。
求大神们帮忙看看！

那我等着毕业加入你们团队@东岳

谢谢哈!

esi最新发布的isoAdvector方法碉堡了

如果你说的是这个，那么在OpenFOAM-1712里面有。