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

EclDefaultMaterial: do not add the connate water saturation to Sw_ow

Browse Source
This commit is contained in:
Andreas Lauser 2015-08-14 16:27:41 +02:00
parent 8e51155f24
commit 14e4692386
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -296,7 +296,7 @@ public:
Evaluation Sg = FsToolbox::template toLhs<Evaluation>(fluidState.saturation(gasPhaseIdx));
Evaluation Sw_ow = Sg + Sw;
Evaluation So_go = 1.0 + Swco - Sw_ow;
Evaluation So_go = 1.0 - Sw_ow;
const Evaluation& kro_ow = OilWaterMaterialLaw::twoPhaseSatKrn(params.oilWaterParams(), Sw_ow);
const Evaluation& kro_go = GasOilMaterialLaw::twoPhaseSatKrw(params.gasOilParams(), So_go);
@ -348,7 +348,7 @@ public:
Sg = std::min(1.0, std::max(0.0, Sg));
Scalar Sw_ow = Sg + std::max(Swco, Sw);
Scalar So_go = 1 + Swco - Sw_ow;
Scalar So_go = 1 + Sw_ow;
params.oilWaterParams().update(/*pcSw=*/Sw, /*krwSw=*/1 - Sg, /*krnSw=*/Sw_ow);
params.gasOilParams().update(/*pcSw=*/1 - Sg, /*krwSw=*/So_go, /*krnSw=*/1 - Sg);