From 66610a3150ff37deb197ca901484b46f60e8f5eb Mon Sep 17 00:00:00 2001
From: Kilian <kilian.weishaupt@iws.uni-stuttgart.de>
Date: Wed, 27 Feb 2019 09:53:43 +0100
Subject: [PATCH] [ff][fluxoversurface] Consider extrusion factor for volume
flux
---
dumux/freeflow/navierstokes/staggered/fluxoversurface.hh | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/dumux/freeflow/navierstokes/staggered/fluxoversurface.hh b/dumux/freeflow/navierstokes/staggered/fluxoversurface.hh
index 5191502f46..c26b6cb634 100644
--- a/dumux/freeflow/navierstokes/staggered/fluxoversurface.hh
+++ b/dumux/freeflow/navierstokes/staggered/fluxoversurface.hh
@@ -250,7 +250,10 @@ public:
const auto& elemFluxVarsCache)
{
CellCenterPrimaryVariables result(0.0);
- result[0] = elemFaceVars[scvf].velocitySelf() * scvf.area() * scvf.directionSign();
+ const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()];
+ const auto& outsideVolVars = elemVolVars[scvf.outsideScvIdx()];
+ const Scalar extrusionFactor = harmonicMean(insideVolVars.extrusionFactor(), outsideVolVars.extrusionFactor());
+ result[0] = elemFaceVars[scvf].velocitySelf() * scvf.area() * extrusionFactor * scvf.directionSign();
return result;
};
--
GitLab