Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    owner neighbour 网格寻址

    OpenFOAM
    2
    3
    2092
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • 李东岳
      李东岳 管理员 last edited by 李东岳

          const labelUList& owner = mesh_.owner();
          const labelUList& neighbour = mesh_.neighbour();
      
          forAll(weis_, kth)
          {
              forAll(sumWeisF[kth], cellI)
              {
      
                  Info << "Current cell index is " << cellI << nl;
      
                  //- Loop the faces of cellI
                  forAll(mesh_.cells()[cellI], faceI)
                  {
                      //- Check if it is an internal face
                      if (mesh_.isInternalFace(mesh_.cells()[cellI][faceI]))   
                      {
                          label faceIndex = mesh_.cells()[cellI][faceI];
                          Info << " its faceIndex is " << faceIndex;
                          Info << ", and this face's owner is "
                              << owner[faceIndex] << ", its neighbour is "
                              << neighbour[faceIndex] << nl;
                      }
                  }
                  
                  Info << nl;
              } 
          }
      

      输出如下:

      Current cell index is 851
       its faceIndex is 1684, and this face's owner is 851, its neighbour is 852
       its faceIndex is 1685, and this face's owner is 851, its neighbour is 896
       its faceIndex is 1596, and this face's owner is 806, its neighbour is 851
       its faceIndex is 1682, and this face's owner is 850, its neighbour is 851
      
      Current cell index is 852
       its faceIndex is 1686, and this face's owner is 852, its neighbour is 853
       its faceIndex is 1687, and this face's owner is 852, its neighbour is 897
       its faceIndex is 1598, and this face's owner is 807, its neighbour is 852
       its faceIndex is 1684, and this face's owner is 851, its neighbour is 852
      
      

      希望有用

      CFD课程 改成线上了 http://dyfluid.com/class.html
      CFD高性能服务器 http://dyfluid.com/servers.html

      1 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by

        注意在变量赋值的时候要使用faceIndex而不是faceI,因为p[faceI]只是一直在1 2 3 4 5...中循环,p[faceIndex]是真正的面址

        CFD课程 改成线上了 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        1 Reply Last reply Reply Quote
        • 李东岳
          李东岳 管理员 last edited by

          4年过去了,当时发完这个贴就出去玩了。:jingya:

          CFD课程 改成线上了 http://dyfluid.com/class.html
          CFD高性能服务器 http://dyfluid.com/servers.html

          1 Reply Last reply Reply Quote
          • First post
            Last post

          CFD中文网 | 东岳流体 | 京ICP备15017992号-2