From 1bc2550541590fb73bb80aa80731b7dcac5363ff Mon Sep 17 00:00:00 2001 From: Tor Harald Sandve Date: Fri, 23 Jun 2017 11:13:00 +0200 Subject: [PATCH] Fix polymer in well model for producers. --- opm/autodiff/BlackoilModelEbos.hpp | 1 - opm/autodiff/StandardWellsDense_impl.hpp | 2 +- 2 files changed, 1 insertion(+), 2 deletions(-) diff --git a/opm/autodiff/BlackoilModelEbos.hpp b/opm/autodiff/BlackoilModelEbos.hpp index 740e116fa..6c84f8d70 100644 --- a/opm/autodiff/BlackoilModelEbos.hpp +++ b/opm/autodiff/BlackoilModelEbos.hpp @@ -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(); diff --git a/opm/autodiff/StandardWellsDense_impl.hpp b/opm/autodiff/StandardWellsDense_impl.hpp index e651902de..c7f50afe8 100644 --- a/opm/autodiff/StandardWellsDense_impl.hpp +++ b/opm/autodiff/StandardWellsDense_impl.hpp @@ -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) {