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

Use well index when iterating over polymer perforation rates

Browse Source
This commit is contained in:
Joakim Hove 2021-05-21 07:48:50 +02:00
parent 5437cfa37a
commit f9e142dc8d
2 changed files with 6 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -794,7 +794,8 @@ namespace Opm
cq_s_poly *= extendEval(intQuants.polymerConcentration() * intQuants.polymerViscosityCorrection()); cq_s_poly *= extendEval(intQuants.polymerConcentration() * intQuants.polymerViscosityCorrection());
} }
// Note. Efficiency factor is handled in the output layer // Note. Efficiency factor is handled in the output layer
well_state.perfRatePolymer()[first_perf_ + perf] = cq_s_poly.value(); auto * perf_rate_polymer = &well_state.perfRatePolymer()[first_perf_];
perf_rate_polymer[perf] = cq_s_poly.value();
cq_s_poly *= well_efficiency_factor_; cq_s_poly *= well_efficiency_factor_;
connectionRates[perf][contiPolymerEqIdx] = Base::restrictEval(cq_s_poly); connectionRates[perf][contiPolymerEqIdx] = Base::restrictEval(cq_s_poly);

7
opm/simulators/wells/WellState.cpp
View File

@ -788,7 +788,8 @@ void WellState::reportConnections(data::Well& well,
comp.rates.set( pi [ i ], *(connPI + i) ); comp.rates.set( pi [ i ], *(connPI + i) );
} }
if ( pu.has_polymer ) { if ( pu.has_polymer ) {
comp.rates.set( rt::polymer, this->perfRatePolymer()[wt.second[1] + local_comp_index]); const auto * perf_polymer_rate = &this->perfRatePolymer()[wt.second[1]];
comp.rates.set( rt::polymer, perf_polymer_rate[local_comp_index]);
} }
if ( pu.has_brine ) { if ( pu.has_brine ) {
comp.rates.set( rt::brine, this->perfRateBrine()[wt.second[1] + local_comp_index]); comp.rates.set( rt::brine, this->perfRateBrine()[wt.second[1] + local_comp_index]);
@ -1017,8 +1018,8 @@ double WellState::solventWellRate(const int w) const
double WellState::polymerWellRate(const int w) const double WellState::polymerWellRate(const int w) const
{ {
return parallel_well_info_[w]->sumPerfValues(&perfRatePolymer_[0] + first_perf_index_[w], const auto * perf_rates_polymer = &perfRatePolymer_[first_perf_index_[w]];
&perfRatePolymer_[0] + first_perf_index_[w] + num_perf_[w]); return parallel_well_info_[w]->sumPerfValues(perf_rates_polymer, perf_rates_polymer + this->num_perf_[w]);
} }
double WellState::brineWellRate(const int w) const double WellState::brineWellRate(const int w) const