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

Merge pull request #3873 from totto82/aprox_thp

Browse Source 
tighten approximated solution bhpFromThpLimit from 3 to 0.1 bar
This commit is contained in:
Tor Harald Sandve 2022-04-21 08:56:57 +02:00 committed by GitHub
commit 03c1884533
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/simulators/wells/WellInterfaceGeneric.cpp
View File

@ -565,9 +565,9 @@ bisectBracket(const std::function<double(const double)>& eq,
} }
} else { // eq_low * eq_high > 0.0 } else { // eq_low * eq_high > 0.0
// Still failed bracketing! // Still failed bracketing!
const double limit = 3.0 * unit::barsa; const double limit = 0.1 * unit::barsa;
if (std::min(abs_low, abs_high) < limit) { if (std::min(abs_low, abs_high) < limit) {
// Return the least bad solution if less off than 3 bar. // Return the least bad solution if less off than 0.1 bar.
deferred_logger.warning("FAILED_ROBUST_BHP_THP_SOLVE_BRACKETING_FAILURE", deferred_logger.warning("FAILED_ROBUST_BHP_THP_SOLVE_BRACKETING_FAILURE",
"Robust bhp(thp) not solved precisely for well " + this->name()); "Robust bhp(thp) not solved precisely for well " + this->name());
approximate_solution = abs_low < abs_high ? low : high; approximate_solution = abs_low < abs_high ? low : high;