@东岳 好的，谢谢东岳大神

感谢分享！

buoyantPimpleFoam大体上可以参考rhoPimpleFoam 不过既然你们都说要buoyantPimpleFoam 我可以写一写 http://dyfluid.com/rhoPimpleFoam.html 网上也有一些其他的资料应该，不过自从我自己写自己的推导之后，就没看过别人的 :qichuang: 这个求解器只不过就是可压缩求解器+传热，没啥大的特殊

@东岳 好的，谢谢东岳大神

@xiaolin 是的

@宝丁 请问这个边界条件中指定的压力值是总压吗？

@nuaa_zp
自己的计算模型及温度边界条件设置：
sheet3D1.png
T1.png
T2.png
报错信息：
error.png

@东岳 问题已解决，感谢东岳老师回答，是我自己代码输入有问题，不好意思。

问题描述清楚的好处，除了“使人理解”（描述不清无法使人理解），还有“使更多的人理解”。有价值的问题值得转化成文本与世界讨论。所以下一步才是体现“描述清楚问题”价值的一步：问题翻译成英文，去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

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

@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有什么算例能参照下，正在初学动网格，好些不懂。

@bestucan 好的，我试一下，谢谢前辈

@bestucan 好的谢谢。感觉你说的挺有道理的。之前尾部网格是很扁的矩形（因为存在边界层），结果微小的幅度旋转都会使网格拉扯出现负体积，所以我才逐步减小尾部网格的横纵比，让网格看起来不会那么扁，成功实现了8°范围的俯仰，然后为了做更大的俯仰角度，继续减小网格横纵比，试图可以让他多转一些角度，当网格的横纵比减小到现在的情况时，并没有出现网格畸形，也出现了发散，就不知道是为什么了。

@星星星星晴 可以留个联系方式

@cccrrryyy 非常感谢~

@东岳 好的，谢谢李老师。今天早上破案了，OpenFOAM7成功编译了。

@东岳 感谢李老师！！是我边界条件设的不合理加上网格密度取值的问题。已经更正。谢谢李老师耐心解答！！不可压的alpha就是groMetricOneField，我想太简单了，以为稳流模型把空气自动过滤掉了。

@东岳 好的，谢谢老师～～那主要还是从模拟结果上分析最为有效吧～～

@东岳 好的东岳老师，谢谢东岳老师，实在不好意思，上次断网没发出去刚发现。

@东岳 对的，我也正在琢磨在cfd-online和stackoverflow上推一下的，此插件支持双语版本，根据用户的vscode界面语言选择使用中文还是英文显示

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

@东岳 顶部#include已加

我找着我第一次发的方程改得，就是第二个block的代码，算出来的结果和MPPICInterfoam的基本一样。当然有可能我的粒子直径比较小，对计算结果影响不大？

我没明白您下面写的这个是什么意思。如果不看mrf的话， 你这个和interfoam相比就是减上了一项

fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)

这一项就是65293cb6-a527-49eb-bdc4-2355ec82c507-image.png
我不知道怎么推出的这一项？

https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/OpenFOAM/fields/GeometricFields/GeometricField/GeometricBoundaryField.C

这是边界场

@星星星星晴 谢谢！！

@红豆沙 太感谢了 啦 确实是这个命令，，半年没用OF 都给忘了 命令了 这错误犯的，，害 丢人列，，:mianmo:

@cccrrryyy 哦哦，好的。感谢大神的解答！

@pcxu 在 关于interFoam dambreak算例的一点疑问 中说：

defaultFieldValues
(
volScalarFieldValue alpha.water 0
);

看引用的部分
在上面setFieldsDict里，设置的默认场就是 alpha.water 0
所以water对应的就是0

你也可以把水设置成 alpha.water 1
这样水的参数对应的就都是1

@东岳 果然，linux版本的和windows版本的不一样。
linux在skip zero点掉以后，可以显示alpha.water的初始场了。
我之前用的windows版本的，0时刻就显示不出来。
至于paraFoam能不能编译出来，随缘吧