OilfieldToolsNet/opm-simulators
4
0 Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-10-24 14:47:44 -05:00
Code Issues Packages Projects Releases Wiki Activity

Missing permFactTransMultiplier by MICP

Browse Source
This commit is contained in:
David Landa Marban committed 2025-04-09 14:58:57 +02:00
1 parent d201ecfc84
commit 7685984058
1 file changed
+9 -12
Download Patch File Download Diff File Expand all files Collapse all files
opm/simulators/flow/FlowProblemBlackoil.hpp
+9 -12
View File
@@ -688,23 +688,20 @@ public:
LhsEval permFactTransMultiplier(const IntensiveQuantities& intQuants, unsigned elementIdx) const
{
OPM_TIMEBLOCK_LOCAL(permFactTransMultiplier);
if (!enableSaltPrecipitation && !enableMICP)
return 1.0;
const auto& fs = intQuants.fluidState();
unsigned tableIdx = this->simulator().problem().satnumRegionIndex(elementIdx);
LhsEval porosityFactor = decay<LhsEval>(1. - fs.saltSaturation());
porosityFactor = min(porosityFactor, 1.0);
if constexpr (enableSaltPrecipitation) {
const auto& fs = intQuants.fluidState();
unsigned tableIdx = this->simulator().problem().satnumRegionIndex(elementIdx);
LhsEval porosityFactor = decay<LhsEval>(1. - fs.saltSaturation());
porosityFactor = min(porosityFactor, 1.0);
const auto& permfactTable = BrineModule::permfactTable(tableIdx);
return permfactTable.eval(porosityFactor, /*extrapolation=*/true);
}
if constexpr (enableMICP){
const auto& permfactTable = MICPModule::permfactTable(tableIdx);
return permfactTable.eval(porosityFactor, /*extrapolation=*/true);
else if constexpr (enableMICP) {
return intQuants.permFactor().value();
}
else {
return 1.0;
}
return 1.0;
}
// temporary solution to facilitate output of initial state from flow