#19895 [CEA] Overlapping volume cells with viscous layer on shared faces

diff --git a/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx b/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
index 64825c0501f0fe4b4cb0ce43506cd8510ad0f71f..1f11791ccbc8c43f11ee6f68292f8b17d9e5d965 100644 (file)
--- a/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
+++ b/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
@@ -263,9 +263,15 @@ static int checkTmpFace(const SMDS_MeshNode* node1,
         bool isDomainInPyramid = ( i3base != i3 );
         return isDomainInPyramid ? HOLE_ID : vol->getshapeId();
       }
-      else
+      else // nbNodes == 6
       {
-        return vol->getshapeId(); // triangle is a prism top
+        int i1 = vol->GetNodeIndex( node1 );
+        if (( i1 == 5 && i2 == 4 && i3 == 3 ) ||
+            ( i1 == 4 && i2 == 3 && i3 == 5 ) ||
+            ( i1 == 3 && i2 == 5 && i3 == 4 ))
+          return HOLE_ID;
+        else
+          return vol->getshapeId(); // triangle is a prism top
       }
     }
   }