@supersoldier OpenFOAM大量使用了C++面向对象的设计思想，不能用一般面向过程的编程思维去阅读代码。了解面向过程的设计思想有助于去理解OpenFOAM的框架和具体实现。

方程实际要求的是：
$$\begin{array}{}\text{(1)}& \frac{\mathrm{\partial }\rho \mathbf{U}}{\mathrm{\partial }t}+\mathrm{\nabla }\cdot (\rho \mathbf{U}\mathbf{U})=0\end{array}$$
对应：

fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U)

连续性方程：
$$\begin{array}{}\text{(2)}& \frac{\mathrm{\partial }\rho }{\mathrm{\partial }t}+\mathrm{\nabla }\cdot (\rho \mathbf{U})=0\end{array}$$
对应：

- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)

因此实际求解的就是：

fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) + turbulence->divDevRhoReff(rho, U) );

把- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)包含进来的原因为这有利于使求解变量有界以更好地收敛。

是不是可以设置一下relTol

@cfd-china 就是能否给定初始场，像fluent每次在计算之前都要进行initialize那样:confused:

nevermind,用gMax代替max就可以了

@cfd-china 我自己写的边界条件里面的，单核计算没有问题，如果不是昨天停电都快算完了，结果也还可以。

@cfd-china 哦，直接后处理的话要挨个打开有点麻烦呀

@cfd-china 厉害，知道怎么配置吗？刚学

你用什么做的网格？如果是icem的话改一下scale factor就可以。

:kiss:

厉害了我的哥哥姐姐们

谢谢东岳大神的回复，看了您推荐的，感觉有点明白了，还在挣扎中。

对了，这个问题我在2015年曾经汇报过bug，Henry直接表示建议采用高版本的OpenFOAM:sunglasses:

bubbleFoam已经消失很久了。然而作为双流体模型的首个求解器，并引以为傲的引入了Henry的速度Phase intensive equation。但是在后续版本中被惨痛抛弃。其中的Phase intensive equation在开发之后解决了一些相分离问题，但是被一些sci文章指出在某些情况下非守恒形式的方程会引起误差。并且在模拟工业规模的鼓泡床、萃取塔的时候，bubbleFoam没有考虑可压缩性（实际上床底和床顶压差很大），最终bubbleFoam被compressibleTwoPhaseEulerFoam（附加颗粒动理学模拟气固）完全取代。

然而bubbleFoam依然是一个经典的求解器，或许身份等同于icoFoam。没人用，但他的意义不在于此。

@xpqiu 谢谢:happy:

gcc在Windows下的表现不如Linux。

最近看了一下，你可以加我QQ435569485，交流一下。

个人感觉区别不大，不过个人习惯在OPENFOAM里面就用scalar,vector等等吧

我想可能的话，可以让slave patch 以任意角度拷贝 master patch

@李东岳 多谢东岳老师，我去尝试一下