From 71ac8676e5127c478af1cad118d0087e94c7c365 Mon Sep 17 00:00:00 2001 From: goncalvesmachadoc Date: Mon, 12 Sep 2022 17:12:44 +0200 Subject: [PATCH] corrrection --- opm/models/blackoil/blackoilbrinemodules.hh | 3 ++- opm/models/blackoil/blackoilnewtonmethod.hh | 2 +- 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/opm/models/blackoil/blackoilbrinemodules.hh b/opm/models/blackoil/blackoilbrinemodules.hh index 4b2098428..1cae8ba08 100644 --- a/opm/models/blackoil/blackoilbrinemodules.hh +++ b/opm/models/blackoil/blackoilbrinemodules.hh @@ -239,6 +239,7 @@ public: double saltDensity = 2170; // Solid salt density kg/m3 const LhsEval solidSalt = Toolbox::template decay(intQuants.porosity()) + / (1.0 - Toolbox::template decay(intQuants.saltSaturation()) + 1.e-8) * saltDensity * Toolbox::template decay(intQuants.saltSaturation()); @@ -459,7 +460,7 @@ public: const auto& permfactTable = BrineModule::permfactTable(elemCtx, dofIdx, timeIdx); - permFactor_ = permfactTable.eval(scalarValue(porosityFactor)); + permFactor_ = permfactTable.eval(porosityFactor); for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) { if (!FluidSystem::phaseIsActive(phaseIdx)) continue; diff --git a/opm/models/blackoil/blackoilnewtonmethod.hh b/opm/models/blackoil/blackoilnewtonmethod.hh index aa0a07774..e1325cc37 100644 --- a/opm/models/blackoil/blackoilnewtonmethod.hh +++ b/opm/models/blackoil/blackoilnewtonmethod.hh @@ -391,7 +391,7 @@ protected: nextValue[pvIdx] = std::max(nextValue[pvIdx], 0.0); // keep the salt saturation below upperlimit if ((enableSaltPrecipitation && currentValue.primaryVarsMeaningBrine() == PrimaryVariables::Sp)) - nextValue[pvIdx] = std::min(nextValue[pvIdx], 0.9); + nextValue[pvIdx] = std::min(nextValue[pvIdx], 1.0-1.e-8); } // keep the temperature within given values