From d91f6fb5fefba2762ae7dc1c6acebaaf777645d9 Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Wed, 21 Jul 2021 13:51:20 +0200
Subject: [PATCH] [THM/localAssembler] Replaced the hard coded solid density
---
.../ThermoHydroMechanicsFEM-impl.h | 12 ++----------
1 file changed, 2 insertions(+), 10 deletions(-)
diff --git a/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsFEM-impl.h b/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsFEM-impl.h
index 637c510ff92..125cf46ae59 100644
--- a/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsFEM-impl.h
+++ b/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsFEM-impl.h
@@ -319,15 +319,6 @@ void ThermoHydroMechanicsLocalAssembler<ShapeFunctionDisplacement,
dthermal_strain =
solid_linear_thermal_expansion_coefficient * dT_int_pt;
- double const T_ref =
- _process_data.reference_temperature(t, x_position)[0];
- double const rho_s =
- solid_density *
- (1 -
- Invariants::trace(solid_linear_thermal_expansion_coefficient) *
- T_int_pt -
- T_ref);
-
auto const K_pT_thermal_osmosis =
(solid_phase.hasProperty(
MaterialPropertyLib::PropertyType::thermal_osmosis_coefficient)
@@ -362,7 +353,8 @@ void ThermoHydroMechanicsLocalAssembler<ShapeFunctionDisplacement,
displacement_index, displacement_index)
.noalias() += B.transpose() * C * B * w;
- auto const rho = rho_s * (1 - porosity) + porosity * fluid_density;
+ auto const rho =
+ solid_density * (1 - porosity) + porosity * fluid_density;
local_rhs.template segment<displacement_size>(displacement_index)
.noalias() -=
(B.transpose() * sigma_eff - N_u_op.transpose() * rho * b) * w;
--
GitLab