From ed965d366053193b22ff9fcbc2867f146ff70a15 Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Mon, 9 Mar 2020 14:41:02 +0100
Subject: [PATCH] [PL/LiquidFlow] Mv MPL porosity to MPL solid phase.
---
ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h | 2 +-
ProcessLib/LiquidFlow/LiquidFlowProcess.cpp | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
index da90cbdc39e..421bc2331ee 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
+++ b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
@@ -176,7 +176,7 @@ void LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>::
t, dt);
auto const porosity =
- medium->property(MaterialPropertyLib::PropertyType::porosity)
+ solid_phase.property(MaterialPropertyLib::PropertyType::porosity)
.template value<double>(vars, pos, t, dt);
auto const storage =
solid_phase.property(MaterialPropertyLib::PropertyType::storage)
diff --git a/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp b/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
index 15aa35995ed..d8afe847f51 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
+++ b/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
@@ -32,7 +32,6 @@ void checkMPLProperties(
std::array const requiredPropertyMedium = {
MaterialPropertyLib::reference_temperature,
- MaterialPropertyLib::PropertyType::porosity,
MaterialPropertyLib::PropertyType::permeability};
std::array const requiredPropertyLiquidPhase = {
@@ -40,6 +39,7 @@ void checkMPLProperties(
MaterialPropertyLib::PropertyType::density};
std::array const requiredPropertySolidPhase = {
+ MaterialPropertyLib::PropertyType::porosity,
MaterialPropertyLib::PropertyType::storage};
for (auto const& element : mesh.getElements())
--
GitLab