OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-28 18:21:00 -06:00
Code Issues Packages Projects Releases Wiki Activity

iterate in computeBhpAtThpLimitProdWithAlq if no solution or potentials are negative

Browse Source
This commit is contained in:
Tor Harald Sandve 2022-02-09 10:59:24 +01:00
parent ec08f80405
commit e0573e99b1
1 changed files with 20 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -2328,11 +2328,30 @@ namespace Opm
double pressure_cell = this->getPerfCellPressure(fs).value(); double pressure_cell = this->getPerfCellPressure(fs).value();
max_pressure = std::max(max_pressure, pressure_cell); max_pressure = std::max(max_pressure, pressure_cell);
} }
return this->StandardWellGeneric<Scalar>::computeBhpAtThpLimitProdWithAlq(frates, auto bhpAtLimit = this->StandardWellGeneric<Scalar>::computeBhpAtThpLimitProdWithAlq(frates,
summary_state, summary_state,
deferred_logger, deferred_logger,
max_pressure, max_pressure,
alq_value); alq_value);
auto v = frates(*bhpAtLimit);
if(bhpAtLimit && std::all_of(v.cbegin(), v.cend(), [](double i){ return i <= 0; }))
return bhpAtLimit;
auto fratesIter = [this, &ebos_simulator, &deferred_logger](const double bhp) {
// Solver the well iterations to see if we are
// able to get a solution with an update
// solution
std::vector<double> rates(3);
computeWellRatesWithBhpIterations(ebos_simulator, bhp, rates, deferred_logger);
return rates;
};
return this->StandardWellGeneric<Scalar>::computeBhpAtThpLimitProdWithAlq(fratesIter,
summary_state,
deferred_logger,
max_pressure,
alq_value);
} }