From 3fa8abcf44b02fc71b27fd3127b743a6048c514b Mon Sep 17 00:00:00 2001
From: Boyan Meng <meng.boyan@ufz.de>
Date: Tue, 24 Sep 2019 17:06:02 +0200
Subject: [PATCH] read in the fluid properties: density, heat capacity and
Darcy velocity
---
.../HeatTransportBHELocalAssemblerSoil-impl.h | 16 ++++++++++++++++
1 file changed, 16 insertions(+)
diff --git a/ProcessLib/HeatTransportBHE/LocalAssemblers/HeatTransportBHELocalAssemblerSoil-impl.h b/ProcessLib/HeatTransportBHE/LocalAssemblers/HeatTransportBHELocalAssemblerSoil-impl.h
index a94a1a1c4a7..875430b76b3 100644
--- a/ProcessLib/HeatTransportBHE/LocalAssemblers/HeatTransportBHELocalAssemblerSoil-impl.h
+++ b/ProcessLib/HeatTransportBHE/LocalAssemblers/HeatTransportBHELocalAssemblerSoil-impl.h
@@ -89,6 +89,7 @@ void HeatTransportBHELocalAssemblerSoil<ShapeFunction, IntegrationMethod>::
auto const& medium = *_process_data.media_map->getMedium(_element_id);
auto const& solid_phase = medium.phase("Solid");
+ auto const& liquid_phase = medium.phase("AqueousLiquid");
MaterialPropertyLib::VariableArray vars;
@@ -114,11 +115,26 @@ void HeatTransportBHELocalAssemblerSoil<ShapeFunction, IntegrationMethod>::
MaterialPropertyLib::PropertyType::specific_heat_capacity)
.template value<double>(vars, pos, t);
+ auto const heat_capacity_f =
+ liquid_phase
+ .property(
+ MaterialPropertyLib::PropertyType::specific_heat_capacity)
+ .template value<double>(vars, pos, t);
+
auto const density_s =
solid_phase.property(MaterialPropertyLib::PropertyType::density)
.template value<double>(vars, pos, t);
+
+ auto const density_f =
+ liquid_phase.property(MaterialPropertyLib::PropertyType::density)
+ .template value<double>(vars, pos, t);
// for now only using the solid phase parameters
+ auto const velocity =
+ liquid_phase
+ .property(MaterialPropertyLib::PropertyType::phase_velocity)
+ .value(vars, pos, t);
+
// assemble Conductance matrix
local_K.noalias() += K_ip * k_s;
--
GitLab