CFD中文网

extend版本应该是类似“A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries”这篇文献的实现方法 ，力源项应该是“An immersed boundary method with direct forcing for the simulation of particulate flows”文章的思想，虽然本质上差不多，但是实现边界条件的方式还是不一样的

@OItoCFD extend版本的IBM是采用插值重构的方式植入的，我们本来想植入直接在动量方程中施加体积力源项的形式的IBM方法，不过最后还是使用了extend版本带的IBM方法

@东岳 谢谢

感谢，最近正好在用它画网格！还想问一下李老师平时用什么软件画网格！

还有这种情况会造成一个owner（或者neighbour） 对应两个neighbour（或者owner）那他们共享的通量有特殊处理么？

谢谢东岳老师及时回复，那关于第二点呢？

最近要使用refineMesh，使用后用paraview打开发现粗细网格交接处有四面体网格过渡，如下图

在网上查找发现和paraview中properties中选项有关，若选择默认的interpolated Volume Fields则会出现上图情况；
若修改该选项为use vtk polyhedorn过渡的四面体消失，如下图

我想知道这是paraview显示的问题，还是网格本身就存在这个过渡四面体？
还有这种情况会造成一个owner（或者neighbour） 对应两个neighbour（或者owner）那他们共享的通量有特殊处理么？

非常感谢李老师，cellLimited Gauss linear 1 格式比Gauss linear好很多，运行程序之前运行renumberMesh -overwrite改善网格质量也有作用

最近发现在openfoam计算压力梯度时，用blockMesh生成的均匀网格，给定压力场为均匀初始场，当压力初始场为零时，全场压力梯度为零，但是当压力初始场为101325时，梯度在x,y（2D）方向上绝对值大小为1e-6(由于网格质量引起的数值误差？)，虽然这个值量级比较小，但是目前计算的系统网格长度量级为1e-5，会对最终的结果产生误差，我分别在openfoam2.2.x以及openfoam6.0上测试都存在该问题，压力梯度计算方法为Gauss linear，不知道有没有人碰到类似的情况，或者有没有推荐的绘制网格质量较高的软件，我试试在不同网格下计算看看
下图为压力初始值为101325的压力梯度分布
（另附测试算例链接https://pan.baidu.com/s/1uFbx_xqbO-5QAbcmTSr-fA）

@东岳 老师，DEM这块不是我编写的，工作主要是在CFD部分，据我了解DEM是计算单个颗粒的曳力，不用指定气泡直径大小吧；那我再看一下CFD部分曳力计算吧，谢谢您

@东岳 确实有聚并破碎

@东岳 CFD-DEM不过DEM部分是自己写的软件！

图中是鼓泡床横截面孔隙率图，我想识别出黑色圈中白色的两个区域气泡以及气泡直径，在网上查了一些资料是用matlab处理，但是自己跟着教程发现不能识别，不知道有没有人有相关资料推荐，万分感谢！

@东岳 在 关于分步法和PISO算法 中说：

Peric

好的 谢谢李老师，确实没有认真看过computational methods for fluid dynamic这本书 我再仔细看看

想请教一下有没有对分步法fractional-step method 了解比较多的，和PISO算法区别，有没有推荐看的资料和文献，还有PISO算法的速度预测步是不是可有可无的一步？

好的。谢谢李老师，我设置成zeroGradient也能算，目前先这样，我再查阅一下相关资料。

想设置一个2D圆柱绕流算例

文献中速度出口设置为对流边界条件，想请教一下在openfoam中如何设置这种类型边界条件！

@cyw 可以在场变量更新的地方添加correctBoundaryConditions();如U.correctBoundaryConditions();目的是可以让进程间边界面上的值更新；具体可以参考 https://www.cfd-online.com/Forums/openfoam-programming-development/88050-how-handel-fields-parallel-computation.html 。希望有用

自己定义了一个类，如图，通过该类求解出来的数据在进程间不连续，我应该怎么写进程间的通信，有没有大神碰见过和我一样的问题，或者openfoam内有没有可以参考的MPI案例，求解答

再更贴一下，openfoam中矢量数组可以新建为矢量列表形式，初始化方法为List<vector> A(n,vector::zero),希望对和我用同样疑问的人有帮助

好的 谢谢李老师啦 我自己再看看