From f46ef54302bcabc960f671d222b103d94e1b2621 Mon Sep 17 00:00:00 2001
From: tnagel <nagelt@tcd.ie>
Date: Thu, 14 May 2020 07:32:01 +0200
Subject: [PATCH] Correction to bulk modulus computation
---
ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index 4799cfc6a32..84cba972d67 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -668,7 +668,8 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
t, x_position, dt, temperature);
auto const K_SR =
- _ip_data[ip].solid_material.getBulkModulus(t, x_position, &C_el);
+ _ip_data[ip].solid_material.getBulkModulus(t, x_position, &C_el) /
+ (1 - alpha);
auto const K_LR =
liquid_phase.property(MPL::PropertyType::bulk_modulus)
--
GitLab