Skip to content
Snippets Groups Projects
Commit cbc9f712 authored by Steffen Beese's avatar Steffen Beese Committed by Dmitri Naumov
Browse files

[RM/DSM] bugfix for transport porosity computation in case of DSM

parent 0b204bc3
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+ 24
− 18
View file @ cbc9f712
...@@ -941,16 +941,19 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement, ...@@ -941,16 +941,19 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
if (medium->hasProperty(MPL::PropertyType::transport_porosity)) if (medium->hasProperty(MPL::PropertyType::transport_porosity))
{ {
variables_prev[static_cast<int>( if (!medium->hasProperty(MPL::PropertyType::saturation_micro))
MPL::Variable::transport_porosity)] = {
_ip_data[ip].transport_porosity_prev; variables_prev[static_cast<int>(
MPL::Variable::transport_porosity)] =
_ip_data[ip].transport_porosity_prev;
_ip_data[ip].transport_porosity = _ip_data[ip].transport_porosity =
medium->property(MPL::PropertyType::transport_porosity) medium->property(MPL::PropertyType::transport_porosity)
.template value<double>(variables, variables_prev, .template value<double>(variables, variables_prev,
x_position, t, dt); x_position, t, dt);
variables[static_cast<int>(MPL::Variable::transport_porosity)] = variables[static_cast<int>(MPL::Variable::transport_porosity)] =
_ip_data[ip].transport_porosity; _ip_data[ip].transport_porosity;
}
} }
else else
{ {
...@@ -1722,16 +1725,19 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement, ...@@ -1722,16 +1725,19 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
if (medium->hasProperty(MPL::PropertyType::transport_porosity)) if (medium->hasProperty(MPL::PropertyType::transport_porosity))
{ {
variables_prev[static_cast<int>( if (!medium->hasProperty(MPL::PropertyType::saturation_micro))
MPL::Variable::transport_porosity)] = {
_ip_data[ip].transport_porosity_prev; variables_prev[static_cast<int>(
MPL::Variable::transport_porosity)] =
_ip_data[ip].transport_porosity_prev;
_ip_data[ip].transport_porosity = _ip_data[ip].transport_porosity =
medium->property(MPL::PropertyType::transport_porosity) medium->property(MPL::PropertyType::transport_porosity)
.template value<double>(variables, variables_prev, .template value<double>(variables, variables_prev,
x_position, t, dt); x_position, t, dt);
variables[static_cast<int>(MPL::Variable::transport_porosity)] = variables[static_cast<int>(MPL::Variable::transport_porosity)] =
_ip_data[ip].transport_porosity; _ip_data[ip].transport_porosity;
}
} }
else else
{ {
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment