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

keeping putting group control in.

Browse Source
This commit is contained in:
Kai Bao 2016-09-30 17:53:55 +02:00
parent f62d40fc05
commit e0dd9bfabd
2 changed files with 43 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/autodiff/StandardWells.hpp
View File

@ -71,7 +71,7 @@ namespace Opm {
Eigen::Dynamic,
Eigen::RowMajor>;
// --------- Public methods ---------
StandardWells(const Wells* wells_arg, const WellCollection* well_collection);
StandardWells(const Wells* wells_arg, WellCollection* well_collection);
void init(const BlackoilPropsAdInterface* fluid_arg,
const std::vector<bool>* active_arg,
@ -196,7 +196,7 @@ namespace Opm {
const Wells* wells_;
const WellOps wops_;
// TODO: It will probably need to be updated during running time.
const WellCollection* well_collection_;
WellCollection* well_collection_;
const BlackoilPropsAdInterface* fluid_;
const std::vector<bool>* active_;

42
opm/autodiff/StandardWells_impl.hpp
View File

@ -71,7 +71,7 @@ namespace Opm
StandardWells::StandardWells(const Wells* wells_arg, const WellCollection* well_collection)
StandardWells::StandardWells(const Wells* wells_arg, WellCollection* well_collection)
: wells_active_(wells_arg!=nullptr)
, wells_(wells_arg)
, wops_(wells_arg)
@ -711,6 +711,17 @@ namespace Opm
if( !localWellsActive() ) return ;
if (well_collection_->needUpdateWellTargets() ) {
well_collection_->updateWellTargets(xw);
std::cout << " well_collection_ need to update well targets " << std::endl;
std::cin.ignore();
for (size_t i = 0; i < well_collection_->numNode(); ++i) {
well_collection_->getNode(i)->setShouldUpdateWellTargets(false);
}
return;
}
std::string modestring[4] = { "BHP", "THP", "RESERVOIR_RATE", "SURFACE_RATE" };
// Find, for each well, if any constraints are broken. If so,
// switch control to first broken constraint.
const int np = wells().number_of_phases;
@ -757,8 +768,37 @@ namespace Opm
logger.wellSwitched(wells().name[w],
well_controls_iget_type(wc, current),
well_controls_iget_type(wc, ctrl_index));
std::ostringstream ss;
ss << "Switching control mode for well " << wells().name[w]
<< " from " << modestring[well_controls_iget_type(wc, current)]
<< " to " << modestring[well_controls_iget_type(wc, ctrl_index)] << std::endl;
const double* distr_before = well_controls_iget_distr(wc, current);
const double* distr_after = well_controls_iget_distr(wc, ctrl_index);
std::cout << "distr_before " << distr_before[0] << " " << distr_before[1] << " " << distr_before[2] << std::endl;
std::cout << "distr_after " << distr_after[0] << " " << distr_after[1] << " " << distr_after[2] << std::endl;
OpmLog::info(ss.str());
xw.currentControls()[w] = ctrl_index;
current = xw.currentControls()[w];
// not good practice, revising the interface for the better implementation later.
auto* well_node = dynamic_cast<Opm::WellNode *>(well_collection_->findNode(std::string(wells().name[w])));
// where to initialize this variable?
if (well_node->individualControl()) {
if (ctrl_index == well_node->groupControlIndex()) {
std::cout << "well " << well_node->name() << " is switching from individual control to group control " << std::endl;
well_node->setIndividualControl(false);
well_node->setShouldUpdateWellTargets(true);
}
} else {
if (ctrl_index != well_node->groupControlIndex()) {
well_node->setIndividualControl(true);
std::cout << "well " << well_node->name() << " is switching from group control to individual control " << std::endl;
well_node->setShouldUpdateWellTargets(true);
}
}
}
// Updating well state and primary variables.