非均匀温度边界条件可以通过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

顶一下

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办法限制相关量的范围也没有取得很好的效果
各位前辈有没有好的建议

@星星星星晴 谢谢大佬回复，我后面也看看这个头文件如何使用

@星星星星晴 在 Openfoam全系统在vscode下的版本切换问题 中说：

把你的默认of改一下 因为每开一个terminal 就要读~/.bashrc文件

感谢前辈指点！

@李东岳 感谢老师回复，我试了一下，把算例文件夹从虚拟机拷贝到主机里，用Windows的paraView读取算例就可以了，应该就是虚拟机调用显卡的问题。

ubuntu 24.04 OF13 编译 找不到scotch.h
解决办法：https://bugs.openfoam.org/view.php?id=4270
安装bison

找到之前的一篇帖子解决这个问题了，用Surface Vectors filter之后在用glyph
https://www.cfd-china.com/topic/1640/paraview画矢量图-如何将某速度分量设为0

已经在求解器中植入了温度方程求解T，是想要根据每一步的T的变化更新物性参数之后准确的算两相传热，现在大概是这样的思路，求教

两个问题：

你的图里面为什么 wallShearStressMean 这个场在体网格也有值？正常情况下，wallShearStress 场只在壁面上才有值，所以你对比的数据是否确定是 wallShearStress 如果对比的数据没有问题，LES 计算得到的 WSS 显著小于 RANS，这就是典型的 Modeled stress depletion，原因是你的网格太粗，尤其是近壁面的网格太粗，远达不到 LES 的要求。由于没有足够小的网格来解析近壁面附近的小尺度脉动，同时数值上计算出来的 nut 又非常小，表现出来的现象就是WSS远小于实际值。改善的办法一个是加密壁面附近网格，达到 WRLES 的标准；或者，使用 WMLES（参考libwmles库）。

发错了 非常抱歉 请问怎么删除这个帖子呢