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

Bugfix in impes transport for well perforations, function signature change.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2011-04-18 13:51:54 +02:00
parent 40c90929a6
commit e2e5ef8caf
2 changed files with 19 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
src/cfs_tpfa.c
View File

@ -1198,13 +1198,14 @@ cfs_tpfa_impes_maxtime(grid_t *G,
/* ---------------------------------------------------------------------- */ /* ---------------------------------------------------------------------- */
void void
cfs_tpfa_expl_mass_transport(grid_t *G, cfs_tpfa_expl_mass_transport(grid_t *G,
well_t *W, well_t *W,
int np, struct completion_data *wdata,
double dt, int np,
const double *porevol, double dt,
struct cfs_tpfa_data *h, const double *porevol,
double *surf_vol) struct cfs_tpfa_data *h,
double *surf_vol)
/* ---------------------------------------------------------------------- */ /* ---------------------------------------------------------------------- */
{ {
int c, i, f, c2, p, w; int c, i, f, c2, p, w;
@ -1251,10 +1252,10 @@ cfs_tpfa_expl_mass_transport(grid_t *G,
for (w = i = 0; w < W->number_of_wells; w++) { for (w = i = 0; w < W->number_of_wells; 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];
/* Get pressure difference between cell and well perforation */ /* Get difference between cell pressure and well bhp */
dp = wpress[i] - cpress[c]; dp = wpress[w] - cpress[c];
for (p = 0; p < np; p++) { for (p = 0; p < np; p++) {
dz = masstrans_p[i*np + p]* dp; dz = masstrans_p[i*np + p] * (dp + wdata->gpot[i*np + p]);
/* A positive dz means flow from the well perforation /* A positive dz means flow from the well perforation
i into the cell c */ i into the cell c */
surf_vol[c*np + p] += dz * dt / porevol[c]; surf_vol[c*np + p] += dz * dt / porevol[c];

15
src/cfs_tpfa.h
View File

@ -108,13 +108,14 @@ cfs_tpfa_impes_maxtime(grid_t *G,
const double *gravity); const double *gravity);
void void
cfs_tpfa_expl_mass_transport(grid_t *G, cfs_tpfa_expl_mass_transport(grid_t *G,
well_t *W, well_t *W,
int np, struct completion_data *wdata,
double dt, int np,
const double *porevol, double dt,
struct cfs_tpfa_data *h, const double *porevol,
double *surf_vol); struct cfs_tpfa_data *h,
double *surf_vol);
void void
cfs_tpfa_destroy(struct cfs_tpfa_data *h); cfs_tpfa_destroy(struct cfs_tpfa_data *h);