Fix polymer in well model for producers.

This commit is contained in:
Tor Harald Sandve
2017-06-23 11:13:00 +02:00
parent e9a1aa2a83
commit 1bc2550541
2 changed files with 1 additions and 2 deletions

View File

@@ -835,7 +835,6 @@ namespace Opm {
const double tol_cnv = param_.tolerance_cnv_;
const double tol_wells = param_.tolerance_wells_;
const int nc = Opm::AutoDiffGrid::numCells(grid_);
const int np = numPhases();
const int numComp = numComponents();

View File

@@ -257,7 +257,7 @@ namespace Opm {
if (wells().type[w] == INJECTOR) {
cq_s_poly *= wpolymer(w);
} else {
cq_s_poly *= extendEval(intQuants.polymerConcentration());
cq_s_poly *= extendEval(intQuants.polymerConcentration() * intQuants.polymerViscosityCorrection());
}
if (!only_wells) {
for (int pvIdx = 0; pvIdx < numEq; ++pvIdx) {