From 70c697fdcb0ec6a8f9e5d4e8d301631f9ee4e040 Mon Sep 17 00:00:00 2001
From: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
Date: Fri, 2 Jul 2021 15:40:40 +0100
Subject: [PATCH] BUG: momentum: boundary field mismatch on "empty" patches
(fixes #2121)
---
src/functionObjects/field/momentum/momentum.C | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/src/functionObjects/field/momentum/momentum.C b/src/functionObjects/field/momentum/momentum.C
index 618690f2a58..0162732793b 100644
--- a/src/functionObjects/field/momentum/momentum.C
+++ b/src/functionObjects/field/momentum/momentum.C
@@ -30,6 +30,7 @@ License
#include "volFields.H"
#include "cellSet.H"
#include "cylindricalRotation.H"
+#include "emptyPolyPatch.H"
#include "addToRunTimeSelectionTable.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
@@ -556,6 +557,11 @@ bool Foam::functionObjects::momentum::write()
forAll(pbm, patchi)
{
+ if (isA<emptyPolyPatch>(pbm[patchi]))
+ {
+ continue;
+ }
+
const auto& pts = pbm[patchi].faceCentres();
const label len = pts.size();
--
GitLab