OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-24 16:30:02 -06:00
Code Issues Packages Projects Releases Wiki Activity

Make boundary condition data accessor function.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2022-09-16 16:17:37 +02:00
parent 68e1479caf
commit 5e08a9fc12
1 changed files with 37 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

37
ebos/eclproblem.hh
View File

@ -2061,6 +2061,43 @@ public:
return this->rockCompTransMultWc_[tableIdx].eval(effectiveOilPressure, SwDeltaMax, /*extrapolation=*/true);
}
std::pair<bool, RateVector> boundaryCondition(const unsigned int globalSpaceIdx, const int directionId)
{
if (!nonTrivialBoundaryConditions_ || directionId == -1) {
return std::make_pair(false, RateVector(0.0));
}
const std::vector<bool>* flag = nullptr;
const std::vector<RateVector>* massrate = nullptr;
switch (directionId) {
case 0:
flag = &freebcXMinus_;
massrate = &massratebcXMinus_;
break;
case 1:
flag = &freebcX_;
massrate = &massratebcX_;
break;
case 2:
flag = &freebcYMinus_;
massrate = &massratebcYMinus_;
break;
case 3:
flag = &freebcY_;
massrate = &massratebcY_;
break;
case 4:
flag = &freebcZMinus_;
massrate = &massratebcZMinus_;
break;
case 5:
flag = &freebcZ_;
massrate = &massratebcZ_;
break;
}
const bool free = (*flag)[globalSpaceIdx];
return { free, (*massrate)[globalSpaceIdx] };
}
private:
// update the parameters needed for DRSDT and DRVDT
void updateCompositionChangeLimits_()