From d716f535135720a3be2a59e1567207f7d6c9ed60 Mon Sep 17 00:00:00 2001
From: heck <katharina.heck@iws.uni-stuttgart.de>
Date: Wed, 28 Apr 2021 16:53:20 +0200
Subject: [PATCH] [fix][richards] add extrusion factor in storage derivative
---
dumux/porousmediumflow/richards/localresidual.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dumux/porousmediumflow/richards/localresidual.hh b/dumux/porousmediumflow/richards/localresidual.hh
index dac5b3d933..6e58d391e2 100644
--- a/dumux/porousmediumflow/richards/localresidual.hh
+++ b/dumux/porousmediumflow/richards/localresidual.hh
@@ -196,7 +196,7 @@ public:
static_assert(!FluidSystem::isCompressible(0),
"richards/localresidual.hh: Analytic Jacobian only supports incompressible fluids!");
- const auto poreVolume = Extrusion::volume(scv)*curVolVars.porosity();
+ const auto poreVolume = Extrusion::volume(scv)*curVolVars.porosity()*curVolVars.extrusionFactor();
static const auto rho = curVolVars.density(0);
// partial derivative of storage term w.r.t. p_w
--
GitLab