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

moved transCompFactor to intensive quantities

Browse Source
This commit is contained in:
hnil 2022-06-19 19:50:38 +02:00 committed by Atgeirr Flø Rasmussen
parent 07a7070dc0
commit d9c59efcf1
1 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
ebos/eclfluxmoduletpfa.hh
View File

@ -498,8 +498,9 @@ public:
// TODO: should the rock compaction transmissibility multiplier be upstreamed // TODO: should the rock compaction transmissibility multiplier be upstreamed
// or averaged? all fluids should see the same compaction?! // or averaged? all fluids should see the same compaction?!
//const auto& globalIndex = stencil.globalSpaceIndex(upstreamIdx); //const auto& globalIndex = stencil.globalSpaceIndex(upstreamIdx);
const Evaluation& transMult = const Evaluation& transMult = up.rockCompTransMultiplier();
problem.template rockCompTransMultiplier<Evaluation>(up, globalIndex); //const Evaluation& transMult =
// problem.template rockCompTransMultiplier<Evaluation>(up, globalIndex);
if (upIdx[phaseIdx] == interiorDofIdx) if (upIdx[phaseIdx] == interiorDofIdx)
volumeFlux[phaseIdx] = volumeFlux[phaseIdx] =
@ -593,8 +594,9 @@ protected:
// or averaged? all fluids should see the same compaction?! // or averaged? all fluids should see the same compaction?!
//const auto& globalIndex = stencil.globalSpaceIndex(upstreamIdx); //const auto& globalIndex = stencil.globalSpaceIndex(upstreamIdx);
//NB as long as this is upwinded it could be an intensive quantity //NB as long as this is upwinded it could be an intensive quantity
const Evaluation& transMult = const Evaluation& transMult = up.rockCompTransMultiplier();
problem.template rockCompTransMultiplier<Evaluation>(up, globalIndex); // const Evaluation& transMult =
// problem.template rockCompTransMultiplier<Evaluation>(up, globalIndex);
if (upIdx_[phaseIdx] == interiorDofIdx_) if (upIdx_[phaseIdx] == interiorDofIdx_)
volumeFlux_[phaseIdx] = volumeFlux_[phaseIdx] =
@ -690,7 +692,9 @@ protected:
const auto& up = elemCtx.intensiveQuantities(upstreamIdx, timeIdx); const auto& up = elemCtx.intensiveQuantities(upstreamIdx, timeIdx);
// deal with water induced rock compaction // deal with water induced rock compaction
transModified *= problem.template rockCompTransMultiplier<double>(up, stencil.globalSpaceIndex(upstreamIdx)); const double transMult = Toolbox::value(up.rockCompTransMultiplier());
transModified *= transMult;
//problem.template rockCompTransMultiplier<double>(up, stencil.globalSpaceIndex(upstreamIdx));
volumeFlux_[phaseIdx] = volumeFlux_[phaseIdx] =
pressureDifference_[phaseIdx]*up.mobility(phaseIdx)*(-transModified/faceArea); pressureDifference_[phaseIdx]*up.mobility(phaseIdx)*(-transModified/faceArea);