CFD中文网

@李东岳 感谢老师回复，我试了一下，把算例文件夹从虚拟机拷贝到主机里，用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

@chenboyao 非常感谢您答疑解惑，我马上去openfoam里去尝试一下。

已经在求解器中植入了温度方程求解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完成后再更新密度，这是出于什么考虑？如果简化为正压流体的话，忽略能量方程，顺序会有变化吗？

对的。第二张图中的总流量就是单个喷嘴的流量（20kg/s 除以 10）

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

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

没有低版本的openfoam，你可以试下下面的方法。
理论上来说在applications/solvers/lagrangian/DPMFoam/DPMTurbulenceModels/DPMTurbulenceModels.C中加上其他湍流模型应该就可以，比如：

#define makeLESModel(Type) \ makeTemplatedTurbulenceModel \ (singlePhaseTransportModelPhaseIncompressibleTurbulenceModel, LES, Type) #include "kEpsilon.H" makeRASModel(kEpsilon); #include "realizableKE.H" makeRASModel(realizableKE); #include "Smagorinsky.H" makeLESModel(Smagorinsky);

@雨罢_清宵 你好，我在模拟论文中的微通道泊肃叶流动时想使用并行，
blockMesh
createCellZones
createFaceZones
dsmcInitialise
decomposePar
然后开并行的时候就计算不下去了，请问你遇到过这种情况吗

我看的一些OpenFOAM的文章，他们使用的ILES好像是湍流模型直接使用laminar，相较于q-DNS，网格较粗。