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 phase index for WELLPI for MSW

Browse Source
This commit is contained in:
Tor Harald Sandve
2021-01-21 11:58:24 +01:00
parent 036a021c63
commit 91ab0f4d93
2 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
opm/simulators/wells/MultisegmentWell.hpp
View File

@@ -240,6 +240,7 @@ namespace Opm
using Base::phaseUsage; using Base::phaseUsage;
using Base::name; using Base::name;
using Base::flowPhaseToEbosCompIdx; using Base::flowPhaseToEbosCompIdx;
using Base::flowPhaseToEbosPhaseIdx;
using Base::ebosCompIdxToFlowCompIdx; using Base::ebosCompIdxToFlowCompIdx;
using Base::getAllowCrossFlow; using Base::getAllowCrossFlow;
using Base::scalingFactor; using Base::scalingFactor;

4
opm/simulators/wells/MultisegmentWell_impl.hpp
View File

@@ -4336,7 +4336,7 @@ namespace Opm
// the reciprocal FVF. // the reciprocal FVF.
const auto connMob = const auto connMob =
mobility[ flowPhaseToEbosCompIdx(p) ].value() mobility[ flowPhaseToEbosCompIdx(p) ].value()
* fs.invB(p).value(); * fs.invB(flowPhaseToEbosPhaseIdx(p)).value();
connPI[p] = connPICalc(connMob); connPI[p] = connPICalc(connMob);
} }
@@ -4393,6 +4393,6 @@ namespace Opm
const auto zero = EvalWell { 0.0 }; const auto zero = EvalWell { 0.0 };
const auto mt = std::accumulate(mobility.begin(), mobility.end(), zero); const auto mt = std::accumulate(mobility.begin(), mobility.end(), zero);
connII[phase_pos] = connIICalc(mt.value() * fs.invB(phase_pos).value()); connII[phase_pos] = connIICalc(mt.value() * fs.invB(flowPhaseToEbosPhaseIdx(phase_pos)).value());
} }
} // namespace Opm } // namespace Opm