QQ截图20190318163023.jpg
压力云图是这样的。

@东岳 好吧

@东岳 在 OpenFOAM多组分返回热物性 中说：

实际就是rho_

这样！明白了明白了，谢谢。

@东岳 好，谢谢，得啃求解器了，

给大家分享一个编译安装的例子，在Ubuntu18.04上安装OpenFoam2.3.1
https://openfoamwiki.net/index.php/Installation/Linux/OpenFOAM-2.3.1/Ubuntu#Ubuntu_18.04

@cccrrryyy 非常感谢

还有 轴缩放

@LiuZhuang 在 twoPhaseEulerFoam下的bubbleColumn算例的问题 中说：

在outlet边界上气体自由逸出，液体则为壁面（zeroGradient）

好长没太细看:zoule: 不过我觉得你可以设置液体壁面速度为slip试试

1.jpg
图片是CES中 作者用openform 模拟的单气泡上升时的 质量传递 过程。 云图为氧气浓度分布云图。
想问下老师们，这个能用Fluent模拟嘛？ Fluent能做出类似的氧气浓度分布云图吗？

@gogo 请问你现在有在边界层生成功吗？我看了一下你的文件，觉得你的边界层网格设置的不好，finalThickness太小，这样即便生出来长宽比也会很大

@队长别开枪 谢谢指导。
这么看来，OpenFOAM对可压缩跟不可压缩的问题实现了代码的复用。
如果是不可压的，收敛后∇⋅U=0，2/3这项自动为0，也就省略了，不知道这么理解对不对。

那么对于动量方程2.png
这种形式的动量方程，没有体现2/3μ这一项，是因为不够严谨，还是因为这一项因为∇⋅U=0而省略掉？

如果是后者的话，那么这就变成不可压缩的动量方程？那为何还有密度这一项？

还是说这是一个统称的写法，具体还要根据问题判断。谢谢。

在做VoF 2d模拟的时候发现z方向网格的厚度对alpha.water的值居然有影响，厚度太小，alpha值会略大于1，不知道大家有相同的经历没。有点想不通啊，方程在z方向应该是不进行求解的啊，厚度怎么会有影响呢。

@东岳 哦，明白了，谢谢李老师。

@队长别开枪 多谢老师

东岳老师您的邮箱多少？关于twoPhaseEulerFoam的问题我有几点想跟您请教

我在网上见到的代码，你试试：

const cellList& cells = mesh.cells(); const vectorField& faceAreas = mesh.faceAreas(); List<bool> outsideFaces(faceAreas.size(), false); forAll(mesh.cellZones(), zoneI) { const labelList& cellLabels = mesh_.cellZones()[zoneI]; outsideFaces = false; // mark all faces that are NOT internal to the cellZone: forAll(cellLabels, i) { const cell& c = cells[cellLabels[i]]; forAll(c, cFaceI) { const label faceI = c[cFaceI]; // xor operation // internal faces get marked twice, outside faces get marked once if (outsideFaces[faceI]) { outsideFaces[faceI] = false; } else { outsideFaces[faceI] = true; } } } // now calculate the area scalar zoneOutsideArea = 0; label zoneOutsideNFaces = 0; forAll(outsideFaces, faceI) { if (outsideFaces[faceI]) { zoneSurfaceArea += mag(faceAreas[faceI]); zoneOutsideNFaces++; } } Info<<"zone:" << zoneI << " nFaces:" << zoneOutsideNFaces << " area:" << zoneOutsideArea << endl; }

https://www.cfd-online.com/Forums/openfoam-programming-development/69661-how-compute-cellzone-volume.html

inlet { .. code #{ forAll(patch().Cf(),i) //遍历整个边界 { Aera+=patch().magSf()[i]; } info<<Aera<<nl; //面积输出到界面 }; .. }

上面是在0/U中，对入口的面积进行计算，利用simpleFoam非定常计算时，time=1, Aera输出了3次，后面的时间Aera都输出了两次，请问这是为什么啊？

@conghongchuan 找到了这个，但是具体的内容我还没研究
https://blog.csdn.net/hedongya/article/details/78976283

想对物面边界的面积进行计算，在Openfoam中应该用什么命令？如何实现？

对称轴可以不设任何边界条件，因为轴上没有网格面。