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

EclTransmissibility: reduce division-by-zero cutoff to 10^-30

Browse Source 
transmissibilities of 10^-20 can occur in real decks. (i.e., Norne;
although the face which belongs to the transmissibility is tiny it
nonetheless affects convergence for some reason.)
This commit is contained in:
Andreas Lauser 2016-08-02 17:15:02 +02:00
parent 64a2c9df51
commit 22ecf54d41
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
applications/ebos/ecltransmissibility.hh
View File

@ -206,7 +206,7 @@ public:
// convert half transmissibilities to full face // convert half transmissibilities to full face
// transmissibilities using the harmonic mean // transmissibilities using the harmonic mean
Scalar trans; Scalar trans;
if (std::abs(halfTrans1) < 1e-20 || std::abs(halfTrans2) < 1e-20) if (std::abs(halfTrans1) < 1e-30 || std::abs(halfTrans2) < 1e-30)
// avoid division by zero // avoid division by zero
trans = 0.0; trans = 0.0;
else else