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

correcting the index of phase rate of the top segment

Browse Source 
When handling the rate control equations.
This commit is contained in:
Kai Bao 2015-10-08 15:14:23 +02:00
parent 70e1ccb736
commit 4fe0ae4a58
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/BlackoilMultiSegmentModel_impl.hpp
View File

@ -1135,7 +1135,7 @@ namespace Opm {
bhp_targets(w) = well_controls_iget_target(wc, current); bhp_targets(w) = well_controls_iget_target(wc, current);
rate_targets(w) = -1e100; rate_targets(w) = -1e100;
for (int p = 0; p < np; ++p) { for (int p = 0; p < np; ++p) {
rate_top_phase_elems.push_back(start_segment + p * nseg_total); rate_top_phase_elems.push_back(np * start_segment + p);
} }
} }
break; break;
@ -1152,7 +1152,7 @@ namespace Opm {
rate_well_elems.push_back(w); rate_well_elems.push_back(w);
rate_top_elems.push_back(start_segment); rate_top_elems.push_back(start_segment);
for (int p = 0; p < np; ++p) { for (int p = 0; p < np; ++p) {
rate_top_phase_elems.push_back(start_segment + p * nseg_total); rate_top_phase_elems.push_back(np * start_segment + p);
} }
// RESERVOIR and SURFACE rates look the same, from a // RESERVOIR and SURFACE rates look the same, from a
// high-level point of view, in the system of // high-level point of view, in the system of
@ -1218,6 +1218,7 @@ namespace Opm {
// 2, if the segment is not the top segment, then the pressure equation // 2, if the segment is not the top segment, then the pressure equation
const ADB bhp_residual = subset(state.segp, bhp_top_elems) - subset(bhp_targets, bhp_well_elems); const ADB bhp_residual = subset(state.segp, bhp_top_elems) - subset(bhp_targets, bhp_well_elems);
const ADB rate_residual = subset(rate_distr * subset(state.segqs, rate_top_phase_elems) - rate_targets, rate_well_elems); const ADB rate_residual = subset(rate_distr * subset(state.segqs, rate_top_phase_elems) - rate_targets, rate_well_elems);
// const ADB rate_residual = subset(rate_distr * state.segqs - rate_targets, rate_well_elems);
ADB others_residual = ADB::constant(V::Zero(nseg_total)); ADB others_residual = ADB::constant(V::Zero(nseg_total));
start_segment = 0; start_segment = 0;