@史浩 我不是大佬。最近在做一个VAR的项目，考虑温度场-电磁场-流场-组织耦合计算，需要自编程序，有程序验证的问题。对OpenFoam，我没做过这方面的案例，但我想考察2个方面：一是看一下OpenFoam与自编程序的结果对比，二是看一下OpenFoam与自编程序的计算速度对比。
ESR和VAR比较类似，有相通的地方，你如果做了OpenFoam算例的话可以交流一下。

@史浩 xiexie

在这里讨论buoyantPimpleFoam，rhoPimpleFoam和buoyantPimpleFoam大同小异

首先，buoyantPimpleFoam解析中的方程(3)中的$A_{\mathrm{P}}$包含密度$\rho$；

1/Ap也是除掉密度$\rho$的

不可压缩算法和可压缩算法的压力方程不一样，可压缩算法的压力方程需要乘以$\rho$，不可压缩不需要，buoyantPimpleFoam解析中的方程(15)都需要乘以$\rho$，你看那个$\frac{\rho^{n}}{A}$，所以，

你会发现surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU))中出现了rho*rAU

@桎梏 我没有报错:quwan:

@Dingcy 谢谢

@夜阑烟寒 您好，我是把这个曲面展开成这个样子的，之前的在tecplot里面变换了坐标，还是曲面，想变成这种平面
QQ截图20190423005626.png

@Dingcy 哈哈我是在用低雷诺数湍流模型算超临界CO2流动传热的例子，物性变化很剧烈，不好收敛。

@东岳 就是一个二维方腔的自然对流，但是不太清楚这个温度边界条件该如何设定

入门不久，尝试着答一下。
我认为openFoam中的边界条件大致分两类：

强制性边界条件，直接规定数值，如fixedValue，fixedMean，这类边界条件的信息会一直存在于求解所得的结果文件，算很多步以后查看文件就会发现边界条件并没有变化； 参考性边界条件，这类边界条件大多是第二类边界条件，如出口定义inletOutlet， 对于流出的气体而言，该边界其实是zeroGradient，这时定义的value只是初始迭代的时候给的一个初值，该边界的数值由流场内部的点插值得到

1/3次幂咯

@bestucan 好的，谢谢你的指点~:ok2:

@东岳 我的边界是没有网格的，但是默认的显示就是在紧挨边界的网格上，就是你说的内网格呗，但是同样的算例，Tutorials里就能把那个数据给放到壁面上，应该是经过处理了吧。

@OF初学者周 最后没有用gdb

感谢分享 :xiexie:

@hurricane007 这个思路挺不不错。我没见过手动分区的例子，不懂了 :zoule: 加油。

最近想做一个简化的模拟，需要一种边界条件，在这个边界上，某个组分的流量等于min(最大可能的流量，一个常数）。考虑准稳态情况下，菲克扩散定律 J =-D gradient (C)*n. 在使得边界网格内这个组分不为零的情况下，从这个网格向边界扩散的流量应该小于这个网格接受其他网格扩散的量。
那么问题来了…… 这个边界上的网格，除了边界面，其他面上的扩散流怎么求……求助一下大神们……

@xpqiu 想请教一下，如果我想找使用了某个模型的算例该怎么找。例如我现在想找个使用了LangmuirHinshelwoodReactionRate 模型的算例能不能找到呢，这个模型在：
/src/thermophysicalModels/specie/reaction/reactionRate/LangmuirHinshelwood/LangmuirHinshelwoodReactionRate.H

可以算一下compressibleInterFoam下的传热:duang:

@东岳
点赞了，自己写出来的第一个OpenFOAM代码。。。

对的，实质是求气泡的质心。