From 09431e03be51f559c3b165f306b01ec29bf6ac31 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Wed, 29 Sep 2021 21:25:41 +0200
Subject: [PATCH] [PL/TH2M] Extend dfW_4_LWpGa with k_over_mu_X'.
---
ProcessLib/TH2M/TH2MFEM-impl.h | 12 ++++++++++--
1 file changed, 10 insertions(+), 2 deletions(-)
diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index c8280cf7e1b..017ea2b8986 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -650,10 +650,18 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
ip_cv.ds_L_dp_cap * rho_W_LR_dot - s_L * c.drho_W_LR_dp_LR / dt;
ip_cv.dfW_3a_dT = s_G * c.drho_W_GR_dT / dt + s_L * c.drho_W_LR_dT / dt;
- ip_cv.dfW_4_LWpG_a_dp_GR = c.drho_W_GR_dp_GR * k_over_mu_G;
+ ip_cv.dfW_4_LWpG_a_dp_GR = c.drho_W_GR_dp_GR * k_over_mu_G
+ // + rhoWGR * (dk_over_mu_G_dp_GR = 0)
+ + 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_dT =
- c.drho_W_GR_dT * k_over_mu_G + c.drho_W_LR_dT * k_over_mu_L;
+ c.drho_W_GR_dT * k_over_mu_G
+ //+ rhoWGR * (dk_over_mu_G_dT != 0 TODO for mu_G(T))
+ + c.drho_W_LR_dT * k_over_mu_L
+ //+ rhoWLR * (dk_over_mu_L_dT != 0 TODO for mu_G(T))
+ ;
// TODO (naumov) for dxmW*/d* != 0
ip_cv.dfW_4_LWpG_d_dp_GR =
--
GitLab