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

Correct blunders concerning presence of prescribed pressure values.

Browse Source 
Specifically, split detection of prescribed reservoir pressure values
from prescribed well (i.e., BHP) pressure values.  Previously, we
would not even assemble any well contributions if there were any
prescribed reservoir pressure values.  Moreover, the return value from
assemble_cell_contrib() was exactly opposite of its intended value...
This commit is contained in:
Bård Skaflestad 2010-12-08 14:50:19 +01:00
parent 30087dd925
commit 2777a6d8d9
1 changed files with 8 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/cfs_tpfa.c
View File

@ -471,7 +471,7 @@ assemble_cell_contrib(grid_t *G,
const double *ctrans = h->pimpl->dd->ctrans;
is_neumann = 0;
is_neumann = 1;
for (c = i = 0; c < G->number_of_cells; c++) {
j1 = csrmatrix_elm_index(c, c, h->A);
@ -734,22 +734,24 @@ cfs_tpfa_assemble(grid_t *G,
struct cfs_tpfa_data *h)
/* ---------------------------------------------------------------------- */
{
int is_neumann;
int res_is_neumann, well_is_neumann;
csrmatrix_zero( h->A);
vector_zero (h->A->m, h->b);
compute_psys_contrib(G, cq, dt, trans, gravcap_f, cpress0, porevol, h);
is_neumann = assemble_cell_contrib(G, bc, src, h);
res_is_neumann = assemble_cell_contrib(G, bc, src, h);
if ((W != NULL) && (wctrl != NULL) &&
(WI != NULL) && (wdp != NULL)) {
is_neumann = is_neumann &&
assemble_well_contrib(G->number_of_cells, W, wctrl, WI, wdp, h);
well_is_neumann = assemble_well_contrib(G->number_of_cells,
W, wctrl, WI, wdp, h);
} else {
well_is_neumann = 1;
}
if (is_neumann) {
if (res_is_neumann && well_is_neumann) {
h->A->sa[0] *= 2;
}
}