Merge pull request #445 from bska/prevent-neg-kro-stone2

Explicitly Prevent Negative 3-p Kro Values in Stone 2 Model
This commit is contained in:
Atgeirr Flø Rasmussen 2021-03-22 11:27:56 +01:00 committed by GitHub
commit 3218ae4286

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});
}
/*!