mirror of
https://github.com/OPM/ResInsight.git
synced 2026-10-24 22:27:55 -05:00
#5747 Completion Export: Avoid nan values for zero permeability
This commit is contained in:
1 parent
787999f936
commit
1272ee92d0
1 file changed
+7
@@ -35,6 +35,13 @@ double RigTransmissibilityEquations::wellBoreTransmissibilityComponent( double c
|
||||
{
|
||||
double K = cvf::Math::sqrt( permeabilityNormalDirection1 * permeabilityNormalDirection2 );
|
||||
|
||||
const double lowerLimit = 1.0e-9;
|
||||
if ( std::fabs( permeabilityNormalDirection1 * permeabilityNormalDirection2 ) < lowerLimit )
|
||||
{
|
||||
// Guard further computations to avoid nan values
|
||||
return 0.0;
|
||||
}
|
||||
|
||||
double nominator = cDarcyForRelevantUnit * 2 * cvf::PI_D * K * cellPerforationVectorComponent;
|
||||
|
||||
double peaceManRad = peacemanRadius( permeabilityNormalDirection1,
|
||||
|
||||
Reference in new issue
Block a user