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

Improved ntg lookup

Browse Source
This commit is contained in:
Antonella Ritorto 2023-11-09 17:36:31 +01:00
parent 9ae887681e
commit 1f4583f674
1 changed files with 3 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
ebos/ecltransmissibility_impl.hh
View File

@ -156,10 +156,10 @@ update(bool global, const std::function<unsigned int(unsigned int)>& map, const
const auto& comm = gridView_.comm(); const auto& comm = gridView_.comm();
ElementMapper elemMapper(gridView_, Dune::mcmgElementLayout()); ElementMapper elemMapper(gridView_, Dune::mcmgElementLayout());
// get the ntg values, the ntg values are modified for the cells merged with minpv
const std::vector<double>& ntg = eclState_.fieldProps().get_double("NTG");
const bool updateDiffusivity = eclState_.getSimulationConfig().isDiffusive() && enableDiffusivity_;
unsigned numElements = elemMapper.size(); unsigned numElements = elemMapper.size();
// get the ntg values, the ntg values are modified for the cells merged with minpv
const std::vector<double>& ntg = this->lookUpData_.assignFieldPropsDoubleOnLeaf(eclState_.fieldProps(), "NTG", numElements);
const bool updateDiffusivity = eclState_.getSimulationConfig().isDiffusive() && enableDiffusivity_;
if (map) if (map)
extractPermeability_(map); extractPermeability_(map);
@ -474,9 +474,6 @@ update(bool global, const std::function<unsigned int(unsigned int)>& map, const
axisCentroids), axisCentroids),
porosity_[outsideElemIdx]); porosity_[outsideElemIdx]);
auto coarseElemIdx = this->lookUpData_.template getFieldPropIdx<Grid>(elemIdx);
auto coarseOutsideElemIdx = this-> lookUpData_.template getFieldPropIdx<Grid>(outsideElemIdx);
applyNtg_(halfDiffusivity1, insideFaceIdx, coarseElemIdx, ntg); applyNtg_(halfDiffusivity1, insideFaceIdx, coarseElemIdx, ntg);
applyNtg_(halfDiffusivity2, outsideFaceIdx, coarseOutsideElemIdx, ntg); applyNtg_(halfDiffusivity2, outsideFaceIdx, coarseOutsideElemIdx, ntg);