@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

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

@bestucan 谢谢回复。我看了一下，OpenFOAM7中是指定一个点的，通过location；OpenFOAM1912中有一个option可选的是locations，可指定多个点

@bestucan

intermediateNew$ wmake libso wmake libso . g++ -std=c++11 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -O3 -DNoRepository -ftemplate-depth-100 -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/lagrangian/basic/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/lagrangian/distributionModels/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/transportModels/compressible/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/specie/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/basic/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/thermophysicalProperties/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/reactionThermo/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/SLGThermo/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/radiation/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/transportModels -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/regionModels/regionModel/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/regionModels/surfaceFilmModels/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/dynamicFvMesh/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/sampling/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/finiteVolume/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/meshTools/lnInclude -IlnInclude -I. -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/OpenFOAM/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/OSspecific/POSIX/lnInclude -fPIC -c parcels/derived/basicKinematicParcel/defineBasicKinematicParcel.C -o Make/linux64GccDPInt32Opt/parcels/derived/basicKinematicParcel/defineBasicKinematicParcel.o In file included from lnInclude/KinematicParcel.H:692:0, from parcels/derived/basicKinematicParcel/basicKinematicParcel.H:40, from parcels/derived/basicKinematicParcel/defineBasicKinematicParcel.C:26: lnInclude/KinematicParcelI.H:127:8: error: prototype for ‘Foam::KinematicParcel<ParcelType>::KinematicParcel(const Foam::polyMesh&, const barycentric&, Foam::label, Foam::label, Foam::label, Foam::label, Foam::scalar, Foam::scalar, Foam::scalar, const vector&, const vector&, const Foam::KinematicParcel<ParcelType>::constantProperties&)’ does not match any in class ‘Foam::KinematicParcel<ParcelType>’ inline Foam::KinematicParcel<ParcelType>::KinematicParcel

sorry我看眼花了，多加了个点这下输出真正的编译错误了

@东岳 谢谢东岳老师！:xiexie:

想在cavity算例模型内部加入一单位大小的矩形悬浮块，想观察它随着周围流态变化的运动过程，请问动网格这块应该设置什么类型，或者of有什么算例能参照下，正在初学动网格，好些不懂。