颗粒的重力也会导致颗粒的运动，颗粒有运动就会导致流体的速度改变，就会发生动量耦合。

@Tens 谢谢老师的指点！我学习一下这篇文献

@李东岳 对了老师，请问一下reactingFoam是不可压缩求解器（论坛上有位大佬好像提过）吗？但是我看它源代码里面能量、动力、连续性方程里面都包含密度。我算的是一个扩散燃烧的算例，有两个速度入口（甲烷空气预混+仅空气），但是如果我把其中的空气入口设为壁面，这样变成一个预混燃烧问题，残差就能降到10^(-3)了，感觉很疑惑

用ICEM画了一个粗网格一个细网格，两套网格质量都在0.8以上，第一层网格高度一样，检查网格没有问题，但就是粗网格向细网格映射的时候出现图中错误。请问这是什么原因呢!粗细网格映射不成功.png

您好，已经2024年了，我在使用DPMFoam时候也遇到了同样问题，请问您已经解决了么。我目前用的of300，版本太老了，而且昨天使用时候都没有这个报错，今天就改了一下粒子数和初速度就有这个报错。是不是of新版本的拉格朗日的求解器稍微好用一点

就是这个问题，参照这里的已解决~
https://www.cfd-online.com/Forums/openfoam-installation/229950-icofoam-not-installed.html

openfoam3太老了，之前拉格朗日粒子总有这个bug，新版本好了一些。

@李东岳 我先将网格细化做一下测试，排除或者确认一下网格的原因。

半圆部分做两个圆o切，然后把里面部分删掉就好了？？可以试一试，多分几个块

@李东岳 好的，多谢李老师，那我尝试联系一下原作者吧

可以考虑合并后再将前面的流场映射到新网格上去，就是那个mapFields。

各位前辈好，我在计算波浪反射透射的时候，尝试写了一个用两个浪高仪计算的程序，但是计算出来的结果反射波都大于入射波，请问哪位前辈写过类似的程序，能跟我分享一下经验吗

感觉就是不同的研究者对这个的规定不一样，论文里面推导过程正确就行。
当然如果按照张量的知识，向量并积计算的话，我觉得应该是
$$
\nabla {\bf{U}} = \frac{{\partial {U_j}}}{{\partial {x_i}}}
$$

感谢分享

@李东岳 在 各位老师，请问我是算出棋盘型压强场了吗？ 中说：

知识盲区了老铁。运动大体上没问题就问题不大。不过看起来有一些小bug没处理好，也不能说肯定是压力期盼分布

运动大体上差不多，但是肯定不完全对:xinlei:

@李东岳 太好了！学习中！:142: :142: :142: :ok3:

好的，谢谢李老师

@zhaoym 老铁，后来这个做的怎么样了，有进展么

interPlicFoam求解器已经并入扩展包geometricVofExt，仓库地址：https://github.com/daidezhi/geometricVofExt

并行计算时候，每个进程掌握的网格就固定那些。
举例b点只存在于processor x1，不在processor x2当中，processor x2自然无法获取b点的信息。
要想processor x2获取b点信息，可以用reduce函数，也就是进程之间通讯一下。
首先初始化intU为零向量，然后赋值，再reduce。

point b (1,1,1); label cellID = mesh.findCell(b); vector intU(0,0,0); intU = U[cellID]; reduce(U,sumOp<vector>());

具体代码你再看一下，我这个是思路，代码可能哪里不对。