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

scale diffusion by saturation

Browse Source
This commit is contained in:
Tor Harald Sandve 2021-01-18 14:01:08 +01:00
parent 87abc2c36b
commit 0b0fae607f
1 changed files with 10 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -122,16 +122,20 @@ public:
continue;
}
// arithmetic mean of the phase's b factor
Evaluation bAvg = fluidStateI.invB(phaseIdx);
bAvg += Toolbox::value(fluidStateJ.invB(phaseIdx));
bAvg /= 2;
// arithmetic mean of the phase's b factor weighed by saturation
Evaluation bSAvg = fluidStateI.saturation(phaseIdx) * fluidStateI.invB(phaseIdx);
bSAvg += Toolbox::value(fluidStateJ.saturation(phaseIdx)) * Toolbox::value(fluidStateJ.invB(phaseIdx));
bSAvg /= 2;
// phase not present, skip
if(bSAvg < 1.0e-6)
continue;
// mass flux of solvent component (oil in oil or gas in gas)
unsigned solventCompIdx = FluidSystem::solventComponentIndex(phaseIdx);
unsigned activeSolventCompIdx = Indices::canonicalToActiveComponentIndex(solventCompIdx);
flux[conti0EqIdx + activeSolventCompIdx] +=
-bAvg
-bSAvg
* extQuants.moleFractionGradientNormal(phaseIdx, solventCompIdx)
* extQuants.effectiveDiffusionCoefficient(phaseIdx, solventCompIdx);
@ -139,7 +143,7 @@ public:
unsigned soluteCompIdx = FluidSystem::soluteComponentIndex(phaseIdx);
unsigned activeSoluteCompIdx = Indices::canonicalToActiveComponentIndex(soluteCompIdx);
flux[conti0EqIdx + activeSoluteCompIdx] +=
-bAvg
-bSAvg
* extQuants.moleFractionGradientNormal(phaseIdx, soluteCompIdx)
* extQuants.effectiveDiffusionCoefficient(phaseIdx, soluteCompIdx);
}