diff --git a/dumux/discretization/cellcentered/mpfa/fourierslaw.hh b/dumux/discretization/cellcentered/mpfa/fourierslaw.hh index b8072e8289d6f7c20cbe3a3f783b34364872d3f0..27c5d07a3d4759bc3d20fc46f1e54bf0e4d907b1 100644 --- a/dumux/discretization/cellcentered/mpfa/fourierslaw.hh +++ b/dumux/discretization/cellcentered/mpfa/fourierslaw.hh @@ -166,13 +166,14 @@ public: return useTpfaBoundary ? flux : flux + fluxVarsCache.heatNeumannFlux(); // Handle interior boundaries - flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, fluxVarsCache); + flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, elemVolVars, fluxVarsCache); // return overall resulting flux return useTpfaBoundary ? flux : flux + fluxVarsCache.heatNeumannFlux(); } static Scalar computeInteriorBoundaryContribution(const FVElementGeometry& fvGeometry, + const ElementVolumeVariables& elemVolVars, const FluxVariablesCache& fluxVarsCache) { // obtain the transmissibilites associated with all pressures