顶一下

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库）。

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

@xpqiu 谢谢大佬。这个问题困扰我两天了

畸形波无论是在fluent还是在openfoam中模拟方法都未有明确作法，请大家指导一下模拟方法。

好的，谢谢！

@李东岳 李老师我在文章里面没看到cent-offline.tar文件在哪里下载？不知道能否重新发一下

在pEqn.H中对U进行更新时有

const surfaceScalarField phig(-rhorAUf*ghf*fvc::snGrad(rho)*mesh.magSf()); U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf);

我想请问其中的p_rghEqn.flux()代表什么意思，以及为什么phig中要先乘rhorAUf再在U中除以rhorAUf；
以及之后的根据密度差判断收敛性中，我发现compressibleContinuityErrs.H代码如下，只计算了误差值，但没有判断是否符合要求，想请问各位大佬是在哪里判断的呢

dimensionedScalar totalMass = fvc::domainIntegrate(rho); scalar sumLocalContErr = (fvc::domainIntegrate(mag(rho - thermo.rho()))/totalMass).value(); scalar globalContErr = (fvc::domainIntegrate(rho - thermo.rho())/totalMass).value(); cumulativeContErr += globalContErr; Info<< "time step continuity errors : sum local = " << sumLocalContErr << ", global = " << globalContErr << ", cumulative = " << cumulativeContErr << endl; }

谢谢大家！

@李东岳 求教一下东岳老师，可压缩流动分离解中，密度的更新顺序，为什么有的在压力PISO循环中用状态方程更新密度，有的用连续方程更新密度，有的完全不更新密度，再PISO完成后再更新密度，这是出于什么考虑？如果简化为正压流体的话，忽略能量方程，顺序会有变化吗？

@xpqiu @李东岳 我用NASA的建议 感觉会更好点。我把对称边界改成slip之后bounding就解决了。之前计算的时候开动量预测会浮点溢出，然后查了下说开不开影响不大，就一直没开。现在我把动量预测打开了，可以正常计算。谢谢两位老师

我只是有空看下论坛，随缘回复一下我能帮助的问题。平时没什么时间。:135: