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

Add gravity contributions in flux calculations.

Browse Source 
Credit:
  Jostein R. Natvig <jostein.r.natvig@sintef.no>
This commit is contained in:
Bård Skaflestad 2011-09-20 13:18:33 +02:00
parent 82a51172b7
commit 06cb0ab97f
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/ifs_tpfa.c
View File

@ -239,7 +239,8 @@ ifs_tpfa_press_flux(grid_t *G,
double *fflux) double *fflux)
/* ---------------------------------------------------------------------- */ /* ---------------------------------------------------------------------- */
{ {
int c1, c2, f; int c1, c2, f;
double dh;
/* Assign cell pressure directly from solution vector */ /* Assign cell pressure directly from solution vector */
memcpy(cpress, h->x, G->number_of_cells * sizeof *cpress); memcpy(cpress, h->x, G->number_of_cells * sizeof *cpress);
@ -249,7 +250,8 @@ ifs_tpfa_press_flux(grid_t *G,
c2 = G->face_cells[2*f + 1]; c2 = G->face_cells[2*f + 1];
if ((c1 >= 0) && (c2 >= 0)) { if ((c1 >= 0) && (c2 >= 0)) {
fflux[f] = trans[f] * (cpress[c1] - cpress[c2]); dh = cpress[c1] - cpress[c2] + h->pimpl->fgrav[f];
fflux[f] = trans[f] * dh;
} else { } else {
fflux[f] = 0.0; fflux[f] = 0.0;
} }