颗粒的全局ID输出问题?



  • 因为用的是并行计算,颗粒的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?请问怎么修改能使其输出?

    先谢谢各位老师们了:xinxin:



  • //- Cumulative particle counter - used to provode unique ID
    static label particleCount_;
    应该不是

    你可以修改基础particle类,让它生产粒子源时附带一个"字符码",particle在multi processors区域传输的时候把这个字符码也带上。就可以区别开来了。



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



  • 楼主你这个问题解决了吗?



  • @fubianhanshu 并没有在openfoam中解决,不过可以在后处理中根据ID和processor来找到相应的颗粒,我现在就是这样做的



  • @zhangxc0223 我不清楚你怎么做的,不过我刚发现可以在paraview中可以显示粒子ID,随着时间步推移,好像每个粒子的ID号是不变的,这样就可以通过后处理输出某个固定ID的各个时间步的信息了,可以交流下


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2