问题描述清楚的好处，除了“使人理解”（描述不清无法使人理解），还有“使更多的人理解”。有价值的问题值得转化成文本与世界讨论。所以下一步才是体现“描述清楚问题”价值的一步：问题翻译成英文，去cfd online找个合适的板块发帖:xiezuoye:

@bestucan 好的，谢谢

你好，请问这个问题你解决了吗？我也遇到类似问题。如果A程序是稳态问题，B程序是瞬态问题。而且AB还要循环多次求解。请问这个问题该如何解决呢？

找到了，是一个二维湍流下漩涡的精准判据。

@sibo 现在问题解决了吗

@东岳 哈哈，老师，不知道啥情况，本来就一台电脑有这种情况，昨天测试后两台电脑都出现了，最后连虚拟机都变成这样了~幸亏提前找到了病根~要不都不知道啥问题~
alpha.water和alpha.water.orig文件里的object只要都是object alpha.water;在paraview就只有一个alpha.water选项，它实际读取的是alpha.water.orig文件里的object alpha.water;所以显示的是单相；
把alpha.water.orig文件里的object改为object alpha_.water;，然后复制文件去掉orig变成alpha.water文件，这时里面的还是object alpha_.water;但是在终端执行setFields后，alpha.water文件里的object alpha_.water;自动变成了object alpha.water;，并没有人工修改~
这时查看paraview，同时出现了alpha.water和alpha_.water；本来带有.orig的文件不该读取，不知道为什么读取了~这个bug太坑人了；还是要再次感谢 @bestucan ；:xiexie:
814175b86c5411fe631753dbd2548e6.png
1e61b5f81ee34689162e23e7a11d552.png

如题，目前的方法是blockMesh先生成方形单层mesh，再利用extrudeMesh生成360的3D网格。
但是出现个问题，里面的面不会自动的成为interior，不太清楚怎么处理。

以下是extrudeMeshDict

constructFrom patch; sourceCase "$FOAM_CASE"; sourcePatches (side1); exposedPatchName side2; extrudeModel sector; nLayers 10;//90; sectorCoeffs { point (0 0 0); axis (0 1 0); angle 20;//360; } flipNormals false; mergeFaces false;

总感觉缺一个点。
得证明大密度瓜皮可以等效替代低密度瓜瓤。
因为体积、质量都一样，然后质心也在同一点上。

"(U|k|omega)Final"是最终迭代步下的设置，不带Final的是非最终迭代步下的设置

cache用来减少内存调用

以杀死有的时候是内存不够

@东岳 谢谢，我也觉得这个方程很简单，这段代码放到icoFoam里面没有出现过问题。但是放在两种液体混合求解器中就不行了。我感觉边界条件没有问题，流体从一侧流入，另一侧流出。

FoamFile { version 2.0; format ascii; class volScalarField; object T; } dimensions [0 0 0 1 0 0 0]; internalField uniform 293; inlet { type groovyBC; valueExpression "(pos().y>1) ? (295):(293)"; } outlet { type zeroGradient; } freeSurface { type zeroGradient; } bottom { type fixedValue; value 293; } front { type symmetryPlane; } back { type symmetryPlane; } FoamFile { version 2.0; format ascii; class volVectorField; location "0"; object U; } dimensions [0 1 -1 0 0 0 0]; internalField uniform (1 0 0); boundaryField { inlet { type fixedValue; value uniform (1 0 0); } outlet { type inletOutlet; inletValue uniform (0 0 0); value $internalField; } freeSurface { type pressureInletOutletParSlipVelocity; value uniform (1 0 0); } bottom { type fixedValue; value uniform (0 0 0); } front { type symmetryPlane; } back { type symmetryPlane; } }

这种手动去掉网格单元进行计算的，有一些方法可以处理。https://www.cfd-china.com/topic/3138 类似浸没边界条件，但边界需要处理。也可以直接操作矩阵，但是操作的是边界值，比如OpenFOAM中的壁面函数的epsilon场。你可以尝试一下。看下用在内部场效果如何。我之前见过一些人进行这么处理。

封闭方腔自然对流

这个不应该出现问题啊，网格很规则，有实验结果可以对比么？还是自己算着玩的。有实验数据对比我可以给你调调，网格少的话

公式就是那几个，你也知道的

好的 谢谢

inGroups、list可以删掉，没用

@bestucan 感谢回复。看来是我瞎操作对了，运行fluent3DMeshToFoam以后，再在boundary文件和0文件里设置边界条件

虽然我们用 emacs 的和用 vim 的向来不和，但是为了合伙掐死用 vs code 的还是可以站在一起的๑乛◡乛๑ 

emacs插件：疑似OpenFOAM的开发者之一 Henry 的配置文件。
配置里还有很多 OpenFOAM 自带的命令，这下破案了，那些命令根本不是记住然后手动在bash里用的，是写到脚本里，绑到快捷键上调用的。

vim插件：vim三种插件。

躲在角落里擦一擦自己的机械硬盘，还能用emacs再续两年……
:mianmo:

哈哈，之前CN格式那个也有一个笔误

@东岳 好的！谢谢~

最后还是用了手动赋值
a9878440-a6cc-4560-a822-8956f3fbc06d-image.png

@东岳 谢谢老师建议! 回头再遇到类似的问题, 我多看看报错, 这样子更有针对性