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

fix model for efficient diffusion coefficient and add reference

Browse Source
This commit is contained in:
Tor Harald Sandve
2021-02-16 11:07:32 +01:00
parent 59707e0f64
commit 7226209791
2 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -274,16 +274,15 @@ protected:
if (FluidSystem::waterPhaseIdx == phaseIdx) { if (FluidSystem::waterPhaseIdx == phaseIdx) {
continue; continue;
} }
// TODO: let the problem do this (this is a constitutive
// relation of which the model should be free of from the // Based on Millington, R. J., & Quirk, J. P. (1961).
// abstraction POV!)
const Evaluation& base = const Evaluation& base =
Toolbox::max(0.0001, Toolbox::max(0.0001,
intQuants.porosity() intQuants.porosity()
* intQuants.fluidState().saturation(phaseIdx)); * intQuants.fluidState().saturation(phaseIdx));
tortuosity_[phaseIdx] = tortuosity_[phaseIdx] =
1.0 / (intQuants.porosity() * intQuants.porosity()) 1.0 / (intQuants.porosity() * intQuants.porosity())
* Toolbox::pow(base, 7.0/3.0); * Toolbox::pow(base, 10.0/3.0);
for (unsigned compIdx = 0; compIdx < numComponents; ++compIdx) { for (unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
diffusionCoefficient_[phaseIdx][compIdx] = diffusionCoefficient_[phaseIdx][compIdx] =

3
opm/models/common/diffusionmodule.hh
View File

@@ -253,13 +253,14 @@ protected:
// TODO: let the problem do this (this is a constitutive // TODO: let the problem do this (this is a constitutive
// relation of which the model should be free of from the // relation of which the model should be free of from the
// abstraction POV!) // abstraction POV!)
// Based on Millington, R. J., & Quirk, J. P. (1961).
const Evaluation& base = const Evaluation& base =
Toolbox::max(0.0001, Toolbox::max(0.0001,
intQuants.porosity() intQuants.porosity()
* intQuants.fluidState().saturation(phaseIdx)); * intQuants.fluidState().saturation(phaseIdx));
tortuosity_[phaseIdx] = tortuosity_[phaseIdx] =
1.0 / (intQuants.porosity() * intQuants.porosity()) 1.0 / (intQuants.porosity() * intQuants.porosity())
* Toolbox::pow(base, 7.0/3.0); * Toolbox::pow(base, 10.0/3.0);
for (unsigned compIdx = 0; compIdx < numComponents; ++compIdx) { for (unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
diffusionCoefficient_[phaseIdx][compIdx] = diffusionCoefficient_[phaseIdx][compIdx] =