From 66366b9eeeb393eaf9e28cec364a8e381baa08f1 Mon Sep 17 00:00:00 2001
From: Mathis Kelm <mathis.kelm@iws.uni-stuttgart.de>
Date: Mon, 3 Mar 2025 11:25:40 +0100
Subject: [PATCH] [pm] Fix velocity reconstruction from neumann boundary
conditions
---
dumux/porousmediumflow/velocity.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dumux/porousmediumflow/velocity.hh b/dumux/porousmediumflow/velocity.hh
index 00884cbf12..bb2e5a6b30 100644
--- a/dumux/porousmediumflow/velocity.hh
+++ b/dumux/porousmediumflow/velocity.hh
@@ -335,7 +335,7 @@ public:
// We assume that the density at the face equals the one at the cell center and reconstruct a volume flux from the Neumann mass flux.
const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()];
const auto eqIdx = VolumeVariables::Indices::conti0EqIdx + phaseIdx;
- scvfFluxes[fIdx] += neumannFlux[eqIdx] / insideVolVars.density(phaseIdx) * scvf.area() * insideVolVars.extrusionFactor();
+ scvfFluxes[fIdx] += neumannFlux[eqIdx] / insideVolVars.density(phaseIdx) * scvf.area();
}
else
{
--
GitLab