Commit c0aad04d authored by Karin Erbertseder's avatar Karin Erbertseder
Browse files

changes due to "linearelaticity"

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@5160 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 79132307
......@@ -90,6 +90,9 @@ public:
{
scvfIdx_ = scvfIdx;
viscosityAtIP_ = Scalar(0);
molarDensityAtIP_ = Scalar(0);
calculateGradients_(problem, element, elemDat);
calculateK_(problem, element, elemDat);
calculateDiffCoeffPM_(problem, element, elemDat);
......@@ -139,6 +142,12 @@ public:
return diffCoeffPM_;
};
Scalar viscosityAtIP() const
{ return viscosityAtIP_;}
Scalar molarDensityAtIP() const
{ return molarDensityAtIP_; }
/*!
* \brief Given the intrinisc permeability times the pressure
......@@ -202,6 +211,12 @@ protected:
tmp = feGrad;
tmp *= elemDat[idx].concentration(1);
concentrationGrad_ += tmp;
// phase viscosity
viscosityAtIP_ += elemDat[idx].viscosity()*face().shapeValue[idx];
//phase moledensity
molarDensityAtIP_ += elemDat[idx].molarDensity()*face().shapeValue[idx];
}
}
else {
......@@ -334,6 +349,12 @@ protected:
//! the intrinsic permeability tensor
Tensor K_;
//! viscosity of the fluid at the integration point
Scalar viscosityAtIP_;
//! molar densities of the fluid at the integration point
Scalar molarDensityAtIP_;
};
} // end namepace
......
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