我想在openfoam里面引入液体tait状态方程,该怎样实现啊？make文件夹怎么编写？

ee6ad707-d3e6-48f0-b4cf-ab27fb7423f0-image.png
计算虚拟机和本地计算机共享的文件夹里的算例时候卡顿，尤其是并行计算后进行文件合并，打开终端需要等半分钟，有谁有解决办法吗？

@wzlr2023 我看你这个算例里只有两相流，正常用 interFoam足够了。这个相关资料还挺多的。这个多相流的求解器正如我前面所说我没有去验证其是否可以正常求解。就你这个算例来讲，里面涉及到了六自由度运动求解。造成两相交界面求解发散的问题也可能是六自由度求解引起的流场发散。建议可以先设置一个简单的动网格算例，如动边界造波来验证多相流求解是否稳定。如发现是多相流求解出问题了，欢迎进行补充修正。

我使用CFDEM中的解析CFD-DEM模型+RANS+k epsilon计算水平圆管中颗粒沉降的时候(管径115mm，入口流速3m/s，粘度0.01Pa·s，密度1000kg/m3；颗粒密度1100kg/m3，直径15mm；雷诺数几万)，发现颗粒越靠近圆管底部，下沉速度越慢，甚至在距离圆管中心0.4mm的时候，颗粒没有下沉，反而向上移动了。想请教一下各位老师，这种情况是不是不对呢？
印象中层流里面，管壁附近重颗粒才会往管道中间移动，不知道我这个认识对不对。

@luca 在 请教在自开发AMR求解器中应用cyclic边界条件报错 中说：

单核是可以运行的，但是一旦并行运算就会报错。

我在使用并行计算时也遇见了相同的问题，题主解决了吗？

学流体的小明 在 interFoam计算气泡槽道流时的压力问题 中说：

应该是 并行 + 自适应网格 的问题？

更正一下，是 并行+自适应网格+压力参考点 三个要素的共同作用，去掉任何一个都可以算出来好的结果。

一个新的发现是壁面上会有奇异点，如下图，着色是p_rgh，这些奇异点都处于并行分区的各个界面上。更奇怪的是，上壁面并没有这样的问题。Note:参考点放置在流场中间。
01292e97-19c8-45be-889b-189c0ff6fe54-image.png

现在暂时放弃自适应加密了，直接全场加密算，压力就是正确的。
😂

楼上说的对，嗯，这是面法向梯度

可不可以结合codestream实现旋转呢？

@Samuel-Tu 请问你解决了这个问题了吗？我是在controlDict中添加code进行计算得到的一些自定义的变量，但是结果在写入到文件中时，除了时间步数之外，其他的参数就一直是6位数或更少，小数点后边的就更少了。也尝试改了李老师说的这个，但是好像没用。

@zhoushan 简单的几何可以用blockMesh直接生成网格，复杂的就要绘制stl文件，用blockMesh搭配snappyHexMesh生成

fluid@zhoushan:~/$WM_PROJECT_USER_DIR/FluidSolidInteraction/src$ ./Allwmake

wmake libso fluidSolidInteraction
./Allwmake: 4: wmake: not found wmake solvers/fsiFoam
./Allwmake: 6: wmake: not found wmake solvers/ampFsiFoam
./Allwmake: 7: wmake: not found wmake solvers/weakFsiFoam
./Allwmake: 8: wmake: not found wmake solvers/fluidFoam
./Allwmake: 9: wmake: not found wmake solvers/solidFoam
./Allwmake: 10: wmake: not found wmake solvers/thermalSolidFoam
./Allwmake: 11: wmake: not found wmake libso utilities/functionObjects/pointHistory
./Allwmake: 13: wmake: not found wmake libso utilities/functionObjects/patchAvgTractionHistory
./Allwmake: 14: wmake: not found wmake libso utilities/functionObjects/centrifugalBodyForce
./Allwmake: 15: wmake: not found
请问上述的not found怎么解决呀 在网上参考了各种方法 都无法成功

最近在做多相流的仿真验证，但是一直找不到LS方法的多相流案例，想请教以下各位大佬openfoam中有没有关于LS方法的多相流案例

如上所述

这个问题我也遇到过，比如对于同一个流场，同样的空化数，来流速度10m/s就不会产生负压，来流速度20m/s就有负压产生。增加SS空化模型的系数n的值或许能解决问题。

最近在植入一种DES混合模型。期望我植入的DES模型的时间平均速度，与RANS一致。对比后发现不一样。debug几个小时没发现问题。索性采用OpenFOAM原生的Smagorinsky算了一下。发现原生的Smagorinsky模型计算的时间平均的速度跟RANS也不一样。

简述：OpenFOAM的pitzDaily算例，用2DLES计算与2DRANS计算时间结果不一致，我怀疑是不是2D LES的问题？

因为我这跑不动3D LES，所以没法做3D LES与2D RANS的结果对比，大半夜的我感觉啪的一下就起来了，很快啊！问问各位大佬有何高见。

屏幕截图 2023-12-15 080858.png

下图是2D RANS的结果：图片出处
替代文字

下面是我用原生的2D Smagorinsky算的时间平均速度：

屏幕截图 2023-12-14 171725.png

下图是我自己植入的2D DES算的时间平均速度：

屏幕截图 2023-12-14 175456.png

目前没有3D LES的结果。纵然我有万核CFD服务器在闲着，但debug起来太慢了主要。。可以理解为跑（debug）不动，没时间持续研究这个事情。

:136:

今天偶尔看到一个文章。好像意思是通过Optuna + OpenFOAM，可以实现最小的计算时间。感觉类似通过Optuna跑好几次，然后寻找出来最合适的松弛因子。不知道有没有哪个大佬感兴趣玩一玩。

https://ss1.xrea.com/penguinitis.g1.xrea.com/study/OpenFOAM/relax_factor_optimization.html

@wangfei9088 在 paraview 不显示坐标轴数值 中说：

@liujm 假设你用的版本是OpenFOAM-10，其他版本也一样。
路径：OpenFOAM-10/etc/config.sh/
打开paraview文件，将69行取消注释，注释掉72行。
./makeParaView重新编译就是5.0.1版本了。

#export ParaView_VERSION=5.0.1 //69行 #export ParaView_VERSION=5.4.0 #export ParaView_VERSION=5.5.0 export ParaView_VERSION=5.6.3 //72行

感谢！

@LGY 请问一下你这边能够在wave2Foam中生成非定常风了吗

@gooseEast 在 rhoCentralFoam发散，出现Maximum number of iterations exceeded 中说：

在thermoI.H文件第46行标量f是什么

OpenFOAM代码中，这个f是Cp，定压比热容，单位是[J/(kg K)]。这个T函数是通过定压比热容的值迭代求出温度。limit，F，dFdt分别对应的是hConstThermoI.H文件里limit，Cp，dCpdt三个函数。

@李东岳 老师，如果不使用refineMesh的话，结果会好很多，而且和论文中的Cd, Cl指数可以基本对上，但是在overset的边缘处还是会有一点问题。
648bb924-7f1b-4e8c-a39d-48cc7c40906e-image.png

我的oversetInterpolation的method使用的是inverseDistance.