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

Remove unneeded term, fixing oil injecting behaviour.

Browse Source 
The 'comp_term' is supposed to be a total divergence term, which is supposed
to be zero for incompressible flow. It was added for improved robustness in
stagnant areas, but as implemented it would not be computed properly for
oil injection scenarios, due to the convention for two-phase transport
source terms (positive terms are inflow of first phase [water], negative
terms are total outflow).
This commit is contained in:
Atgeirr Flø Rasmussen 2013-10-04 13:16:02 +02:00
parent 13ddbad9a1
commit dcaac93e23
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
opm/core/transport/reorder/TransportSolverTwophaseReorder.cpp
View File

@ -145,15 +145,17 @@ namespace Opm
explicit Residual(const TransportSolverTwophaseReorder& tmodel, int cell_index)
: tm(tmodel)
{
// Initialize [in|out]flux to include effect of transport source terms.
cell = cell_index;
s0 = tm.saturation_[cell];
double src_flux = -tm.source_[cell];
bool src_is_inflow = src_flux < 0.0;
influx = src_is_inflow ? src_flux : 0.0;
outflux = !src_is_inflow ? src_flux : 0.0;
comp_term = tm.source_[cell]; // Note: this assumes that all source flux is water.
dtpv = tm.dt_/tm.porevolume_[cell];
// Compute fluxes over interior edges. Boundary flow is supposed to be
// included in the transport source term, along with well sources.
for (int i = tm.grid_.cell_facepos[cell]; i < tm.grid_.cell_facepos[cell+1]; ++i) {
int f = tm.grid_.cell_faces[i];
double flux;
@ -173,13 +175,13 @@ namespace Opm
} else {
outflux += flux;
}
comp_term -= flux;
}
}
}
double operator()(double s) const
{
return s - s0 + dtpv*(outflux*tm.fracFlow(s, cell) + influx + s*comp_term);
return s - s0 + dtpv*(outflux*tm.fracFlow(s, cell) + influx);
}
};