网格节点编号问题

[piteqiu]

我输出cellPoints，发现六面体单元的节点编号时按照从小到大的顺序输出的：

其顺序与用户指南上的顺序不一致：

请教如何才能输出与上图对应的节点顺序的cellPoints文件 ？

[队长别开枪]

@piteqiu This forms the basis for the input syntax for the blockMesh mesh generator, described in section 5.3.
OpenFOAM的网格数据结构是基于faces的，由faces链接points节点坐标, owner以及neighbour单元信息。单独给出cellPoints是没法确定这个单元的节点分布的。

[piteqiu]

@队长别开枪 请教一个信息遍历问题：cellEdges输出信息如下：

第一行为单元个数，接下来是每个单元含有的edge数及edge的具体编号。
请教一下如何遍历所有单元，获得每个单元包含的edge？
我知道遍历patch获得其开始编号和大小可以按照下述代码：

  const polyMesh& patches = mesh.boundaryMesh();

  forAll(patches, patchi)
    {
      const polyPatch& pp = patches[patchi];

      Info << patchi << '\t' << pp.start() << '\t' << pp.size() << endl;

    }

我自己参照utility写了一个，运行出错：

  const polyMesh& mesh;
  const labelList& cellLabels;
  forAll(cellLabels, i)
    {
        label celli = cellLabels[i];

        const labelList& cEdges = mesh.cellEdges()[celli];

        forAll(cEdges, cEdgeI)
        {
            const edge& e = mesh.edges()[cEdges[cEdgeI]];

            Info<<e[0]<<e[1]<<nl;
        }
    }

请教如何遍历单元获得单元包含edge？

[piteqiu]

没人有类似的需求？