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

OItoCFD

我把代码和测试算例精简后打包发在了网盘，如果有兴趣的老师和朋友可以帮忙看一下，谢谢！blockMesh后运行
链接：https://pan.baidu.com/s/1HEy6w8MgbeXu9EsMcayung
提取码：1111
--来自百度网盘超级会员V4的分享

李东岳

代码上看起来是没问题。

label cellId = searchEngine_.findCell(tmpOffSetPoint);

这句话之前赋值为0然后再赋值呢？

李东岳

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

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

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

OItoCFD

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

代码上看起来是没问题。

label cellId = searchEngine_.findCell(tmpOffSetPoint);

这句话之前赋值为0然后再赋值呢？

老师，这样做也是一样的结果，还是重复的CellID

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 好的好的感谢感谢

CFD中文网

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