@李东岳 谢谢李老师

能否根据一些散点的坐标xyz，即这些散点的风压值，生成paraview可以读取的格式让其能插值显示曲面风压？

@wangwang 谢谢前辈，是的发现问题了

温度越界这个问题在可压缩那面出现的太多了。没有一个普适性的答案。这就是发散了。

使用waves2foam造波，interfoam求解器。运行到后面时间步变得无限小，结构受力无限大。我这里有十多个构件单独输出受力，每个构件的受力全都变得很大。我觉得可能是网格质量不好的原因（网格用snappy生成），但是调整了很多参数，依旧会发散。以下是log文件、snappy文件（由于杆件太多，只展示部分杆件，所有杆件设置一致）以及网格局部图片，请问大家有没有调整的建议或者除了网格导致发散的其他原因。86a8b480-7d79-4850-86d4-6140dfede200-image.png 3db69d5e-0327-426e-a358-54a61e4a5a3a-image.png 25c1ee00-0922-4a62-bf33-b7c3741cc376-image.png f112809c-4b63-43d3-b173-d40aca264e39-image.png 22f99966-bf52-4050-8b3e-271e6bcfcb07-image.png db211ec2-f12a-4897-85af-b21f259fd754-image.png 6601262c-29d1-4f28-aaa7-ffd719e3f306-image.png 1287beb0-7bab-4333-b0af-a96fe9381bb0-image.png b6ae8373-35bf-4f95-b812-caa392da8516-image.png 65fcf0a7-e5cf-4748-8487-0d9cdfbeb344-image.png

@cccrrryyy 大佬，请教个问题，最近这两个求解器也困惑了我许久，拜读这个帖子后仍有疑问。
问题1：按照李老师所说，声速较大会导致rhoReactingFoam中采用的压力方程无法求解，那是不是意味着reactingFoam适用于声速可能无穷大也就是低马赫数时密度不随压力改变的情况？
问题2：同时由于rhoReactingFoam采用rhoPimpleFoam的压力方程，其为压力基求解器，捕捉激波的能力没有rhoCentralFoam强，所以其实际上也不适用于高马赫数流动？
问题3：既然这俩都不能计算高速可压流，那基于psi和基于rho的热物性带来的差异主要体现在哪？

@chpjz0391 请问这个multiNormal的文件在哪个tutorials里面？

绘图1.png
李老师，上面算的算例因为网格问题要重新计算了，我参照您发我的算例在ICEM中画网格。如图所示A的严格对称点应该是B点，您画在了C点，所以要在openfoam里面用cyclicAMI来周期上下两个面。出现了面不匹配的老问题，我尝试在ICEM中设定周期，但是好像这种叶珊上下面应该没法设置节点的周期把？随后我仔细观察了您的算例，发现在boundary文件中的cyclicAMI的容差设置为1，随后我在createPatchDict采用了相同的容差。这一步算是可以通过了，但是在decomposePar中采用scotch方法进行分块时，还是在cyclicAMI这里报错了，麻烦李老师给看看我上述的操作是漏了什么嘛？
另附一下 网格处理命令：
fluentMeshToFoam ICEM.mesh
autoPatch 45 -overwrite
createPatch -overwrite

屏幕截图 2023-12-30 235244.png

@李东岳 在 MULES::correct() 和 MULES::explicitSolve() 有什么区别？ 中说：

MULES要看基金会的版本，这是基金会那面提出的。我记得CMULS已经被删了。我在无痛苦ns方程笔记里面说的很齐全。你看看参考一下

好的，谢谢李老师

导师经费很多，可详谈

按照V8的release说明，除了某些特别的需要依赖时间进行平均的函数，其余的比如sample采样等函数会在0时刻（如果startTime是0）先执行得到初始时刻的值。这样在最后的dat中就会有0时刻的数据了。
commit原文在这。

我想在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位数或更少，小数点后边的就更少了。也尝试改了李老师说的这个，但是好像没用。