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的。。不知道为什么会这样。好奇怪。。


  • 教授

    @Samuel-TupatchInternelField和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 没有分配空间的啊,可能你实际用的时候不是这样写的?



  • @xpqiu 对,nutW2是分配了空间了的,大小就是cellID的大小。
    造成这个的原因我怀疑是不是因为的我边界有转折导致的。。我是提取的roof这个边界上对应的第一层网格的值的。
    121.png
    因为当时我做计算的时候除以了nutW1之后计算报错了,我才发现nutW1[0]=0,不能被除。然后才改成使用了nutW2的


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1