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 explicit flux BCs during flux calculation.

Browse Source
This commit is contained in:
Bård Skaflestad 2012-03-13 23:35:26 +01:00
parent 81f3c31392
commit 7e74da8021
1 changed files with 19 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
opm/core/pressure/tpfa/ifs_tpfa.c
View File

@ -322,7 +322,7 @@ ifs_tpfa_press_flux(struct UnstructuredGrid *G,
{ {
int c1, c2, f; int c1, c2, f;
size_t i, j; size_t i, j;
double dh; double dh, s;
/* 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);
@ -361,7 +361,24 @@ ifs_tpfa_press_flux(struct UnstructuredGrid *G,
} }
} }
/* Other boundary condtions currently not handled */ else if (F->bc->type[ i ] == BC_FLUX_TOTVOL) {
assert (F->bc->cond_pos[i+1] - F->bc->cond_pos[i] == 1);
for (j = F->bc->cond_pos[ i + 0 ];
j < F->bc->cond_pos[ i + 1 ]; j++) {
f = F->bc->face[ j ];
c1 = G->face_cells[2*f + 0];
c2 = G->face_cells[2*f + 1];
assert ((c1 < 0) ^ (c2 < 0));
/* BC flux is positive into reservoir. */
s = 2.0*(c1 < 0) - 1.0;
fflux[f] = s * F->bc->value[ i ];
}
}
} }
} }
} }