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

Only add contribution from positive phase rates to the wellbore mixture

Browse Source 
Instead of adding contribution from all wells with positive total rates,
each phase in each well are checked seperatly.
This commit is contained in:
Tor Harald Sandve 2015-05-04 10:22:43 +02:00
parent a860361217
commit f1d7e0c7ea
1 changed files with 3 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
View File

@ -902,15 +902,7 @@ namespace detail {
selectInjectingPerforations[c] = 1;
else
selectProducingPerforations[c] = 1;
}
// select injectors
V selectInjector = V::Zero(nw);
for (int w = 0; w < nw; ++w) {
if (wells().type[w] == INJECTOR) {
selectInjector[w] = 1;
}
}
}
// HANDLE FLOW INTO WELLBORE
@ -950,8 +942,9 @@ namespace detail {
for (int phase = 0; phase < np; ++phase) {
const ADB& q_ps = wops_.p2w * cq_ps[phase];
const ADB& q_s = subset(state.qs, Span(nw, 1, phase*nw));
Selector<double> injectingPhase_selector(q_s.value(), Selector<double>::GreaterZero);
const int pos = pu.phase_pos[phase];
wbq[phase] = (selectInjector * compi.col(pos)) * q_s - q_ps;
wbq[phase] = (compi.col(pos) * injectingPhase_selector.select(q_s,ADB::constant(V::Zero(nw)))) - q_ps;
wbqt += wbq[phase];
}