diff --git a/dumux/freeflow/navierstokes/staggered/fluxoversurface.hh b/dumux/freeflow/navierstokes/staggered/fluxoversurface.hh index 5191502f46c45a7058b63e393a589a8123230c86..c26b6cb6344219c4e8ecd7f766a055e6f3551857 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; };