我在瞬态投影法的算法里面，没有看到RC插值的操作。

在interIsoFoam基础上加入了拉格朗日粒子，编译成功后在网格不变的情况下运行没有问题，但是加入动网格后，计算在进行一段时间后突然终止，没有报错，就是停在了Solving 3-D cloud kinematicCloud这里：

Courant Number mean: 4.05392e-05 max: 0.899987 Interface Courant Number mean: 3.47277e-06 max: 0.47951 deltaT = 4.6765e-06 Time = 0.102143 PIMPLE: iteration 1 Selected 4 cells for refinement out of 208887. Refined from 208887 to 208915 cells. Selected 0 split points out of a possible 6252. GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.35461e-06, No Iterations 25 GAMGPCG: Solving for pcorr, Initial residual = 0.0776576, Final residual = 7.22028e-06, No Iterations 15 time step continuity errors : sum local = 4.77276e-15, global = 1.316e-17, cumulative = 8.05896e-11 Evolving kinematicCloud Solving 3-D cloud kinematicCloud

重新跑了好多次，每次都是停在这里，有大神知道怎么解决吗？

@李东岳 在 OpenFOAM中移动接触线问题（接触角）的边界条件设置 中说：

我去，竟然来了！！

别人提问的话，我这边能收到邮件呢。工作太忙，顾不上逛帖子了。

另外的话呢，对于像OpenFOAM和fluent这类使用CSF模型来计算气液两相流的话，低速问题时虚拟流动（spurious current）影响较大，做学术研究的话精度恐怕不够。建议用Basilisk去做，精度很高。但是工程问题的话，商软或者OpenFOAM是可以的。其他看到的吧友可以借鉴一下。
我目前已经不做单纯的界面流动问题了，研究方向稍微有了变化。

chatgpt信息量太大了..

看了一篇论文，说是在openfoam上搞了个dem扩展，原文是这个:
OpenHFDIB-DEM.pdf

然后对应的代码库在github是这个：
https://github.com/techMathGroup/openHFDIB-DEM

表面上看起来好像挺有用的，但是毕竟没用过也不大知道这东西的弊端，有没有大佬尝试过这东西，评价一下给小弟开开眼界。

@dabiu 在 关于icoFOAM中通量phi的更新 中说：

icoFOAM中更新通量phi采用的如下语句（注：我理解这个语句）

if (piso.finalNonOrthogonalIter()) { phi = phiHbyA - pEqn.flux(); }

为什么不在下面这条语句的后面，直接用phi = fvc::flux(U)来更新通量，这样更新跟上面那条语句有什么区别吗？用phi=fvc::flux(U)会有什么问题吗？

U = HbyA - rAU*fvc::grad(p); U.correctBoundaryConditions();

初学icoFOAM，求各位高手解答，OpenFOAM用的版本为OpenFOAM 7.

我觉得这里的phi是经过Rhie-Chow插值修正过的面通量。

@李东岳 对哦，谢谢您

看起来是个非常奇怪的现象

你这是个二维算例吗？我记得openfoam二维算例用网格自适应会有问题的啊？

@heike256 帅炸了！

理想气体负压会报错，需要使用pMin限制。求解器里在更新热物性之前对p施加pMin限制（比如mixture.correct()前），之后再把压力恢复原样应该就可以了。不过我认为合理的计算不应该出现负的绝对压力。

@HatoriChise 多谢老铁支持

我之前使用相同网格，
fsiProperties:

couplingReuse 1; couplingScheme Aitken;

时非常好用。
不知道为什么改成了:
fsiProperties:

couplingReuse 1; couplingScheme IQN-ILS;

就出现了以上错误。

是的，类似这个 http://dyfluid.com/openfoam.html#id6 放大尺度就不行了

@xpqiu 太强了:140:

这个bug看起来是湍流模型的，但是基于你的反馈，缺出现在rhoCentralFoam+动网格的情况。看起来像是个bug。

https://cfd-china.com/topic/8018 你看看这个

问题解决！在使用旋转周期性交界面的时候把 rotationAngle 45去掉

你最好把圆孔当做inlet，其他是壁面，这最简单。