From 95bb48d65e101ad1e63cf526eae4e3e0d62b9b5b Mon Sep 17 00:00:00 2001
From: Norbert Grunwald <Norbert.Grunwald@ufz.de>
Date: Wed, 9 Dec 2020 13:04:53 +0100
Subject: [PATCH] fixed bug in dValue computation
---
.../MPL/Properties/RelativePermeability/RelPermUdell.cpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdell.cpp b/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdell.cpp
index 3bbda249441..5e70d82b34a 100644
--- a/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdell.cpp
+++ b/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdell.cpp
@@ -97,7 +97,7 @@ PropertyDataType RelPermUdell::dValue(
auto const dk_rel_LRdse = 3. * s * s;
auto const dk_rel_LRdsL = dk_rel_LRdse * d_se_d_sL;
- auto const dk_rel_GRdse = 3. * (1. - s) * (1. - s);
+ auto const dk_rel_GRdse = -3. * (1. - s) * (1. - s);
auto const dk_rel_GRdsL = dk_rel_GRdse * d_se_d_sL;
return Eigen::Vector2d{dk_rel_LRdsL, dk_rel_GRdsL};
--
GitLab