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

fixing the fraction calculation when initializing perf rates

Browse Source
This commit is contained in:
Kai Bao 2023-06-01 22:43:06 +02:00
parent d1fd5539af
commit 5d950eea46
1 changed files with 7 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
opm/simulators/wells/StandardWellConnections.cpp
View File

@ -513,6 +513,9 @@ computeProperties(const WellState& well_state,
{ {
total_tw = comm.sum(total_tw); total_tw = comm.sum(total_tw);
} }
// for producers where all perforations have zero rates we
// approximate the perforation mixture ration using the (invB * mobility) ratio,
// and weight the perforation rates using the well transmissibility.
for (int perf = 0; perf < nperf; ++perf) { for (int perf = 0; perf < nperf; ++perf) {
const int cell_idx = well_.cells()[perf]; const int cell_idx = well_.cells()[perf];
const double well_tw_fraction = well_.wellIndex()[perf] / total_tw; const double well_tw_fraction = well_.wellIndex()[perf] / total_tw;
@ -526,10 +529,12 @@ computeProperties(const WellState& well_state,
} }
for (int p = 0; p < np; ++p) { for (int p = 0; p < np; ++p) {
int ebosPhaseIdx = well_.flowPhaseToEbosPhaseIdx(p); int ebosPhaseIdx = well_.flowPhaseToEbosPhaseIdx(p);
perfRates[perf * well_.numComponents() + p] = well_tw_fraction * mobility(cell_idx, ebosPhaseIdx) / total_mobility; perfRates[perf * well_.numComponents() + p] = well_tw_fraction *
invB(cell_idx, ebosPhaseIdx) * mobility(cell_idx, ebosPhaseIdx) / total_mobility;
} }
if constexpr (Indices::enableSolvent) { if constexpr (Indices::enableSolvent) {
perfRates[perf * well_.numComponents() + Indices::contiSolventEqIdx] = well_tw_fraction * solventInverseFormationVolumeFactor(cell_idx) / total_mobility; perfRates[perf * well_.numComponents() + Indices::contiSolventEqIdx] = well_tw_fraction *
solventInverseFormationVolumeFactor(cell_idx) * solventMobility(cell_idx) / total_mobility;
} }
} }
} }