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

Implement RESV constraint assembly.

Browse Source
This commit is contained in:
Bård Skaflestad 2012-03-20 23:17:52 +01:00
parent d6a6d24716
commit a3a9e12388
1 changed files with 33 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -235,6 +235,39 @@ assemble_rate_well(int nc, int w,
struct ifs_tpfa_data *h ) struct ifs_tpfa_data *h )
/* ---------------------------------------------------------------------- */ /* ---------------------------------------------------------------------- */
{ {
int c, i, wdof;
size_t jc, jw;
double trans, resv;
struct WellControls *ctrls;
ctrls = W->ctrls[ w ];
wdof = nc + w;
resv = ctrls->target[ ctrls->current ];
for (i = W->well_connpos[w]; i < W->well_connpos[w + 1]; i++) {
c = W->well_cells [ i ];
trans = mt[ c ] * W->WI[ i ];
/* c->w connection */
jc = csrmatrix_elm_index(c, c , h->A);
jw = csrmatrix_elm_index(c, wdof, h->A);
h->A->sa[ jc ] += trans;
h->A->sa[ jw ] -= trans;
h->b [ c ] += trans * wdp[ i ];
/* w->c connection */
jc = csrmatrix_elm_index(wdof, c , h->A);
jw = csrmatrix_elm_index(wdof, wdof, h->A);
h->A->sa[ jc ] -= trans;
h->A->sa[ jw ] += trans;
h->b [ wdof ] -= trans * wdp[ i ];
}
h->b[ wdof ] += resv;
} }