@李东岳 您好东岳老师，我后续发现，这些流失的气体flow rate其实是从壁面流走的，也就是壁面的速度并不为0，尽管我设置的壁面的边界条件是noSlip。
如下图所示，我的圆管壁面遍布着非0的速度值，是因为颗粒进入近壁的网格并改变该网格的气体速度导致的（双向耦合模拟，颗粒比近壁的网格尺寸小）。
所以我还想请问您一下关于这个壁面速度不为0，导致有气体flow rate流失的问题：
@李东岳 flowRatePatch(name=inlet), flowRatePatch(name=outlet). 是通过这个对得到的时间步进行后处理得到的
最近用DPMFoam求解一个不可压缩的颗粒-气相流的算例中发现，因为有颗粒的存在，入口和出口的气体flow rate不守恒。然后我试了一下tutorials里的Goldschmidt算例，发现入口和出口的flow rate也不守恒，出口的flow rate只有入口的60%左右。所以想请问一下大家：
Goldschmidt算例中，因为颗粒占的体积分数是0.2比较大，所以颗粒会明显地减小气体的速度说得过去。但是从流量守恒的角度，应该怎么解释这个问题？算例中，入口比出口多进入的气体流量去哪里了？
我的算例中，颗粒的体积分数只有0.001左右，很小。但是所得到的在出口处的气体flow rate大概只有入口的10%，感觉解释不通。。请问这能是什么原因导致的？
@李东岳 在 能不能设置两组颗粒（一组one-way，一组four-way）同时模拟？ 中说：
我目前有个想法是添加两个不同的拉格朗日cloud，看看能否分开处理，但我现在在出差手头没有openfoam还不确定这么做是否straightforward。
欧拉场这面可以算一下速度的物质导数就是加速度$D \bfU/D t$
现在才反应过来您的这个方法才是最straightforward的。。
还想继续请教东岳老师，在DPMFoam中，从欧拉场计算流体相的加速度，是简单地为
DUcDt = fvc::ddt(Uc)，
还是
DUcDt = fvc::DDt(phic, Uc)？
另外我在pressureGradientForce的C文件中，还看到
fvc::ddt(Uc) + (Uc & fvc::grad(Uc))
请问哪个是颗粒-液相流动中，液相的正确加速度？
非常感谢楼上两位的耐心回复，所以是不是可以认为：
（1）deterministic就是接近于真实的碰撞情况（虽然只考虑一个pair两个球的碰撞，而真实中可能还有三个球四个球的碰撞），只要初始条件一样，不论模拟多少次，出来的结果都是一样的；
（2）而stochastic就是随机性的，哪怕初始条件一样，每次模拟出来的结果都是不一样的？
另外这个事驱和时驱的概念还是不太理解，拿事驱碰撞来说，比如有三个独立的碰撞，分别发生在1.4s，2.0s 和 3.3s，如果我是用固定时间步dt=1.0s进行模拟，那岂不是只能捕捉到第二个发生在2.0s的碰撞，其他两个时间点的碰撞就被错过了？
我原以为硬球和软球的区别就是，硬球只考虑碰撞的速度变化不考虑形变等参数，而软球是考虑的。现在才知道他们更主要的区别是，硬球是事驱，而软球是时驱。然后我的问题是：
在通常的模拟中，都是通过每个时间步来计算颗粒位置和速度，这应该就是所谓的时驱计算吧？这样的话，事驱的模型是怎么应用到这样的时驱模拟中呢？
想问问这句话对不对：颗粒的deterministic碰撞是事驱，stochastic碰撞是时驱？
最后想问问openfoam中的pairCollision, 是属于deterministic碰撞的一种吗？
最近被这几个概念弄得非常迷糊，也找不到很权威的能看懂的概念，所以想问一问论坛里的前辈们，可能我的理解完全不对，请纠正
@lys 有一个叫PatchCollisionDensity的cloudFunctions，你可以试试
最近对两个概念非常混淆，即摩擦系数和切线回弹系数，这两个系数是否就是同一个概念？
比如在以下的颗粒pairCollision参数设置中，法向碰撞的回弹系数根据alpha设定。但是切向的回弹系数，是不是就是这个mu?
pairSpringSliderDashpotCoeffs
{
useEquivalentSize no;
alpha 0.01;
b 1.5;
mu 0.1;
cohesionEnergyDensity 0;
collisionResolutionSteps 1;
};
@东岳 谢谢东岳老师提供idea。
我尝试了比如添加一个kinematicCloud或者添加一组constantProperties，但是好像不修改C或者H文件的话，在计算时只会读取第一组的参数。
欧拉-欧拉模拟的话我没有什么经验，那请问我这样相当于是三相，一个是流相，二是拉格朗日颗粒相，三是欧拉颗粒相，这样是可以模拟的嘛？
@东岳 其实是想通过拉格朗日的tracer颗粒来跟踪流相的颗粒，用以研究流相的加速度。欧拉的话好像不太方便研究加速度吧
确实两组粒径差别太大的颗粒算碰撞太缓慢了，请问有没有办法能让碰撞只针对某一组颗粒？
如题，想在four-way coupling的气固颗粒流中，添加一组one-way coupling的tracer颗粒。添加的tracer颗粒的密度等于气体密度，直径接近于0。
我本想可以通过多个injectionModels来添加一组直径非常小的颗粒。但是这样下来，tracer的直径和气固流中的颗粒直径相差太大，而我全局是four-way的需要计算颗粒碰撞，导致计算非常非常的缓慢。
所以还是想问问，有没有办法在模拟中同时使用两套kinematicCloudProperties？一组是four-way coupling，一组是one-way coupling。我是用的是DPMFoam求解器。
@fubianhanshu 并没有在openfoam中解决，不过可以在后处理中根据ID和processor来找到相应的颗粒，我现在就是这样做的
抱歉问的问题可能太泛了，具体举个例子来说。比如说在直径D=0.01m的空气管流中，速度Uf=15m/s，雷诺数Re=10000。颗粒的斯托克斯数St=10，直径是大概dp=50×10^(-6)m，密度为ρp=1000kg/m^3，颗粒-空气的质量比是Φ=0.1，双向耦合，只考虑颗粒曳力。
在这种情况下，我查看OpenFOAM计算得到的kinamaticCloude:UTrans文件中，动量转移都是在10^(-14)量级（单位是kg×m/s），请问这合理吗，为什么只有这么小。。？
据我所知，这个UTrans文件表征的就是颗粒→流体相的动量转移，也就是颗粒相通过源项Sp，施加在流体相的曳力作用。不知道这个理解对不对？
加了引号好像也不起作用呢，都是显示读取了这个场文件，但是就是没有文件输出出来。
Reading set description:
cloudPlane001
Time = 0.082
Reading fields:
volVectorField: U.air
volScalarField::Internal: kinematicCloud:UCoeff
volVectorField::Internal: kinematicCloud:UTrans
Executing functionObjects
End
结果中的这两个场文件，想用Probe或者sampleDict进行提取输出，但是都没有任何内容，请问这两个场有什么特殊的处理方式嘛？
我试过在probe或者sampleDict中设置如下代码来提取一些点的数据，但是都只能输出U.air和p的数据，kinematicCloud:UCoeff和kinematicCloud:UTrans并没有任何内容，也没有报错。
fields
(
U.air p kinematicCloud:UCoeff kinematicCloud:UTrans
);
另外还想问一问在paraview中也找不到这两个场文件，想问问怎么能可视化这两个场文件？
非常感谢大佬们的对话，受益匪浅啊
@cxzhmg 并没有琢磨明白
creatFields能添加的应该都是流体相的场，离散的颗粒相应该实现不了
@霜染丹枫 入口用的是turbulentInlet，出口是inletOutlet
@qingdong_wang @东岳 谢谢一楼和东岳老师的回复。还想请问东岳老师，这个曳力项，应该写在哪个文件里输出？
我发现求解器(比如DPMFoam)里面的creatFields.H文件，输出的都是流体相的场文件，比如Uc, p, phic之类的，并没有颗粒相的场文件。我也查找了一楼所说的src/lagrangian/basic目录下的文件，也没发现控制输出颗粒场的文件。
所以还想麻烦您进一步指明，这个曳力项应该写在哪里？
RT, 好像现有版本的OP没有相应的utility可以实现这一功能吧。是不是只能在求解器的creatFields.H文件里里添加了？
老实说本人的C++功底太差了。。所以想麻烦各位帮我举个例子怎么修改？比如我使用的是ErgunWenYuDrag，想在每个输出的时间步中添加这个颗粒的ErgunWenYuDragForce数据文件，请问怎么实现？
我试着添加如下，但是报错一大堆。
volScalarField ErgunWenYuDrag
(
IOobject
(
IOobject::groupName("ErgunWenYuDragForce", continuousPhaseName),
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
);