想在集群上安装openfoam，看到东岳老师的公众号链接https://mp.weixin.qq.com/s/0mZdvB280cwVKVA_GZya9g
但是里面没有提供cent-offline.tar下载链接，能否有大佬提供一下这个包

目前网格质量在ICEM中为0.6，在openfoam中求解不可压缩情况时可以收敛，但求解buossinesq假设情况下的方程就会发散，速度取1m/s时，y+值为100，取0.003m/s时y+为1，且不论采用komegasst模型或是湍流模型均会发散，想请问各位大佬是否与y+值过大有关系呢？

我参考HZDR的代码弄了个setbubble小工具，让气泡初始化就是圆的：5be8a963-f837-4652-9c92-a44ffd016afd-image.png
但是我的问题是进行计算之后貌似界面的细化很差，很多界面细化网格级别没达到最低，不知道您怎么解决的这个问题：
2964e223-615d-4d60-81f0-7a70f7db825d-image.png

如果有结构物造成的网格畸变，建议更换为overset。

@random_ran 老师您好 请问您论文题目叫啥呢 我打开链接找不到您的论文

非均匀温度边界条件可以通过codedFixedValue轻松实现，
带壁厚的温度边界条件可以通过externalWallHeatFluxTemperture实现。
如代码所示：

\\非均匀温度边界条件 wall { type codedFixedValue; name dummy; code #{ const vectorField& Cf = patch().Cf(); const scalar a0 = 0.1; forAll(Cf, faceI) { const scalar y = Cf[faceI].y(); (*this)[faceI] = a0*y; } #}; } \\薄壁温度条件 wall { type externalWallHeatFluxTemperature; mode coefficient; h uniform 1e10; Ta 1620; thicknessLayers (0.015); kappaLayers (4); kappaMethod fluidThermo; value $internalField; }

但是，当我想像下面这样编程使Ta为非均匀值时，会报错，即使我使用了codeInclude codeOptions codeLibs之后，仍会报错 ‘patch’ was not declared in this scope，似乎在externalWallHeatFluxTemperature中，无法访问到网格信息。请问各位大佬有解决办法吗？或者有没有别的办法实现带薄壁的非均匀温度边界呢？

‘vectorField’ does not name a type ‘Cf’ was not declared in this scope Ta coded; code #{ const vectorField& Cf = patch().Cf(); const scalar a0 = 0.1; forAll(Cf, faceI) { const scalar y = Cf[faceI].y(); (*this)[faceI] = a0*y; } #};

各位老师好，之前我都在of10中修改求解器，但最近换到了of13，我希望去更改buossinesq假设时的求解器代码，但是对于fluidsolver和isothermalfluid的继承感到困惑，在官网看到fluidsolver继承isothermalfluid，但在fluidsolver中我并未找到相关代码，因此想请教各位在模块化代码下要怎么去修改求解器，谢谢大家！

各位老师好，

最近在使用K-omega SST湍流模型处理两相流的时候发现自由液面处的速度异常，但是不影响波面，如图所示。请问，这种情况是什么引起的，会影响计算结果吗？要如何解决？

万分感谢！

662a54c8-2f33-4bff-ac92-c8c5059e396d-image.png

@李东岳 很奇怪，又画了个很稀的6万网格量的网格，先前会发散的极低雷诺数，现在却不发散了。看来雷诺数极低，网格就不能太密？

通过网盘分享的文件：2D_21cylinders_Re0.0001.zip
链接: https://pan.baidu.com/s/1NdfFZrV2k_T8mAY7fA7ktg?pwd=8888 提取码: 8888

@李东岳 老师，请问我想算空化泡溃灭，但似乎compressibleInterFoam并不包含传质项？是要自己改成compressibleInterPhaseChangeFoam嘛？

@liujm 。。。。。。。。。，除了你画图的时候浮体的重心位置变化了，其他设置都没啥变化的。

simpleFoam求解的定常问题，用的空间离散格式都是二阶的，审稿人让验证网格的二阶收敛精度。不知有没有人用过paraview将稀网格数据插值到密网格上并算过残差E和收敛阶次p？另外对于很复杂的三维网格，网格的非均匀性对收敛阶次公式的计算又该如何考虑呢？
35fb0dd0-05d4-4fcf-b850-3b3e8d317042-image.png

be5eddf7-3dd1-41c7-b320-e9fa9531fffc-image.png

Hi，

我这现在在做欧拉拉格朗日模拟，2-way。然后在模拟中途就爆掉了。。检查了一下flow，发现整个flow都乱掉了，因为parcel的聚集，毕竟两者是相互影响的。

但是神奇的是，同事把PIMPLE改为了PISO就过了。
理解不了为什么。

flow.0068.png

flow.0069.png

flow.0070.png

flow.0071.png

flow.0072.png

flow.0073.png

下面两个图是根据log文件画的，pimple 的拟合造成了温度场和压力场的变化？理解不了？求解惑

WhatsApp Image 2025-09-25 at 12.08.19.jpeg

WhatsApp Image 2025-09-25 at 12.08.07.jpeg

请教一下各位大佬前辈，最近遇到了一个难题，就是有一个从.stl文件中导入的triSurfaceMesh对象（黑色三角形），我需要用他去切割cell单元（蓝色部分），只需要找到triSurfaceMesh和cell的edge的交点然后连线就可以（红色四边形），这样每一个被切割的cell都有一个被切割形成的面，我应该怎么得到这些面，非常感谢。
70c7412a-f1de-4e99-87a1-b2ff6fbcff39-image.png

顶一下

rhoCetralFoam求解器可不可以增加MRF功能。目前，有没有哪个版本做出更新，希望知道的大佬能给点提示

2bfdf92c-d914-442d-ac28-cb584602383e-image.png
在网上查了资料Openfoam残差定义是这个（代码我没找到），和楼主发的一摸一样，这种残差定义的相关文献有吗。
为什么有些openfoam论文里面残差定义是两次迭代值之差初以迭代值呢
58713e7d-182b-458d-8058-defca3c19baa-image.png

目前在用open foam9版本的rheotool求解粘弹性单气泡上升问题，
使用PPTLog模型或者FENE-CRLog模型稳定性都很差，Co数缩小到0.05依然发散，GiesekusLog模型能够收敛但是气液界面看起来很怪，不是光滑的水滴型。
a3619796-3abb-4e40-9a05-b1c5fb886877-image.png
仿真主要设置参考的官方case，即theta的散度离散采用高阶形式，
PPT模型和fene模型主要发散都发生在theta计算特征值或者计算应力时候发散，使用clipping办法限制相关量的范围也没有取得很好的效果
各位前辈有没有好的建议