From 0b92e472725ef7957dd6a2d1b09e4775b2660b6c Mon Sep 17 00:00:00 2001
From: DennisGlaeser <dennis.glaeser@iws.uni-stuttgart.de>
Date: Fri, 10 Feb 2017 18:00:02 +0100
Subject: [PATCH] [mpfa][fickslaw] pass elemVolVars to interior boundary
method()
---
dumux/discretization/cellcentered/mpfa/fickslaw.hh | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/dumux/discretization/cellcentered/mpfa/fickslaw.hh b/dumux/discretization/cellcentered/mpfa/fickslaw.hh
index dc79901c82..be41e0885a 100644
--- a/dumux/discretization/cellcentered/mpfa/fickslaw.hh
+++ b/dumux/discretization/cellcentered/mpfa/fickslaw.hh
@@ -219,7 +219,7 @@ public:
return useTpfaBoundary ? flux*rho*effFactor : flux*rho*effFactor + fluxVarsCache.componentNeumannFlux(compIdx);
// Handle interior boundaries
- flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, fluxVarsCache, getX, phaseIdx, compIdx);
+ flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, elemVolVars, fluxVarsCache, getX, phaseIdx, compIdx);
// return overall resulting flux
return useTpfaBoundary ? flux*rho*effFactor : flux*rho*effFactor + fluxVarsCache.componentNeumannFlux(compIdx);
@@ -314,6 +314,7 @@ public:
template<typename GetXFunction>
static Scalar computeInteriorBoundaryContribution(const FVElementGeometry& fvGeometry,
+ const ElementVolumeVariables& elemVolVars,
const FluxVariablesCache& fluxVarsCache,
const GetXFunction& getX,
unsigned int phaseIdx, unsigned int compIdx)
--
GitLab