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

Don't include connections with negligible flow in thpress caclulations

Browse Source
This commit is contained in:
Tor Harald Sandve 2017-01-02 16:03:19 +01:00
parent ef3de796bc
commit 0dbbfbd8a6
1 changed files with 6 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ebos/eclthresholdpressure.hh
View File

@ -197,6 +197,12 @@ private:
// the current face is not at the boundary between EQUIL regions! // the current face is not at the boundary between EQUIL regions!
continue; continue;
// don't include connections with negligible flow
const Scalar& trans = simulator_.problem().transmissibility(elemCtx, i, j);
const Scalar& faceArea = face.area();
if ( std::abs(faceArea * trans) < 1e-18)
continue;
// determine the maximum difference of the pressure of any phase over the // determine the maximum difference of the pressure of any phase over the
// intersection // intersection
Scalar pth = 0.0; Scalar pth = 0.0;