OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#2456 Non-Darcy Flow : Allow transmissibilities close to zero when computing d-factor

Browse Source
This commit is contained in:
Magne Sjaastad 2018-02-16 11:22:13 +01:00
parent 504e80534b
commit cf11eff319
1 changed files with 7 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
ApplicationCode/Commands/CompletionExportCommands/RicExportFractureCompletionsImpl.cpp
View File

@ -375,22 +375,18 @@ std::vector<RigCompletionData> RicExportFractureCompletionsImpl::generateCompdat
sumOfTransmissibilitiesInFracture += c.transmissibility(); sumOfTransmissibilitiesInFracture += c.transmissibility();
} }
if (fabs(sumOfTransmissibilitiesInFracture) > 1e-10)
{
for (auto& c : allCompletionsForOneFracture) for (auto& c : allCompletionsForOneFracture)
{ {
if (fabs(c.transmissibility()) > 1e-10) // NOTE : What is supposed to happen when the transmissibility is close to zero?
{
double dFactorForOneConnection = dFactorForFracture * sumOfTransmissibilitiesInFracture / c.transmissibility(); double dFactorForOneConnection = dFactorForFracture * sumOfTransmissibilitiesInFracture / c.transmissibility();
c.setDFactor(dFactorForOneConnection); c.setDFactor(dFactorForOneConnection);
}
double khForOneConnection = khForFracture * c.transmissibility() / sumOfTransmissibilitiesInFracture; double khForOneConnection = khForFracture * c.transmissibility() / sumOfTransmissibilitiesInFracture;
c.setKh(khForOneConnection); c.setKh(khForOneConnection);
} }
} }
}
std::copy(allCompletionsForOneFracture.begin(), allCompletionsForOneFracture.end(), std::back_inserter(fractureCompletions)); std::copy(allCompletionsForOneFracture.begin(), allCompletionsForOneFracture.end(), std::back_inserter(fractureCompletions));