OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #716 from goncalvesmachadoc/patch-2

Browse Source 
Bug fix on solid salt calculation
This commit is contained in:
Bård Skaflestad 2022-09-15 14:12:46 +02:00 committed by GitHub
commit a5e9af95fc
2 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/models/blackoil/blackoilbrinemodules.hh
View File

@ -239,6 +239,7 @@ public:
double saltDensity = 2170; // Solid salt density kg/m3
const LhsEval solidSalt =
Toolbox::template decay<LhsEval>(intQuants.porosity())
/ (1.0 - Toolbox::template decay<LhsEval>(intQuants.saltSaturation()) + 1.e-8)
* saltDensity
* Toolbox::template decay<LhsEval>(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;

2
opm/models/blackoil/blackoilnewtonmethod.hh
View File

@ -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