From e8a7cfef9ce2180dcb9977f5c06769d3d44c5b1a Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Tue, 5 Apr 2022 17:29:31 +0200
Subject: [PATCH] [PL/TH2M] Move up evaluation of T_dot and T_prev
---
ProcessLib/TH2M/TH2MFEM-impl.h | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index ce77e8a5eb9..451253b5085 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -263,6 +263,8 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
ip_data.beta_T_SR = Invariants::trace(ip_data.alpha_T_SR);
double const T_dot = NT.dot(temperature_dot);
+ double const T_prev = T - T_dot * dt;
+
MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const
dthermal_strain = ip_data.alpha_T_SR * T_dot * dt;
@@ -312,7 +314,6 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
auto const rhoSR = rho_ref_SR;
#endif // NON_CONSTANT_SOLID_PHASE_VOLUME_FRACTION
- auto const T_prev = T - T_dot * dt;
ip_cv.C = ip_data.updateConstitutiveRelation(vars, t, pos, dt, T_prev);
// constitutive model object as specified in process creation
--
GitLab