Commit b5355b50 authored by Bernd Flemisch's avatar Bernd Flemisch
Browse files

2p2c computeOutflowValues: use boundaryFace() instead of face()

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@7391 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 25cd3b78
......@@ -405,7 +405,7 @@ class TwoPTwoCLocalResidual: public GET_PROP_TYPE(TypeTag, BaseLocalResidual)
}
// add diffusive flux of gas component in liquid phase
Scalar tmp = boundaryVars.molarConcGrad(lPhaseIdx)*boundaryVars.face().normal;
Scalar tmp = boundaryVars.molarConcGrad(lPhaseIdx)*boundaryVars.boundaryFace().normal;
tmp *= -1;
tmp *= boundaryVars.porousDiffCoeff(lPhaseIdx) *
boundaryVars.molarDensityAtIP(lPhaseIdx);
......@@ -416,7 +416,7 @@ class TwoPTwoCLocalResidual: public GET_PROP_TYPE(TypeTag, BaseLocalResidual)
flux[contiLEqIdx] -= tmp * FluidSystem::molarMass(lCompIdx);
// add diffusive flux of liquid component in gas phase
tmp = boundaryVars.molarConcGrad(gPhaseIdx)*boundaryVars.face().normal;
tmp = boundaryVars.molarConcGrad(gPhaseIdx)*boundaryVars.boundaryFace().normal;
tmp *= -1;
tmp *= boundaryVars.porousDiffCoeff(gPhaseIdx) *
boundaryVars.molarDensityAtIP(gPhaseIdx);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment