我想用MPPICInterFoam做泥沙冲淤，水槽底部需要铺5厘米左右的沙子，我就先在静水中加入颗粒，让颗粒在重力下自然沉降，但是算的过程中会经常出现算到某一步时间步库朗数和颗粒动量都突然增大，然后计算就卡住不再进行了，有时候关掉从上一个存储文件接着算还能再算下去，但是时不时就出现这种卡住的情况是什么原因呢？该怎么解决呢？

@李东岳 我在水槽底部5公分范围内加了颗粒，然后开始计算，水槽内是静水没有流速。开始颗粒正常沉降，最大相分数也随之增大，但很快就超过0.63了，之后颗粒的动量值还在一直增大，才算了0.2秒的过程就崩了 ，怎么都已经达到packing limit了颗粒动量和最大相分数还在持续增大，而且大得离谱，而且我流场是静态的啊，我不清楚是哪儿出的问题

@李东岳 李老师，您好！请问该怎么判断它是否到达packing limit呢？

我在距水槽底部5公分范围内加入了41万个初始颗粒想进行冲淤模拟，结果模拟了一秒钟这些颗粒就都沉底铺开了，河床高度一下就没了，有人知道是怎么回事儿吗？

看官网介绍说这个driftFluxFoam也就是老版本的settlingFoam是用来计算离散相在连续相中的沉积的，想尝试用它计算泥沙冲淤。我用setFields在计算域中加了底沙，但是算起来发现没一会儿底沙就都冲没啦 想知道这个求解器适不适合做泥沙冲淤类的问题，或者还需要做哪些改动。

请问有人知道沿着弯道分布的一排圆柱的拖曳力该怎么求吗？

如题，有没有人在MPPICFoam中添加过多孔介质源项的？怎样实现多孔介质阻挡颗粒啊？

@xpqiu 您好，我现在发现了一个新的问题。我在论坛里看到在paraview里求解拖曳力的方法，可以提取多孔介质区域的表面，然后生成表面法向量Normals，用计算器功能计算p*Normals的分量，再用Integrate Variables 积分就可以得到多孔介质区域的拖曳力啦。但是不知道为什么，我用中国方法计算得到的拖曳力和用forces function后处理得到的拖曳力差很多，甚至数量级都不一样。想请教一下您这是怎么回事。谢谢！

@xpqiu 好的，非常感谢！

@xpqiu 您好，想再请教您一下。这个函数求得结果在postProcess文件夹里的是整个计算域内多孔介质源项的porous force的合力吗？我的计算域是弯道的，多孔介质区域是好多根并排的圆柱体，那这个合力作用位置和方向该怎么确定呢？