From 2ba619e3ccdbfbb63097eebe174786de43eecc6e Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Fri, 28 Feb 2020 18:15:33 +0100
Subject: [PATCH] [PL/RM] Fix forgotten dkrel/dSL derivative.
---
.../RichardsMechanics/RichardsMechanicsFEM-impl.h | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index c1368c241a0..4643675858b 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -733,11 +733,11 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
div_u_dot * N_p * w;
*/
- double const dk_rel_dS_l = 0;
- /*
- _process_data.flow_material->getRelativePermeabilityDerivative(
- t, x_position, -p_cap_ip, temperature, S_L);
- */
+ double const dk_rel_dS_l =
+ medium->property(MPL::PropertyType::relative_permeability)
+ .template dValue<double>(variables,
+ MPL::Variable::liquid_saturation,
+ x_position, t, dt);
typename ShapeMatricesTypeDisplacement::GlobalDimVectorType const
grad_p_cap = -dNdx_p * p_L;
local_Jac
--
GitLab