pairCollision中,A和B分别指网格还是粒子?



  •     // Loop over all Parcels in cell A (a)
        forAll(cellOccupancy[realCelli], a)
        {
            pA_ptr = cellOccupancy[realCelli][a];
    
            forAll(dil[realCelli], interactingCells)
            {
                List<typename CloudType::parcelType*> cellBParcels =
                    cellOccupancy[dil[realCelli][interactingCells]];
    
                // Loop over all Parcels in cell B (b)
                forAll(cellBParcels, b)
                {
                    pB_ptr = cellBParcels[b];
    
                    evaluatePair(*pA_ptr, *pB_ptr);
                }
            }
    
            // Loop over the other Parcels in cell A (aO)
            forAll(cellOccupancy[realCelli], aO)
            {
                pB_ptr = cellOccupancy[realCelli][aO];
    
                // Do not double-evaluate, compare pointers, arbitrary
                // order
                if (pB_ptr > pA_ptr)
                {
                    evaluatePair(*pA_ptr, *pB_ptr);
                }
            }
        }
    

    上述一段code是来自paircollision的,原本一直一位A 和B 是分别表示两个parcel,但是看上面的意思是指代两个网格吗?希望有人能帮忙解释一下,谢谢。


  • 版主 副教授

    这段一看 A B 就是网格。但是把 A B 理解成 parcel 相关也不算没有依据,parcel 好像就是从网格上划分出来的。比如,进口面有10个网格,那么每个网格喷出来1个 parcel。一次喷出来10个parcel。

    如果把 parcel 理解成单位时间内单位空间上的粒子物理量。那理解成网格还是parcel的区别不就是欧拉和拉格朗日的区别。
    我猜的:mianmo:



  • @bestucan 谢谢回复。但是我还是没有理解您的意思。A和B从code上来看,确实是代表的网格。那么对于parcel A 和 parcel B之间的作用,怎么就能当作从网格里提取出来的呢?那是不是说明一个网格就只能代表一个parcel了?感觉这不是他们想表达的意思。希望您可以帮忙再想想。我也会再研究一下,希望还能再多沟通,谢谢🙏


  • 版主 副教授

    @zhe 原来你是疑惑怎么从网格里提取 parcel 啊。

    https://www.openfoam.com/documentation/guides/latest/api/PairCollision_8C_source.html
    你放的代码是这个89行吧,84行创建了cellOccupancy,

    那个创建容积率的函数在多个模型下有同名的,我猜这里这个是
    https://www.openfoam.com/documentation/guides/latest/api/KinematicCloudI_8H_source.html
    这个的371行,

    这个函数就是判断是否存在,否则创建。就是这个函数buildCellOccupancy(),在
    https://www.openfoam.com/documentation/guides/latest/api/KinematicCloud_8C_source.html
    这里140行

    你给的代码里 pA_ptr, pB_ptr 都是通过 cellOccupancy 的返回值得到 parcel 的指针,就是 那个371行,如果存在就返回
    return *cellOccupancyPtr_;

    而这个 cellOccupancyPtr_ 也有多个定义,在这里应该是这个
    https://www.openfoam.com/documentation/guides/latest/api/KinematicCloud_8H_source.html
    173行。

    173行的注释大概解释了cell 和parcel怎么联系起来的。

    好了
    模板继承用的太多,我不太熟悉这部分,同名函数定义太多,剩下的我也没招了:136:



  • @bestucan 抱歉回复的晚了,已经很感谢您的解答了,帮我捋清了很多了。我按着您的这个思路走了一遍,大概明白了我之前给出的的确是A和B都是网格名称。而同时,在分别中也代表着其中包含的粒子a和b,也就是您说的也可以理解为A 和B 也是粒子。

    我的理解是A 和B 是有相关的两个cells,而最后的

    evaluatePair(*pA_ptr, *pB_ptr);

    是建立两个网格之间(包括网格里所有粒子)的一种配对机制,不知道对不对?

            // Loop over the other Parcels in cell A (aO)
            forAll(cellOccupancy[realCelli], aO)
            {
                pB_ptr = cellOccupancy[realCelli][aO];
    
                // Do not double-evaluate, compare pointers, arbitrary
                // order
                if (pB_ptr > pA_ptr)
                {
                    evaluatePair(*pA_ptr, *pB_ptr);
                }
            }
    

    在之后的loop里,对于相同cell里的粒子相对配对,也就是说把一个粒子单拎出来a,然后其它的所有粒子相当于a0(而这个其实与之前B网格的指示相同而以)。这样理解不知道对不对?

    最后一点儿避免两次被征用,所以有了一个对比。

    知道您肯定也花时间翻了一遍了,实在是感谢!

    如果有其他朋友一起讨论,感激不尽。


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1