OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-11 15:15:42 -06:00
Code Issues Packages Projects Releases Wiki Activity

divide diffusivity by face area since opm-models expects pr area flux

Browse Source
This commit is contained in:
Tor Harald Sandve 2021-01-20 09:01:28 +01:00
parent 0b0fae607f
commit 59707e0f64
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
opm/models/blackoil/blackoildiffusionmodule.hh
View File

@ -138,7 +138,6 @@ public:
-bSAvg
* extQuants.moleFractionGradientNormal(phaseIdx, solventCompIdx)
* extQuants.effectiveDiffusionCoefficient(phaseIdx, solventCompIdx);
// mass flux of solute component (gas in oil or oil in gas)
unsigned soluteCompIdx = FluidSystem::soluteComponentIndex(phaseIdx);
unsigned activeSoluteCompIdx = Indices::canonicalToActiveComponentIndex(soluteCompIdx);
@ -398,7 +397,9 @@ protected:
const auto& intQuantsInside = elemCtx.intensiveQuantities(extQuants.interiorIndex(), timeIdx);
const auto& intQuantsOutside = elemCtx.intensiveQuantities(extQuants.exteriorIndex(), timeIdx);
Scalar diffusivity = elemCtx.problem().diffusivity(elemCtx, face.interiorIndex(), face.exteriorIndex());
const Scalar diffusivity = elemCtx.problem().diffusivity(elemCtx, face.interiorIndex(), face.exteriorIndex());
const Scalar faceArea = face.area();
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
if (!FluidSystem::phaseIsActive(phaseIdx)) {
continue;
@ -412,7 +413,7 @@ protected:
(intQuantsOutside.fluidState().moleFraction(phaseIdx, compIdx)
-
intQuantsInside.fluidState().moleFraction(phaseIdx, compIdx))
* diffusivity;
* diffusivity / faceArea; //opm-models expects pr area flux
Opm::Valgrind::CheckDefined(moleFractionGradientNormal_[phaseIdx][compIdx]);
// use the arithmetic average for the effective