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

Add bool return value to forceShutWellByNameIfPredictionMode().

Browse Source 
Returns true if the well requested was actually shut.
This commit is contained in:
Atgeirr Flø Rasmussen
2018-12-03 13:24:34 +01:00
parent 7b550ba06c
commit 6b896d178e
2 changed files with 12 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
opm/autodiff/BlackoilWellModel.hpp
View File

@@ -239,7 +239,8 @@ namespace Opm {
bool hasTHPConstraints() const;
/// Shut down any single well, but only if it is in prediction mode.
void forceShutWellByNameIfPredictionMode(const std::string& wellname, const double simulation_time);
/// Returns true if the well was actually found and shut.
bool forceShutWellByNameIfPredictionMode(const std::string& wellname, const double simulation_time);
protected:

13
opm/autodiff/BlackoilWellModel_impl.hpp
View File

@@ -133,30 +133,37 @@ namespace Opm {
/// Return true if any well has a THP constraint.
/// Return true if the well was found and shut.
template<typename TypeTag>
void
bool
BlackoilWellModel<TypeTag>::
forceShutWellByNameIfPredictionMode(const std::string& wellname,
const double simulation_time)
{
// Only add the well to the closed list on the
// process that owns it.
int well_was_shut = 0;
for (const auto& well : well_container_) {
if (well->name() == wellname) {
if (well->underPredictionMode()) {
wellTestState_.addClosedWell(wellname, WellTestConfig::Reason::PHYSICAL, simulation_time);
well_was_shut = 1;
}
break;
}
}
// Communicate across processes if a well was shut.
well_was_shut = ebosSimulator_.vanguard().grid().comm().max(well_was_shut);
// Only log a message on the output rank.
if (terminal_output_) {
if (terminal_output_ && well_was_shut) {
const std::string msg = "Well " + wellname
+ " will be shut because it cannot get converged.";
OpmLog::info(msg);
}
return (well_was_shut == 1);
}