Skip to content
Snippets Groups Projects
Commit 20e4ff18 authored by Katherina Baber's avatar Katherina Baber
Browse files

added porosity to spatial parameters, corrected calculation of

Kmvpnormal_


git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@6693 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent d7845600
No related branches found
No related tags found
No related merge requests found
......@@ -266,7 +266,7 @@ protected:
K_.mv(potentialGrad_, Kmvp);
KmvpNormal_ = 0;
for (int i = 0; i < dim; ++i)
KmvpNormal_ = Kmvp[i]*boundaryFace_->normal[i];
KmvpNormal_ += Kmvp[i]*boundaryFace_->normal[i];
KmvpNormal_ *= -1;
const VolumeVariables &vertDat = elemDat[scvIdx_];
......
......@@ -67,12 +67,20 @@ public:
int scvIdx) const
{ return 1e-10; }
/* \brief Define the porosity.
*
* \param element The finite element
* \param fvElemGeom The finite volume geometry
* \param scvIdx The local index of the sub-control volume where
*/
Scalar porosity(const Element &element,
const FVElementGeometry &fvElemGeom,
int scvIdx) const
{ return 0.4; }
};
} // end namespace
#endif
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment