距壁面第二个网格中心速度代码求助

[yangniannian]

壁面函数代码里面距壁面第一个网格中心速度是Up=turbModel.U().boundaryField()[patchi].patchInternalField()
我还需要第二个网格中心的速度，请问大家代码应该怎么写？

[李东岳]

下面这种情况红点是第一个网格中心，哪个绿点是第二个网格呢？

[yangniannian]

抱歉问题描述不严谨
是结构六面体网格的中心点

[李东岳]

即使六面体，边边角角的第二层网格也不好判断是哪个，
目前OpenFOAM没有现成的可以调用其他层网格的代码，虽然寻找某个网格周围的网格是可以实现的，

[马乔]

这个问题主要是找到边界面正对的面，然后通过neighbor-owner找到单元。可以通过边界面的点找到共点的edge，然后找到另一顶点，再来查找面。呃，好像有点繁琐。第二个可以通过particle类，将边界面法相取反作为位移方向，edge长度作为位移，trackToFace找到面。记得用labelPair保存相邻单元。