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

Added well handling to cfs_tpfa_expl_mass_transport(), simplified interface.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2011-02-03 12:51:53 +01:00
parent a779a061cc
commit a0c54cb6de
3 changed files with 56 additions and 27 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
src/TPFACompressiblePressureSolver.hpp
View File

@ -344,7 +344,6 @@ public:
/// @brief
/// Explicit IMPES transport.
void explicitTransport(const double dt,
const double* cell_pressures,
double* cell_surfvols)
{
if (wells_.number_of_wells != 0) {
@ -352,13 +351,13 @@ public:
"This function does not work with wells yet.");
}
int np = 3; // Number of phases.
std::vector<double> masstrans_f(np*grid_.c_grid()->number_of_faces);
std::vector<double> gravtrans_f(np*grid_.c_grid()->number_of_faces);
cfs_tpfa_retrieve_masstrans(grid_.c_grid(), np, data_, &masstrans_f[0]);
cfs_tpfa_retrieve_gravtrans(grid_.c_grid(), np, data_, &gravtrans_f[0]);
cfs_tpfa_expl_mass_transport(grid_.c_grid(), np, dt, &porevol_[0],
&masstrans_f[0], &gravtrans_f[0],
cell_pressures, cell_surfvols);
well_t* wells = NULL;
if (wells_.number_of_wells != 0) {
wells = &wells_;
}
cfs_tpfa_expl_mass_transport(grid_.c_grid(), wells, np, dt, &porevol_[0],
data_, cell_surfvols);
}

43
src/cfs_tpfa.c
View File

@ -1051,24 +1051,36 @@ cfs_tpfa_retrieve_gravtrans(grid_t *G,
/* ---------------------------------------------------------------------- */
void
cfs_tpfa_expl_mass_transport(grid_t *G,
well_t *W,
int np,
double dt,
const double *porevol,
const double *masstrans_f,
const double *gravtrans_f,
const double *cpress,
struct cfs_tpfa_data *h,
double *surf_vol)
/* ---------------------------------------------------------------------- */
{
int c, i, f, c2, p;
double dp, dz;
int c, i, f, c2, p, w;
double dp, dz, gsgn;
const double *masstrans_f, *gravtrans_f, *masstrans_p;
const double *cpress, *wpress;
/* Set up convenience pointers */
masstrans_f = h->pimpl->masstrans_f;
gravtrans_f = h->pimpl->gravtrans_f;
masstrans_p = h->pimpl->masstrans_p;
cpress = h->x;
wpress = h->x + G->number_of_cells;
/* Transport through interiour faces */
for (c = i = 0; c < G->number_of_cells; c++) {
for (; i < G->cell_facepos[c + 1]; i++) {
f = G->cell_faces[i];
if ((c2 = G->face_cells[2*f + 0]) == c) {
gsgn = 1.0;
c2 = G->face_cells[2*f + 1];
} else {
gsgn = -1.0;
}
if (c2 >= 0) {
@ -1076,13 +1088,32 @@ cfs_tpfa_expl_mass_transport(grid_t *G,
for (p = 0; p < np; p++) {
dz = masstrans_f[f*np + p] * dp;
dz += gravtrans_f[f*np + p] ;
dz += gravtrans_f[f*np + p] * gsgn;
/* A positive dz means flow from the cell c into
the cell c2. */
surf_vol[c*np + p] -= dz * dt / porevol[c];
}
}
}
}
/* Transport through well perforations */
if (W != NULL) {
for (w = i = 0; w < W->number_of_wells; w++) {
for (; i < W->well_connpos[w + 1]; i++) {
c = W->well_cells[i];
/* Get pressure difference between cell and well perforation */
dp = wpress[i] - cpress[c];
for (p = 0; p < np; p++) {
dz = masstrans_p[i*np + p]* dp;
/* A positive dz means flow from the well perforation
i into the cell c */
surf_vol[c*np + p] += dz * dt / porevol[c];
}
}
}
}
}

5
src/cfs_tpfa.h
View File

@ -99,12 +99,11 @@ cfs_tpfa_retrieve_gravtrans(grid_t *G,
void
cfs_tpfa_expl_mass_transport(grid_t *G,
well_t *W,
int np,
double dt,
const double *porevol,
const double *masstrans_f,
const double *gravtrans_f,
const double *cpress,
struct cfs_tpfa_data *h,
double *surf_vol);
void