@东岳 谢谢东岳老师提供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