OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

ebos: slightly improve the transmissibility calculation

Browse Source
This commit is contained in:
Andreas Lauser 2016-03-02 13:25:38 +01:00
parent 567adcab6f
commit 86daf25aa6
1 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
applications/ebos/ecltransmissibility.hh
View File

@ -279,18 +279,20 @@ private:
const DimVector& distance, const DimVector& distance,
const DimMatrix& perm) const const DimMatrix& perm) const
{ {
Scalar isArea = is.geometry().volume();
DimVector n = is.centerUnitOuterNormal();
unsigned dimIdx = faceIdx/2; unsigned dimIdx = faceIdx/2;
assert(dimIdx < dimWorld); assert(dimIdx < dimWorld);
halfTrans = perm[dimIdx][dimIdx]; halfTrans = perm[dimIdx][dimIdx];
halfTrans *= is.geometry().volume(); halfTrans *= isArea;
const auto &normal = is.centerUnitOuterNormal();
Scalar val = 0; Scalar val = 0;
for (unsigned i = 0; i < normal.size(); ++i) for (unsigned i = 0; i < n.size(); ++i)
val += is.centerUnitOuterNormal()[i]*distance[i]; val += n[i]*distance[i];
halfTrans *= std::abs<Scalar>(val); halfTrans *= std::abs<Scalar>(val);
halfTrans /= distance*distance; halfTrans /= distance.two_norm2();
} }
DimVector distanceVector_(const Intersection& is, DimVector distanceVector_(const Intersection& is,