OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

EclHysteresisTwoPhaseLaw: use '<=' instead of '<' to decide whether to use the imbibition or the drainage curve

Browse Source 
this does not change the results in any way, but for the Norne deck,
it improves performance (i.e., It Must Be Right (TM))...
This commit is contained in:
Andreas Lauser
2015-08-18 13:53:18 +02:00
parent 57948012a5
commit ef28c31368
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
opm/material/fluidmatrixinteractions/EclHysteresisTwoPhaseLaw.hpp
View File

@@ -229,7 +229,7 @@ public:
// if it is enabled, use either the drainage or the imbibition curve. if the
// imbibition curve is used, the saturation must be shifted.
if (Sw < params.krwSwMdc())
if (Sw <= params.krwSwMdc())
return EffectiveLaw::twoPhaseSatKrw(params.drainageParams(), Sw);
return EffectiveLaw::twoPhaseSatKrw(params.imbibitionParams(),
@@ -255,7 +255,7 @@ public:
// if it is enabled, use either the drainage or the imbibition curve. if the
// imbibition curve is used, the saturation must be shifted.
if (Sw < params.krnSwMdc())
if (Sw <= params.krnSwMdc())
return EffectiveLaw::twoPhaseSatKrn(params.drainageParams(), Sw);
return EffectiveLaw::twoPhaseSatKrn(params.imbibitionParams(),