From 6e60228d77c3f6b93c20ec238b3e8234d0fe8148 Mon Sep 17 00:00:00 2001
From: renchao_lu <renchao.lu@gmail.com>
Date: Sun, 2 May 2021 22:22:22 +0200
Subject: [PATCH] [PL/RCT] get material properties in local assembler.
---
.../RichardsComponentTransportFEM-impl.h | 9 ++++++++-
1 file changed, 8 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM-impl.h b/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM-impl.h
index a970aa99c95..aead9ac5e18 100644
--- a/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM-impl.h
+++ b/ProcessLib/RichardsComponentTransport/RichardsComponentTransportFEM-impl.h
@@ -10,6 +10,9 @@
#include "RichardsComponentTransportFEM.h"
+#include "MaterialLib/MPL/Medium.h"
+#include "MaterialLib/MPL/Property.h"
+
namespace ProcessLib
{
namespace RichardsComponentTransport
@@ -56,7 +59,7 @@ LocalAssemblerData<ShapeFunction, IntegrationMethod, GlobalDim>::
template <typename ShapeFunction, typename IntegrationMethod,
unsigned GlobalDim>
void LocalAssemblerData<ShapeFunction, IntegrationMethod, GlobalDim>::assemble(
- double const t, double const /*dt*/, std::vector<double> const& local_x,
+ double const t, double const dt, std::vector<double> const& local_x,
std::vector<double> const& /*local_xdot*/,
std::vector<double>& local_M_data, std::vector<double>& local_K_data,
std::vector<double>& local_b_data)
@@ -89,6 +92,10 @@ void LocalAssemblerData<ShapeFunction, IntegrationMethod, GlobalDim>::assemble(
GlobalDimMatrixType const& I(
GlobalDimMatrixType::Identity(GlobalDim, GlobalDim));
+ // Get material properties
+ auto const& medium = *_process_data.media_map->getMedium(_element_id);
+ auto const& phase = medium.phase("AqueousLiquid");
+
auto KCC = local_K.template block<concentration_size, concentration_size>(
concentration_index, concentration_index);
auto MCC = local_M.template block<concentration_size, concentration_size>(
--
GitLab