@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这个我没找到，其他的都能找到。

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

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

可以使用inletOutlet或者outletInlet，满足你的要求。

非常感谢 @程迪 的讨论，我重新回顾一下看看能不能分享一些看法。

有关ddtPhi()这个函数的形式大家已经明白了，更重要的是明白为什么要这么做。MULES限制器和ddtPhi都是Henry Weller发明的，MULES限制器的植入以及原因目前都清楚了。但是ddtPhi的原因尚不明了。

@dyj19901127
宇老师大作被翻出来了，哈哈

我大略看了一下：

#include "createFluidMeshes.H" #include "createSolidMeshes.H"

划分两个网格的时候就存在了边界场，在求解fluid之后，边界场信息会更新，你不需要显性的传递边界场。

fvMesh& mesh = solidRegions[i]; solidThermo& thermo = thermos[i]; tmp<volScalarField> trho = thermo.rho(); const volScalarField& rho = trho(); tmp<volScalarField> tcp = thermo.Cp(); const volScalarField& cp = tcp();

这几行代码自动包含边界场，其中给固体的是网格共享的边界场信息。

应该是这个思路， 你可以确认下。