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

Now accepts both BHP and SURFACE_RATE wells.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2013-05-22 00:09:57 +02:00
parent 9aea3b4589
commit edc08c0679
1 changed files with 19 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
opm/autodiff/ImpesTPFAAD.hpp
View File

@ -314,16 +314,30 @@ namespace Opm {
well_flow_residual_[phase] = well_rates - subset(qs, well_flow_res_phase_idx); well_flow_residual_[phase] = well_rates - subset(qs, well_flow_res_phase_idx);
} }
cell_residual_ = cell_residual_ + divcontrib_sum; cell_residual_ = cell_residual_ + divcontrib_sum;
// Assuming bhp wells. // Handling BHP and SURFACE_RATE wells.
V bhp_targets(nw,1); V bhp_targets(nw,1);
V rate_targets(nw,1);
M rate_distr(nw, np*nw);
for (int w = 0; w < nw; ++w) { for (int w = 0; w < nw; ++w) {
const WellControls* wc = wells_.ctrls[w]; const WellControls* wc = wells_.ctrls[w];
if (wc->type[wc->current] != BHP) { if (wc->type[wc->current] == BHP) {
THROW("Only BHP wells for now, please!"); bhp_targets[w] = wc->target[wc->current];
rate_targets[w] = -1e100;
} else if (wc->type[wc->current] == SURFACE_RATE) {
bhp_targets[w] = -1e100;
rate_targets[w] = wc->target[wc->current];
for (int phase = 0; phase < np; ++phase) {
rate_distr.insert(w, phase*nw + w) = wc->distr[phase];
}
} else {
THROW("Can only handle BHP and SURFACE_RATE type controls.");
} }
bhp_targets[w] = wc->target[wc->current];
} }
well_residual_ = bhp - bhp_targets; const ADB bhp_residual = bhp - bhp_targets;
const ADB rate_residual = rate_distr * qs - rate_targets;
// Choose bhp residual for positive bhp targets.
Selector<double> bhp_selector(bhp_targets);
well_residual_ = bhp_selector.select(bhp_residual, rate_residual);
ASSERT(np == 2); ASSERT(np == 2);
const ADB well_flow_res = vertcat(well_flow_residual_[0], well_flow_residual_[1]); const ADB well_flow_res = vertcat(well_flow_residual_[0], well_flow_residual_[1]);