diff --git a/MaterialLib/MPL/Properties/PermeabilityOrthotropicPowerLaw.cpp b/MaterialLib/MPL/Properties/PermeabilityOrthotropicPowerLaw.cpp
index f03c06e630f6ab2049e3741d1bc3d3e1c0d3b01a..4974f4544d681b6dc045ceba40fa85ba4ebf98ed 100644
--- a/MaterialLib/MPL/Properties/PermeabilityOrthotropicPowerLaw.cpp
+++ b/MaterialLib/MPL/Properties/PermeabilityOrthotropicPowerLaw.cpp
@@ -37,20 +37,11 @@ PermeabilityOrthotropicPowerLaw<DisplacementDim>::
template <int DisplacementDim>
void PermeabilityOrthotropicPowerLaw<DisplacementDim>::checkScale() const
{
- if (!std::holds_alternative<Phase*>(scale_))
+ if (!std::holds_alternative<Medium*>(scale_))
{
OGS_FATAL(
"The property 'PermeabilityOrthotropicPowerLaw' is "
- "implemented on the 'phase' scales only.");
- }
-
- auto const phase = std::get<Phase*>(scale_);
- if (phase->name != "Solid")
- {
- OGS_FATAL(
- "The property 'PermeabilityOrthotropicPowerLaw' must be given in "
- "the 'Solid' phase, not in '{:s}' phase.",
- phase->name);
+ "implemented on the 'medium' scales only.");
}
}
template <int DisplacementDim>
@@ -64,13 +55,13 @@ PropertyDataType PermeabilityOrthotropicPowerLaw<DisplacementDim>::value(
// TODO (naumov) The phi0 must be evaluated once upon
// creation/initialization and be stored in a local state.
// For now assume porosity's initial value does not change with time.
- auto const phase = std::get<Phase*>(scale_);
+ auto const medium = std::get<Medium*>(scale_);
auto const phi_0 =
- phase->hasProperty(PropertyType::transport_porosity)
- ? phase->property(PropertyType::transport_porosity)
+ medium->hasProperty(PropertyType::transport_porosity)
+ ? medium->property(PropertyType::transport_porosity)
.template initialValue<double>(
pos, std::numeric_limits<double>::quiet_NaN())
- : phase->property(PropertyType::porosity)
+ : medium->property(PropertyType::porosity)
.template initialValue<double>(
pos, std::numeric_limits<double>::quiet_NaN());
diff --git a/MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp b/MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp
index da5296b23fc2337574cc8ed5e8af236dbe856f09..c58c63b5e4b8ee0bfcadc24c08c9cd03893bacb4 100644
--- a/MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp
+++ b/MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp
@@ -18,19 +18,11 @@ namespace MaterialPropertyLib
{
void PorosityFromMassBalance::checkScale() const
{
- if (!std::holds_alternative<Phase*>(scale_))
+ if (!std::holds_alternative<Medium*>(scale_))
{
OGS_FATAL(
"The property 'PorosityFromMassBalance' is "
- "implemented on the 'phase' scales only.");
- }
- auto const phase = std::get<Phase*>(scale_);
- if (phase->name != "Solid")
- {
- OGS_FATAL(
- "The property 'PorosityFromMassBalance' must be given in the "
- "'Solid' phase, not in '{:s}' phase.",
- phase->name);
+ "implemented on the 'medium' scales only.");
}
}
@@ -53,7 +45,7 @@ PropertyDataType PorosityFromMassBalance::value(
double const beta_SR = std::get<double>(
variable_array[static_cast<int>(Variable::grain_compressibility)]);
auto const alpha_b =
- std::get<Phase*>(scale_)
+ std::get<Medium*>(scale_)
->property(PropertyType::biot_coefficient)
.template value<double>(variable_array, pos, t, dt);
diff --git a/MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp b/MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp
index f53492e289c13fc3f6fe5537b575fe75cd78b531..09f6e2f68abac4d1d30523c004d0f83291bdbb08 100644
--- a/MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp
+++ b/MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp
@@ -18,19 +18,11 @@ namespace MaterialPropertyLib
{
void TransportPorosityFromMassBalance::checkScale() const
{
- if (!std::holds_alternative<Phase*>(scale_))
+ if (!std::holds_alternative<Medium*>(scale_))
{
OGS_FATAL(
"The property 'TransportPorosityFromMassBalance' is "
- "implemented on the 'phase' scales only.");
- }
- auto const phase = std::get<Phase*>(scale_);
- if (phase->name != "Solid")
- {
- OGS_FATAL(
- "The property 'TransportPorosityFromMassBalance' must be given in "
- "the 'Solid' phase, not in '{:s}' phase.",
- phase->name);
+ "implemented on the 'medium' scales only.");
}
}
@@ -43,7 +35,7 @@ PropertyDataType TransportPorosityFromMassBalance::value(
double const beta_SR = std::get<double>(
variable_array[static_cast<int>(Variable::grain_compressibility)]);
auto const alpha_b =
- std::get<Phase*>(scale_)
+ std::get<Medium*>(scale_)
->property(PropertyType::biot_coefficient)
.template value<double>(variable_array, pos, t, dt);
diff --git a/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp b/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
index 91b50a3c999efe77b43c64ccaf6ea0a9bee9f2df..9095d1b4cf933a28be4bd57c8b6e3ae8f8ed6064 100644
--- a/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
+++ b/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
@@ -35,11 +35,12 @@ void checkMPLProperties(
MaterialPropertyLib::reference_temperature,
MaterialPropertyLib::bishops_effective_stress,
MaterialPropertyLib::relative_permeability,
- MaterialPropertyLib::saturation};
+ MaterialPropertyLib::saturation,
+ MaterialPropertyLib::porosity, MaterialPropertyLib::biot_coefficient
+ };
std::array const required_liquid_properties = {
MaterialPropertyLib::viscosity, MaterialPropertyLib::density};
std::array const required_solid_properties = {
- MaterialPropertyLib::porosity, MaterialPropertyLib::biot_coefficient,
MaterialPropertyLib::density};
for (auto const& m : media)
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index 6823983d15c023a551ca82c7d988d50504a01dac..ad200b0860d8f48a5fb9c47af140c7d5c0908517 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -66,7 +66,6 @@ RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
e.getID());
auto const& medium = _process_data.media_map->getMedium(_element.getID());
- auto const& solid_phase = medium->phase("Solid");
ParameterLib::SpatialPosition x_position;
x_position.setElementID(_element.getID());
@@ -99,17 +98,17 @@ RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
// Initial porosity. Could be read from integration point data or mesh.
ip_data.porosity =
- solid_phase.property(MPL::porosity)
+ medium->property(MPL::porosity)
.template initialValue<double>(
x_position,
std::numeric_limits<
double>::quiet_NaN() /* t independent */);
ip_data.transport_porosity = ip_data.porosity;
- if (solid_phase.hasProperty(MPL::PropertyType::transport_porosity))
+ if (medium->hasProperty(MPL::PropertyType::transport_porosity))
{
ip_data.transport_porosity =
- solid_phase.property(MPL::transport_porosity)
+ medium->property(MPL::transport_porosity)
.template initialValue<double>(
x_position,
std::numeric_limits<
@@ -370,7 +369,7 @@ void RichardsMechanicsLocalAssembler<
variables[static_cast<int>(MPL::Variable::temperature)] = temperature;
auto const alpha =
- solid_phase.property(MPL::PropertyType::biot_coefficient)
+ medium->property(MPL::PropertyType::biot_coefficient)
.template value<double>(variables, x_position, t, dt);
auto const C_el = _ip_data[ip].computeElasticTangentStiffness(
t, x_position, dt, temperature);
@@ -430,7 +429,7 @@ void RichardsMechanicsLocalAssembler<
{ // Porosity update
variables_prev[static_cast<int>(MPL::Variable::porosity)] =
_ip_data[ip].porosity_prev;
- phi = solid_phase.property(MPL::PropertyType::porosity)
+ phi = medium->property(MPL::PropertyType::porosity)
.template value<double>(variables, variables_prev,
x_position, t, dt);
variables[static_cast<int>(MPL::Variable::porosity)] = phi;
@@ -474,14 +473,14 @@ void RichardsMechanicsLocalAssembler<
identity2.transpose() * C_el.inverse() * sigma_sw_prev;
}
- if (solid_phase.hasProperty(MPL::PropertyType::transport_porosity))
+ if (medium->hasProperty(MPL::PropertyType::transport_porosity))
{
variables_prev[static_cast<int>(
MPL::Variable::transport_porosity)] =
_ip_data[ip].transport_porosity_prev;
_ip_data[ip].transport_porosity =
- solid_phase.property(MPL::PropertyType::transport_porosity)
+ medium->property(MPL::PropertyType::transport_porosity)
.template value<double>(variables, variables_prev,
x_position, t, dt);
variables[static_cast<int>(MPL::Variable::transport_porosity)] =
@@ -522,7 +521,7 @@ void RichardsMechanicsLocalAssembler<
_ip_data[ip].material_state_variables->getEquivalentPlasticStrain();
auto const K_intrinsic = MPL::formEigenTensor<DisplacementDim>(
- solid_phase.property(MPL::PropertyType::permeability)
+ medium->property(MPL::PropertyType::permeability)
.value(variables, x_position, t, dt));
GlobalDimMatrixType const rho_K_over_mu =
@@ -738,7 +737,7 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
auto& S_L = _ip_data[ip].saturation;
auto const S_L_prev = _ip_data[ip].saturation_prev;
auto const alpha =
- solid_phase.property(MPL::PropertyType::biot_coefficient)
+ medium->property(MPL::PropertyType::biot_coefficient)
.template value<double>(variables, x_position, t, dt);
auto const C_el = _ip_data[ip].computeElasticTangentStiffness(
@@ -799,7 +798,7 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
variables_prev[static_cast<int>(MPL::Variable::porosity)] =
_ip_data[ip].porosity_prev;
- phi = solid_phase.property(MPL::PropertyType::porosity)
+ phi = medium->property(MPL::PropertyType::porosity)
.template value<double>(variables, variables_prev,
x_position, t, dt);
variables[static_cast<int>(MPL::Variable::porosity)] = phi;
@@ -894,14 +893,14 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
sigma_sw = sigma_sw_prev + delta_sigma_sw;
}
- if (solid_phase.hasProperty(MPL::PropertyType::transport_porosity))
+ if (medium->hasProperty(MPL::PropertyType::transport_porosity))
{
variables_prev[static_cast<int>(
MPL::Variable::transport_porosity)] =
_ip_data[ip].transport_porosity_prev;
_ip_data[ip].transport_porosity =
- solid_phase.property(MPL::PropertyType::transport_porosity)
+ medium->property(MPL::PropertyType::transport_porosity)
.template value<double>(variables, variables_prev,
x_position, t, dt);
variables[static_cast<int>(MPL::Variable::transport_porosity)] =
@@ -934,7 +933,7 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
_ip_data[ip].material_state_variables->getEquivalentPlasticStrain();
auto const K_intrinsic = MPL::formEigenTensor<DisplacementDim>(
- solid_phase.property(MPL::PropertyType::permeability)
+ medium->property(MPL::PropertyType::permeability)
.value(variables, x_position, t, dt));
GlobalDimMatrixType const rho_Ki_over_mu = K_intrinsic * rho_LR / mu;
@@ -1612,7 +1611,7 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
double const chi_S_L_prev = chi(S_L_prev);
auto const alpha =
- solid_phase.property(MPL::PropertyType::biot_coefficient)
+ medium->property(MPL::PropertyType::biot_coefficient)
.template value<double>(variables, x_position, t, dt);
auto const C_el = _ip_data[ip].computeElasticTangentStiffness(
@@ -1640,7 +1639,7 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
{ // Porosity update
variables_prev[static_cast<int>(MPL::Variable::porosity)] =
_ip_data[ip].porosity_prev;
- phi = solid_phase.property(MPL::PropertyType::porosity)
+ phi = medium->property(MPL::PropertyType::porosity)
.template value<double>(variables, variables_prev,
x_position, t, dt);
variables[static_cast<int>(MPL::Variable::porosity)] = phi;
@@ -1676,14 +1675,14 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
identity2.transpose() * C_el.inverse() * sigma_sw_prev;
}
- if (solid_phase.hasProperty(MPL::PropertyType::transport_porosity))
+ if (medium->hasProperty(MPL::PropertyType::transport_porosity))
{
variables_prev[static_cast<int>(
MPL::Variable::transport_porosity)] =
_ip_data[ip].transport_porosity_prev;
_ip_data[ip].transport_porosity =
- solid_phase.property(MPL::PropertyType::transport_porosity)
+ medium->property(MPL::PropertyType::transport_porosity)
.template value<double>(variables, variables_prev,
x_position, t, dt);
variables[static_cast<int>(MPL::Variable::transport_porosity)] =
@@ -1720,7 +1719,7 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
_ip_data[ip].material_state_variables->getEquivalentPlasticStrain();
auto const K_intrinsic = MPL::formEigenTensor<DisplacementDim>(
- solid_phase.property(MPL::PropertyType::permeability)
+ medium->property(MPL::PropertyType::permeability)
.value(variables, x_position, t, dt));
double const k_rel =
diff --git a/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos.prj b/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos.prj
index 3ab898901f72c74f90204467ec4ad998e2edc849..ab932b49e06cfcaf700b23d0208ec129d19eaf09 100644
--- a/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos.prj
+++ b/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<meshes>
<mesh axially_symmetric="true">liakopoulos.vtu</mesh>
@@ -69,25 +69,25 @@
<type>Constant</type>
<value>1e12</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
- <property>
- <name>porosity</name>
- <type>Constant</type>
- <value>0.2975</value>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>4.5e-13</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.2975</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>4.5e-13</value>
+ </property>
<property>
<name>saturation</name>
<type>SaturationLiakopoulos</type>
diff --git a/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj b/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj
index d650a131fe8fb83c718f22429104e1312d721f0e..f2c6c325561f61a0d05b58df684cb13fb4a8a72f 100644
--- a/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj
+++ b/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj
@@ -89,25 +89,25 @@
<type>Constant</type>
<value>1e12</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
- <property>
- <name>porosity</name>
- <type>Constant</type>
- <value>0.2975</value>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>4.5e-13</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.2975</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>4.5e-13</value>
+ </property>
<property>
<name>saturation</name>
<type>SaturationLiakopoulos</type>
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_quasinewton.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_quasinewton.prj
index 5fc7f7f50384f83b3fcabd45b48ee3ce46396054..d2401a89aa873c1f22ee2cdd0e35a3f726fef5a7 100644
--- a/Tests/Data/RichardsMechanics/RichardsFlow_2d_quasinewton.prj
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_quasinewton.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>Richards_2d.vtu</mesh>
<geometry>Richards_2d.gml</geometry>
@@ -57,37 +57,37 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.38000000000000006</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>4.46e-13</value>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value> 0.0 </value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.38000000000000006</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>4.46e-13</value>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value> 0.0 </value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_richardsflow.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_richardsflow.prj
index c84e569c8df916827c415f47091d44ba120565d4..48892f6632b10d5613df37fe4e25923af38cc6e5 100644
--- a/Tests/Data/RichardsMechanics/RichardsFlow_2d_richardsflow.prj
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_richardsflow.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>Richards_2d_linear.vtu</mesh>
<geometry>Richards_2d.gml</geometry>
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small.prj
index e48a44a5ad27d9956924f382bfd486efb52c6158..d725516efce5e67051fb8da3be3427385fadf322 100644
--- a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small.prj
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>RichardsFlow_2d_small.vtu</mesh>
<geometry>RichardsFlow_2d_small.gml</geometry>
@@ -52,37 +52,37 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.38</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>4.46e-13</value>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value> 0.0 </value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.38</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>4.46e-13</value>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value> 0.0 </value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj
index cffc79f4541b4f6e63ef37c5a689945333c8c43b..219071745313779b109bcfea3ef0022aa13b4517 100644
--- a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>Richards_2d.vtu</mesh>
<geometry>Richards_2d.gml</geometry>
@@ -58,32 +58,32 @@
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.38</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value> 0.0 </value>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>4.46e-13</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.38</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value> 0.0 </value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>4.46e-13</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/alternative_mass_balance_anzInterval_10.prj b/Tests/Data/RichardsMechanics/alternative_mass_balance_anzInterval_10.prj
index 173a64e13fab978b3c6f2700a671b68b8cfc7a9f..a5345576e634fc2c8e42f67ffcfea8a7fe505901 100644
--- a/Tests/Data/RichardsMechanics/alternative_mass_balance_anzInterval_10.prj
+++ b/Tests/Data/RichardsMechanics/alternative_mass_balance_anzInterval_10.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad_1e0.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -57,27 +57,27 @@
<type>Constant</type>
<value>2.6e3</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.72222222222222222222</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1e-12</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.72222222222222222222</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-12</value>
+ </property>
<property>
<name>saturation</name>
<type>SaturationVanGenuchten</type>
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj
index 49082707f50d50fa69848ece8e222bc1ac8b55fb..534e4245ac9e56626de3e38e4e0280c587d9b3bb 100644
--- a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>cube_1x1x1_hex20_1e0.vtu</mesh>
<geometry>cube_1x1x1.gml</geometry>
@@ -60,32 +60,32 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj
index c3fc61049a49558d1a4c885204d396c7f0f6b58c..266d4b4efe9033b02b89aa63d7f6fc4e4b39485c 100644
--- a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>cube_1x1x1_hex20_1e0.vtu</mesh>
<geometry>cube_1x1x1.gml</geometry>
@@ -60,32 +60,32 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/confined_compression_fully_saturated.prj b/Tests/Data/RichardsMechanics/confined_compression_fully_saturated.prj
index 767d566c67655de2ce8d4098ae8508185611e4ac..7cc12f8749aab9e9d9276748ac91372d290043d9 100644
--- a/Tests/Data/RichardsMechanics/confined_compression_fully_saturated.prj
+++ b/Tests/Data/RichardsMechanics/confined_compression_fully_saturated.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -64,32 +64,32 @@
<type>Constant</type>
<value>1.2e-6</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1e-12</value>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value>1e+300</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-12</value>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value>1e+300</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/confined_compression_fully_saturated_restart.prj b/Tests/Data/RichardsMechanics/confined_compression_fully_saturated_restart.prj
index 151f3b8e7bfc131bb2a0a5404df61c8b9a39e906..c26b568d822ca0e3377a3049e221bd6c3a61aebe 100644
--- a/Tests/Data/RichardsMechanics/confined_compression_fully_saturated_restart.prj
+++ b/Tests/Data/RichardsMechanics/confined_compression_fully_saturated_restart.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>confined_compression_fully_saturated_ts_20_t_100.000000.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -64,32 +64,32 @@
<type>Constant</type>
<value>1.2e-6</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1e-12</value>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value>1e+300</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-12</value>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value>1e+300</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/deformation_dependent_porosity.prj b/Tests/Data/RichardsMechanics/deformation_dependent_porosity.prj
index f4af3f63b4cca5d2ec3281fc61cf7e890f721a64..d3edf8356999b20873ce08840f2618f36e9eb440 100644
--- a/Tests/Data/RichardsMechanics/deformation_dependent_porosity.prj
+++ b/Tests/Data/RichardsMechanics/deformation_dependent_porosity.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -60,32 +60,32 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.9166666666666666666</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.9166666666666666666</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/deformation_dependent_porosity_swelling.prj b/Tests/Data/RichardsMechanics/deformation_dependent_porosity_swelling.prj
index b89a4fc64d6df7aef0a4b0676026b92650a9e682..8e1cca41544274a2a65d054bde84efe7983c8509 100644
--- a/Tests/Data/RichardsMechanics/deformation_dependent_porosity_swelling.prj
+++ b/Tests/Data/RichardsMechanics/deformation_dependent_porosity_swelling.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -63,11 +63,6 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.9166666666666666666</value>
- </property>
<property>
<name>density</name>
<type>Exponential</type>
@@ -82,25 +77,6 @@
<factor>0.1</factor>
</exponent>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>transport_porosity</name>
- <type>TransportPorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1</value>
- </property>
<property>
<name>swelling_stress_rate</name>
<type>SaturationDependentSwelling</type>
@@ -113,6 +89,30 @@
</phase>
</phases>
<properties>
+ <property>
+ <name>transport_porosity</name>
+ <type>TransportPorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.9166666666666666666</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/double_porosity_swelling.prj b/Tests/Data/RichardsMechanics/double_porosity_swelling.prj
index ddd75f69dc2a2982c5bd755dd015d84e65d0bc2a..1b80e4b2fb3d77c2d4139fce434269899c101c05 100644
--- a/Tests/Data/RichardsMechanics/double_porosity_swelling.prj
+++ b/Tests/Data/RichardsMechanics/double_porosity_swelling.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -68,30 +68,6 @@
<type>Constant</type>
<value>2e3</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.6</value>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1e-21</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>transport_porosity</name>
- <type>TransportPorosityFromMassBalance</type>
- <initial_porosity>phi_tr0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
<property>
<name>swelling_stress_rate</name>
<type>SaturationDependentSwelling</type>
@@ -104,6 +80,30 @@
</phase>
</phases>
<properties>
+ <property>
+ <name>transport_porosity</name>
+ <type>TransportPorosityFromMassBalance</type>
+ <initial_porosity>phi_tr0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.6</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-21</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/flow_fully_saturated.prj b/Tests/Data/RichardsMechanics/flow_fully_saturated.prj
index 1e0c4f33d2b4b551b60ac15b2d84cf3e5a398423..f6399282c2afc99163796af2c5ba02958767a35f 100644
--- a/Tests/Data/RichardsMechanics/flow_fully_saturated.prj
+++ b/Tests/Data/RichardsMechanics/flow_fully_saturated.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -57,32 +57,32 @@
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value>1e+300</value>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1e-12</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value>1e+300</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-12</value>
+ </property>
<property>
<name>saturation</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/flow_fully_saturated_anisotropic.prj b/Tests/Data/RichardsMechanics/flow_fully_saturated_anisotropic.prj
index b94b8598c69e53da0566fdb6e695735d8e0120ee..83bbe072e3f606c32eb4068bfa4b8f85f6dc2fb4 100644
--- a/Tests/Data/RichardsMechanics/flow_fully_saturated_anisotropic.prj
+++ b/Tests/Data/RichardsMechanics/flow_fully_saturated_anisotropic.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -57,32 +57,32 @@
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>7.75e-13 3.8971143170299739104e-13 3.8971143170299739104e-13 3.25e-13</value>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value>1e+300</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>7.75e-13 3.8971143170299739104e-13 3.8971143170299739104e-13 3.25e-13</value>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value>1e+300</value>
+ </property>
<property>
<name>saturation</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/flow_fully_saturated_coordinate_system.prj b/Tests/Data/RichardsMechanics/flow_fully_saturated_coordinate_system.prj
index a8718407e652c538a90f094810e2df9d2df5bb9e..9054472c1e6563be37e0ef060d4640ab37dba7ac 100644
--- a/Tests/Data/RichardsMechanics/flow_fully_saturated_coordinate_system.prj
+++ b/Tests/Data/RichardsMechanics/flow_fully_saturated_coordinate_system.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -52,37 +52,37 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Parameter</type>
- <parameter_name>k</parameter_name>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value>1e+300</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Parameter</type>
+ <parameter_name>k</parameter_name>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value>1e+300</value>
+ </property>
<property>
<name>saturation</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/flow_fully_saturated_linear.prj b/Tests/Data/RichardsMechanics/flow_fully_saturated_linear.prj
index af25126e129f51fe4e28bf1ec55e26e9bb3ba089..f1e346bdd9354ee0de3f137621c6385cf16e4493 100644
--- a/Tests/Data/RichardsMechanics/flow_fully_saturated_linear.prj
+++ b/Tests/Data/RichardsMechanics/flow_fully_saturated_linear.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -57,32 +57,32 @@
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value>1e+300</value>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1e-12</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value>1e+300</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-12</value>
+ </property>
<property>
<name>saturation</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/gravity.prj b/Tests/Data/RichardsMechanics/gravity.prj
index f58db23e3e007e94e5fedd2da48feeba0fbb6ddf..5852cfe24ffed067d01871dfb33588836cc5be46 100644
--- a/Tests/Data/RichardsMechanics/gravity.prj
+++ b/Tests/Data/RichardsMechanics/gravity.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -59,37 +59,37 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.5</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1e3</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1e-12</value>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value> 0.0 </value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-12</value>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value> 0.0 </value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/mechanics_linear.prj b/Tests/Data/RichardsMechanics/mechanics_linear.prj
index e6b12323ca6bfae296acb408d5234cb6a859b6bc..9ed864d8803135fbc7330b05cf42f3cdcfc2edb0 100644
--- a/Tests/Data/RichardsMechanics/mechanics_linear.prj
+++ b/Tests/Data/RichardsMechanics/mechanics_linear.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<!--
Mechanics setup identical to Mechanics/Linear/square_1e2.prj.
No flow and no coupling between pressure and displacement.
@@ -65,32 +65,32 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.5</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1</value>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
<property>
<name>saturation</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/orthotropic_power_law_permeability_xyz.prj b/Tests/Data/RichardsMechanics/orthotropic_power_law_permeability_xyz.prj
index 613b6e754ec613f0ac67999c67a552d5cea91ce9..ddd221c60194db248bc83e22c605a1c835fd5b0b 100644
--- a/Tests/Data/RichardsMechanics/orthotropic_power_law_permeability_xyz.prj
+++ b/Tests/Data/RichardsMechanics/orthotropic_power_law_permeability_xyz.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -61,40 +61,40 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.94444444444444444444</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>transport_porosity</name>
- <type>TransportPorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>PermeabilityOrthotropicPowerLaw</type>
- <intrinsic_permeabilities>0 1</intrinsic_permeabilities>
- <exponents>1 1</exponents>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>transport_porosity</name>
+ <type>TransportPorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.94444444444444444444</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>PermeabilityOrthotropicPowerLaw</type>
+ <intrinsic_permeabilities>0 1</intrinsic_permeabilities>
+ <exponents>1 1</exponents>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/orthotropic_swelling_xy.prj b/Tests/Data/RichardsMechanics/orthotropic_swelling_xy.prj
index 7a252dfe3bf3ebb4369ce92e1bd3b4ee7268271b..a45669255e7cac57acc5b67e28ee6de0b3b4212f 100644
--- a/Tests/Data/RichardsMechanics/orthotropic_swelling_xy.prj
+++ b/Tests/Data/RichardsMechanics/orthotropic_swelling_xy.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>square_1x1_quad8_1e0.vtu</mesh>
<geometry>square_1x1.gml</geometry>
@@ -61,28 +61,11 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.5</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1</value>
- </property>
<property>
<name>swelling_stress_rate</name>
<type>SaturationDependentSwelling</type>
@@ -95,6 +78,23 @@
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/orthotropic_swelling_xyz.prj b/Tests/Data/RichardsMechanics/orthotropic_swelling_xyz.prj
index 39dcbb6a7ab037b5f1b80ead1243745fa1ffa3e5..4fa60dc19ab243114ad03a1162f8c75720a3834b 100644
--- a/Tests/Data/RichardsMechanics/orthotropic_swelling_xyz.prj
+++ b/Tests/Data/RichardsMechanics/orthotropic_swelling_xyz.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<OpenGeoSysProject>
<mesh>cube_1x1x1_hex20_1e0.vtu</mesh>
<geometry>cube_1x1x1.gml</geometry>
@@ -64,28 +64,11 @@
<phase>
<type>Solid</type>
<properties>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>0.5</value>
- </property>
<property>
<name>density</name>
<type>Constant</type>
<value>1</value>
</property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
- <property>
- <name>permeability</name>
- <type>Constant</type>
- <value>1</value>
- </property>
<property>
<name>swelling_stress_rate</name>
<type>SaturationDependentSwelling</type>
@@ -98,6 +81,23 @@
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
diff --git a/Tests/Data/RichardsMechanics/rotated_consolidation.prj b/Tests/Data/RichardsMechanics/rotated_consolidation.prj
index 515d67cc279843cb6993a02b81c90774ec89849b..196e685253c844c1e711d6faf2356ad2df0fc9e1 100644
--- a/Tests/Data/RichardsMechanics/rotated_consolidation.prj
+++ b/Tests/Data/RichardsMechanics/rotated_consolidation.prj
@@ -1,4 +1,4 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml version='1.0' encoding='ISO-8859-1'?>
<!--
Modified consolidation test
- Column geometry rotated at -45 degrees
@@ -85,32 +85,32 @@
<type>Constant</type>
<value>2780.0</value>
</property>
- <property>
- <name>biot_coefficient</name>
- <type>Constant</type>
- <value>1</value>
- </property>
- <property>
- <name>permeability</name>
- <type>Parameter</type>
- <parameter_name>k_i</parameter_name>
- </property>
- <property>
- <name>storage</name>
- <type>Constant</type>
- <value>1e-10</value>
- </property>
- <property>
- <name>porosity</name>
- <type>PorosityFromMassBalance</type>
- <initial_porosity>phi0</initial_porosity>
- <minimal_porosity>0</minimal_porosity>
- <maximal_porosity>1</maximal_porosity>
- </property>
</properties>
</phase>
</phases>
<properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Parameter</type>
+ <parameter_name>k_i</parameter_name>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value>1e-10</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>PorosityFromMassBalance</type>
+ <initial_porosity>phi0</initial_porosity>
+ <minimal_porosity>0</minimal_porosity>
+ <maximal_porosity>1</maximal_porosity>
+ </property>
<property>
<name>reference_temperature</name>
<type>Constant</type>