diff --git a/opm/simulators/wells/BlackoilWellModelGeneric.cpp b/opm/simulators/wells/BlackoilWellModelGeneric.cpp index f269c2b76..3071849c6 100644 --- a/opm/simulators/wells/BlackoilWellModelGeneric.cpp +++ b/opm/simulators/wells/BlackoilWellModelGeneric.cpp @@ -1841,6 +1841,17 @@ updateFiltrationModelsPostStep(const double dt, } } +template +void BlackoilWellModelGeneric:: +updateInjMult(DeferredLogger& deferred_logger) +{ + for (const auto& well : this->well_container_generic_) { + if (well->isInjector() && well->wellEcl().getInjMultMode() != Well::InjMultMode::NONE) { + well->updateInjMult(this->prev_inj_multipliers_[well->name()], deferred_logger); + } + } +} + template void BlackoilWellModelGeneric:: updateFiltrationModelsPreStep(DeferredLogger& deferred_logger) @@ -1849,26 +1860,14 @@ updateFiltrationModelsPreStep(DeferredLogger& deferred_logger) if (well->isInjector()) { const auto it = filter_cake_.find(well->name()); if (it != filter_cake_.end()) { - well->updateFilterCakeMultipliers(it->second.multipliers()); it->second.updatePreStep(*well, deferred_logger); + well->updateFilterCakeMultipliers(it->second.multipliers()); } } } } -template -void BlackoilWellModelGeneric:: -updateInjMult(DeferredLogger& deferred_logger) -{ - for (const auto& well : this->well_container_generic_) { - if (well->isInjector() && well->wellEcl().getInjMultMode() != Well::InjMultMode::NONE) { - well->updateInjMult(this->prev_inj_multipliers_[well->name()], deferred_logger); - } - } -} - - template class BlackoilWellModelGeneric; #if FLOW_INSTANTIATE_FLOAT diff --git a/opm/simulators/wells/WellFilterCake.cpp b/opm/simulators/wells/WellFilterCake.cpp index d0d572f12..a5a71210e 100644 --- a/opm/simulators/wells/WellFilterCake.cpp +++ b/opm/simulators/wells/WellFilterCake.cpp @@ -92,6 +92,8 @@ applyCleaning(const WellInterfaceGeneric& well, filter_cake.sf_multiplier = 1.0; skin_factor_[perf] *= factor; + const auto denom = connection.ctfProperties().peaceman_denom; + inj_fc_multiplier_[perf] = denom / (denom + skin_factor_[perf]); const Scalar rw = connection.getFilterCakeRadius(); switch (filter_cake.geometry) { case FilterCake::FilterCakeGeometry::LINEAR: {