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

Merge pull request #4995 from aritorto/removeSmallNonCart

Browse Source 
Lgrs neighbors considered in removeSmallNonCart
This commit is contained in:
Markus Blatt 2023-11-14 10:00:34 +01:00 committed by GitHub
commit a70df8cf0a
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
ebos/ecltransmissibility_impl.hh
View File

@ -628,7 +628,9 @@ removeSmallNonCartesianTransmissibilities_()
int gc2 = std::max(cartMapper_.cartesianIndex(elements.first), cartMapper_.cartesianIndex(elements.second)); int gc2 = std::max(cartMapper_.cartesianIndex(elements.first), cartMapper_.cartesianIndex(elements.second));
// only adjust the NNCs // only adjust the NNCs
if (gc2 - gc1 == 1 || gc2 - gc1 == cartDims[0] || gc2 - gc1 == cartDims[0]*cartDims[1]) // When LGRs, all neighbors in the LGR are cartesian neighbours on the level grid representing the LGR.
// When elements on the leaf grid view have the same parent cell, gc1 and gc2 coincide.
if (gc2 - gc1 == 1 || gc2 - gc1 == cartDims[0] || gc2 - gc1 == cartDims[0]*cartDims[1] || gc2 - gc1 == 0)
continue; continue;
//remove transmissibilities less than the threshold (by default 1e-6 in the deck's unit system) //remove transmissibilities less than the threshold (by default 1e-6 in the deck's unit system)