From a41c5269dba127ca36ed31226416951f5426ff2f Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Mon, 17 Feb 2020 20:06:08 +0100
Subject: [PATCH] [PL/RM] Fix permeability computation in v output.
Now the permeability computation in RM requires the porosity to be set
in the .
---
ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h | 3 +++
1 file changed, 3 insertions(+)
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index 05f99990569..4d1b32c3ad6 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -777,6 +777,9 @@ std::vector<double> const& RichardsMechanicsLocalAssembler<
.template value<double>(variables, x_position, t, dt);
variables[static_cast<int>(MPL::Variable::temperature)] = temperature;
+ variables[static_cast<int>(MPL::Variable::porosity)] =
+ _ip_data[ip].porosity;
+
auto const mu = liquid_phase.property(MPL::PropertyType::viscosity)
.template value<double>(variables, x_position, t, dt);
auto const rho_LR =
--
GitLab