From 9ff13b714450e9bb7181935967c71f7cc001e5d1 Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Tue, 20 Mar 2018 08:25:00 +0100
Subject: [PATCH] [navierstokesnc] Use an effective diffusion coefficient
---
dumux/discretization/staggered/freeflow/fickslaw.hh | 4 ++--
dumux/freeflow/navierstokesnc/volumevariables.hh | 8 ++++++++
2 files changed, 10 insertions(+), 2 deletions(-)
diff --git a/dumux/discretization/staggered/freeflow/fickslaw.hh b/dumux/discretization/staggered/freeflow/fickslaw.hh
index 80247a5d63..7589cc0537 100644
--- a/dumux/discretization/staggered/freeflow/fickslaw.hh
+++ b/dumux/discretization/staggered/freeflow/fickslaw.hh
@@ -161,7 +161,7 @@ public:
const auto& insideVolVars = elemVolVars[insideScv];
const Scalar insideDistance = (insideScv.dofPosition() - scvf.ipGlobal()).two_norm();
- const Scalar insideD = insideVolVars.diffusionCoefficient(phaseIdx, compIdx);
+ const Scalar insideD = insideVolVars.effectiveDiffusivity(phaseIdx, compIdx);
const Scalar ti = calculateOmega_(insideDistance, insideD, 1.0);
if(scvf.boundary())
@@ -171,7 +171,7 @@ public:
const auto& outsideScv = fvGeometry.scv(scvf.outsideScvIdx());
const auto& outsideVolVars = elemVolVars[scvf.outsideScvIdx()];
const Scalar outsideDistance = (outsideScv.dofPosition() - scvf.ipGlobal()).two_norm();
- const Scalar outsideD = outsideVolVars.diffusionCoefficient(phaseIdx, compIdx);
+ const Scalar outsideD = outsideVolVars.effectiveDiffusivity(phaseIdx, compIdx);
const Scalar tj = calculateOmega_(outsideDistance, outsideD, 1.0);
tij = scvf.area()*(ti * tj)/(ti + tj);
diff --git a/dumux/freeflow/navierstokesnc/volumevariables.hh b/dumux/freeflow/navierstokesnc/volumevariables.hh
index 101073af6f..0e47ecd39e 100644
--- a/dumux/freeflow/navierstokesnc/volumevariables.hh
+++ b/dumux/freeflow/navierstokesnc/volumevariables.hh
@@ -204,6 +204,14 @@ public:
DUNE_THROW(Dune::InvalidStateException, "Diffusion coefficient called for phaseIdx = compIdx");
}
+ /*!
+ * \brief Returns the effective diffusion coefficient \f$\mathrm{[m^2/s]}\f$
+ */
+ Scalar effectiveDiffusivity(int pIdx, int compIdx) const
+ {
+ return diffusionCoefficient(pIdx, compIdx);
+ }
+
protected:
Implementation &asImp_()
--
GitLab