OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

effectiveRadius: Fix numerator formula

Browse Source 
We must take the square root of the permeability ratios, not use them
as is.  This was lost in the refactorisation commit 96cf137.

Pointy hat: @bska.
This commit is contained in:
Bård Skaflestad 2014-09-02 15:12:50 +02:00
parent a8bce79f15
commit 9591aec757
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/core/wells/WellsManager.cpp
View File

@ -241,7 +241,8 @@ namespace WellsManagerDetail
const double D0_sq = D[0] * D[0]; const double D0_sq = D[0] * D[0];
const double D1_sq = D[1] * D[1]; const double D1_sq = D[1] * D[1];
const double num = std::sqrt(K10*D0_sq + K01*D1_sq); const double num = std::sqrt((std::sqrt(K10) * D0_sq) +
(std::sqrt(K01) * D1_sq));
const double den = std::pow(K01, 0.25) + std::pow(K10, 0.25); const double den = std::pow(K01, 0.25) + std::pow(K10, 0.25);
// Note: Analytic constant 0.28 derived for infintely sized // Note: Analytic constant 0.28 derived for infintely sized