OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-25 18:50:19 -06:00
Code Issues Packages Projects Releases Wiki Activity

Fix bug in saturation initialisation.

Browse Source 
We shall only use gas-water capillary to initialise when we would get
unphysical saturations otherwise.
This commit is contained in:
Atgeirr Flø Rasmussen 2014-02-21 14:47:14 +01:00 committed by Andreas Lauser
parent e51b161c89
commit 734046e16e
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/core/simulator/initStateEquil.hpp
View File

@ -796,9 +796,10 @@ namespace Opm
sg = satFromPc(props, gaspos, cell, pcog, increasing);
phase_saturations[gaspos][local_index] = sg;
}
if (gas && water && sg > smin[gaspos] && sw > smin[waterpos]) {
if (gas && water && (sg + sw > 1.0)) {
// Overlapping gas-oil and oil-water transition
// zones. Must recalculate using gas-water
// zones can lead to unphysical saturations when
// treated as above. Must recalculate using gas-water
// capillary pressure.
const double pcgw = phase_pressures[gaspos][local_index] - phase_pressures[waterpos][local_index];
sw = satFromSumOfPcs(props, waterpos, gaspos, cell, pcgw);