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

get bhp from thp without solving

Browse Source
This commit is contained in:
Tor Harald Sandve
2021-04-14 15:58:58 +02:00
parent 481c72a0fc
commit ebcaea3faa
1 changed files with 5 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -1547,17 +1547,12 @@ namespace Opm
}
case Well::ProducerCMode::THP:
{
well_state.thp()[well_index] = this->getTHPConstraint(summaryState);
gliftDebug(
"computing BHP from THP to update well state",
deferred_logger);
auto bhp = computeBhpAtThpLimitProd(well_state, ebos_simulator, summaryState, deferred_logger);
if (bhp) {
well_state.bhp()[well_index] = *bhp;
} else {
deferred_logger.warning("FAILURE_GETTING_CONVERGED_BHP",
"Failed to find BHP when switching to THP control for well " + name());
std::vector<double> rates(3, 0.0);
for (int p = 0; p<np; ++p) {
rates[p] = well_state.wellRates()[well_index*np + p];
}
double bhp = calculateBhpFromThp(well_state, rates, well, summaryState, deferred_logger);
well_state.bhp()[well_index] = bhp;
break;
}
case Well::ProducerCMode::GRUP: