diff --git a/dumux/geomechanics/poroelastic/localresidual.hh b/dumux/geomechanics/poroelastic/localresidual.hh
index 4a7f6544f8cdc53b51872ab425d8aa187c4bd9d6..6e0b8e1c27aa44a2ef0eceabd93f8fec81e1ae8b 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());