@东岳 对，后来查看了一下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是能做无粘的，但相似的现象仍然有。不知道是哪里的问题。