diff --git a/dumux/multidomain/2cnistokes2p2cni/2cnistokes2p2cnilocaloperator.hh b/dumux/multidomain/2cnistokes2p2cni/2cnistokes2p2cnilocaloperator.hh
index 7672bed00633b85159630467ccf6de6b0a145377..90c217942a5446b4242cb73f33e18da0208aaa04 100644
--- a/dumux/multidomain/2cnistokes2p2cni/2cnistokes2p2cnilocaloperator.hh
+++ b/dumux/multidomain/2cnistokes2p2cni/2cnistokes2p2cnilocaloperator.hh
@@ -119,7 +119,7 @@ public:
                 const Scalar conductiveFlux =
                     bfNormal1 *
                     boundaryVars1.temperatureGrad() *
-                    boundaryVars1.thermalConductivity();
+                    (boundaryVars1.thermalConductivity() + boundaryVars1.thermalEddyConductivity());
 
                 couplingRes2.accumulate(lfsu_n.child(energyEqIdx2), vertInElem2,
                                         -(convectiveFlux - conductiveFlux));
diff --git a/dumux/multidomain/2cstokes2p2c/2cstokes2p2clocaloperator.hh b/dumux/multidomain/2cstokes2p2c/2cstokes2p2clocaloperator.hh
index 4a27349c026992069a31a1aaf8d903d196f68256..1f6e962c86bc66b34572035bf78545362d4e8d52 100644
--- a/dumux/multidomain/2cstokes2p2c/2cstokes2p2clocaloperator.hh
+++ b/dumux/multidomain/2cstokes2p2c/2cstokes2p2clocaloperator.hh
@@ -388,7 +388,7 @@ class TwoCStokesTwoPTwoCLocalOperator :
                 const Scalar diffusiveFlux =
                     bfNormal1.two_norm() *
                     normalMoleFracGrad *
-                    boundaryVars1.diffusionCoeff(transportCompIdx1) *
+                    (boundaryVars1.diffusionCoeff(transportCompIdx1) + boundaryVars1.eddyDiffusivity()) *
                     boundaryVars1.molarDensity() *
                     FluidSystem::molarMass(transportCompIdx1);
 
@@ -412,7 +412,7 @@ class TwoCStokesTwoPTwoCLocalOperator :
                         const Scalar diffusiveFluxAtCorner =
                             bfNormal1 *
                             boundaryVars1.moleFractionGrad(transportCompIdx1) *
-                            boundaryVars1.diffusionCoeff(transportCompIdx1) *
+                            (boundaryVars1.diffusionCoeff(transportCompIdx1) + boundaryVars1.eddyDiffusivity()) *
                             boundaryVars1.molarDensity() *
                             FluidSystem::molarMass(transportCompIdx1);
 
@@ -439,7 +439,7 @@ class TwoCStokesTwoPTwoCLocalOperator :
                     const Scalar diffusiveFlux =
                         bfNormal1 *
                         boundaryVars1.moleFractionGrad(transportCompIdx1) *
-                        boundaryVars1.diffusionCoeff(transportCompIdx1) *
+                        (boundaryVars1.diffusionCoeff(transportCompIdx1) + boundaryVars1.eddyDiffusivity()) *
                         boundaryVars1.molarDensity() *
                         FluidSystem::molarMass(transportCompIdx1);
 
diff --git a/dumux/multidomain/couplinglocalresiduals/stokesnccouplinglocalresidual.hh b/dumux/multidomain/couplinglocalresiduals/stokesnccouplinglocalresidual.hh
index 33853f1f85d6d91e2b0fc768b5b9b73e183229da..52ae8ee49deb155a1b455e8094aab8864920bf58 100644
--- a/dumux/multidomain/couplinglocalresiduals/stokesnccouplinglocalresidual.hh
+++ b/dumux/multidomain/couplinglocalresiduals/stokesnccouplinglocalresidual.hh
@@ -159,7 +159,8 @@ public:
                                 boundaryVars.face().normal;
 
                         boundaryVars.velocityGrad().umv(boundaryFaceNormal, muGradVelNormal);
-                        muGradVelNormal *= boundaryVars.dynamicViscosity();
+                        muGradVelNormal *= (boundaryVars.dynamicViscosity()
+                                            + boundaryVars.dynamicEddyViscosity());
 
                         for (int i=0; i < this->residual_.size(); i++)
                             Valgrind::CheckDefined(this->residual_[i]);
@@ -352,7 +353,8 @@ protected:
                     this->residual_[scvIdx][dimIdx] += beaversJosephCoeff *
                                             boundaryFaceArea *
                                             tangentialV[dimIdx] *
-                                            boundaryVars.dynamicViscosity();
+                                            (boundaryVars.dynamicViscosity()
+                                             + boundaryVars.dynamicEddyViscosity());
 
                 ////////////////////////////////////////////////////////////////////////////////////
                 //normal component has only to be set if no coupling conditions are defined
diff --git a/dumux/multidomain/couplinglocalresiduals/stokesncnicouplinglocalresidual.hh b/dumux/multidomain/couplinglocalresiduals/stokesncnicouplinglocalresidual.hh
index 04178fd731502098e5f92a54293e728a93e7e66a..5bfeb58f894f1e6fba74a6b4423e40d8f39255f7 100644
--- a/dumux/multidomain/couplinglocalresiduals/stokesncnicouplinglocalresidual.hh
+++ b/dumux/multidomain/couplinglocalresiduals/stokesncnicouplinglocalresidual.hh
@@ -162,7 +162,8 @@ namespace Dumux
 							boundaryVars.face().normal;
 							
 							boundaryVars.velocityGrad().umv(boundaryFaceNormal, muGradVelNormal);
-							muGradVelNormal *= boundaryVars.dynamicViscosity();
+							muGradVelNormal *= (boundaryVars.dynamicViscosity()
+                                                + boundaryVars.dynamicEddyViscosity());
 							
 							for (int i=0; i < this->residual_.size(); i++)
 								Valgrind::CheckDefined(this->residual_[i]);
@@ -351,7 +352,7 @@ namespace Dumux
 						this->residual_[scvIdx][dimIdx] += beaversJosephCoeff *
 						boundaryFaceArea *
 						tangentialV[dimIdx] *
-						boundaryVars.dynamicViscosity();
+						(boundaryVars.dynamicViscosity() + boundaryVars.dynamicEddyViscosity());
 					
 					////////////////////////////////////////////////////////////////////////////////////
 					//normal component has only to be set if no coupling conditions are defined