patchInternelField和faceCells对应的值竟然可能不一致

Samuel-Tu

今天发现一件奇怪的事，我用patchInternalField找出nut的一个边界紧邻的内部网格上的值，和用faceCells先找到这个紧邻网格的编号，再再nut场中取出这一个值，发现大多数网格两个值是一样的，但是存在一个网格值不一样，代码如下：

scalarField nutW1=nut.boundaryField()[patchID].patchInternalField();

labelList cellID=mesh.boundary()[patchID].faceCells();
scalarField nutW2;
forAll(cellID,i)
{
    nutW2[i]=nut[cellID[i]];
}

最后发现nutW1和nutW2几乎所有值一样，但是第一个值，nutW1[0]=0，nutW2[0]是一个正常值。我们知道nut是湍流粘度，一定是大于0的。。不知道为什么会这样。好奇怪。。

xpqiu

@Samuel-Tu 在 patchInternelField和faceCells对应的值竟然可能不一致中说：

今天发现一件奇怪的事，我用patchInternalField找出nut的一个边界紧邻的内部网格上的值，和用faceCells先找到这个紧邻网格的编号，再再nut场中取出这一个值，发现大多数网格两个值是一样的，但是存在一个网格值不一样，代码如下：
scalarField nutW1=nut.boundaryField()[patchID].patchInternalField();

labelList cellID=mesh.boundary()[patchID].faceCells();
scalarField nutW2;
forAll(cellID,i)
{
    nutW2[i]=nut[cellID[i]];
}
最后发现nutW1和nutW2几乎所有值一样，但是第一个值，nutW1[0]=0，nutW2[0]是一个正常值。我们知道nut是湍流粘度，一定是大于0的。。不知道为什么会这样。好奇怪。。

试了一下，没法复现这个情况，我这边用这样方式算出来 nutw1和nutw2 是一样的。
PS: 你代码里面先定义 scalarField nutw2, 然后直接用 nutw2[i] 这样来赋值，会导致段错误的吧，nutw2 没有分配空间的啊，可能你实际用的时候不是这样写的？

Samuel-Tu

@xpqiu 对，nutW2是分配了空间了的，大小就是cellID的大小。
造成这个的原因我怀疑是不是因为的我边界有转折导致的。。我是提取的roof这个边界上对应的第一层网格的值的。

因为当时我做计算的时候除以了nutW1之后计算报错了，我才发现nutW1[0]=0，不能被除。然后才改成使用了nutW2的