From cf11eff31993b8e008e1f53d419454fd1e39d876 Mon Sep 17 00:00:00 2001 From: Magne Sjaastad Date: Fri, 16 Feb 2018 11:22:13 +0100 Subject: [PATCH] #2456 Non-Darcy Flow : Allow transmissibilities close to zero when computing d-factor --- .../RicExportFractureCompletionsImpl.cpp | 18 +++++++----------- 1 file changed, 7 insertions(+), 11 deletions(-) diff --git a/ApplicationCode/Commands/CompletionExportCommands/RicExportFractureCompletionsImpl.cpp b/ApplicationCode/Commands/CompletionExportCommands/RicExportFractureCompletionsImpl.cpp index 16a2633334..8e3384a2a1 100644 --- a/ApplicationCode/Commands/CompletionExportCommands/RicExportFractureCompletionsImpl.cpp +++ b/ApplicationCode/Commands/CompletionExportCommands/RicExportFractureCompletionsImpl.cpp @@ -375,20 +375,16 @@ std::vector RicExportFractureCompletionsImpl::generateCompdat sumOfTransmissibilitiesInFracture += c.transmissibility(); } - if (fabs(sumOfTransmissibilitiesInFracture) > 1e-10) + for (auto& c : allCompletionsForOneFracture) { - for (auto& c : allCompletionsForOneFracture) - { - if (fabs(c.transmissibility()) > 1e-10) - { - double dFactorForOneConnection = dFactorForFracture * sumOfTransmissibilitiesInFracture / c.transmissibility(); - c.setDFactor(dFactorForOneConnection); - } + // NOTE : What is supposed to happen when the transmissibility is close to zero? - double khForOneConnection = khForFracture * c.transmissibility() / sumOfTransmissibilitiesInFracture; + double dFactorForOneConnection = dFactorForFracture * sumOfTransmissibilitiesInFracture / c.transmissibility(); + c.setDFactor(dFactorForOneConnection); - c.setKh(khForOneConnection); - } + double khForOneConnection = khForFracture * c.transmissibility() / sumOfTransmissibilitiesInFracture; + + c.setKh(khForOneConnection); } }