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

Minor simplification

Browse Source 
Simplify code using the unknown facedir value.
This commit is contained in:
Håkon Hægland
2022-09-14 15:33:59 +02:00
parent 7d0265dd44
commit 9ed974fda3
1 changed files with 4 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
ebos/eclfluxmodule.hh
View File

@@ -286,21 +286,17 @@ public:
const Evaluation& transMult = up.rockCompTransMultiplier();
const auto& materialLawManager = problem.materialLawManager();
const Evaluation *mob = nullptr;
FaceDir::DirEnum facedir = FaceDir::DirEnum::Unknown;
if (materialLawManager->hasDirectionalRelperms()) {
auto facedir = scvf.dirId(); // direction (X, Y, or Z) of the face
mob = &(up.mobility(phaseIdx, facedir));
}
else {
mob = &(up.mobility(phaseIdx));
facedir = scvf.dirId(); // direction (X, Y, or Z) of the face
}
if (upwindIsInterior)
volumeFlux[phaseIdx] =
pressureDifferences[phaseIdx]*(*mob)*transMult*(-trans/faceArea);
pressureDifferences[phaseIdx]*up.mobility(phaseIdx, facedir)*transMult*(-trans/faceArea);
else
volumeFlux[phaseIdx] =
pressureDifferences[phaseIdx]*
(Toolbox::value(*mob)*Toolbox::value(transMult)*(-trans/faceArea));
(Toolbox::value(up.mobility(phaseIdx, facedir))*Toolbox::value(transMult)*(-trans/faceArea));
}
}