From 7572ff8b15a2d1f935f1a14ba68bffe483876681 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Mon, 11 Oct 2021 17:33:31 +0200
Subject: [PATCH] [PL/TH2M] Correct sign of LWpG Jacobian contributions.
---
ProcessLib/TH2M/TH2MFEM-impl.h | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index 23d6388e81a..41ca5c41ece 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -1565,18 +1565,18 @@ void TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
LWpG.noalias() += gradNpT * (advection_W + diffusion_W_p) * gradNp * w;
// fW_4 LWpG' parts; LWpG = \int grad (a + d) grad
- local_Jac.template block<W_size, C_size>(W_index, C_index).noalias() -=
+ local_Jac.template block<W_size, C_size>(W_index, C_index).noalias() +=
gradNpT * (ip_cv.dfW_4_LWpG_a_dp_GR + ip_cv.dfW_4_LWpG_d_dp_GR) *
gradpGR * Np * w;
- local_Jac.template block<W_size, W_size>(W_index, W_index).noalias() -=
+ local_Jac.template block<W_size, W_size>(W_index, W_index).noalias() +=
gradNpT * (ip_cv.dfW_4_LWpG_a_dp_cap + ip_cv.dfW_4_LWpG_d_dp_cap) *
gradpGR * Np * w;
local_Jac
.template block<W_size, temperature_size>(W_index,
temperature_index)
- .noalias() -= gradNpT *
+ .noalias() += gradNpT *
(ip_cv.dfW_4_LWpG_a_dT + ip_cv.dfW_4_LWpG_d_dT) *
gradpGR * NT * w;
--
GitLab