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

Explicitly Prevent Negative 3-p Kro Values in Stone 2 Model

Browse Source 
This makes the material balance equations converge more smoothly in
corner cases and makes a real field case run at all.
This commit is contained in:
Bård Skaflestad 2021-02-02 16:24:46 +01:00
parent 372bb122a6
commit 1479adf6af
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/material/fluidmatrixinteractions/EclStone2Material.hpp
View File

@ -362,7 +362,7 @@ public:
const Evaluation krg = GasOilMaterialLaw::twoPhaseSatKrn(params.gasOilParams(), 1 - Swco - Sg);
const Evaluation krog = GasOilMaterialLaw::twoPhaseSatKrw(params.gasOilParams(), 1 - Swco - Sg);
return krocw * ((krow/krocw + krw) * (krog/krocw + krg) - krw - krg);
return Opm::max(krocw * ((krow/krocw + krw) * (krog/krocw + krg) - krw - krg), Evaluation{0});
}
/*!