From b20645786aa8a5edd0dc6099ece5f88420d13bc9 Mon Sep 17 00:00:00 2001 From: Kilian Weishaupt Date: Thu, 8 Sep 2016 14:16:44 +0200 Subject: [PATCH 1/2] [stokes][localResidual] Account for extrusion factor in outflow *Flux leaving the domain needs to be scaled by extrusion factor --- dumux/freeflow/stokes/localresidual.hh | 1 + 1 file changed, 1 insertion(+) diff --git a/dumux/freeflow/stokes/localresidual.hh b/dumux/freeflow/stokes/localresidual.hh index 93b9e4d833..7951945fc9 100644 --- a/dumux/freeflow/stokes/localresidual.hh +++ b/dumux/freeflow/stokes/localresidual.hh @@ -519,6 +519,7 @@ protected: PrimaryVariables values(0.0); asImp_()->computeFlux(values, boundaryFaceIdx, /*onBoundary=*/true); + values *= this->curVolVars_(scvIdx).extrusionFactor(); Valgrind::CheckDefined(values); for (int eqIdx = 0; eqIdx < numEq; ++eqIdx) -- GitLab From 1d797b8b93cd69a9234acb8b76fdf5a1d8885fe7 Mon Sep 17 00:00:00 2001 From: Kilian Weishaupt Date: Thu, 8 Sep 2016 14:17:07 +0200 Subject: [PATCH 2/2] [box][localResidual] Account for extrusion factor in outflow *Flux leaving the domain needs to be scaled be extrusion factor --- dumux/implicit/box/localresidual.hh | 1 + 1 file changed, 1 insertion(+) diff --git a/dumux/implicit/box/localresidual.hh b/dumux/implicit/box/localresidual.hh index 4b69761d83..1eacceb3b1 100644 --- a/dumux/implicit/box/localresidual.hh +++ b/dumux/implicit/box/localresidual.hh @@ -204,6 +204,7 @@ protected: //calculate outflow fluxes PrimaryVariables values(0.0); this->asImp_().computeFlux(values, boundaryFaceIdx, true); + values *= this->curVolVars_(scvIdx).extrusionFactor(); Valgrind::CheckDefined(values); for (int equationIdx = 0; equationIdx < numEq; ++equationIdx) -- GitLab