diff --git a/opm/polymer/TransportModelPolymer.cpp b/opm/polymer/TransportModelPolymer.cpp index b5a5018a8..2b7af5bb5 100644 --- a/opm/polymer/TransportModelPolymer.cpp +++ b/opm/polymer/TransportModelPolymer.cpp @@ -570,7 +570,7 @@ namespace Opm { ResidualC res(*this, cell); const double a = 0.0; - const double b = polyprops_.cMax(); + const double b = polyprops_.cMax()*1.1; // Add 10% to account for possible non-monotonicity of hyperbolic system. int iters_used; // Check if current state is an acceptable solution.