fvMeshStitcher: Prevent addition of small couplings from interfering with stabilisation

Resolves bug report https://bugs.openfoam.org/view.php?id=3965
This commit is contained in:
Will Bainbridge
2023-03-21 16:12:03 +00:00
parent 00e25bd828
commit 182490e0b2

View File

@ -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;
}
}
}
}