Commit a457dcac authored by Thomas Fetzer's avatar Thomas Fetzer
Browse files

[multidomain] added eddy viscosity, diffusivity and conductivity to

couplinglocalresiduals. this is done in consistency with the freeflow
stokes models. the eddy coefficient are zero for the stokes case.

approved by gruenich


git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@12738 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent e3df4179
......@@ -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));
......
......@@ -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);
......
......@@ -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
......
......@@ -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
......
Markdown is supported
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