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

Account for gravity in well connection flux calculation

Browse Source 
We previously ignored effects of gravity in the calculation of the well
connection fluxes (i.e., perforation fluxes).  This commit includes
those effects where appropriate.
This commit is contained in:
Bård Skaflestad
2012-11-01 13:34:27 +01:00
parent 7406847201
commit a0c8028ea8
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
opm/core/pressure/tpfa/cfs_tpfa_residual.c
View File

@@ -1090,7 +1090,7 @@ compute_wflux(int np ,
{ {
int w, c, i, p; int w, c, i, p;
double pw, dp, t; double pw, dp, t;
const double *pmob; const double *pmob, *wdp;
struct Wells *W; struct Wells *W;
struct CompletionData *cdata; struct CompletionData *cdata;
@@ -1101,13 +1101,14 @@ compute_wflux(int np ,
W = wells->W; W = wells->W;
cdata = wells->data; cdata = wells->data;
wdp = cdata->wdp;
for (w = i = 0; w < W->number_of_wells; w++) { for (w = i = 0; w < W->number_of_wells; w++) {
pw = wpress[w]; pw = wpress[w];
for (; i < W->well_connpos[w + 1]; i++) { for (; i < W->well_connpos[w + 1]; i++) {
c = W->well_cells[ i ]; c = W->well_cells[ i ];
dp = pw - cpress[c]; dp = pw + wdp[ i ] - cpress[c];
if (dp > 0) { pmob = cdata->phasemob + (i * np); } /* w->c */ if (dp > 0) { pmob = cdata->phasemob + (i * np); } /* w->c */
else { pmob = pmobc + (c * np); } /* c->w */ else { pmob = pmobc + (c * np); } /* c->w */