From 2c2237b2a4e8e22a0bb7d72adaccdd70f3be5d1c Mon Sep 17 00:00:00 2001
From: Jan Thiedau <jan.thiedau@bgr.de>
Date: Thu, 18 Jul 2019 13:39:10 +0200
Subject: [PATCH] RM: Add missing saturation in bulk density computation
---
ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index f64ec015285..6cfa2b4c248 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -243,7 +243,7 @@ void RichardsMechanicsLocalAssembler<
displacement_index, displacement_index)
.noalias() += B.transpose() * C * B * w;
- double const rho = rho_SR * (1 - porosity) + porosity * rho_LR;
+ double const rho = rho_SR * (1 - porosity) + S_L * porosity * rho_LR;
rhs.template segment<displacement_size>(displacement_index).noalias() +=
N_u_op.transpose() * rho * b * w;
@@ -452,7 +452,7 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
displacement_index, displacement_index)
.noalias() += B.transpose() * C * B * w;
- double const rho = rho_SR * (1 - porosity) + porosity * rho_LR;
+ double const rho = rho_SR * (1 - porosity) + S_L * porosity * rho_LR;
local_rhs.template segment<displacement_size>(displacement_index)
.noalias() -=
(B.transpose() * sigma_eff - N_u_op.transpose() * rho * b) * w;
--
GitLab