diff --git a/dumux/implicit/box/fvelementgeometry.hh b/dumux/implicit/box/fvelementgeometry.hh index 012e72d2cdf5529441a61db363822a139391b0b6..4f0f4cadb86d870a441dfd876df2eb497df67106 100644 --- a/dumux/implicit/box/fvelementgeometry.hh +++ b/dumux/implicit/box/fvelementgeometry.hh @@ -943,6 +943,10 @@ public: this->crossProduct(elemNormal, elemVec1, elemVec2); this->crossProduct(scvFace.normal, faceVec, elemNormal); + // scale normal vector by the area + scvFace.normal /= scvFace.normal.two_norm(); + scvFace.normal *= faceVec.two_norm(); + const auto diffVec = subContVol[j].global - subContVol[i].global; // make sure the normal points to the right direction if (scvFace.normal * diffVec < 0)