@肖恩曼 看了下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是浓度。第二项是菲克扩散定律的结果。
求大神们帮忙看看！

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

谢谢哈!