From 1644c79232c0605ded1aab1319e96d26c940d58c Mon Sep 17 00:00:00 2001
From: Martin Schneider <martin.schneider@iws.uni-stuttgart.de>
Date: Wed, 26 Apr 2023 19:59:03 +0200
Subject: [PATCH] [disc][pq1bubble] Allow to specify if boundary scvf is
overlapping
---
dumux/discretization/pq1bubble/fvgridgeometry.hh | 3 ++-
dumux/discretization/pq1bubble/geometryhelper.hh | 5 +++++
2 files changed, 7 insertions(+), 1 deletion(-)
diff --git a/dumux/discretization/pq1bubble/fvgridgeometry.hh b/dumux/discretization/pq1bubble/fvgridgeometry.hh
index c8bd0df627..3f694fc5a7 100644
--- a/dumux/discretization/pq1bubble/fvgridgeometry.hh
+++ b/dumux/discretization/pq1bubble/fvgridgeometry.hh
@@ -349,7 +349,8 @@ private:
intersection.centerUnitOuterNormal(),
std::move(scvPair),
scvfLocalIdx,
- typename SubControlVolumeFace::Traits::BoundaryFlag{ intersection }
+ typename SubControlVolumeFace::Traits::BoundaryFlag{ intersection },
+ geometryHelper.isOverlappingBoundaryScvf(localFacetIndex)
);
// store look-up map to construct boundary scvf geometries
diff --git a/dumux/discretization/pq1bubble/geometryhelper.hh b/dumux/discretization/pq1bubble/geometryhelper.hh
index 966e6366f1..843c93b4f8 100644
--- a/dumux/discretization/pq1bubble/geometryhelper.hh
+++ b/dumux/discretization/pq1bubble/geometryhelper.hh
@@ -406,6 +406,11 @@ public:
return true;
}
+ bool isOverlappingBoundaryScvf(unsigned int localFacetIndex) const
+ {
+ return false;
+ }
+
bool isOverlappingScv(unsigned int localScvIndex) const
{
if (localScvIndex < boxHelper_.numScv())
--
GitLab