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

ECL transmissibility: implement region multipliers

Browse Source 
a.k.a. support the MULTREGT keyword. This has not really been tested
and may contain bugs...
This commit is contained in:
Andreas Lauser 2015-01-27 13:17:01 +01:00
parent 019c272a11
commit 0e973f7108
1 changed files with 32 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

33
applications/ebos/ecltransmissibility.hh
View File

@ -124,7 +124,9 @@ public:
} }
} }
Opm::EclipseStateConstPtr eclState = simulator_.gridManager().eclState(); const auto& gridManager = simulator_.gridManager();
Opm::EclipseStateConstPtr eclState = gridManager.eclState();
auto multipliers = eclState->getTransMult();
const std::vector<double>& multx = const std::vector<double>& multx =
eclState->getDoubleGridProperty("MULTX")->getData(); eclState->getDoubleGridProperty("MULTX")->getData();
const std::vector<double>& multy = const std::vector<double>& multy =
@ -180,6 +182,9 @@ public:
if (insideElemIdx > outsideElemIdx) if (insideElemIdx > outsideElemIdx)
continue; continue;
int cartesianElemIdxInside = gridManager.cartesianCellId(insideElemIdx);
int cartesianElemIdxOutside = gridManager.cartesianCellId(outsideElemIdx);
// local indices of the faces of the inside and // local indices of the faces of the inside and
// outside elements which contain the intersection // outside elements which contain the intersection
int insideFaceIdx = intersection.indexInInside(); int insideFaceIdx = intersection.indexInInside();
@ -224,6 +229,32 @@ public:
multy, multyMinus, multy, multyMinus,
multz, multzMinus); multz, multzMinus);
// apply the region multipliers (cf. the MULTREGT keyword)
Opm::FaceDir::DirEnum faceDir;
switch (insideFaceIdx) {
case 0:
case 1:
faceDir = Opm::FaceDir::XPlus;
break;
case 2:
case 3:
faceDir = Opm::FaceDir::YPlus;
break;
case 4:
case 5:
faceDir = Opm::FaceDir::ZPlus;
break;
default:
OPM_THROW(std::logic_error, "Could not determine a face direction");
}
trans *= multipliers->getRegionMultiplier(cartesianElemIdxInside,
cartesianElemIdxOutside,
faceDir);
trans_[isId_(insideElemIdx, outsideElemIdx)] = trans; trans_[isId_(insideElemIdx, outsideElemIdx)] = trans;
} }
} }