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

Limit access to face direction.

Browse Source 
Only access the grid face direction information if directional relperms
have been enabled.
This commit is contained in:
Håkon Hægland 2022-09-06 19:27:20 +02:00
parent 8b1bcb36bc
commit 7d0265dd44
1 changed files with 11 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
ebos/eclfluxmodule.hh
View File

@ -230,7 +230,6 @@ public:
unsigned I = stencil.globalSpaceIndex(interiorDofIdx);
unsigned J = stencil.globalSpaceIndex(exteriorDofIdx);
auto facedir = scvf.dirId(); // direction (X, Y, or Z) of the face
Scalar trans = problem.transmissibility(elemCtx, interiorDofIdx, exteriorDofIdx);
Scalar faceArea = scvf.area();
Scalar thpres = problem.thresholdPressure(I, J);
@ -286,13 +285,22 @@ public:
// or averaged? all fluids should see the same compaction?!
const Evaluation& transMult = up.rockCompTransMultiplier();
const auto& materialLawManager = problem.materialLawManager();
const Evaluation *mob = nullptr;
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));
}
if (upwindIsInterior)
volumeFlux[phaseIdx] =
pressureDifferences[phaseIdx]*up.mobility(phaseIdx, facedir)*transMult*(-trans/faceArea);
pressureDifferences[phaseIdx]*(*mob)*transMult*(-trans/faceArea);
else
volumeFlux[phaseIdx] =
pressureDifferences[phaseIdx]*
(Toolbox::value(up.mobility(phaseIdx, facedir))*Toolbox::value(transMult)*(-trans/faceArea));
(Toolbox::value(*mob)*Toolbox::value(transMult)*(-trans/faceArea));
}
}