pairCollision中，A和B分别指网格还是粒子？

zhe

    // 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，但是看上面的意思是指代两个网格吗？希望有人能帮忙解释一下，谢谢。

bestucan

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

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

zhe

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

bestucan

@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怎么联系起来的。

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

zhe

@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网格的指示相同而以）。这样理解不知道对不对？

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

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

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