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

Fix updateState() saturation update.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2017-08-09 15:47:39 +02:00
parent 608a674858
commit 1e0facec4a
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/autodiff/BlackoilReorderingTransportModel.hpp
View File

@ -855,13 +855,14 @@ namespace Opm {
// Get saturation updates.
const double dsw = dx[0];
double dso = -dsw;
double dsg = 0.0;
auto& hcstate = state_.reservoir_state.hydroCarbonState()[cell];
if (hcstate == HydroCarbonState::GasAndOil) {
dsg = dx[1];
dso -= dsg;
} else if (hcstate == HydroCarbonState::GasOnly) {
dsg = -dsw;
}
const double dso = -(dsw + dsg);
// Handle too large saturation changes.
const double maxval = std::max(std::fabs(dsw), std::max(std::fabs(dso), std::fabs(dsg)));