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

Now computeTransportSource() ignores crossflow.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen
2012-04-24 15:18:19 +02:00
parent ce98195001
commit f0aed255b0
1 changed files with 5 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
opm/core/utility/miscUtilities.cpp
View File

@@ -327,10 +327,12 @@ namespace Opm
if (perf_rate > 0.0) {
// perf_rate is a total inflow rate, we want a water rate.
if (wells->type[w] != INJECTOR) {
std::cout << "**** Warning: crossflow in well with index " << w << std::endl;
std::cout << "**** Warning: crossflow in well with index " << w << " ignored." << std::endl;
perf_rate = 0.0;
} else {
ASSERT(std::fabs(zfrac[WATER] + zfrac[OIL] - 1.0) < 1e-6);
perf_rate *= zfrac[WATER];
}
ASSERT(std::fabs(zfrac[WATER] + zfrac[OIL] - 1.0) < 1e-6);
perf_rate *= zfrac[WATER];
}
transport_src[perf_cell] += perf_rate;
}