@李东岳 谢谢老师，我会试试看

已解决，解决方法如下：（有效但不是唯一）
Step1:
cd .config/
Step2:
cd ParaView
Step3:
gedit ParaView-UserSettings.json
Step4:
把图中的 "ColorSpace" : 0, 改成 "ColorSpace" : 1,(把数字改成1就OK了）0_1516159750930_2018-01-17 11-26-32屏幕截图.png

@李东岳 哦哦，是这样。顺带问大哥个问题，我用另一个翼型照wingmotion例子的前处理用simplefoam算，结果time step continuity errors 总在十几步发散，想问下这个时间步连续误差是哪个参数影响的？
替代文字
我做的一些尝试：deltT放到1e-8也是到那一步发散，应该不是时间步的问题。湍流模型也改为层流也是到那一步停止，网格大小也没问题，fvscheme用的迎风，应该也没问题。例子里的网格正交性很好，我用的三角网格，把非正交修正也改成5,10也没什么效果。

@yhdthu 呵呵，OF的精度问题不止如此，湍流模型如果是后加的，算子分裂可能会导致精度下降到1阶和稳态无法达到的问题。所以Jasak在FE中玩耦合求解是另外的套路。

fvc::flux(phi, alpha1, alphaScheme) = fv::convectionScheme<Type>::New ( mesh(), phi, alphaScheme )().flux(phi, alpha1) = phi*interpolate(phi, alpha1) = phi*alpha1f

因此，代码fvc::flux(phi, alpha1, alphaScheme)即为$\nabla\cdot\left(\alpha \mathbf{U}\right)$

https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/finiteVolume/finiteVolume/convectionSchemes/gaussConvectionScheme/gaussConvectionScheme.C

https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/finiteVolume/finiteVolume/fvc/fvcFluxTemplates.C

在 coalChemistryFoam 的 Make/options 文件里，包含了下面一行

-I$(FOAM_SOLVERS)/combustion/reactingFoam \

我对比了 reactingFoam 文件夹下的头文件和 coalChemistryFoam 文件夹下的文件，发现并没有需要调用的。但是去掉这一行之后编译时会给警告

In file included from /home/user/OpenFOAM/OpenFOAM-5.x/src/OpenFOAM/lnInclude/postProcess.H:129:0, from coalChemistryFoam.C:50: ./createFields.H: In function ‘int main(int, char**)’: ./createFields.H:30:23: warning: unused variable ‘T’ [-Wunused-variable] const volScalarField& T = thermo.T(); ^ ./createFields.H:31:23: warning: unused variable ‘psi’ [-Wunused-variable] const volScalarField& psi = thermo.psi(); ^

可是这两个变量在 EEqn.H 和 pEqn.H 中有用到。不知道这个警告是否可以忽略？

@李东岳 好的 谢谢:happy:

虚假流的产生原因是气泡的数值受力不平衡，由于将表面张力&重力差值到后并不能保证面上满足动量守恒，我记得alamos lab有过一篇paper讲这个问题怎么处理，貌似是在staggered grid上做的，可以找找看看

这个问题可以到stackoverflow上搜或者问，纯代码、linux系统问题。segmentation fault是程序访问地址超出范围，然后被限制了，所以出现很多问号。一般这种问题出现在嵌入式开发交叉编译中，就是硬件A（比如电脑）为硬件B（比如安卓）编译可执行文件。当然编出来的在A上是没法运行的。

自己编的gcc有c++库，编译求解器应该用的是自己编的。运行的时候用了系统自带的c库，系统自带的c库又去呼喊glibc-2.11，调用start.S（程序运行通用的汇编代码）就错误了。大概是这个样子。

具体的解决方法，猜不出来:confused: 得摸着电脑慢慢试

@mohui 那我感觉应该可以用，我只是知道，没用过，期待你的结果

大四小朋友 :cheeky:

如果还没决定要用OpenFOAM，如果只是做本科毕业设计，ANSYS Fluent那面上手要快的多。

如果可以读硕士，使用OpenFOAM自带的求解器算一算还是可以的，发文章这个东西影响因素太多了，导师是很重要的一部分。

如果可以读博士，首选OpenFOAM。通过使用OpenFOAM做CFD模拟，会比用Fluent了解更多关于CFD的东西，并且时间足够，也可以出东西。

OpenFOAM超音速求解器那面一直没有赞助，读博士还是有的搞的。

行业和就业前景之类的

这个和学历、期望待遇、是否同意轻度的转行关系很大，很难说。

想了解一下这个坑有多深，难度怎么样，还需要什么先修知识和最注重哪些知识，

总之，个人建议本科毕设不要入坑，很难毕业。硕士基本只能把OpenFOAM当做商软用用。博士才有的搞。并且还得用心搞。

但是总体来说OpenFOAM值得去玩，我这7年来，一直在用OpenFOAM，我学到了非常多有关CFD的东西，虽然刚开始的时候，多次有放弃的念头，但走到现在，我庆幸入行OpenFOAM。

@李东岳 谢谢岳哥~

我的16.04也有这个问题

0_1515589858123_捕获.JPG

采用相同的求解器，上面是并行计算，下面是单核计算，计算结果有偏差。

感谢关注，过年的时候应该可以 :cheeky:

在自己的笔记本上跑算例的时候，如果打开Chrome，经常由于内存不足导致网页自动崩溃。回国赶紧换工作站。

我觉得可能是win10系统占用的内存更高了，或者是新版本Chrome内存占用过高。

@李东岳 我最后还是把这个打开了。我模拟的是潮流能水轮机，开之前功率系数不过30%出头，开了之后到了45%左右（TSR=6），功率系数-叶尖速比曲线基本上与实验对上了。而且开了之后对时间步长就不敏感了，时间步从0.0001s到0.002s，CoN从40多一直到700多，功率系数变化在5%以内。

而且，OpenFOAM/tutorials/multiphase/interPhaseChangeDyMFoam/propeller这个算例也是开启了动量预测，我考虑是不是固定转速的旋转机械都需要开这个选项。

@流卜 的确是接口差别很大，目前准备把extend上的部分内容移到v1706，感觉崩溃了快

@李东岳 0_1515462127425_2018-01-09 09-40-58屏幕截图.png
我按照您说的，我把整个算例搬过来，只是改了一下0文件里面各个变量，在system文件加了几个字典文件，结果又出现了上图的问题。我怕是我哪个地方，出错了，又重新操作了一遍，还是这个问题。我又去运行了一下原来的算例，没有出现任何报错。我看了一下，除了data这个我没找到，其他的都能找到。

对于单气泡的曳力，升力更倾向于是一种面力。但是由于这个面不是网格单元的面。处理起来存在一定困境。我们把这些力当做体积力？有一些模拟结果表明存在震荡。插值做面力（类似基金会这种处理）？但是结果是错的。

需要深入研究力的处理方式。