diff --git a/dumux/multidomain/2cnistokes2p2cni/stokesncnicouplinglocalresidual.hh b/dumux/multidomain/2cnistokes2p2cni/stokesncnicouplinglocalresidual.hh index 948499aa359c902e163b35a59a5869a4e9c0d97c..5ca274a9b13d423bbd811c589bbc9bbffeb33df0 100644 --- a/dumux/multidomain/2cnistokes2p2cni/stokesncnicouplinglocalresidual.hh +++ b/dumux/multidomain/2cnistokes2p2cni/stokesncnicouplinglocalresidual.hh @@ -147,10 +147,8 @@ public: // if p.n comes from the pm if (bcTypes.isCouplingNeumann(momentumYIdx) || bcTypes.isCouplingMortar(momentumYIdx)) { - DimVector pressureCorrection(intersection.centerUnitOuterNormal()); - pressureCorrection *= volVars.pressure(); - this->residual_[scvIdx][momentumYIdx] += pressureCorrection[momentumYIdx] - * boundaryVars.face().area; + this->residual_[scvIdx][momentumYIdx] += volVars.pressure() + * boundaryVars.face().normal[momentumYIdx]; Valgrind::CheckDefined(this->residual_[scvIdx][momentumYIdx]); } diff --git a/dumux/multidomain/2cstokes2p2c/stokesnccouplinglocalresidual.hh b/dumux/multidomain/2cstokes2p2c/stokesnccouplinglocalresidual.hh index ec2f7bcf18c1effb516c2ce4c55d3ff2fb6d877e..915db02897db946687621202d600e9774e0930aa 100644 --- a/dumux/multidomain/2cstokes2p2c/stokesnccouplinglocalresidual.hh +++ b/dumux/multidomain/2cstokes2p2c/stokesnccouplinglocalresidual.hh @@ -157,10 +157,8 @@ public: // if p.n comes from the pm if (bcTypes.isCouplingNeumann(momentumYIdx) || bcTypes.isCouplingMortar(momentumYIdx)) { - DimVector pressureCorrection(intersection.centerUnitOuterNormal()); - pressureCorrection *= volVars.pressure(); - this->residual_[scvIdx][momentumYIdx] += pressureCorrection[momentumYIdx] - * boundaryVars.face().area; + this->residual_[scvIdx][momentumYIdx] += volVars.pressure() + * boundaryVars.face().normal[momentumYIdx]; Valgrind::CheckDefined(this->residual_[scvIdx][momentumYIdx]); }