OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-18 21:02:57 -06:00
Code Issues Packages Projects Releases Wiki Activity

Corrected bug (residual should not be divided by dt).

Browse Source
This commit is contained in:
Xavier Raynaud 2012-08-21 17:12:07 +02:00
parent 52ab67201d
commit 660f14df3b
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/core/pressure/tpfa/cfs_tpfa_residual.c
View File

@ -1242,7 +1242,7 @@ cfs_tpfa_res_comprock_assemble(
int c, w, well_is_neumann, rock_is_incomp;
size_t j;
double dpvdt;
double dpv;
const struct Wells* W;
/* Assemble usual system (without rock compressibility). */
@ -1270,13 +1270,13 @@ cfs_tpfa_res_comprock_assemble(
for (c = 0; c < G->number_of_cells; c++) {
j = csrmatrix_elm_index(c, c, h->J);
dpvdt = (porevol[c] - porevol0[c]) / dt;
if (dpvdt != 0.0 || rock_comp[c] != 0.0) {
dpv = (porevol[c] - porevol0[c]);
if (dpv != 0.0 || rock_comp[c] != 0.0) {
rock_is_incomp = 0;
}
h->J->sa[j] += porevol[c] * rock_comp[c] / dt;
h->F[c] += dpvdt;
h->J->sa[j] += porevol[c] * rock_comp[c];
h->F[c] += dpv;
}
/* Re-do the singularity-removing adjustment if necessary */