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

implement MINPV treatment

Browse Source 
this is done my modifying the NTG values. it is just a
straight-forward port of the corresponding functionality from
opm-simulators.
This commit is contained in:
Andreas Lauser 2016-08-09 12:40:33 +02:00
parent 461f3bba1a
commit e8b460fe10
1 changed files with 41 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

42
applications/ebos/ecltransmissibility.hh
View File

@ -100,9 +100,11 @@ public:
const auto eclGrid = eclState->getInputGrid();
const auto& transMult = eclState->getTransMult();
const std::vector<double>& ntg =
std::vector<double> ntg =
eclState->get3DProperties().getDoubleGridProperty("NTG").getData();
applyMinPv_(ntg);
unsigned numElements = elementMapper.size();
// this code assumes that the DOFs are the elements. (i.e., an
@ -329,6 +331,44 @@ private:
}
}
void applyMinPv_(std::vector<double> &ntg)
{
const auto& model = simulator_.model();
const auto& gridManager = simulator_.gridManager();
const auto& eclState = gridManager.eclState();
const Opm::EclipseGridConstPtr eclGrid = gridManager.eclGrid();
const auto& porv = eclState->get3DProperties().getDoubleGridProperty("PORV").getData();
Scalar minpvValue = eclGrid->getMinpvValue();
if (minpvValue <= 0.0)
return;
int numCells = model.numGridDof();
for (int cellIdx = 0; cellIdx < numCells; ++cellIdx) {
const int nx = eclGrid->getNX();
const int ny = eclGrid->getNY();
const int cartesianCellIdx = gridManager.cartesianIndex(cellIdx);
const Scalar cellVolume = eclGrid->getCellVolume(cartesianCellIdx);
ntg[cartesianCellIdx] *= cellVolume;
Scalar totalCellVolume = cellVolume;
// Average properties as long as there exist cells above which exhibit a pore
// volume less than the MINPV threshold
int cartesianCellIdxAbove = cartesianCellIdx - nx*ny;
while (cartesianCellIdxAbove >= 0
&& porv[cartesianCellIdxAbove] > 0.0
&& porv[cartesianCellIdxAbove] < minpvValue)
{
const double cellVolumeAbove = eclGrid->getCellVolume(cartesianCellIdxAbove);
totalCellVolume += cellVolumeAbove;
ntg[cartesianCellIdx] += ntg[cartesianCellIdxAbove]*cellVolumeAbove;
cartesianCellIdxAbove -= nx*ny;
}
ntg[cartesianCellIdx] /= totalCellVolume;
}
}
void applyNtg_(Scalar &trans, unsigned faceIdx, unsigned cartElemIdx,
const std::vector<double>& ntg) const
{