OpenFOAM编程findCell的诡异问题，对同一坐标寻找cell结果不一致

[OItoCFD]

@李东岳 在 OpenFOAM编程findCell的诡异问题，对同一坐标寻找cell结果不一致 中说：

另外你没用过primitiveMesh的findCell()试试么，这种

point a(1,2,3);
mesh.findCell(a);

我几年前记得在某个SCI里面看到过说meshSearch好像效率低，不过具体细节我忘了

老师，mesh search和mesh的findCell我都是用了的，一样的结果，都是循环外单独找那两个点是对的，循环内就不对，是重复的CellID

[马乔]

@oitocfd 可能是精度问题？手动写入的有截断误差？刚好这个点位于face上呢？还有不管是primitiveMesh和meshSearch的findcell效率都挺低的，都要去遍历所有的cell，不知道写什么需求要这么写？

[OItoCFD]

@马乔 您好，这么讲的话，也许是可能是截断误差影响。然后我之所以用findCell，是因为我用到一个delta function，做一个点的物理量分布到周围网格范围内流体网格的distribution的操作，和点周围流体网格物理量插值到点的interpolation操作。就是grid to point and point to grid。需要找到点所在cell及其周围一圈两倍网格间距范围内的cell，不知道想实现这种插值分布函数，还有没有其他更简单的方法吗？请问您有什么看法吗？谢谢！

[马乔]

@oitocfd 我错了，meshSearch的查找效率很高，是primitiveMesh不行。你这是在写pic吗？

[上级]

@oitocfd 之前写过一个小函数，具体是在dpmfoam里用到的，是寻找一个颗粒所在的cell label以及这个cell周围的一圈cell label，不知道能不能帮到你。

template<class Type>
const Foam::labelList& Foam::interpolationCellAround<Type>::cellpointCells
(
    const label celli,
    DynamicList<label>& storage
) const
{
        const labelList& cPoints = this->psi_.mesh().cellPoints()[celli];
 
        storage.clear();
 
        forAll(cPoints, i)
        {
            const label pointi = cPoints[i];
 
            const labelList& pointcellsi = this->psi_.mesh().pointCells()[pointi];
 
            forAll(pointcellsi, j)
            {
                storage.append(pointcellsi[j]);
            }
        }
 
        // Filter duplicates
        if (storage.size() > 1)
        {
            sort(storage);
 
            label n = 1;
            for (label i = 1; i < storage.size(); i++)
            {
                if (storage[i-1] != storage[i])
                {
                    storage[n++] = storage[i];
                }
            }
 
            // truncate addressed list
            storage.setSize(n);
        }
 
        return storage;
}
 
template<class Type>
const Foam::labelList& Foam::interpolationCellAround<Type>::cellpointCells(const label celli) const
{
    return cellpointCells(celli, labels_);
}
 

interpolationCellAround是自定义的类，主要是函数的实现。

[OItoCFD]

@马乔 您好，算是吧 就是每个时间步 我要知道我每个拉格朗日点在哪个cell，然后把力按照cell距离的远近按照delta函数分布到周围所有cell里 我的固体是拉格朗日点代表的 固体可变形 所以每个时间步都要所有点来一次findCell

[OItoCFD]

@上级 在 OpenFOAM编程findCell的诡异问题，对同一坐标寻找cell结果不一致 中说：

mesh().pointCells()

好的 我学习一下 谢谢

[马乔]

@oitocfd 如果粒子追踪的话可以尝试下lagrangian库，当前网格周围一圈网格可以直接用primitiveMesh.cellCells()

[Tong]

应该是精度的问题，polyMesh 的findCell 的原型为：findCell(const point & p,const cellDecomposition decompMode = CELL_TETS) 默认是使用polyMesh里的八叉树进行空间搜索的，vector每个点是一个双精度值，有效数字大约有11位置，而Info只输出了部分有效数字，感觉你直接这样写出来，就会损失精度

vector tmpKnotPoint(1.5,0.72,1.359);
vector point1(1.47727,0.694583,1.30181);
vector point2(1.5,0.694583,1.30181);

如果网格尺寸很小的话可能造成差异，
还有就是，我的记得polyMesh 网格索引树默认的精度是1e-3 如果你的网格尺寸小的话你可以调整一下，具体的你可以去看indexedOctree<>这个类的实现

[OItoCFD]

@tong 好的好的 感谢感谢