From c0aad04d1980309e19edc1c5bc55a7f0343fe064 Mon Sep 17 00:00:00 2001
From: Karin Erbertseder <k.erbertseder@web.de>
Date: Wed, 2 Feb 2011 13:57:50 +0000
Subject: [PATCH] changes due to "linearelaticity"
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@5160 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
dumux/boxmodels/1p2c/1p2cfluxvariables.hh | 21 +++++++++++++++++++++
1 file changed, 21 insertions(+)
diff --git a/dumux/boxmodels/1p2c/1p2cfluxvariables.hh b/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
index 1490d4a2ba..6a56063ce0 100644
--- a/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
+++ b/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
@@ -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
--
GitLab