@李东岳 多谢多谢 嘿嘿嘿

压力泊松方程的离散问题

of中构建的压力泊松方程似乎是这样的形式：
$\nabla\cdot(\frac{1}{a_P}\nabla p)=\nabla\cdot(HbyA)$
(icoFoam解析中的式27、33）

按我的理解，对上面的方程似乎应该写成下面这种离散形式？

fvm::laplacian(rAU, p) == fvc::div(HbyA)

为什么实际上icoFoam使用了phiHbyA作为中间变量呢？

phiHybA的问题：

phiHybA是一个surfaceScalarField 对象，即面上标量场；随后求解器执行了对这个面上标量的散度：

fvc::div(phiHbyA)

为什么可以对一个面上标量计算散度呢？是否div函数在对surfaceScalarField作用的时候采用了不同的方式？我查了div的代码，确实看到div对surfaceScalarField似乎有不同的处理，但是我没看懂。请问这是否有关系呢？

期待各位大佬指点

@linhan-ge 在 多相流VOF 中说：

VOF-DEM Study of Solid Distribution Characteristics in Slurry Taylor Flow-Based Multiphase Microreactors
强啊前辈，我这进展还缓慢着呢
刚看了您“VOF-DEM Study of Solid Distribution Characteristics in Slurry Taylor Flow-Based Multiphase Microreactors”的模型部分，请问您是在interFoam上加了Laplacian filter就可以抑制假速度、算出液膜了吗？因为我现在就做气液流动，速度也跟您差不多（1cm/s~10cm/s），所以想请教一下这个方案可行性。我直接用interFoam算的液膜厚度非常厚，不知道问题在哪里。
另外我看到您是设置接触角90°，那最开始设置气相区域的时候需要留出一个液膜厚度吗？还是直接把整个管径宽度都设置成气体都可以？我比较担心如果把整段管都设置成气体会让气体“粘”在壁面上。
谢谢！

话说什么叫较大体积力？考虑重力的气液两相流动算吗？反正我看到文献力一般处理这种都是PRESTO！

请问李老师，计算域内网格尺寸不一样的话，库朗数是按照最小网格计算的吗？也就是说，计算的速度是被最小单元控制的？

@东岳 感觉还是自适应网格的问题，我去掉了dynamicMeshDict，就能算了

@东岳 李老师好。我改了pcg，还是发散了。log最后是这样的：

我改的fvSolution文件是这样的：

运行了of-6自带算例 damBreakWithObstacle ，在自己的小破笔记本运行可以收敛，但是在服务器上运行就会很快发散。请问有人遇到过类似问题么？
以下是发散时候的报错：

谢谢各位大佬前辈

@linhan-ge 哇太棒了 到时候搜您的名字ge linhan 可以搜到么？

使用的软件不限，主要做算法还是实现也不限。课题组自己一个人搞这东西真的头大。要是有了解的拜托拉一下，多谢啦！