From 0283ebb5f10e4eac7633c84a1aea620cdccadc99 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Thu, 7 Oct 2021 21:12:41 +0200
Subject: [PATCH] [PL/TH2M] Extend dfW_4_LWpG_a/dp_cap with k/mu derivatives.
---
ProcessLib/TH2M/TH2MFEM-impl.h | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index 11e98a7cdf5..fe85628cce3 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -662,7 +662,11 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
+ c.drho_W_LR_dp_GR * k_over_mu_L
// + rhoWLR * (dk_over_mu_L_dp_GR = 0)
;
- ip_cv.dfW_4_LWpG_a_dp_cap = -c.drho_W_LR_dp_LR * k_over_mu_L;
+ ip_cv.dfW_4_LWpG_a_dp_cap = c.drho_W_GR_dp_cap * k_over_mu_G +
+ ip_data.rhoWGR * ip_cv.dk_over_mu_G_dp_cap +
+ -c.drho_W_LR_dp_LR * k_over_mu_L +
+ ip_data.rhoWLR * ip_cv.dk_over_mu_L_dp_cap;
+
ip_cv.dfW_4_LWpG_a_dT =
c.drho_W_GR_dT * k_over_mu_G
//+ rhoWGR * (dk_over_mu_G_dT != 0 TODO for mu_G(T))
--
GitLab