OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-14 01:13:52 -06:00
Code Issues Packages Projects Releases Wiki Activity

#3627 Ensure volume weights for cells in 2d element extrusion doesn't exceed volume of extrusion.

Browse Source
This commit is contained in:
Gaute Lindkvist 2018-11-08 12:44:38 +01:00
parent ad9f86a517
commit 6ed598e40b
1 changed files with 22 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
ApplicationCode/ProjectDataModel/Rim2dGridProjection.cpp
View File

@ -671,6 +671,9 @@ std::vector<std::pair<size_t, double>> Rim2dGridProjection::visibleCellsAndWeigh
std::vector<std::tuple<size_t, double, double>> matchingVisibleCellsWeightAndHeight;
std::array<cvf::Vec3d, 8> hexCorners;
double sumOverlapVolumes = 0.0;
double maxHeight = -std::numeric_limits<double>::infinity();
double minHeight = std::numeric_limits<double>::infinity();
for (size_t globalCellIdx : allCellIndices)
{
if ((*m_cellGridIdxVisibility)[globalCellIdx])
@ -689,10 +692,28 @@ std::vector<std::pair<size_t, double>> Rim2dGridProjection::visibleCellsAndWeigh
{
double height = overlapBBox.max().z();
matchingVisibleCellsWeightAndHeight.push_back(std::make_tuple(globalCellIdx, overlapVolume, height));
sumOverlapVolumes += overlapVolume;
maxHeight = std::max(maxHeight, height);
minHeight = std::min(minHeight, overlapBBox.min().z());
}
}
}
cvf::Vec3d improvedBottomSwCorner = bottomSWCorner; improvedBottomSwCorner.z() = minHeight;
cvf::Vec3d improvedTopNECorner = topNECorner; improvedTopNECorner.z() = maxHeight;
cvf::BoundingBox improvedBbox2dElement = cvf::BoundingBox(improvedBottomSwCorner, improvedTopNECorner);
cvf::Vec3d improvedBboxExtent = improvedBbox2dElement.extent();
double improvedBboxVolume = improvedBboxExtent.x() * improvedBboxExtent.y() * improvedBboxExtent.z();
double volAdjustmentFactor = 1.0;
if (sumOverlapVolumes > improvedBboxVolume)
{
// Total volume weights for 2d Element should never be larger than the volume of the extruded 2d element.
volAdjustmentFactor = improvedBboxVolume / sumOverlapVolumes;
}
std::vector<std::pair<size_t, double>> matchingVisibleCellsAndWeight;
if (!matchingVisibleCellsWeightAndHeight.empty())
{
@ -704,7 +725,7 @@ std::vector<std::pair<size_t, double>> Rim2dGridProjection::visibleCellsAndWeigh
for (const auto& visWeightHeight : matchingVisibleCellsWeightAndHeight)
{
matchingVisibleCellsAndWeight.push_back(std::make_pair(std::get<0>(visWeightHeight), std::get<1>(visWeightHeight)));
matchingVisibleCellsAndWeight.push_back(std::make_pair(std::get<0>(visWeightHeight), std::get<1>(visWeightHeight) * volAdjustmentFactor));
}
}