From 2c5fd082584c95d20a075ed2fe9a86e487e24f0e Mon Sep 17 00:00:00 2001
From: Klaus Mosthaf <klmos@env.dtu.dk>
Date: Mon, 5 Mar 2012 13:27:01 +0000
Subject: [PATCH] The averaged viscosity is used for the computation of mu grad
v (part of the momentum flux) instead of upwinding
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@7948 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
dumux/freeflow/stokes/stokeslocalresidual.hh | 6 +-----
1 file changed, 1 insertion(+), 5 deletions(-)
diff --git a/dumux/freeflow/stokes/stokeslocalresidual.hh b/dumux/freeflow/stokes/stokeslocalresidual.hh
index 5a983e145e..fdc3c85761 100644
--- a/dumux/freeflow/stokes/stokeslocalresidual.hh
+++ b/dumux/freeflow/stokes/stokeslocalresidual.hh
@@ -227,11 +227,7 @@ protected:
// TODO: dilatation term has to be accounted for in outflow, coupling, neumann
// velGradComp[velIdx] += 2./3*fluxVars.velocityDivAtIP;
- if (massUpwindWeight_ == 1.0) // fully upwind
- velGradComp *= up.viscosity();
- else
- velGradComp *= massUpwindWeight_ * up.viscosity() +
- (1.0 - massUpwindWeight_) * dn.viscosity();
+ velGradComp *= fluxVars.viscosityAtIP();
flux[momentumXIdx + velIdx] -=
velGradComp*fluxVars.face().normal;
--
GitLab