From fc3df5de6caec5030459bb9186e3048327f24c27 Mon Sep 17 00:00:00 2001
From: Yue Wang <yue.wang@iws.uni-stuttgart.de>
Date: Thu, 16 Mar 2023 18:23:41 +0100
Subject: [PATCH] [poroelastic] fix the average density in localresidual
---
dumux/geomechanics/poroelastic/localresidual.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dumux/geomechanics/poroelastic/localresidual.hh b/dumux/geomechanics/poroelastic/localresidual.hh
index 4a7f6544f8..6e0b8e1c27 100644
--- a/dumux/geomechanics/poroelastic/localresidual.hh
+++ b/dumux/geomechanics/poroelastic/localresidual.hh
@@ -92,7 +92,7 @@ public:
const auto& vv = elemVolVars[scv];
const auto phi = vv.porosity();
const auto rhoFluid = problem.spatialParams().effectiveFluidDensity(element, scv);
- const auto rhoAverage = phi*rhoFluid + (1.0 - phi*vv.solidDensity());
+ const auto rhoAverage = phi*rhoFluid + (1.0 - phi)*vv.solidDensity();
// add body force
const auto& g = problem.spatialParams().gravity(scv.center());
--
GitLab