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

ECL peaceman well: throw an exception if the BHP cannot determined

Browse Source 
this used to be an assert(), but it is not a programmer error but
something which can happen during valid simulation runs.
This commit is contained in:
Andreas Lauser 2015-12-31 13:20:39 +01:00
parent c44ad1ca3c
commit 2c8e1da7f3
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
applications/ebos/eclpeacemanwell.hh
View File

@ -1358,7 +1358,10 @@ protected:
Scalar fStar = wellResidual_(bhp + eps);
Scalar fPrime = (fStar - f)/eps;
assert(std::abs(fPrime) > 1e-20);
if (std::abs(fPrime) < 1e-20)
OPM_THROW(Opm::NumericalProblem,
"Cannot determine the bottom hole pressure for well " << name()
<< ": Derivative of the well residual is too small");
Scalar delta = f/fPrime;
bhp -= delta;