@霜染丹枫 入口用的是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,
);

@Ant 请问是不是考虑了颗粒-颗粒碰撞之后，就不会有这个颗粒堆积使用的问题，因为颗粒颗粒之间的间距考虑进去了？

@刘雄国 非常感谢你的回复，请问你说的这个字符码，是在particle.C里修改之后，再用creatField输出的意思嘛？

因为用的是并行计算，颗粒的ID被输出到origProcId文件（计算颗粒的processor）和origId文件（颗粒在当前processor中的ID）。 比如我有4000个颗粒用4个processors计算，我的origProcId值是0-3，origId可能是0-999。这种情况下需要两个文件中的origProcId和origId值完全对应，才能找到指定的颗粒，对于大量颗粒时间步多的情况来说，效率有点低。

所以想请问一下各位老师们，怎么让openfoam运行的时候就输出一个全局的颗粒ID？按上述例子能得到一个颗粒的ID文件是0-3999，而不是4个0-999文件。

我注意到particle.H文件中有一段是：

        //- Cumulative particle counter - used to provode unique ID
        static label particleCount_;

        //- Fraction of distance to tet centre to move a particle to
        // 'rescue' it from a tracking problem
        static const scalar trackingCorrectionTol;

        //- Fraction of the cell volume to use in determining tolerance values
        //  for the denominator and numerator of lambda
        static const scalar lambdaDistanceToleranceCoeff;

        //- Minimum stepFraction tolerance
        static const scalar minStepFractionTol;

这个particleCount是不是我要的unique particle ID？请问怎么修改能使其输出？

先谢谢各位老师们了

RT，在颗粒流求解器中，好像没有现成的utility能统计某时间段中颗粒碰撞的次数吧？请问如果想实现这一功能，有什么思路嘛？

@东岳 对，竖直管道，重力和就是流向方向，出入口是周期边界，然后用添加了一个Ubar控制总体速度，初始速度场是从充分发展湍流的管流映射过来的。颗粒-气体体积分数是1o^-4量级，St在5左右。这个截面也是充分发展之后的结果。

@星星星星晴 哈哈这个不夸张，确实就是这个理，谢谢你的建议。

@东岳 确实求出的结果就是这样，我也纳闷。然后我试了很多组模拟，发现在总体速度稍低（Ub≈3）的颗粒流中，颗粒速度总是会大于气体速度，然而在总体速度稍高（Ub≈12）的颗粒流中，颗粒速度只有在近壁区域会大于气体速度，其他区域都是小于气体速度。如下两图所示，其中纵坐标是颗粒或者气体的绝对速度，横坐标是距离中心线的距离，r/R=1表示管壁。我的模拟都是重颗粒（颗粒空气密度比为1000）

我一开始以为是网格的问题，然后根据速度的不同，按相同比例细化之后，发现仍然是这个趋势。不知道是什么原因导致的？会是DPMFoam中对颗粒速度处理的问题嘛？

我查找过相关文献，但是确实在稍高速度的条件下的研究很少，所以也无法确定是模拟的问题，还是颗粒流的真实情况就是这样？但是我的理解是颗粒收重力的作用，任何情况下都会比气体快的，这个理解没错吧？