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 compressed BC representation during flux calculation.

Browse Source
This commit is contained in:
Bård Skaflestad 2012-03-13 23:28:43 +01:00
parent 11671f0cad
commit 81f3c31392
1 changed files with 15 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -321,7 +321,7 @@ ifs_tpfa_press_flux(struct UnstructuredGrid *G,
/* ---------------------------------------------------------------------- */
{
int c1, c2, f;
size_t i;
size_t i, j;
double dh;
/* Assign cell pressure directly from solution vector */
@ -342,20 +342,23 @@ ifs_tpfa_press_flux(struct UnstructuredGrid *G,
if ((F != NULL) && (F->bc != NULL)) {
for (i = 0; i < F->bc->nbc; i++) {
if (F->bc->type[ i ] == BC_PRESSURE) {
f = F->bc->face[ i ];
c1 = G->face_cells[2*f + 0];
c2 = G->face_cells[2*f + 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));
assert ((c1 < 0) ^ (c2 < 0));
if (c1 < 0) { /* Environment -> c2 */
dh = F->bc->value[ i ] - cpress[c2];
if (c1 < 0) { /* Environment -> c2 */
dh = F->bc->value[ i ] - cpress[c2];
}
else { /* c1 -> environment */
dh = cpress[c1] - F->bc->value[ i ];
}
fflux[f] = trans[f] * (dh + h->pimpl->fgrav[f]);
}
else { /* c1 -> environment */
dh = cpress[c1] - F->bc->value[ i ];
}
fflux[f] = trans[f] * (dh + h->pimpl->fgrav[f]);
}
/* Other boundary condtions currently not handled */