@李东岳 在 Simple算法中通量phi的计算 中说：

也就是说在边界处phi是不正确的。不知道是我没有找到phi边界修正，

边界处的phi不需要计算，是用户给的。加入速度1m/s，乘以网格边界面积，就是phi啊，并且不变。

如果是pressureInletVelocity 边界呢？ 代码里可以看出pressureInletVelocity 需要用边界上的phi来更新出新的边界速度。如果是这种边界边界上的phi又应该如何确认呢？

@李东岳 李老师，您好，我按照您的建议进行了相关的测试。结果发现当更换为GAMG时计算比PCG慢非常多，如果我在计算一开始取残差与相对残差值较大时，完全不收敛的情况基本上到500步左右就改善了。但是我又发现一个特别奇怪的现象：
首先通过pimpleFoam (nOuterCorrectors=1) 计算到一定步时计算结果收敛很好，残差基本下降至1e-5左右；但是当我在运行时更改nOuterCorrectors=2时，初始残差立即上升至0.4左右，计算异常缓慢，并且继续运行算例残差也下降不了。请问这个可能是什么原因导致的呢？

替代文字

你说上面这个图。这个算力对应的不是这个图。所以不一样。不过如果你增加近气速度的话，应该也是可以的。你加大10倍试试。

您好，请问这个问题您是否有好的解决办法了吗？
强制设为0可能不太对，反映出来的可能不是压实效果。。

@李东岳 东岳老师，这个主要是因为我的是粘弹性流体，在动量方程里边加入了弹性相，不只有雷诺数，还有魏森贝格数，可能是这两个无量纲数导致的:135:

4年过去了，当时发完这个贴就出去玩了。:jingya:

@ruibin_tongji 你看你的log是不是有些迭代次数太多了

@是嘿嘿牛啊 密度是被动地，是根据流体种类、温度和压力变化的，你可以设置上下两层温差或不同的流体种类

@李东岳 噢噢，谢谢东岳老师。

@shaneheee 那就不用管了。

有时候加上一行啥用没有的代码，算例就能用。那就不用管了，放在这。

@xpqiu 太感谢您了 您教的好细， 非常感谢 ，对我帮助非常大 ，谢谢您:140:

各位老师好！
报错信息：
Screenshot from 2021-11-20 14-31-59.png

solidParticleCloud.C的代码：

solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); Cloud<solidParticle>::move(*this, td, mesh_.time().deltaTValue()); if(mesh_.time().value()> td.cloud().tInjStart_ && // mesh_.time().value()< td.cloud().tInjEnd_) {this->inject(td);} void Foam::solidParticleCloud::inject(solidParticle::trackingData &td) //. { //. label cellI = -1; //. label tetFaceI = -1; label tetPtI = -1; mesh_.findCellFacePt(td.cloud().posB_, cellI, tetFaceI, tetPtI); if(cellI > 0){ solidParticle* ptr1 = new solidParticle(mesh_, td.cloud().posB_, cellI, tetFaceI, tetPtI,td.cloud().dB_, td.cloud().UB_); Cloud<solidParticle>::addParticle(ptr1); } }

@bestucan 嗯嗯，谢谢老师，我试试:146:

@bestucan 你们都溜的很啊！

@是嘿嘿牛啊 按照我的理解，这样分的size不应该都是10吗？为什么size会变来变去呢？求大佬指点一下，拜托啦！

好像我问的有问题，$A=0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})$这个东西看起来是个矢量。但是面积不是矢量啊。我回去对对OpenFOAM的代码。

我看了下，OpenFOAM里面mesh.Sf()对应的是这个：$0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})$，mesh.magSf()对应的是这个：$\left|0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})\right|$，后面这个应该就是面积了。

@evensun $\left|0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})\right|$ 这个面积可以理解不

大家好，我想问下在rhoCentralFoam中想植入一种新的重构格式，文献中以标量守恒律作为控制方程，对其的未知量进行了界面左右值的重构，这个应该在rhoCentralFoam中哪一块进行植入呢，我不太清楚，望解答一下 :yes: ，谢谢！

@冠竹 有可能是入口速度太小了

@tens 好的！非常感谢您

@李东岳 可以把日常的流体都看成稀薄流体呀。从各种benchmark来看的话，精度还不错，处于同一个水平。
现在主要的问题是，显式LBM并行效率特别的高，然后我就寻思着有限体积也可以显式呀，谁怕谁。结果性能卡在了压力泊松方程的求解上:136: