fvMeshStitcher: Prevent addition of small couplings from interfering …

…with stabilisation Resolves bug report https://bugs.openfoam.org/view.php?id=3965
OpenFOAM · Mar 21, 2023 · 182490e · 182490e
1 parent 00e25bd
commit 182490e
Showing 1 changed file with 12 additions and 8 deletions.
diff --git a/src/finiteVolume/fvMesh/fvMeshStitchers/fvMeshStitcher/fvMeshStitcher.C b/src/finiteVolume/fvMesh/fvMeshStitchers/fvMeshStitcher/fvMeshStitcher.C
@@ -359,6 +359,7 @@ void Foam::fvMeshStitcher::intersectNonConformalCyclic
                                 )
                             );
                     }
+
                     if (!owner)
                     {
                         c.nbr = c;
@@ -367,15 +368,18 @@ void Foam::fvMeshStitcher::intersectNonConformalCyclic
                     SfBf[patchi][patchFacei] = c.nbr.area;
                     CfBf[patchi][patchFacei] = c.nbr.centre;
 
-                    part origP
-                    (
-                        SfBf[origPp.index()][origFacei],
-                        CfBf[origPp.index()][origFacei]
-                    );
-                    origP -= c;
+                    if (i != -1)
+                    {
+                        part origP
+                        (
+                            SfBf[origPp.index()][origFacei],
+                            CfBf[origPp.index()][origFacei]
+                        );
+                        origP -= c;
 
-                    SfBf[origPp.index()][origFacei] = origP.area;
-                    CfBf[origPp.index()][origFacei] = origP.centre;
+                        SfBf[origPp.index()][origFacei] = origP.area;
+                        CfBf[origPp.index()][origFacei] = origP.centre;
+                    }
                 }
             }
         }