From a6719a03bdec8c6925454e4584a3020ea74ac7fe Mon Sep 17 00:00:00 2001
From: Mathis Kelm <mathis.kelm@iws.uni-stuttgart.de>
Date: Wed, 7 Jun 2023 16:44:20 +0200
Subject: [PATCH] [flux][staggered] ficks use average density at boundary
---
dumux/flux/staggered/freeflow/fickslaw.hh | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/dumux/flux/staggered/freeflow/fickslaw.hh b/dumux/flux/staggered/freeflow/fickslaw.hh
index 9f08dd723d..2ff4f468e9 100644
--- a/dumux/flux/staggered/freeflow/fickslaw.hh
+++ b/dumux/flux/staggered/freeflow/fickslaw.hh
@@ -115,7 +115,10 @@ public:
if (scvf.boundary())
{
- flux[compIdx] = insideDensity * insideDiffCoeff
+ // When the face lies on a boundary, use the average density
+ const Scalar outsideDensity = massOrMolarDensity(outsideVolVars, referenceSystem, phaseIdx);
+ const Scalar avgDensity = 0.5*(insideDensity + outsideDensity);
+ flux[compIdx] = avgDensity * insideDiffCoeff
* (massOrMoleFractionInside - massOrMoleFractionOutside) / insideDistance;
}
else
--
GitLab