@队长别开枪 谢谢指导。
这么看来，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中应该用什么命令？如何实现？

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

@我是河滩 是这样的，代码看懂了，谢谢指导。

@liailei0627 恩恩 谢谢大佬！！

@东岳 最近根据东岳先生的指导，对这个求解器不那么陌生了，也简单做了一个自己的小例子。目前程序能跑起来。拜谢东岳先生！！！
这是我做出来的一点东西：

52da8e03-f286-4265-af05-d6e9fc489eb9-image.png

这个是两相流的含气率分布。但是我看了一些论文在顶部设置的边界条件叫做液相壁面采用无滑移边界条件。气泡达到顶部以上浮速度自由逸出。不太理解这句话对应的OF边界条件是什么？
PS：目前顶部采用的边界条件为：U.air

4b6d193c-82a4-4b02-bcff-35f1c7175f10-image.png

U.water

6476ca4f-42b8-4589-9b2b-cd630e777ae0-image.png

p文件：

1bb31ddf-45fb-4d4e-a0d3-1e02293e4b32-image.png

目前还是依据算例进行自己改动，仍然没达到气泡以上浮速度自由逸出的效果，请问东岳先生对此有什么好的建议么？
叨扰叨扰，再次拜谢！！！

@就是个胖子 嗯，是，可以，只是感觉略微有点麻烦。

@就是个胖子
一些太基础以至于很难搜到的点：

#是注释符号， #!/bin/bash 是指定运行改脚本的shell， $1 是脚本运行的第一个参数，同理$2 $3是第二第三个参数。

下面这个脚本功能是，编译某个文件夹下，所有以.tex结尾的文件。
脚本文件的名称叫 latexall.sh。 执行的时候输入命令"latexall ./"
这时候对于脚本来说$1就变成了“./”即，当前文件夹。
循环.png
这是个根据前后缀特征提取字幕文件里的字幕信息的脚本。这些乱七八糟的符号叫“正则表达式”。
处理字幕.png

以DTChull算例为例，该船6m长，速度是1.668 m/s 自己根据网上的找的公式算出来的k和omega和算例中给定的差别很大，自己用的公式是6ac30593-c076-46d9-8b9b-af132f89a9f5-图片.png
3d1cc183-2b26-40b4-9c6c-e71332a0b616-图片.png

@东岳 李老师太客气啦!:shangxue:

请问大侠，一个对称圆柱的二维计算，怎样计算其出口质量流量，还是说二维计算没有质量流量一说，在preview中计算质量流量，或者怎样转化为三维情况下的结果？

@东岳 如果是像fluent中DPM里颗粒都是mass point的话就不存在这个问题了。

externalWallHeatFluxTemperature，将q设为0.

不过我用这个BC加恒定热流密度，后处理用wallheatflux计算，发现有个min/max区间，integral的值比应该注入的能量少0.3W（48.8W VS 48.5W），不知道这是否正常。