From 09c4d061f4102f9f5afa16eb482b9a8e374b2190 Mon Sep 17 00:00:00 2001
From: Kilian Weishaupt <kilian.weishaupt@iws.uni-stuttgart.de>
Date: Mon, 24 Jul 2017 15:21:38 +0200
Subject: [PATCH] [velocityoutput][box1d] Include extrusion factor
---
dumux/porousmediumflow/implicit/velocityoutput.hh | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/dumux/porousmediumflow/implicit/velocityoutput.hh b/dumux/porousmediumflow/implicit/velocityoutput.hh
index 11b5bf0601..69d526dfa5 100644
--- a/dumux/porousmediumflow/implicit/velocityoutput.hh
+++ b/dumux/porousmediumflow/implicit/velocityoutput.hh
@@ -161,9 +161,11 @@ public:
// get the volume flux divided by the area of the
// subcontrolvolume face in the reference element
- // TODO: Divide by extrusion factor!!?
Scalar localArea = scvfReferenceArea_(geomType, scvf.index());
Scalar flux = fluxVars.advectiveFlux(phaseIdx, upwindTerm) / localArea;
+ flux /= problem_.extrusionFactor(element,
+ fvGeometry.scv(scvf.insideScvIdx()),
+ problem_.model().elementSolution(element, problem_.model().curSol()));
tmpVelocity *= flux;
const int eIdxGlobal = problem_.elementMapper().index(element);
--
GitLab