OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-27 17:50:59 -06:00
Code Issues Packages Projects Releases Wiki Activity

Compute perforation fluxes (no gravity adjustments).

Browse Source 
Requires additional parameter 'pmobc' (phase mobilities per cell) in
public function cfs_tpfa_res_flux().
This commit is contained in:
Bård Skaflestad 2011-12-20 02:13:50 +01:00
parent 87989fd12b
commit e7e5d94628
2 changed files with 40 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -987,6 +987,38 @@ compute_flux(grid_t *G,
} }
static void
compute_wflux(int np ,
struct cfs_tpfa_res_wells *W ,
const double *pmobc ,
const double *cpress,
const double *wpress,
double *wflux )
{
int w, c, i, p;
double pw, dp, t;
const double *pmob;
for (w = i = 0; w < W->conn->number_of_wells; w++) {
pw = wpress[w];
for (; i < W->conn->well_connpos[w + 1]; i++) {
c = W->conn->well_cells[ i ];
dp = pw - cpress[c];
if (dp > 0) { pmob = W->data->phasemob + (i * np); } /* w->c */
else { pmob = pmobc + (c * np); } /* c->w */
for (p = 0, t = 0.0; p < np; p++) {
t += pmob[ p ];
}
wflux[ i ] = W->data->WI[i] * t * dp;
}
}
}
static size_t static size_t
maxconn(grid_t *G) maxconn(grid_t *G)
{ {
@ -1141,6 +1173,7 @@ cfs_tpfa_res_flux(grid_t *G ,
struct cfs_tpfa_res_forces *forces , struct cfs_tpfa_res_forces *forces ,
int np , int np ,
const double *trans , const double *trans ,
const double *pmobc ,
const double *pmobf , const double *pmobf ,
const double *gravcap_f, const double *gravcap_f,
const double *cpress , const double *cpress ,
@ -1149,11 +1182,13 @@ cfs_tpfa_res_flux(grid_t *G ,
double *wflux ) double *wflux )
/* ---------------------------------------------------------------------- */ /* ---------------------------------------------------------------------- */
{ {
(void) forces;
(void) wpress;
(void) wflux;
compute_flux(G, np, trans, pmobf, gravcap_f, cpress, fflux); compute_flux(G, np, trans, pmobf, gravcap_f, cpress, fflux);
if ((forces != NULL) && (forces->W != NULL) &&
(wpress != NULL) && (wflux != NULL)) {
compute_wflux(np, forces->W, pmobc, cpress, wpress, wflux);
}
} }

1
opm/core/pressure/tpfa/cfs_tpfa_residual.h
View File

@ -79,6 +79,7 @@ cfs_tpfa_res_flux(grid_t *G ,
struct cfs_tpfa_res_forces *forces , struct cfs_tpfa_res_forces *forces ,
int np , int np ,
const double *trans , const double *trans ,
const double *pmobc ,
const double *pmobf , const double *pmobf ,
const double *gravcap_f, const double *gravcap_f,
const double *cpress , const double *cpress ,