diff --git a/MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.cpp b/MaterialLib/MPL/Utils/FormKelvinVector.cpp
similarity index 75%
rename from MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.cpp
rename to MaterialLib/MPL/Utils/FormKelvinVector.cpp
index 869581d222085d2f8c8e1fe5d804420b0528d777..e1fd045d23f1c0837883d0e4af1c6f047ddd6759 100644
--- a/MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.cpp
+++ b/MaterialLib/MPL/Utils/FormKelvinVector.cpp
@@ -9,19 +9,19 @@
* Created on February 12, 2021, 4:34 PM
*/
-#include "FormKelvinVectorFromThermalExpansivity.h"
+#include "FormKelvinVector.h"
#include "BaseLib/Error.h"
namespace MaterialPropertyLib
{
static constexpr const char error_info[] =
- "The thermal expansivity can only be either a scalar number for isotropic "
- "thermal expansion or a three element array for anisotropic thermal "
- "expansion.";
+ "The conversion to a Kelvin vector of correct dimensionality is ambiguous."
+ "Please use a scalar number for isotropic properties, a three element "
+ "array or a 3 x 3 matrix for anisotropic properties.";
template <int GlobalDim>
-struct FormKelvinVectorFromThermalExpansivity
+struct FormKelvinVector
{
MathLib::KelvinVector::KelvinVectorType<GlobalDim> operator()(
double const& value) const
@@ -54,9 +54,9 @@ struct FormKelvinVectorFromThermalExpansivity
}
MathLib::KelvinVector::KelvinVectorType<GlobalDim> operator()(
- Eigen::Matrix<double, 3, 3> const& /*values*/) const
+ Eigen::Matrix<double, 3, 3> const& values) const
{
- OGS_FATAL(error_info);
+ return MathLib::KelvinVector::tensorToKelvin<GlobalDim>(values);
}
MathLib::KelvinVector::KelvinVectorType<GlobalDim> operator()(
@@ -79,20 +79,16 @@ struct FormKelvinVectorFromThermalExpansivity
};
template <int GlobalDim>
-MathLib::KelvinVector::KelvinVectorType<GlobalDim>
-formKelvinVectorFromThermalExpansivity(
+MathLib::KelvinVector::KelvinVectorType<GlobalDim> formKelvinVector(
MaterialPropertyLib::PropertyDataType const& values)
{
- return std::visit(FormKelvinVectorFromThermalExpansivity<GlobalDim>(),
- values);
+ return std::visit(FormKelvinVector<GlobalDim>(), values);
}
-template MathLib::KelvinVector::KelvinVectorType<2>
-formKelvinVectorFromThermalExpansivity<2>(
+template MathLib::KelvinVector::KelvinVectorType<2> formKelvinVector<2>(
MaterialPropertyLib::PropertyDataType const& values);
-template MathLib::KelvinVector::KelvinVectorType<3>
-formKelvinVectorFromThermalExpansivity<3>(
+template MathLib::KelvinVector::KelvinVectorType<3> formKelvinVector<3>(
MaterialPropertyLib::PropertyDataType const& values);
} // namespace MaterialPropertyLib
diff --git a/MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.h b/MaterialLib/MPL/Utils/FormKelvinVector.h
similarity index 51%
rename from MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.h
rename to MaterialLib/MPL/Utils/FormKelvinVector.h
index 910cfa23ea48a2af868e5bb37827cd4d2b1d2cc5..a7db805b7fe2fbd4aa3cd542ad93e6458bfb5ee1 100644
--- a/MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.h
+++ b/MaterialLib/MPL/Utils/FormKelvinVector.h
@@ -20,20 +20,19 @@
namespace MaterialPropertyLib
{
/**
- * \brief A function to form a Kelvin vector from thermal expansivity for
- * thermal strain.
+ * \brief A function to form a Kelvin vector from strain or stress alike
+ * property like thermal expansivity for thermal strain.
*
- * It takes the thermal expansivity, either a scalar number for isotropic
- * thermal expansion or a three element vector for anisotropic thermal
- * expansion, to get a Kelvin vector for thermal strain.
+ * It takes either a scalar number for isotropic thermal expansion or a
+ * three element vector or a 3 x 3 matrix for anisotropic properties, to
+ * get a Kelvin vector for strain or stress.
*
- * @param values: Thermal expansivity, which can be scalar number or a three
- * element vector.
- * @return Thermal expansivity in Kelvin vector type.
+ * @param values: e.g., Thermal expansivity, which can be scalar number, a three
+ * element vector or a 3 x 3 matrix.
+ * @return Kelvin vector type property.
*/
template <int GlobalDim>
-MathLib::KelvinVector::KelvinVectorType<GlobalDim>
-formKelvinVectorFromThermalExpansivity(
+MathLib::KelvinVector::KelvinVectorType<GlobalDim> formKelvinVector(
MaterialPropertyLib::PropertyDataType const& values);
} // namespace MaterialPropertyLib
diff --git a/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsFEM-impl.h b/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsFEM-impl.h
index 910522822eba7b67bdd92677044784eb1f3d2a2a..c9597713ae7a4387a4b99e5099f14b231a0d3cad 100644
--- a/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsFEM-impl.h
+++ b/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsFEM-impl.h
@@ -15,7 +15,7 @@
#include "MaterialLib/MPL/Medium.h"
#include "MaterialLib/MPL/Property.h"
#include "MaterialLib/MPL/Utils/FormEigenTensor.h"
-#include "MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.h"
+#include "MaterialLib/MPL/Utils/FormKelvinVector.h"
#include "MaterialLib/MPL/Utils/GetLiquidThermalExpansivity.h"
#include "MaterialLib/SolidModels/SelectSolidConstitutiveRelation.h"
#include "MathLib/KelvinVector.h"
@@ -363,7 +363,7 @@ void ThermoHydroMechanicsLocalAssembler<
// Consider also anisotropic thermal expansion.
MathLib::KelvinVector::KelvinVectorType<
DisplacementDim> const solid_linear_thermal_expansion_coefficient =
- MPL::formKelvinVectorFromThermalExpansivity<DisplacementDim>(
+ MPL::formKelvinVector<DisplacementDim>(
solid_phase
.property(
MaterialPropertyLib::PropertyType::thermal_expansivity)
@@ -834,7 +834,7 @@ void ThermoHydroMechanicsLocalAssembler<
// Consider also anisotropic thermal expansion.
MathLib::KelvinVector::KelvinVectorType<
DisplacementDim> const solid_linear_thermal_expansion_coefficient =
- MPL::formKelvinVectorFromThermalExpansivity<DisplacementDim>(
+ MPL::formKelvinVector<DisplacementDim>(
solid_phase
.property(
MaterialPropertyLib::PropertyType::thermal_expansivity)
diff --git a/ProcessLib/ThermoMechanics/ThermoMechanicsFEM-impl.h b/ProcessLib/ThermoMechanics/ThermoMechanicsFEM-impl.h
index a44b9d7b433ffa77765c457800d3385e0ffcf0f3..01e3a9edfca0b50d554a3230cd0003b3e8f77065 100644
--- a/ProcessLib/ThermoMechanics/ThermoMechanicsFEM-impl.h
+++ b/ProcessLib/ThermoMechanics/ThermoMechanicsFEM-impl.h
@@ -15,7 +15,7 @@
#include "MaterialLib/MPL/Medium.h"
#include "MaterialLib/MPL/Property.h"
#include "MaterialLib/MPL/Utils/FormEigenTensor.h"
-#include "MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.h"
+#include "MaterialLib/MPL/Utils/FormKelvinVector.h"
#include "ProcessLib/Utils/SetOrGetIntegrationPointData.h"
#include "ProcessLib/Utils/TransposeInPlace.h"
@@ -197,7 +197,7 @@ void ThermoMechanicsLocalAssembler<ShapeFunction, DisplacementDim>::
// Consider also anisotropic thermal expansion.
auto const solid_linear_thermal_expansivity_vector =
- MPL::formKelvinVectorFromThermalExpansivity<DisplacementDim>(
+ MPL::formKelvinVector<DisplacementDim>(
solid_phase
.property(
MaterialPropertyLib::PropertyType::thermal_expansivity)
@@ -429,7 +429,7 @@ void ThermoMechanicsLocalAssembler<ShapeFunction, DisplacementDim>::
// Consider also anisotropic thermal expansion.
auto const solid_linear_thermal_expansivity_vector =
- MPL::formKelvinVectorFromThermalExpansivity<DisplacementDim>(
+ MPL::formKelvinVector<DisplacementDim>(
solid_phase
.property(
MaterialPropertyLib::PropertyType::thermal_expansivity)
diff --git a/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/SolidThermalExpansion.cpp b/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/SolidThermalExpansion.cpp
index e7cc364a94034ffb24d3d01ad41b09e3dd8d6113..70879c949403d159719396590cc3f867099c1384 100644
--- a/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/SolidThermalExpansion.cpp
+++ b/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/SolidThermalExpansion.cpp
@@ -10,7 +10,7 @@
#include "SolidThermalExpansion.h"
-#include "MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.h"
+#include "MaterialLib/MPL/Utils/FormKelvinVector.h"
namespace ProcessLib::ThermoRichardsMechanics
{
@@ -24,7 +24,7 @@ void SolidThermalExpansionModel<DisplacementDim>::eval(
// Consider also anisotropic thermal expansion.
out.solid_linear_thermal_expansivity_vector =
- MPL::formKelvinVectorFromThermalExpansivity<DisplacementDim>(
+ MPL::formKelvinVector<DisplacementDim>(
media_data.solid.property(MPL::PropertyType::thermal_expansivity)
.value(variables, x_t.x, x_t.t, x_t.dt));
}
diff --git a/ProcessLib/ThermoRichardsMechanics/Tests.cmake b/ProcessLib/ThermoRichardsMechanics/Tests.cmake
index 2d37f1658d1dd6c62b66a1bdf9047ce00cec0d8d..ca778a279cf17cbc34e8a3518ab84a0b64b605b5 100644
--- a/ProcessLib/ThermoRichardsMechanics/Tests.cmake
+++ b/ProcessLib/ThermoRichardsMechanics/Tests.cmake
@@ -7,6 +7,96 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE ThermoRichardsMechanics/OrthotropicSwelling/orthotropic_swelling_xyz.xml)
endif()
+AddTest(
+ NAME ThermoRichardsMechanics_anisotropic_thermal_expansion_vector
+ PATH ThermoRichardsMechanics/anisotropic_thermal_expansion
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS aniso_expansion.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ RUNTIME 1
+ DIFF_DATA
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu sigma sigma 1e-12 1e-12
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu displacement displacement 1e-14 1e-14
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu saturation saturation 1e-12 1e-12
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu pressure pressure 1e-10 1e-10
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu epsilon epsilon 1e-14 1e-14
+)
+
+AddTest(
+ NAME ThermoRichardsMechanics_anisotropic_thermal_expansion_expansivity_matrix
+ PATH ThermoRichardsMechanics/anisotropic_thermal_expansion
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS aniso_expansion_expansivity_matrix.xml
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ RUNTIME 1
+ DIFF_DATA
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_ts_1_t_1000000.000000.vtu sigma sigma 1e-12 1e-12
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_ts_1_t_1000000.000000.vtu displacement displacement 1e-14 1e-14
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_ts_1_t_1000000.000000.vtu saturation saturation 1e-14 1e-14
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_ts_1_t_1000000.000000.vtu pressure pressure 1e-10 1e-10
+ expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_ts_1_t_1000000.000000.vtu epsilon epsilon 1e-12 1e-12
+)
+
+AddTest(
+ NAME ThermoRichardsMechanics_anisotropic_thermal_expansion_expansivity_matrix_z90
+ PATH ThermoRichardsMechanics/anisotropic_thermal_expansion
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS aniso_expansion_expansivity_matrix_z90.xml
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ RUNTIME 1
+ DIFF_DATA
+ expected_anisotropic_thermal_expansion_z90_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_z90_ts_1_t_1000000.000000.vtu sigma sigma 1e-12 1e-12
+ expected_anisotropic_thermal_expansion_z90_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_z90_ts_1_t_1000000.000000.vtu displacement displacement 1e-14 1e-14
+ expected_anisotropic_thermal_expansion_z90_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_z90_ts_1_t_1000000.000000.vtu saturation saturation 1e-12 1e-12
+ expected_anisotropic_thermal_expansion_z90_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_z90_ts_1_t_1000000.000000.vtu pressure pressure 1e-10 1e-10
+ expected_anisotropic_thermal_expansion_z90_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_expansivity_matrix_z90_ts_1_t_1000000.000000.vtu epsilon epsilon 1e-12 1e-12
+)
+
+AddTest(
+ NAME ThermoRichardsMechanics_anisotropic_thermal_expansion_x45
+ PATH ThermoRichardsMechanics/anisotropic_thermal_expansion
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS aniso_expansion_x45.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ RUNTIME 1
+ DIFF_DATA
+ expected_anisotropic_thermal_expansion_x45_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_x45_ts_1_t_1000000.000000.vtu epsilon epsilon 1e-12 1e-12
+)
+
+AddTest(
+ NAME ThermoRichardsMechanics_anisotropic_thermal_expansion_y45
+ PATH ThermoRichardsMechanics/anisotropic_thermal_expansion
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS aniso_expansion_y45.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ RUNTIME 1
+ DIFF_DATA
+ expected_anisotropic_thermal_expansion_y45_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_y45_ts_1_t_1000000.000000.vtu epsilon epsilon 1e-12 1e-12
+)
+
+AddTest(
+ NAME ThermoRichardsMechanics_anisotropic_thermal_expansion_z45
+ PATH ThermoRichardsMechanics/anisotropic_thermal_expansion
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS aniso_expansion_z45.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ RUNTIME 1
+ DIFF_DATA
+ expected_anisotropic_thermal_expansion_z45_ts_1_t_1000000.000000.vtu anisotropic_thermal_expansion_z45_ts_1_t_1000000.000000.vtu epsilon epsilon 1e-12 1e-12
+)
+
AddTest(
NAME ThermoRichardsMechanics_liakopoulosHM
PATH ThermoRichardsMechanics/LiakopoulosHM
diff --git a/ProcessLib/ThermoRichardsMechanics/ThermoRichardsMechanicsFEM-impl.h b/ProcessLib/ThermoRichardsMechanics/ThermoRichardsMechanicsFEM-impl.h
index 19e618a14d47d88041cc8fcfcb64c8c053f5c523..a2c741b908051634781467748951e73d1b524fc4 100644
--- a/ProcessLib/ThermoRichardsMechanics/ThermoRichardsMechanicsFEM-impl.h
+++ b/ProcessLib/ThermoRichardsMechanics/ThermoRichardsMechanicsFEM-impl.h
@@ -18,7 +18,6 @@
#include "MaterialLib/MPL/Medium.h"
#include "MaterialLib/MPL/Utils/FormEigenTensor.h"
-#include "MaterialLib/MPL/Utils/FormKelvinVectorFromThermalExpansivity.h"
#include "MaterialLib/MPL/Utils/GetLiquidThermalExpansivity.h"
#include "MaterialLib/PhysicalConstant.h"
#include "MathLib/KelvinVector.h"
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion.prj b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion.prj
new file mode 100644
index 0000000000000000000000000000000000000000..47ce3430c0dbfc37d99c947c311934c834c0ea9f
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion.prj
@@ -0,0 +1,340 @@
+<?xml version='1.0' encoding='ISO-8859-1'?>
+<OpenGeoSysProject>
+ <mesh>cube_1x1x1_hex_1e0.vtu</mesh>
+ <geometry>cube_1x1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>TRM</name>
+ <type>THERMO_RICHARDS_MECHANICS</type>
+ <integration_order>3</integration_order>
+ <constitutive_relation id="0">
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <temperature>temperature</temperature>
+ <pressure>pressure</pressure>
+ <displacement>displacement</displacement>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="velocity" output_name="velocity"/>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>600.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2600</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>600.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Parameter</type>
+ <parameter_name>expansivity_vector</parameter_name>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>relative_permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>saturation</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>bishops_effective_stress</name>
+ <type>BishopsSaturationCutoff</type>
+ <cutoff_value>0.1</cutoff_value>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-17</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>EffectiveThermalConductivityPorosityMixing</type>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="TRM">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-10 5e-12 1e-12 1e-12 1e-12</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1000000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1000000</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>anisotropic_thermal_expansion</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1000000</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>temperature</variable>
+ <variable>displacement</variable>
+ <variable>pressure</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity</variable>
+ <variable>saturation</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <local_coordinate_system>
+ <basis_vector_0>e0</basis_vector_0>
+ <basis_vector_1>e1</basis_vector_1>
+ <basis_vector_2>e2</basis_vector_2>
+ </local_coordinate_system>
+ <parameters>
+ <parameter>
+ <name>e0</name>
+ <type>Constant</type>
+ <values>0 1 0</values>
+ </parameter>
+ <parameter>
+ <name>e1</name>
+ <type>Constant</type>
+ <values>-1 0 0</values>
+ </parameter>
+ <parameter>
+ <name>e2</name>
+ <type>Constant</type>
+ <values>0 0 1</values>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e6</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.3</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>pressure_ic</name>
+ <type>Constant</type>
+ <values>1e6</values>
+ </parameter>
+ <parameter>
+ <name>0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>expansivity_tensor</name>
+ <type>Constant</type>
+ <values>1e-4 0 0 0 2e-4 0 0 0 3e-4</values>
+ <use_local_coordinate_system>false</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>expansivity_vector</name>
+ <type>Constant</type>
+ <values>1e-4 2e-4 3e-4</values>
+ <use_local_coordinate_system>false</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>T_ref</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>CurveScaled</type>
+ <curve>temp_c</curve>
+ <parameter>T_ref</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>temp_c</name>
+ <coords>1 1000000</coords>
+ <values>1 1.003411222923418</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T_ref</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>back</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>pressure_ic</initial_condition>
+ </process_variable>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</type>
+ <max_iter>50</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_expansivity_matrix.xml b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_expansivity_matrix.xml
new file mode 100644
index 0000000000000000000000000000000000000000..b2504ef7a269831f59e84bf0bad7ae62e381e3e0
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_expansivity_matrix.xml
@@ -0,0 +1,5 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProjectDiff base_file="aniso_expansion.prj">
+ <replace sel="/*/media/medium/phases/phase/properties/property/name[text()="thermal_expansivity"]/../parameter_name/text()">expansivity_tensor</replace>
+ <replace sel="/*/time_loop/output/prefix/text()">anisotropic_thermal_expansion_expansivity_matrix</replace>
+</OpenGeoSysProjectDiff>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_expansivity_matrix_z90.xml b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_expansivity_matrix_z90.xml
new file mode 100644
index 0000000000000000000000000000000000000000..230e5ecd9fb5897a17d94dd277f07265cbcbba1f
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_expansivity_matrix_z90.xml
@@ -0,0 +1,6 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProjectDiff base_file="aniso_expansion.prj">
+ <replace sel="/*/media/medium/phases/phase/properties/property/name[text()="thermal_expansivity"]/../parameter_name/text()">expansivity_tensor</replace>
+ <replace sel="/*/parameters/parameter/name[text()="expansivity_tensor"]/../use_local_coordinate_system/text()">true</replace>
+ <replace sel="/*/time_loop/output/prefix/text()">anisotropic_thermal_expansion_expansivity_matrix_z90</replace>
+</OpenGeoSysProjectDiff>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_x45.prj b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_x45.prj
new file mode 100644
index 0000000000000000000000000000000000000000..600a835f70aa90dd1dd351287e004bfd51cdc218
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_x45.prj
@@ -0,0 +1,340 @@
+<?xml version='1.0' encoding='ISO-8859-1'?>
+<OpenGeoSysProject>
+ <mesh>cube_1x1x1_hex_1e0.vtu</mesh>
+ <geometry>cube_1x1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>TRM</name>
+ <type>THERMO_RICHARDS_MECHANICS</type>
+ <integration_order>3</integration_order>
+ <constitutive_relation id="0">
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <temperature>temperature</temperature>
+ <pressure>pressure</pressure>
+ <displacement>displacement</displacement>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="velocity" output_name="velocity"/>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>600.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2600</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>600.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Parameter</type>
+ <parameter_name>expansivity_vector</parameter_name>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>relative_permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>saturation</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>bishops_effective_stress</name>
+ <type>BishopsSaturationCutoff</type>
+ <cutoff_value>0.1</cutoff_value>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-17</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>EffectiveThermalConductivityPorosityMixing</type>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="TRM">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-10 5e-4 1e-12 1e-12 1e-12</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1000000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1000000</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>anisotropic_thermal_expansion_x45</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1000000</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>temperature</variable>
+ <variable>displacement</variable>
+ <variable>pressure</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity</variable>
+ <variable>saturation</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <local_coordinate_system>
+ <basis_vector_0>e0</basis_vector_0>
+ <basis_vector_1>e1</basis_vector_1>
+ <basis_vector_2>e2</basis_vector_2>
+ </local_coordinate_system>
+ <parameters>
+ <parameter>
+ <name>e0</name>
+ <type>Constant</type>
+ <values>1 0 0</values>
+ </parameter>
+ <parameter>
+ <name>e1</name>
+ <type>Constant</type>
+ <values>0 0.7071067811865476 0.7071067811865476</values>
+ </parameter>
+ <parameter>
+ <name>e2</name>
+ <type>Constant</type>
+ <values>0 -0.7071067811865476 0.7071067811865476</values>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e60</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.0</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>pressure_ic</name>
+ <type>Constant</type>
+ <values>1e6</values>
+ </parameter>
+ <parameter>
+ <name>0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>expansivity_tensor</name>
+ <type>Constant</type>
+ <values>0 0 0 0 2.0e-4 0 0 0 0</values>
+ <use_local_coordinate_system>true</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>expansivity_vector</name>
+ <type>Constant</type>
+ <values>0 2e-4 0</values>
+ <use_local_coordinate_system>true</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>T_ref</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>CurveScaled</type>
+ <curve>temp_c</curve>
+ <parameter>T_ref</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>temp_c</name>
+ <coords>1 1000000</coords>
+ <values>1 1.003411222923418</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T_ref</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>back</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>pressure_ic</initial_condition>
+ </process_variable>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ <!--boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>0</parameter>
+ </boundary_condition-->
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</type>
+ <max_iter>50</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_y45.prj b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_y45.prj
new file mode 100644
index 0000000000000000000000000000000000000000..7b588aefdf28f704f3900975a90ad057c2fda2f2
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_y45.prj
@@ -0,0 +1,340 @@
+<?xml version='1.0' encoding='ISO-8859-1'?>
+<OpenGeoSysProject>
+ <mesh>cube_1x1x1_hex_1e0.vtu</mesh>
+ <geometry>cube_1x1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>TRM</name>
+ <type>THERMO_RICHARDS_MECHANICS</type>
+ <integration_order>3</integration_order>
+ <constitutive_relation id="0">
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <temperature>temperature</temperature>
+ <pressure>pressure</pressure>
+ <displacement>displacement</displacement>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="velocity" output_name="velocity"/>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>600.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2600</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>600.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Parameter</type>
+ <parameter_name>expansivity_vector</parameter_name>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>relative_permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>saturation</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>bishops_effective_stress</name>
+ <type>BishopsSaturationCutoff</type>
+ <cutoff_value>0.1</cutoff_value>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-17</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>EffectiveThermalConductivityPorosityMixing</type>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="TRM">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-10 5e-4 1e-12 1e-12 1e-12</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1000000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1000000</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>anisotropic_thermal_expansion_y45</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1000000</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>temperature</variable>
+ <variable>displacement</variable>
+ <variable>pressure</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity</variable>
+ <variable>saturation</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <local_coordinate_system>
+ <basis_vector_0>e0</basis_vector_0>
+ <basis_vector_1>e1</basis_vector_1>
+ <basis_vector_2>e2</basis_vector_2>
+ </local_coordinate_system>
+ <parameters>
+ <parameter>
+ <name>e0</name>
+ <type>Constant</type>
+ <values>0.7071067811865476 0.0 -0.7071067811865476</values>
+ </parameter>
+ <parameter>
+ <name>e1</name>
+ <type>Constant</type>
+ <values>0 1 0</values>
+ </parameter>
+ <parameter>
+ <name>e2</name>
+ <type>Constant</type>
+ <values>0.7071067811865476 0.0 0.7071067811865476</values>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e60</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.0</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>pressure_ic</name>
+ <type>Constant</type>
+ <values>1e6</values>
+ </parameter>
+ <parameter>
+ <name>0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>expansivity_tensor</name>
+ <type>Constant</type>
+ <values>2.0e-4 0 0 0 0 0 0 0 0</values>
+ <use_local_coordinate_system>true</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>expansivity_vector</name>
+ <type>Constant</type>
+ <values>2e-4 0 0</values>
+ <use_local_coordinate_system>true</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>T_ref</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>CurveScaled</type>
+ <curve>temp_c</curve>
+ <parameter>T_ref</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>temp_c</name>
+ <coords>1 1000000</coords>
+ <values>1 1.003411222923418</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T_ref</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>back</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>pressure_ic</initial_condition>
+ </process_variable>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ <!--boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>0</parameter>
+ </boundary_condition-->
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</type>
+ <max_iter>50</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_z45.prj b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_z45.prj
new file mode 100644
index 0000000000000000000000000000000000000000..e1e7f91cdddeceb8bbe989023d034f5b10260b1e
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/aniso_expansion_z45.prj
@@ -0,0 +1,340 @@
+<?xml version='1.0' encoding='ISO-8859-1'?>
+<OpenGeoSysProject>
+ <mesh>cube_1x1x1_hex_1e0.vtu</mesh>
+ <geometry>cube_1x1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>TRM</name>
+ <type>THERMO_RICHARDS_MECHANICS</type>
+ <integration_order>3</integration_order>
+ <constitutive_relation id="0">
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <temperature>temperature</temperature>
+ <pressure>pressure</pressure>
+ <displacement>displacement</displacement>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="velocity" output_name="velocity"/>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>600.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2600</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>600.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Parameter</type>
+ <parameter_name>expansivity_vector</parameter_name>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>relative_permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>saturation</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>bishops_effective_stress</name>
+ <type>BishopsSaturationCutoff</type>
+ <cutoff_value>0.1</cutoff_value>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1e-17</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>EffectiveThermalConductivityPorosityMixing</type>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="TRM">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-10 5e-4 1e-12 1e-12 1e-12</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1000000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1000000</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>anisotropic_thermal_expansion_z45</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1000000</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>temperature</variable>
+ <variable>displacement</variable>
+ <variable>pressure</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity</variable>
+ <variable>saturation</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <local_coordinate_system>
+ <basis_vector_0>e0</basis_vector_0>
+ <basis_vector_1>e1</basis_vector_1>
+ <basis_vector_2>e2</basis_vector_2>
+ </local_coordinate_system>
+ <parameters>
+ <parameter>
+ <name>e0</name>
+ <type>Constant</type>
+ <values>0.7071067811865476 0.7071067811865476 0</values>
+ </parameter>
+ <parameter>
+ <name>e1</name>
+ <type>Constant</type>
+ <values>-0.7071067811865476 0.7071067811865476 0</values>
+ </parameter>
+ <parameter>
+ <name>e2</name>
+ <type>Constant</type>
+ <values>0 0 1</values>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e60</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.0</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>pressure_ic</name>
+ <type>Constant</type>
+ <values>1e6</values>
+ </parameter>
+ <parameter>
+ <name>0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>expansivity_tensor</name>
+ <type>Constant</type>
+ <values>0 0 0 0 2.0e-4 0 0 0 0</values>
+ <use_local_coordinate_system>true</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>expansivity_vector</name>
+ <type>Constant</type>
+ <values>0 2e-4 0</values>
+ <use_local_coordinate_system>true</use_local_coordinate_system>
+ </parameter>
+ <parameter>
+ <name>T_ref</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>CurveScaled</type>
+ <curve>temp_c</curve>
+ <parameter>T_ref</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>temp_c</name>
+ <coords>1 1000000</coords>
+ <values>1 1.003411222923418</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T_ref</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>back</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>pressure_ic</initial_condition>
+ </process_variable>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ <!--boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>0</parameter>
+ </boundary_condition-->
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</type>
+ <max_iter>50</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/cube_1x1x1.gml b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/cube_1x1x1.gml
new file mode 100644
index 0000000000000000000000000000000000000000..6ff4684bcbc304c16240037249ddeff6c7b692ea
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/cube_1x1x1.gml
@@ -0,0 +1,44 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml-stylesheet type="text/xsl" href="OpenGeoSysGLI.xsl"?>
+
+<OpenGeoSysGLI xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ogs="http://www.opengeosys.org">
+ <name>cube_1x1x1_geometry</name>
+ <points>
+ <point id="0" x="0" y="0" z="0"/>
+ <point id="1" x="0" y="0" z="1"/>
+ <point id="2" x="0" y="1" z="1"/>
+ <point id="3" x="0" y="1" z="0"/>
+
+ <point id="4" x="1" y="0" z="0"/>
+ <point id="5" x="1" y="0" z="1"/>
+ <point id="6" x="1" y="1" z="1"/>
+ <point id="7" x="1" y="1" z="0"/>
+ </points>
+
+ <surfaces>
+ <surface id="0" name="left"><!-- x=0 -->
+ <element p1="0" p2="1" p3="2"/>
+ <element p1="0" p2="3" p3="2"/>
+ </surface>
+ <surface id="1" name="right"><!-- x=1 -->
+ <element p1="4" p2="6" p3="5"/>
+ <element p1="4" p2="6" p3="7"/>
+ </surface>
+ <surface id="2" name="top"><!-- z=1 -->
+ <element p1="1" p2="2" p3="5"/>
+ <element p1="5" p2="2" p3="6"/>
+ </surface>
+ <surface id="3" name="bottom"><!-- z=0 -->
+ <element p1="0" p2="3" p3="4"/>
+ <element p1="4" p2="3" p3="7"/>
+ </surface>
+ <surface id="4" name="front"><!-- y=0 -->
+ <element p1="0" p2="1" p3="4"/>
+ <element p1="4" p2="1" p3="5"/>
+ </surface>
+ <surface id="5" name="back"><!-- y=1 -->
+ <element p1="2" p2="3" p3="6"/>
+ <element p1="6" p2="3" p3="7"/>
+ </surface>
+ </surfaces>
+</OpenGeoSysGLI>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/cube_1x1x1_hex_1e0.vtu b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/cube_1x1x1_hex_1e0.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..3af224289fd471786c2d3136ee6d1bc4248323d7
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/cube_1x1x1_hex_1e0.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="7" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="96" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.7320508076" offset="120" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="388" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="484" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="508" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _QAAAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAAAAAAAAAAAA==wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAAAAAAAADwPwAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA8D8AAAAAAADwPwAAAAAAAAAAAAAAAAAA8D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8=QAAAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAwAAAAAAAAACAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAHAAAAAAAAAAYAAAAAAAAACAAAAAAAAAAIAAAAAAAAAA==AQAAAAAAAAAM
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..28a53034255f6741da7a61d6539d7991c72d6dc1
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_ts_1_t_1000000.000000.vtu
@@ -0,0 +1,55 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="761" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45" RangeMax="103" offset="272" />
+ <DataArray type="Float64" Name="dry_density_solid_ip" NumberOfTuples="27" format="appended" RangeMin="2600" RangeMax="2600" offset="352" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0.00037416573868" RangeMax="0.00037416573868" offset="420" />
+ <DataArray type="Float64" Name="liquid_density_ip" NumberOfTuples="27" format="appended" RangeMin="1000" RangeMax="1000" offset="1232" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="1300" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="27" format="appended" RangeMin="1" RangeMax="1" offset="1360" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="1.2346729037e-13" RangeMax="5.7323930869e-13" offset="1428" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="2724" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="2796" />
+ <DataArray type="Float64" Name="velocity_ip" NumberOfComponents="3" NumberOfTuples="27" format="appended" RangeMin="6.7764574201e-27" RangeMax="3.9959091286e-26" offset="2856" />
+ <DataArray type="Float64" Name="viscosity_ip" NumberOfTuples="27" format="appended" RangeMin="0.001" RangeMax="0.001" offset="3352" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0.195" RangeMax="0.195" offset="3424" />
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-0" RangeMax="2.6469779602e-23" offset="3520" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="3" format="appended" RangeMin="5.7477486955e-15" RangeMax="1.7048364009e-14" offset="3596" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="7" offset="3912" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.00037416573868" offset="3996" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="6" format="appended" RangeMin="0.00037416573868" RangeMax="0.00037416573868" offset="4132" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-1.2035251336e-12" RangeMax="2.29234456e-12" offset="4568" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="-1.2035251336e-12" RangeMax="2.29234456e-12" offset="4712" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="1" RangeMax="1" offset="4856" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="6" format="appended" RangeMin="1.7901926447e-13" RangeMax="8.0606285363e-13" offset="4944" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="294.15" RangeMax="294.15" offset="5516" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="294.15" RangeMax="294.15" offset="5584" />
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="3" format="appended" RangeMin="8.4948758498e-27" RangeMax="5.7419719217e-26" offset="5652" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="5896" />
+ <DataArray type="Float64" Name="liquid_density_avg" format="appended" RangeMin="1000" RangeMax="1000" offset="5956" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0" RangeMax="0" offset="6020" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="1" RangeMax="1" offset="6080" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="6" format="appended" RangeMin="1.8352047754e-13" RangeMax="1.8352047754e-13" offset="6144" />
+ <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001" RangeMax="0.001" offset="6268" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.7320508076" offset="6336" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="6420" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="6504" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="6564" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAqwAAAAAAAAA=eF6t0TEKwzAMBdC7eM5SAh1ylVKMG6tBYFuupKSEkLvXQyl0iQdnlASPL2kzmBQmdoqUbKZSWcfsVjHDbfsbEntgM/SdSS6CGYw4nb8zzKa05/gAtvS0I8VMCZIW5bJ3FScTk6CubYqyS1Ioted4kAXD8WrXKiJvCAHTZEUZRBoxnKJrIxYINFZO01eVgK8ZvfWQ2s+8oIwnfMvz+ssjFEq6Q+6+fwDcAhd3AQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9a1MDLXTTcxSjVKNTYxAwApJwRbAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYACCgCUODMOMBgB2oyB9AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAPgIAAAAAAAA=eF574Czz+pGZlP0jMK1l/8oo1EBrpbG9m7K/8aaSl7vPWNenzVvwzRrKt34CVf8ETX3Nv7Zn5zf9tYaof2Z9M+f+mbIfB3a/gKp/jKZeLqmqzN5UyIYBCn4Xrb5Z6Pxh9xM097yEqt/PJc9wzOMN3D0Q/h/rR2jugamHqYO5Jzhk/cyuM5+tH6K5B6b+Am9rm/M7bhuI+mPW7xy9i2e+fmr9HM09L6DqIf77sBui7qv1ffUXNe/ecNg8RHPPC9TwhLtnpc+GkxfWsdjcQHMPTL1ToKcwrzoL3D0OGW7NVg2MNg9wuB/ingu7ITSDDcy/sPB8iKZeRbZZbbvmL2tYOEHUX7aGxdd9NPW/O1vnpcXz2ED9DU0PN+Hxhe5+iLnX4PGVrintsHb/D3h8PURTvz1N+ED12o/QeH1jDU0P8Pi6j6a+77KohfZjZqh7HlhD0t9deHyhuweqbndJVPWyVvHP1vqpD/7sucUIjy9098DMnZbgvrkh6ZH1J+HAyXaWv61v4HAPg96WzpfmP6DuP2Ot4F3LJnTni/UDNPOfQ9Xze4g3TDE7Cg1/NptMIbkL9Wm/4fF1D039f5eNP8UmfoXGE7NNmNrWnfpSL+HxdRtN/b6uDaFJz1lsUNPDZoz08AzFvwvg+Qrm/kdo7oGph4THa2vU+H0Bj6/baOqvcLzYtvTCD2tY+oGF03Mc7oGkr527LT+9O3NT5r01JLy+ws1Hdw9E3Xlo+DyGl1c3cLgHFq/QdApLz9YARCNZKQ==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYAACh34HhmFGAwDXLxyWAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEQAAAAAAAAA=eF5jYACBD/YMw4wGADZoH/Y=AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAqwMAAAAAAAA=eF5VU21MU2cURoS6INNtgFIyQ62dLYwWWspn+5RaOysfLQUFmXZRa0YWlkXmHH/AbSbbD4bLRkSiTMQt4CbTMbZF9sErZikxDMMm/hCMkWBUklnWWOeC0MWd2/vem3j/nHtunnPu8zznnC+mWrs7pjdD2+TXfrmtDP2U27qq4Xzx7tenzqnZQETV2XgiR8rR1j3VOu6ogicY8Ey+vx1mqktctgN95StH943lQcSrsKLZOKu7tzDcQ/2+L68Bhan767ejk+pH0upwdk238cGFMsTwx/LmH9fbb2rZcnqPDJXKfB5SXUptJW5Yzv56ZUQj8xFzM9ZaBpOyI1uj+J9MVfARr4sBLyScxMf0y8+KjxL1UXxzgRPv+rVN6owKfEbRb/Vgbsj7uP6Qk+NjMPbP/uf7klXIpL5r4t0IkOC/i0qhp/ql17ZyfTom4nJQV7snXlHqQDX9PznOgb3UNnSnAm9Q3pG5RfZP4hNyH5w8XGeDkvpvzNKz+/SDxAsW7KL8Yjrp+q2hzdsOmc+MwqLbeb4IN8kP/ZlNSCIe+1eXw0BxQO3mfJ5hYiyU9doGkyzHPnSBwuAK4ysyPlPTqW8w5so+ifhncZoGFfG7cZKivc8W1Svgr/tKsur/dEq6uZ4XmNRf4jNOLy6ni/d9Tp7X7mSz5prWhBbye+e/lqifZ/RAC/mU1G/H1Y+P/jUSypDmDL4PuEN6SxOzcIL41I0ooRr1BDXb8vFy5O7cENWLfFL5/qUgRI1bCitRTLxafR6EqX5ve5mEY/X9O/6bj81CYFXw3EumYsTTvl26lo8n9JQsFiNXyFUlct+TXY8e1j5OQ+x3PQc/X8hFF/HoTbUz4Q4+YZmsZ1lMjF0Th3c8Ax8UB4ycfxxuX73SMddlwALxSb68Hm7ag5ZwNvSVwcCtURUuTx9S/pD+xCr6Y8fSg/GaCaUZQ+RfW1EOdMLebUhlr3uDgVWzq/HtrOHUe73ZfE5A9fEj9wob1AhT/yOfFqCI9I4lWJlrhj703xieCfvjD39jw9P7MG/dRLhXT5tRwfm8TYVLEwl8b6as0l1J/B3CnfyuYMJ9GZDLakh3yiMzE/3YgKfnq4aPgK6cPCbcu1JVxnxUf8DgYOvO9xZuMZvk/ZF8quL3ItzjgaP5iMto8n8Vm8D3KzzcuOuticXJjdyvbCzSYJtD6YzKgh2NeSxN0KErYCJOwXWuk+9tD/GZn/ZG+f+4vJpNUt2llAomzZXvqbTP+B9YNDVWAQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAEwAAAAAAAAA=eF5jYBgFo2AUjIJRAAIABRAAAQ==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAAIgCAAAAAAAAUgEAAAAAAAA=eF7LdvsbyBPZb5Wv9blNjHuq1dMVSbN/Ns2xMtquzsjSVLGrl5lxOWdswS6Y+Nm6SZodb2fvMllTHDXbcu6uJ1Bxk9sm+T+tA6zyoObMeVMZysqabZW0ya6571od3ByY+JNlk26w5U+AmwMTn8tv8I3Xq3sXzJx9uixuOw5O3XXy2qqmnqlNcHNg4q+tveXEDhXDzYGJZ0P9tTl4/mX7zkqr7viMEyevV8H9teefuOaUjWW7YOIwf3lfKeaU2j0RLg7zF8ycR5/8TLiSS3fB/AUzByYO8xfMHJg4zF8wc0InOvaUp0+E+wtmDkwc5i+YOTBxmL90GN2t3vY07LL9FJOnqde3C+av9z/1Vsux1cLFYf5q/Xj+YdKnCrg4zF8wczjz1nZ6q/TA/QUzByYO8xfMHJg4zF8wc668/LskV6oL7i+YOTBxmL9g5sDEAQR2QAU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAFAAAAAAAAAA=eF77s/LjJd+kAPs/w4wGALBrfYM=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAJgAAAAAAAAA=eF5T0InpP/T1hP1CXQjNoAXlQ8XfqUNoBSh/gTaEFoXSADi0IcY=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF5jYAABB2swxWCARjM0YOcjaABqhAO9AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAywAAAAAAAAA=eF4BwAA//yGOMAt6IBG9iSKGoEA0Ab3Q85+XneLQvIwe8uyJ0uG8/yFwyPECAT0yejjVeb7uPEE4s89WJhA9XRN8qIAD5bzU0H/jYq6tPL4dyQBP6O28IiwBOWv07jxsY/itpRARPUXaHsNqHQI9b4BG/DRv4Dze/cD/2aLTPAc1Pg341uq8pcC8qNLt8jyS6BvZ6PLmvHo4enLy5Pc8OJ91jCzp8bycNveVe2P9vPTcRY1S7vK88NkBZTCk67x5cI7oMjEGvQhWZ24=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNBKWZoTQLlGaF0mxQmh1KAwAC4AAdAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAQwAAAAAAAAA=eF5jYMAOrjjLvH5kJmWPQ5rhJlheCy7/Bqr+KZo4OnhmFGqgtdLY/jWa+S+g4jD+Y6g5z6Hi96DqYeKvoeIATO0otg==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAJgEAAAAAAAA=eF676Czz+pGZlP01MK1l/8Qo1EBrpbF9iNjiL6pu83eXt69td1i7Z7dektkFWZ2W3Veg6p9C1b+Aqhe3OnjWLo/Dhi/irdks2ynW2md6ZecpsNu8haq/BVX/FKr+297mo2uj2W3mJkh89klhsQncpNZg1tC0+w1U/TOo+rdQ9U8/GG66+1Dcxjn547n2/ke7b8YzGO1q/7f7OVT9A6j6R1D13MUzC/ceWWB9JNp3/7Gcx9aL9qoldux/b/0Sqv4+VP1LqPofX9iu7rP6vfvyv1nPNZb+tF5/XzsnIIPP5h5U/SOo+mdQ9bLCb7bezJtjfdXaYknYUS4bydm1uScXMdjcQlP/Cqp+pb/GlC/Lf+8O1jKQPXL3q7XdFq8NSQz3rAEV7K+0AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAASwAAAAAAAAA=eF4BQAC//yZ7ELMz0XC9vmD3s6Dscb3ON0rgsmdtvdjnaqvmKYQ9aOkVeoWEQL2Ynd5Q0MJ6PZhg+8WRwHE9PAwEAjIsdb17liFVAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAASwAAAAAAAAA=eF4BQAC//yZ7ELMz0XC9vmD3s6Dscb3ON0rgsmdtvdjnaqvmKYQ9aOkVeoWEQL2Ynd5Q0MJ6PZhg+8WRwHE9PAwEAjIsdb17liFVAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAIQAAAAAAAAA=eF778x8E3tszgMEHe0YoDeP/gsr/g9Iw+f9QPgDaISFWAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAiwEAAAAAAAA=eF4BgAF//m6ik6Xzg1m9s7CRco4iW734v4bDsW5UveZ+PH5G8da8FzYSq4IJ9bzBLXQbv5i9vPO59vruADm9KeOhmFJRRD0HgdSmmPQ5PS5jfkcW7EE9+UsP2GayzTzensXWVvdAPe97hX2Gsj495TiWZo+WR70U7oOs2eIjvYe5dI2qIkE9VmpD6y3rPT252dsd7iS7vJMrwSNam149fV4hqfAqXT32LaETP39gPY4LZvkahVE9FuLUbGc7G707m2pZlWU2vSt9wfWBl049ExXeulpmRT2c3a/xskopPckz8kZFmtg8k4O+lAOAHD16OPuNUkonPWc8VLr2uFU9zEPBFeAJUj3HZQ8krD5RPR6bNli1+TO9v4nsKSzQMj2Aj5/wiU5GPflIyzlAPge9MAV0S8LMMz3tg4FOiaWxPCvDUJy729Q8G4goQZSjQz1eLVxtdGg4PTz4XbaIvEI9UiZy0cfpTz3oA8SXwKhWPYhtdwbTSTS97J/EKH8KMD3LOQEXTAIWPWMXtpU=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAlAAAAAAAAAA=eF47fzm6Qu+dh1W1YNc0sfPhuzyz08+oHG7ZdRgqrvfKwfVG3ZJdNy2DJp8qnL8rWmv1tWUxi3dNEoCoX/3g+er2a9N2RUDF9aHql+b1W+QuXmo1dWn0zmMeU3apzbCMO/u4e5cz1PxZUPH0hd9tqu4usLoFNd87c0O+SMpsKwuo+g1Q872g4ilQ9auh5gMAc/NmTA==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBD/YAAicBMA==AQAAAAAAAAAAgAAAAAAAADAAAAAAAAAAOwAAAAAAAAA=eF4BMADP/z5rszZrFjo93/It33IFOT0goh7q4dI9PToemQYeEik99MFPWz3SKT1kIpDQ514oPdYiEz0=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77v/LjJd+kAHsAGZgEqg==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAHAAAAAAAAAA=eF5jYMAHPtjjlcaQh/GJ1YdLPy51mDQAp2EONQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQLlGaF0uxQmg1KAwAC8AAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jAQAADQAN
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_x45_ts_1_t_1000000.000000.vtu b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_x45_ts_1_t_1000000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b6ceb60cbbce434c8ef218b58bc341c12184dd0e
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_x45_ts_1_t_1000000.000000.vtu
@@ -0,0 +1,55 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="761" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45" RangeMax="103" offset="272" />
+ <DataArray type="Float64" Name="dry_density_solid_ip" NumberOfTuples="27" format="appended" RangeMin="2600" RangeMax="2600" offset="352" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0.00017320508076" RangeMax="0.00017320508076" offset="420" />
+ <DataArray type="Float64" Name="liquid_density_ip" NumberOfTuples="27" format="appended" RangeMin="1000" RangeMax="1000" offset="1276" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="1344" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="27" format="appended" RangeMin="1" RangeMax="1" offset="1404" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="9.8776550476e+39" RangeMax="5.5224861524e+40" offset="1472" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="2420" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="2492" />
+ <DataArray type="Float64" Name="velocity_ip" NumberOfComponents="3" NumberOfTuples="27" format="appended" RangeMin="3.5579952679e-35" RangeMax="5.7967338587e-34" offset="2552" />
+ <DataArray type="Float64" Name="viscosity_ip" NumberOfTuples="27" format="appended" RangeMin="0.001" RangeMax="0.001" offset="2948" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0.195" RangeMax="0.195" offset="3020" />
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-4.9568952744e-31" RangeMax="1.2270242447e-30" offset="3116" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="3" format="appended" RangeMin="1.9060767508e+39" RangeMax="7.0368358581e+39" offset="3236" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="7" offset="3552" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="3" format="appended" RangeMin="0.00014884019815" RangeMax="0.00042309589734" offset="3636" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="6" format="appended" RangeMin="0.00017320508076" RangeMax="0.00017320508076" offset="3848" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="8.2807366724e-06" RangeMax="8.2807366724e-06" offset="4212" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="8.2807366724e-06" RangeMax="8.2807366724e-06" offset="4288" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="1" RangeMax="1" offset="4364" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="6" format="appended" RangeMin="1.9901061074e+40" RangeMax="6.3422077943e+40" offset="4452" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="294.15" RangeMax="294.15" offset="4908" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="294.15" RangeMax="294.15" offset="4976" />
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="3" format="appended" RangeMin="6.9553239456e-34" RangeMax="9.390011918e-34" offset="5044" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="5268" />
+ <DataArray type="Float64" Name="liquid_density_avg" format="appended" RangeMin="1000" RangeMax="1000" offset="5328" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0" RangeMax="0" offset="5392" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="1" RangeMax="1" offset="5452" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="6" format="appended" RangeMin="2.3896067146e+40" RangeMax="2.3896067146e+40" offset="5516" />
+ <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001" RangeMax="0.001" offset="5640" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.7320508076" offset="5708" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5792" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5876" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5936" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAqwAAAAAAAAA=eF6t0TEKwzAMBdC7eM5SAh1ylVKMG6tBYFuupKSEkLvXQyl0iQdnlASPL2kzmBQmdoqUbKZSWcfsVjHDbfsbEntgM/SdSS6CGYw4nb8zzKa05/gAtvS0I8VMCZIW5bJ3FScTk6CubYqyS1Ioted4kAXD8WrXKiJvCAHTZEUZRBoxnKJrIxYINFZO01eVgK8ZvfWQ2s+8oIwnfMvz+ssjFEq6Q+6+fwDcAhd3AQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9a1MLLQTTcySTW3TEtOBAApuwTCAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYACCgCUODMOMBgB2oyB9AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAYAIAAAAAAAA=eF4L8IoP9nu9ePc9Z5nXj8yk7GF0VtD8VNOQOdZ3ofwTe3b6h996svvAxh92W05G776Ppt57dtiN3Uo7d8PU/+WYJKQr93i36f58/s16i6zRzZddZ+gTWXMfrv7SnSq3iWcf7GYAg5IdMHUPoLTILPnTS+J7rO9A+cZ/pQTXHb8GVZ9hcR9N/UyZ062tlod2w9Q/e/do1rGq6xD1Du0W6ObbRPQ+yOhaBFcfHlJsWMB1abcPJHys0dV3ekifXSmyGa7+uP+eA1Vmq3fLbAKHjzW6eza/45Yu+r8R7v7eP4GvFdYt2a0JCZ/d6OYbl86c33dhAVx9p0f9iobQDdYB0PiCqXsEpasyd/d8VlwBd+eiBj2v946I+EJXL2bobDmhaRfc3mtn2HiPtj6Cx9d9NPXTdkdckol+AFe/ZsF7B01TRHzBzIfRp5Zt5jnEvxDuHvT4Qldfl+64Lf/lIbj56PGFHp67F3QvyY5cA1cfhhZf6OY/9P+qc/HNBnh6E+bc0VMbtwkeX+jqn0fI75o7fwM8/b/ax+qy0GQdPL7Q3XNmU06ePlsPXP3mc14/GlSWY8QXzL33pwQbP+PZCA+fGpfFLu0rH8Pj6yGa+oz78WcD+HbCxTdOqGaPffsQHl/o5sstmPGP5wcivnhPnihrsryPM7740+c/m8W+Gu5+9Ph6iKY+ZdP2Z5EzD8P9hR5f6ObP++JiYg7MLzDzCcVXl0D8s4WWiHzh+Epw517F7fD4QnePWVvqpkn16+Diiw/EnJtwdzM8vtDNz7gn4bw9uwxu/joJPXtGm9nWAMp1WY8=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYAACh34HhmFGAwDXLxyWAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEQAAAAAAAAA=eF5jYACBD/YMw4wGADZoH/Y=AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAApgIAAAAAAAA=eF7jtmifO+eUxIl2/65ei8n+cDqpf3E1/0JBj6mSH1POMgWd2GC4fLrsl8gTEbdnFN7zuni8fUdf4PIgc7j6Ho/dZRs89E/A1B8XOVrqdzr8xEKZkPMOC8U80M0vC21a9n16MFz988YbFftehpwQ3fBZNu7JiyMwdYx638uVGg08OqzOmeiuYPQQfbEi7Th3wonEmT55ClE+J57lqIRfv3fXBeYemPoJ/3vf71hndgKm3u943Ep2S98T9Z7Zcx68/+OCbv70hVHxqcLicPWvky2C5UvcTrBBwscDXX3nvd7m8v3qcPVtu9+bObxXODEJEj7u6O7pzzd9XpCmBnf/z7nX1ovtlTwxCRI+J9DN/zvjMlv5WxG4esvN86f6h6h4cEPjC6bu6L4DTbpcER4cqSvuSXDKefwRb34SX6B/Iu7IZ8U8M0R8oavnMX7x2rPe4ARMfclyQYb5vEjxBXU/TH39ky2c3tEhcPVTRVSZ+XUQ8QUzH0a//pT9V95GFO4e9PhCVx/WJ7FgVac53HyM+EILT4GTjid9lyjC1b9Ciy9086MXKussFlCFpzcpte4nu++oweMLXT1zkEJb7wYVePqfsLONlf2MEiK+0NzzOCB0eWMBI1y9z0SOC8vPymDEFyyew+deU5hooQoPH55whzX7KiLg8eXatU+hTMsHrv5v6Nl7bvr6J3S4DK+UWUd4KFnFdb33D4PHF7r57cXPc+bJhcLDx3l5YW98XjDO+PJMsdggckkB7n70+IK5B6ZeRnZxwJwUixM6yxdcz/J18UCPL3Tzwx2vPXqxTx0e/oTi66igkOOddWJw96+ZoDbJ/JQOPL7Q3VNuwq266peSBzx8VJQ2N+3UgMcXuvmXHT79XvX89XGY+ep7LpZs+CTgAQC+Zk/RAQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAEwAAAAAAAAA=eF5jYBgFo2AUjIJRAAIABRAAAQ==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAAIgCAAAAAAAABgEAAAAAAAA=eF47+v2t35bJv3bMORQh8mXZrx0nku+vr9j7a0fWpbwtCQeXW8zWXtg2u2WbBUycSV27bZPVHwvWqLtWAZl/4OIPJczPBcostYCZwwAFEP3H4ObAxDXYzopybLkMNwcm/vhN3Ub2u7/h5jSVeWYcUv1j8Z2PzTIt/jLcHJi4cN3iQvEPH3fAzIGJH4X6C2LPYrg+VH8h3AXzF8xdMHGYv1DNOYaThulHNeeYBcxfqOZcQ/PXMbg4zF+o5lyD+wvmTxgN89eWBz+WlDy+ABeH+WvK2YPB6rc/wvXB/AXjw/TB3AHjw2iYO2DmwMRh/kI158MOmL9QzfmwA+YvVHM+7AAAk9tIMw==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAFAAAAAAAAAA=eF77s/LjJd+kAPs/w4wGALBrfYM=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAJgAAAAAAAAA=eF5boB7Tf+jrCfuF2hC6BEorQuk2LQi9AMpfAFVvrgOhAVj0Imw=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAANwAAAAAAAAA=eF5jYACBqZZgqmHBzofP/E93SC/ZCeEvtWQS2G7SodgGkWfogIgzTNnZs8fll9DjHZYADkISRg==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAywAAAAAAAAA=eF4BwAA//2VMKnTVEApIXIhZ0PJqK0gpVqVCGfvdxww5bbhxsf5HWSWmL6EVEkgry0D+D6ckyHQwFfnvkgVIoJuYRzwILcjQZdlzX/UsSClwhEiGJRBI7oefVCFvDsjjow2o9/roxwdjwPsBzRTI5i+Xyy3iIkhd69M2Ns0XSFf0hKktkg/IhMRbch/GH0guO7DwGAQTSBzb+tYFJ+jHAFGBKI01Kchjco4OPxogyH+IusOXEQDIKb1fT97CEsgsbWxsXzoWyHqGVgQ=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNBKWZoTQLlGaF0mxQmh1KAwAC4AAdAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAfAAAAAAAAAA=eF5jYIAAbcdbehMc1OwzvreqXOdQ3x8h87o8YuG63ejiUOUMfYYmlxjWGdqvn/n4/UcZU/s0iHrrXqj4Oqg4TP3upWGv3E8r7P+w6n3s6VvC+8X8eRxX5W+wRheHqY+Rc1m2ue7vvuWrQMDaXgqifnc4VHwZVBwAKjVPMQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAA8AAAAAAAAAA=eF6LlHldHrFw3e57zjKvH5lJ2d+F0lI3XYIO2h22vgPl9x8rtGyZssk6Fk39Ayg9/+/37e3WZ6xvQfkzon7uDrc6aJ0BUW8NUwfT19F6K6vq0+Pd96H87teF0r9fP9ydBFUPE38Ipdfvuu878enD3TB+AqdTM8OLx7vF/HkcV+VvsIaZC3Mvx0aBqt0eG3bfhvLnSs1OeL11zW5pqHqYO2H2fD/dmuDuun/3NSi/VtzscURwhbUURP1umPth9G1NP6c3cafh7myte3KBp/fZbnmoelg4wmjzxd6xLeKHrW9C+QFztlgIdb7YDQBrh7doAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF7L2OBwYn3sQ7sGNDoBBx9dHABndCKhAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF7L2OBwYn3sQ7sGNDoBBx9dHABndCKhAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAIQAAAAAAAAA=eF778x8E3tszgMEHe0YoDeP/gsr/g9Iw+f9QPgDaISFWAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAANAEAAAAAAAA=eF6r93idvY5f+cT7Y28S7s5x8kgKef37xykLD7HE6EO3Ppt7FCw6ckLsvYXHxmn7pv+UV/dohqqfefbn2tdpASfSoOonvcz8vXmio4ezvlvfhbcpJ3796XFSKDT16ISo9/gLNf9kS0hkyA+lE1bRUpc6L0ecaNq58LJ6pZXHmUA/HnO/sBNtUPUToebvgqpftnn9cmu3sBP7bTuv6TMGnTDZ9/Nl/bGIE3x1p/Pul6h4bNAR0l6m5+0RCXXPnpPhB7JcVE6sh7p/3fL5O+7eVTwhAlVfaXpuIp+Nx4lEqPqgt3p1m0qNT9hD3V+5Z7f3ZLd37lD1JzZDzT8PdY/P3Pnm2RyOHh1Q9ytvErmvLhhzQgKqvgtq/l6o+ri4vGMsh809zkPd3yZ3YE9BZ9wJAAKMwnY=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAhwAAAAAAAAA=eF4r8+fXvnucYaf1gomdZp8Zdm4v8Xbe68uwswIq/nt7/hqnHWyWhcIRzhN82Sx/e1zWZK5is/SDqs9DE2faAVGft+Puf+0Ytp1T3ysUzvRgs2xkfLTtEpDeATV/NlT8Y5nG8fUubDth5neI6t0r9WHb2QNVXwoV74WKf4CqL4WaDwAVz0/KAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBD/YAAicBMA==AQAAAAAAAAAAgAAAAAAAADAAAAAAAAAAOwAAAAAAAAA=eF4BMADP/3Icx3Ecx+FECLqFw8R8Fsj7I4KenZ0eSCaQqfX3RCnIqgL+pWCBR8gEk7+3jdtIyFKRGRk=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77v/LjJd+kAHsAGZgEqg==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAHAAAAAAAAAA=eF5jYMAHPtjjlcaQh/GJ1YdLPy51mDQAp2EONQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQLlGaF0uxQmg1KAwAC8AAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jAQAADQAN
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_y45_ts_1_t_1000000.000000.vtu b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_y45_ts_1_t_1000000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..8b4a0e2877634279fb46800bb94638026a72a198
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_y45_ts_1_t_1000000.000000.vtu
@@ -0,0 +1,55 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="761" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45" RangeMax="103" offset="272" />
+ <DataArray type="Float64" Name="dry_density_solid_ip" NumberOfTuples="27" format="appended" RangeMin="2600" RangeMax="2600" offset="352" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0.00017320508076" RangeMax="0.00017320508076" offset="420" />
+ <DataArray type="Float64" Name="liquid_density_ip" NumberOfTuples="27" format="appended" RangeMin="1000" RangeMax="1000" offset="916" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="984" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="27" format="appended" RangeMin="1" RangeMax="1" offset="1044" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="1.4244373766e+40" RangeMax="5.9939370414e+40" offset="1112" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="1720" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="1792" />
+ <DataArray type="Float64" Name="velocity_ip" NumberOfComponents="3" NumberOfTuples="27" format="appended" RangeMin="3.5579952679e-35" RangeMax="5.7967338587e-34" offset="1852" />
+ <DataArray type="Float64" Name="viscosity_ip" NumberOfTuples="27" format="appended" RangeMin="0.001" RangeMax="0.001" offset="2248" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0.195" RangeMax="0.195" offset="2320" />
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-4.9568952744e-31" RangeMax="1.2270242447e-30" offset="2416" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="3" format="appended" RangeMin="2.8599841533e+39" RangeMax="7.4937003109e+39" offset="2536" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="7" offset="2852" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="3" format="appended" RangeMin="0.00013220274128" RangeMax="0.00056730680619" offset="2936" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="6" format="appended" RangeMin="0.00017320508076" RangeMax="0.00017320508076" offset="3144" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="8.2807366724e-06" RangeMax="8.2807366724e-06" offset="3556" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="8.2807366724e-06" RangeMax="8.2807366724e-06" offset="3632" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="1" RangeMax="1" offset="3708" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="6" format="appended" RangeMin="2.9194269961e+40" RangeMax="5.7626840537e+40" offset="3796" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="294.15" RangeMax="294.15" offset="4332" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="294.15" RangeMax="294.15" offset="4400" />
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="3" format="appended" RangeMin="6.9553239456e-34" RangeMax="9.390011918e-34" offset="4468" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="4692" />
+ <DataArray type="Float64" Name="liquid_density_avg" format="appended" RangeMin="1000" RangeMax="1000" offset="4752" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0" RangeMax="0" offset="4816" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="1" RangeMax="1" offset="4876" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="6" format="appended" RangeMin="1.6298630307e+40" RangeMax="1.6298630307e+40" offset="4940" />
+ <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001" RangeMax="0.001" offset="5064" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.7320508076" offset="5132" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5216" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5300" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5360" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAqwAAAAAAAAA=eF6t0TEKwzAMBdC7eM5SAh1ylVKMG6tBYFuupKSEkLvXQyl0iQdnlASPL2kzmBQmdoqUbKZSWcfsVjHDbfsbEntgM/SdSS6CGYw4nb8zzKa05/gAtvS0I8VMCZIW5bJ3FScTk6CubYqyS1Ioted4kAXD8WrXKiJvCAHTZEUZRBoxnKJrIxYINFZO01eVgK8ZvfWQ2s+8oIwnfMvz+ssjFEq6Q+6+fwDcAhd3AQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9a1MLLQTTcySTW3TEtOBAApuwTCAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYACCgCUODMOMBgB2oyB9AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAUgEAAAAAAAA=eF574Czz+pGZlD0DGGy1fgDln7GuT5u3YNluCH1s912I+P4HKOo7rB9C+Tdz7p8p+7HBGqL+2W6ouv33kdU3LLJ+BOV36j/3WcG1w5oBCh5gNX8LXD3E/AO78avvQFO/wXo/lzzDMY83u+8RcE/Nv7Zn5zfthbr/mDV28zfDw8fy07szN2X2Q8PnGnbzge6BqYe4Z8JuCP0BZ/jA1PN7iDdMMTsKDx+Y+aju2UFifG3ZjeoeWHwdw+Gec2jqJ1hD6Adw9cS6B0f87kZVPw3uHuzxhXAPJF7n7CbkHtT0sGA3XvOB7nmEaj7B8IGph7kf5h7s8XXEGqYf1T3XcMTXgt1o6vG7h+EsXD0kfa4nkJ6PWN9D8e8aAuG5AOYvDPdgD8+zu1HNn0O0e1DTzzMc5iPc46bsb7ypZCWG+3G5Jzhk/cyuM5uh7rkGs3c/AAYKQxw=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYAACh34HhmFGAwDXLxyWAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEQAAAAAAAAA=eF5jYACBD/YMw4wGADZoH/Y=AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAApQEAAAAAAAA=eF5j1PtertRo4OF53z7nj5KWByOUr5178EKTpswJCG1zYqrkx5SzTEFw+Wc5KuHX7/11d+3ap1Cm5ePxQHCijJK8KlRfzAmd5QuuZ/m6nGjf0Re4PMj8hHtm44oCNnGPo/sONOlyRXj8bT77V+qNrgcDFMDUw8wX3fBZNu6JBlw9xHzTE7jUw9yDql7Vw/pJxaGPMxNP/BFvfhJfoO+Byz1L3v+7YVttBHW/jQe6+Q8blp/fv1Qd7n++3y3vO2OMoeHjg2E+zD0w9RD3sJ6A0Kk4wwemPjJs6+JXaVbw8IGZj+pfXaLjCxqeJ1DdA4svGwz33D3PqfXulxOaelYoHXoCe/jjdg+O+D2Bqp4H7h708ER3DyReBU8Qcg9qehA9gct8mHtg6qHmEwwfmHqY+2HuQY8viPkWcPtQ3eODEV8unw89efJS+ASaepzuKVsTkP+d2wnOh6RPZZzpGe4e/65ei8n+JyD+VcQZnnD3QNWjuwc9POHuQTFfkGj3oKafGAzz0d1zYdZNrc3r5DHcj8s9bfYry4KWaEDd4+MBMx8AU1ZAgw==AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAEwAAAAAAAAA=eF5jYBgFo2AUjIJRAAIABRAAAQ==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAAIgCAAAAAAAABgEAAAAAAAA=eF47+v2t35bJv3bMORQh8mXZrx0nku+vr9j7a0fWpbwtCQeXW8zWXtg2u2WbBUycSV27bZPVHwvWqLtWAZl/4OIPJczPBcostYCZwwAFEP3H4ObAxDXYzopybLkMNwcm/vhN3Ub2u7/h5jSVeWYcUv1j8Z2PzTIt/jLcHJi4cN3iQvEPH3fAzIGJH4X6C2LPYrg+VH8h3AXzF8xdMHGYv1DNOYaThulHNeeYBcxfqOZcQ/PXMbg4zF+o5lyD+wvmTxgN89eWBz+WlDy+ABeH+WvK2YPB6rc/wvXB/AXjw/TB3AHjw2iYO2DmwMRh/kI158MOmL9QzfmwA+YvVHM+7AAAk9tIMw==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAFAAAAAAAAAA=eF77s/LjJd+kAPs/w4wGALBrfYM=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAJwAAAAAAAAA=eF5boBbTf+jrCXtFHQjtoQWhF0L5nOoQWgHKXwDlf4CqAwAyeCH4AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAANwAAAAAAAAA=eF5jYACBqZZgqmHBzofP/E93SC/ZCeEvtWQS2G7SodgGkWfogIgzTNnZs8fll9DjHZYADkISRg==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAywAAAAAAAAA=eF4BwAA//+bcWgUWZfDHB/cxGGJgKsg2Yob/IogWyCpZFYGh+RBI5A1t/rYH8kdpwq7SwwQsyA+TaMfa2iVIJlo7jVaSIEhI9aEd0DsHyAU+0PXUjxpIysL0myBtA0ieGR29RygzyCsCTtfaEgxIHRXu27rF9MeL1pmXjBoeSCtQkax58SPISar4jKzp7EfPPXmCwoklSMdoLysfePDHgINXvreoB0gZY8c+Wcw1SFVGDGYslyrIdPPexSZy88dzpv1B7MT5R8mFXGU=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNBKWZoTQLlGaF0mxQmh1KAwAC4AAdAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAewAAAAAAAAA=eF5jYICAztupFgaJEvvLL37adp1Nb/99Z5nXj8yk7GHia7rYA+/K2eyHKmfogoq/dXmyoemXGYb6fGNL4V9cTnD1n0LESte8+WvXvj7FMe21wv4HUPUw8S2VDtNUJprC1X+Fin8zlv12MscAQ/3T0NXZ1h22+wHudU8wAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAEwEAAAAAAAA=eF674yzz+pGZlL0q/+szZ0Rv7L4D5SuYcNfP4yi15o//bnlq2SXr2xDx/feg8l+Xxlhfcqu3fgDlp8Vcn+2+4IT1UT3DJP6pD3ZD1e2/D5UPK2Js2MR6Y/ddKL9MR8D9xffl1tPP9t6+pb3W+iFUPZS2r/D1nxK4vMj6KZRvNJFpdQ7rHus/Xqv5Jly+vvsZVD3MPMFFKk8fKp6DudN+ze+YMJ1P83efFqj04JR9BbMX7h6D0zMfn/m9FGav/fTi1FDdyuO7twhc/DVR8TpcPcx/YeZv7j/dc9YaZp/+hoJfIbtn7uZ4ab3Idf/F3ffRzNesjEhPe7cKHj57T97oUNi8ZLeiwnoH178n4eYDAKhavK0=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF7L2OBwYn3sQ7sGNDoBBx9dHABndCKhAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF7L2OBwYn3sQ7sGNDoBBx9dHABndCKhAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAIQAAAAAAAAA=eF778x8E3tszgMEHe0YoDeP/gsr/g9Iw+f9QPgDaISFWAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAcQEAAAAAAAA=eF5Lvrk0Y/cexxPcxw5a/Pzpe2KhuPc121b3Ew58rYu26b5yf7RA/+Z2M3cP2a9rAhyOanmkQNUrexqYrSj77P656LtWXFigx+MlH9ODl9l5rN94UTdrccgJdaj63OcMtk9DjU7c9uOf8e6p74nbyvP+HvkQcMJ6r8+9I1dkPR7HvnnjqK7kkWRj9tNFPuREPFS9NLOMf4L5c/dpOu84wk8leDT67yvmKDD0OPdm3bmO5b4nYqHqj8VpXq9+pO8xdf2a2JgDTh4Pn79getbuc2I3zw91r4fCJzo36O+Qdkg44Se0KXMKo8OJ81D1pvsLliTESHtc2zXTMsww0mPn9W/uXOa2J/JNjyownvI5EQxVn7r3/f70LgePhLnmc1sanTyebv9+d6a20wnvNOU2P2aBEx6X+H9NFHM7sXPx0zmu7/w8YqHqlQJkvVIPKXjorcxcyvPM0ePuinM+/oclT3S4/lcJ13E4sRWqHgDTLbqzAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAhwAAAAAAAAA=eF4r8+fXvnucYaf1gomdZp8Zdm4v8Xbe68uwswIq/nt7/hqnHWyWhcIRzhN82Sx/e1zWZK5is/SDqs9DE2faAVGft+Puf+0Ytp1T3ysUzvRgs2xkfLTtEpDeATV/NlT8Y5nG8fUubDth5neI6t0r9WHb2QNVXwoV74WKf4CqL4WaDwAVz0/KAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBD/YAAicBMA==AQAAAAAAAAAAgAAAAAAAADAAAAAAAAAAOwAAAAAAAAA=eF4BMADP/ynAYYs6JfPHYdPSV50SDMgGLvd3IoEwSPzaSwPtd/dH6FTPvmXnRcgaqg2riKAiyDw6GAI=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77v/LjJd+kAHsAGZgEqg==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAHAAAAAAAAAA=eF5jYMAHPtjjlcaQh/GJ1YdLPy51mDQAp2EONQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQLlGaF0uxQmg1KAwAC8AAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jAQAADQAN
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_z45_ts_1_t_1000000.000000.vtu b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_z45_ts_1_t_1000000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..da0db902d5ad2d96bf84e268f44e180685facabd
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_z45_ts_1_t_1000000.000000.vtu
@@ -0,0 +1,55 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="761" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45" RangeMax="103" offset="272" />
+ <DataArray type="Float64" Name="dry_density_solid_ip" NumberOfTuples="27" format="appended" RangeMin="2600" RangeMax="2600" offset="352" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0.00017320508076" RangeMax="0.00017320508076" offset="420" />
+ <DataArray type="Float64" Name="liquid_density_ip" NumberOfTuples="27" format="appended" RangeMin="1000" RangeMax="1000" offset="1012" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="1080" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="27" format="appended" RangeMin="1" RangeMax="1" offset="1140" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="1.1481289519e+40" RangeMax="4.5364423482e+40" offset="1208" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="1864" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="1936" />
+ <DataArray type="Float64" Name="velocity_ip" NumberOfComponents="3" NumberOfTuples="27" format="appended" RangeMin="3.5579952679e-35" RangeMax="5.7967338587e-34" offset="1996" />
+ <DataArray type="Float64" Name="viscosity_ip" NumberOfTuples="27" format="appended" RangeMin="0.001" RangeMax="0.001" offset="2392" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0.195" RangeMax="0.195" offset="2464" />
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-4.9568952744e-31" RangeMax="1.2270242447e-30" offset="2560" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="3" format="appended" RangeMin="1.316746706e+39" RangeMax="7.452350766e+39" offset="2680" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="7" offset="2996" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="3" format="appended" RangeMin="7.4223192182e-06" RangeMax="0.00036039450016" offset="3080" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="6" format="appended" RangeMin="0.00017320508076" RangeMax="0.00017320508076" offset="3296" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="8.2807366724e-06" RangeMax="8.2807366724e-06" offset="3672" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="8.2807366724e-06" RangeMax="8.2807366724e-06" offset="3748" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="1" RangeMax="1" offset="3824" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="6" format="appended" RangeMin="1.5362733083e+40" RangeMax="4.3161388225e+40" offset="3912" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="294.15" RangeMax="294.15" offset="4420" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="294.15" RangeMax="294.15" offset="4488" />
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="3" format="appended" RangeMin="6.9553239456e-34" RangeMax="9.390011918e-34" offset="4556" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="4780" />
+ <DataArray type="Float64" Name="liquid_density_avg" format="appended" RangeMin="1000" RangeMax="1000" offset="4840" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0" RangeMax="0" offset="4904" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="1" RangeMax="1" offset="4964" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="6" format="appended" RangeMin="6.4147850379e+39" RangeMax="6.4147850379e+39" offset="5028" />
+ <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001" RangeMax="0.001" offset="5152" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.7320508076" offset="5220" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5304" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5388" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5448" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAqwAAAAAAAAA=eF6t0TEKwzAMBdC7eM5SAh1ylVKMG6tBYFuupKSEkLvXQyl0iQdnlASPL2kzmBQmdoqUbKZSWcfsVjHDbfsbEntgM/SdSS6CGYw4nb8zzKa05/gAtvS0I8VMCZIW5bJ3FScTk6CubYqyS1Ioted4kAXD8WrXKiJvCAHTZEUZRBoxnKJrIxYINFZO01eVgK8ZvfWQ2s+8oIwnfMvz+ssjFEq6Q+6+fwDcAhd3AQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9a1MLLQTTcySTW3TEtOBAApuwTCAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYACCgCUODMOMBgB2oyB9AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAmQEAAAAAAAA=eF677yzz+pGZlP1DKM0ABhOs70H4+89Y16fNW3DNWttnxeQ257W7H0DVPUBRz8DwAE2919ys1XEmG+Dq7+E1/5m176mQBeu0D8PVPyRg/oJs071TD6wi6B4YPyW3RNPEcjkO9zAw3EVTb/Pm/MWLe1fhdA+6+g198XMeL9pgfR+He9DD88BBp7VnSxp2w9Tjcg9M/ZlceRWbtigc5k/Yje6em2FrinJ+z7fGHj4TdqOHD3dV0Jt/G2ZaYw+fCbtR3X/M+ld3sGHd8tk4zEePr2O7iY0vqH93E4ovtPDcjR5fqOYfwAgfbc37unMXL9yNPTwPwP0LU19X/f2Fd/NuHPGFaX7f0n1OM7l24jB/AVz9zZz7Z8p+HLCWNouSmLO1CUd4LoDHF0z9j48z9ViaM+Dq76OpRw2fZdbbHAK166VLCMYXxPwLu9fOmaqcJ7GUyPg6trtw3pV+nokLd2N3DyI9w8z3lDWwWCWKy3zM8MwzLXefbLkAZ3yhu2d+4+Gvbud2wdXfx2E+TP2G1V3+U0J37AYAtKIIfA==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYAACh34HhmFGAwDXLxyWAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEQAAAAAAAAA=eF5jYACBD/YMw4wGADZoH/Y=AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAygEAAAAAAAA=eF5r39EXuDzI/IRr1z6FMi0fD5fPh548ecns8Ue8+Ul8gb6Hdu7BC02aPh5pgQ+f/jFVOsGo971cqdHAA0YzQIHO8gXXs3xdTsDUF8xRfiCkrQJX3+7f1Wsx2f8EdvNjPBz/ft4re8ccrh7mHlzmn87jbnu0VIGge2B8qTh2k+znMhjugclPlfyYcpYpCK5fve7rb80lCjjdg65+7x/f6AVbVDzaoeGJ7h708OywWTbHpfrLcZh6XO6BqY82Trk6qfGCO7r50PA8ge6eTvMDV7mfC3ughw9MPXr4vHMqP7T/J78HevjA1KO638Zj5X5T7ymmghjmw8xDjS+bE8TGF9S/JwjFF1p4nkCPL1T/GmOEj5qN+HzVL6InsIenMdy/MPWGhub3eeUNT6DHFy7zTSv1pizT0sdhvjBc/QPBiTJK8qYe55jrtmp6fnfH7n5heHzB1IusW9f7vPEeXD3MHph61PCR8VjksUJsOttLDPPRwx9ivusJ76Zn4SmOUkTGl82J80zX67nuip5Adw9MPap/XU8wL9zl/NkW03xc4Rmsyl3UayiCM77Q3VOj/eVy+joDRHyhhA/CfJj6M7kraqM89U4AACfO/tY=AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAEwAAAAAAAAA=eF5jYBgFo2AUjIJRAAIABRAAAQ==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAAIgCAAAAAAAABgEAAAAAAAA=eF47+v2t35bJv3bMORQh8mXZrx0nku+vr9j7a0fWpbwtCQeXW8zWXtg2u2WbBUycSV27bZPVHwvWqLtWAZl/4OIPJczPBcostYCZwwAFEP3H4ObAxDXYzopybLkMNwcm/vhN3Ub2u7/h5jSVeWYcUv1j8Z2PzTIt/jLcHJi4cN3iQvEPH3fAzIGJH4X6C2LPYrg+VH8h3AXzF8xdMHGYv1DNOYaThulHNeeYBcxfqOZcQ/PXMbg4zF+o5lyD+wvmTxgN89eWBz+WlDy+ABeH+WvK2YPB6rc/wvXB/AXjw/TB3AHjw2iYO2DmwMRh/kI158MOmL9QzfmwA+YvVHM+7AAAk9tIMw==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAFAAAAAAAAAA=eF77s/LjJd+kAPs/w4wGALBrfYM=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAJgAAAAAAAAA=eF5T0I7pP/T1hL0ilE5Rh9ALoXwWLQitAOUvgMrvhIoDABX2IVo=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAANwAAAAAAAAA=eF5jYACBqZZgqmHBzofP/E93SC/ZCeEvtWQS2G7SodgGkWfogIgzTNnZs8fll9DjHZYADkISRg==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAywAAAAAAAAA=eF4BwAA//4iAVbawtQNIOG3K+Jiepcf/gSJLbN4HyKSJK3UuyB1Io3wg8nl4NMhZGvB6xScCyJUJNi9b1xTINAXVX/3PGUhVL6rJV1D+R5/dOCKxUhPISpNEA8deE0gj2lcCphADSL6La1OtxfxH9+JMBidLCMhC53HOL7r8x1MJI69+IRbIrTR9drTcGsjfTWiBd1kByO5EHo6fIttHTkmilo1nBkhcI1VLTQUGSMe+/cTdiQtIHRxrcTgrMEixyDH8oX4BSIwhVlM=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNBKWZoTQLlGaF0mxQmh1KAwAC4AAdAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAgQAAAAAAAAA=eF5jYICAw8mduTmRkvvTVvxQ+27Luv++s8zrR2ZS9g+qPkw4U2i8PwMqDlXOMIU5ZtP5daft7jlxrp/16LbdA6j6S5eD1qif0Nz/GCoOUy+k9vuqx20luPkw9T+3dD8Q7TTDML9DYEtj7H+2/Xeg5sDcw+jR6b3rmz7cfAD0ElKCAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAA+AAAAAAAAAA=eF677Szz+pGZlP0dKN3Exumc/ffQ7lsQ/n4Np9AoJ5E7u6+opbPuXXNs9z009QWr+xccKK3dDeXv19zhzio1/az1Hbm27302cdb3oepgdMvkM26fp3TBxPcfb+fcnrdw927f7T8TZ8ptt34IVQejM+MWczHY98PVb27Yrnug+L51zMvbfTMPHNp9F6oORjdC3X8bql6hOGHLFpcjux86GM21Lty2G909q9Dcf+xq5nqZ+ZuspzJUJT8xW24N8+8DKN2N5n6DawYLOF9s3L1yXfHKq1yLrGHhAnN/LdT9UHP2X5s5Ny3qx0Xr5vNlX04mrN0NAKETvHg=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF7L2OBwYn3sQ7sGNDoBBx9dHABndCKhAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF7L2OBwYn3sQ7sGNDoBBx9dHABndCKhAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAIQAAAAAAAAA=eF778x8E3tszgMEHe0YoDeP/gsr/g9Iw+f9QPgDaISFWAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAXAEAAAAAAAA=eF5j7v1lfM/Z+UTq5MAXXTnGJ4IZwko6F5mf2NSsv75rs6/Hno4bB50rA0/8bz0QMcPC5gQrVL3drW0VicWhJ0yDaqvtQj8dD8lvllvPFOwx2Xt3wP9MJw81uynS3y5ccm86nC+cUefmEQI1n7Njo8LuLgaPZUf/zeWINz6x93jvuXOGhic6902q8Tum49EMVW8PNV/17127U4XMHvc5ztRXPhA6kVLA/GTupGCP54froq6HmZ+Yt3vKJLub1h6x0U+nPWY39AiEuv+IuXPL9j7VE4+C+4QeWlqemBhje7bjnfaJJVD1zmeil5jfCvEIgbq/KODSK/6pjh5Xlee3a+qqe3yW1lzX9UTGY5HfqstaKaYnIqDmc0Pdz3k0OZX9goOH4vIvS/w81U6Iv2WfO7dEzGMtVL0l1HxTqPsZp8UI/z/remLPBnehubFuHqrCbqeNPJROAABzF7T2AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAhwAAAAAAAAA=eF4r8+fXvnucYaf1gomdZp8Zdm4v8Xbe68uwswIq/nt7/hqnHWyWhcIRzhN82Sx/e1zWZK5is/SDqs9DE2faAVGft+Puf+0Ytp1T3ysUzvRgs2xkfLTtEpDeATV/NlT8Y5nG8fUubDth5neI6t0r9WHb2QNVXwoV74WKf4CqL4WaDwAVz0/KAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBD/YAAicBMA==AQAAAAAAAAAAgAAAAAAAADAAAAAAAAAAOwAAAAAAAAA=eF4BMADP/+2zqftO1AlIJMBhizol88fO5BuwTEgUyPsXg+RfcC9I/Fx1b2/9F0igh/4qDJMXyGryGCE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77v/LjJd+kAHsAGZgEqg==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAHAAAAAAAAAA=eF5jYMAHPtjjlcaQh/GJ1YdLPy51mDQAp2EONQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQLlGaF0uxQmg1KAwAC8AAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jAQAADQAN
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_z90_ts_1_t_1000000.000000.vtu b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_z90_ts_1_t_1000000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..f0ddcc905f44405e822c48cf887af906a278a0f7
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/anisotropic_thermal_expansion/expected_anisotropic_thermal_expansion_z90_ts_1_t_1000000.000000.vtu
@@ -0,0 +1,55 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="761" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45" RangeMax="103" offset="272" />
+ <DataArray type="Float64" Name="dry_density_solid_ip" NumberOfTuples="27" format="appended" RangeMin="2600" RangeMax="2600" offset="352" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0.00037416573868" RangeMax="0.00037416573868" offset="420" />
+ <DataArray type="Float64" Name="liquid_density_ip" NumberOfTuples="27" format="appended" RangeMin="1000" RangeMax="1000" offset="1300" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="1368" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="27" format="appended" RangeMin="1" RangeMax="1" offset="1428" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="4.1941114012e-14" RangeMax="6.38201853e-13" offset="1496" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="6" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="2872" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="27" format="appended" RangeMin="0" RangeMax="0" offset="2944" />
+ <DataArray type="Float64" Name="velocity_ip" NumberOfComponents="3" NumberOfTuples="27" format="appended" RangeMin="6.2942809455e-27" RangeMax="5.1719031468e-26" offset="3004" />
+ <DataArray type="Float64" Name="viscosity_ip" NumberOfTuples="27" format="appended" RangeMin="0.001" RangeMax="0.001" offset="3500" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0.195" RangeMax="0.195" offset="3572" />
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-1.3234889801e-23" RangeMax="2.6469779602e-23" offset="3668" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="3" format="appended" RangeMin="2.9882888951e-15" RangeMax="2.0061766405e-14" offset="3752" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="7" offset="4068" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.00037416573868" offset="4152" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="6" format="appended" RangeMin="0.00037416573868" RangeMax="0.00037416573868" offset="4288" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-2.0395929162e-12" RangeMax="2.7857377172e-12" offset="4724" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="-2.0395929162e-12" RangeMax="2.7857377172e-12" offset="4868" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="1" RangeMax="1" offset="5012" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="6" format="appended" RangeMin="2.9172065238e-13" RangeMax="9.9892981327e-13" offset="5100" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="294.15" RangeMax="294.15" offset="5672" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="294.15" RangeMax="294.15" offset="5740" />
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="3" format="appended" RangeMin="2.2362564883e-26" RangeMax="7.1122240911e-26" offset="5808" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="6048" />
+ <DataArray type="Float64" Name="liquid_density_avg" format="appended" RangeMin="1000" RangeMax="1000" offset="6108" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0" RangeMax="0" offset="6172" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="1" RangeMax="1" offset="6232" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="6" format="appended" RangeMin="7.059853699e-14" RangeMax="7.059853699e-14" offset="6296" />
+ <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001" RangeMax="0.001" offset="6416" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.7320508076" offset="6484" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="6568" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="6652" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="6712" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAqwAAAAAAAAA=eF6t0TEKwzAMBdC7eM5SAh1ylVKMG6tBYFuupKSEkLvXQyl0iQdnlASPL2kzmBQmdoqUbKZSWcfsVjHDbfsbEntgM/SdSS6CGYw4nb8zzKa05/gAtvS0I8VMCZIW5bJ3FScTk6CubYqyS1Ioted4kAXD8WrXKiJvCAHTZEUZRBoxnKJrIxYINFZO01eVgK8ZvfWQ2s+8oIwnfMvz+ssjFEq6Q+6+fwDcAhd3AQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9a1MDLXTTcxSjVKNTYxAwApJwRbAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYACCgCUODMOMBgB2oyB9AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAcgIAAAAAAAA=eF51k19IU1Ecx8dS9odkgbYxUbKSoAUpSer4HdjDRm0lFS3BqAcrEiIjkjSJqEGRWE9GDSMQoYcKs0jDBedS9rIHWWQGo632p417553L1NFIpyvonHNhh7wvlwufe87n9/v+flF71Vyi0WKb//eutIl7WustzxpsAbjRMTQcBIPT5LnfWIryXSOhS/ZNQoLycY53H3v58E4gCfu2H24Yu6xDC5ZDGzrvaQWR8jMcj5vN8Zq8BNal+UCoyoBmcn1G44MCTnI+EuXJuTJ+p9+i8jszEOqMBXp+xyHG+XA8kDpy9J6fEOF8GP/kzKi1oFHTOlWI1P0LJM4nRfk2wxuX+n0eD5W4s1NrGazFYonNZ0S8zyzlrxVuSx/HfkB/Xarlqf4LvvDBtat+sAIFOR/GD6C+rlNzZYjUOQ3NXyu0oqocRSkvcTyrk+SQBbXJZxaieiWv7//l/bSPC7B24sU5R1qt5BXl+Lqz8VUhvIgZn2lzJv2OrJIX82H9Yd7FOUzgGOeTKvIJUk4C9h3hfBhP6gsDmzv2n7SOzx/Hq2XjgIgJ/xmadBmk6dmI1vPxtu8f95xOAJm3AKh2v+6Xm7RKXryP98rW3u7rBWB9HXQdsLt1pUpe3ygvc/PGeLpnSl5JjqfebP5xeWX1sPXWCmZ5yRxv2uyZftShEUj/l7Eps20y+imm5MV8Zot8lpR5ID6ikleS41meKvrQPuEI58N489uDd4WTOcz2mNwzqeTFfNg+kvOew/l8TfvRsgJ0h+Ujcu8K8D5SsQ/UVt/c4duZBl3t1fTo6iIEOR/Gj+wdP/641Ufz8sLFiRastqXgL4zEUeg=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEgAAAAAAAAA=eF5jYAACh34HhmFGAwDXLxyWAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAEQAAAAAAAAA=eF5jYACBD/YMw4wGADZoH/Y=AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAA5QMAAAAAAAA=eF51VHtMW3UUbsyERcAwyzZEuraOtqS83669B9ZSaEnDdbpRVkCirfO1SGd4ZMxMGyYZBDedzsKWlW1u6iQ8xnRzwI5ULWOLy1DcXPfQwabEhIo8ssiWxnpK723iH/afX3+53z3n+77znRshEAjs4gfnVqr7hWk+A0xNu9mqlFTs80kc2w5Ew+j1nY9/IS4A9daxa/tu6fFaq9PTWCyFZk+rc3vhCtzHet2rjAxmDg6ENUeKQBf/+4nObi28NB5eeCRMG8L7/H5/rf1prLIoGuQLxXjU3rj555sS2Fb5+uUH48Xwzkd/KC22XKSynqvaUiA6an/GJnAoGizD9Wqu7pN4Q/350KXhBAh/M2MycSoW2pye1j2oRBPhym1SLKP7ynvZyOEhqCOd6yOHQP2s1A1Y6nWzNRUmdAjV/YoSFs98oLO11ak4nXmc7jQgupaBls2Q2GBR2E+ZIWyE9RoOsjBTnxCbXZWNrsKK9R06GXaLvgp7pGsDNFN/nIuGRTfrjRmVQiW997U4Gz4xRoxYL8pgR+Suw+HHVuDfEYY5No8Fstvb5HoO66hRy2dmlJDuPfISFHU2Vd+dL+J0LoeaV5LihyaNYCTe7/bGw03SkfKpBmoIb/6TCekMziEFopP1p3LOFGI91R08rIX3Aj6p1HCb5nX+Ow2HF3D1E6FO7le1uNR4kvywDbLwPeHk4euW9E5EF4H7UW+3LFOJPL4i9eTwxqJkJBi7/JYBkkinZb0B9DS4hRfzQ7z/O4fZc6+Sn4qZh3BH4LwTA9u/pIZtPobPW/CUAH+np0JrE2A7DW7Cl4k5dCbmJ2FQnxD43PHvBeY7Ld0EjaTbcYOFebr3XCqC3snUzreOi7mcRsGVnVcP+Sd08H98Dh28t2C6HwfBvC2DWrbPrnJrlvJw2mzGHoIVTBnQRr6a84pxd92zgxeqc0O+XimTlZf+kg+BsqvWxqKO/ugFIuTzzOeNx3N7hhQH98y6ErhAZ69cD+8HGjIspzOKzz+KTsz8NVabhSqa/31nOTCBHOw3QTvx2XK2DJ5xatP6lBpuDzLRZF2z8Xz5aqQ19P+k0UI6+Vr542qwcfkJ8lGG8hDkI8bnA7x7ZHCbfDQbn4Iu4jM/oArNU8D9OJ8wYJ/w4UgY5fjn0DOXpADOHqjWHX85nd9LLkeLDJ8fWYy6/0PPIjOyRdFw5/IsE6z3G9MxZt19MSIXOtuP5e3amrWU+5LHnsBAjve64iCK+Hx7Opn7Xs0yygRHymsZa2Hgjd5/5qxKIHjrN2+X4a9k1NC0YGm/rr+gw/3LxuOYNXNc7n5gOu6quz6OksK/M4hOrA==AQAAAAAAAAAAgAAAAAAAABAFAAAAAAAAEwAAAAAAAAA=eF5jYBgFo2AUjIJRAAIABRAAAQ==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAADAAAAAAAAAA=eF5jYBieAAAA2AABAQAAAAAAAAAAgAAAAAAAAIgCAAAAAAAAUQEAAAAAAAA=eF49z00og3EAx/HnsJaXGql5DlIOO2ryMq1/SeKJmFI0OaCd2IlJ1JO0l8aTnsXTw+pJPczhabXETf1RahKWNGpLj3JaNFEulFpc9vtfv4dPffnTYpv1a5H65VzISVVSUVOZa75eoYLia22olclz8KzgM8OsW5ps005jn0iBnoj6NESry537DEbsXVHmnG+Kl9KqQrP+hMOrdzIHXbdkCr/GGnPQSffV9yO3xZwTc2HUmtmlSZG/dxkbFA66uXRwnB6MUjjofPlrUkpot0cy8cw8vKw7JIKv9irOFfubp+j4GrgZeU3l46zjC864R0/W1YvsCw46vuCg4wtOdsLdN9yosS846PiCg46vt5+dFvd7iNhSY3bVu8e+JGEqtvyhUHR8hdO9mjBrsI4vOP3G4UXpLkDwBQcdX3DQ8QWnmC91bM/F2RccdHzBQf8HDI5FWg==AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAFAAAAAAAAAA=eF77s/LjJd+kAPs/w4wGALBrfYM=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAJgAAAAAAAAA=eF5T0InpP/T1hP1CXQjNoAXlQ8XfqUNoBSh/gTaEFoXSADi0IcY=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHAAAAAAAAAA=eF5jYAABB2sGFNpgN5S2RqUZGlBpA2sAf7oEIw==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAywAAAAAAAAA=eF4BwAA//2xhxlQadfW8tyjUcCwjAD3odyBC69nNPG4gvyH+K/U859ws6KJWDr37F9NeZIEEvQPBmKoVxwm9SHvrHFYf17z6tIkIDmkLvdpIU4JxLO08mVIgRsfb1Tyxkrld0z8WPb3ZH+Twt+S81h41PScI/bxUuIsHfCbpPKZBIzDJKcK8aEjgtAPvEz3Pth9LDQ7QvAVdR8SDgAk9rNOn0FPB5Lx5A6bfHSm0PAxltJZy7LG86swnE9o257z2OaR9yt/avDHhX88=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNBKWZoTQLlGaF0mxQmh1KAwAC4AAdAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAARAAAAAAAAAA=eF5jYMAOLjrLvH5kpmWPQ5rhA1heCi7/Aqr+JJo4OnhqFGqgtdLY/hWa+S+h4jD+dag5r6Hid6Hqv0HFn0DFAU0tKKw=AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAJAEAAAAAAAA=eF476yzz+pGZlv1bMC1l/9go1EBrpbF90Cvp5PeBNdaeJW+l23+37P4za9L+s2+m7b4IVX8aqv4FVP1s5RsSf86J2bC1V91dumb77k1/qhJFp0nYvISq/wBV/wqq3utj/7n7tcJ7dnZtuFyuwrmnfMmNqYbzF1vD1J+Bqn8CVb8xKV5hSv5ba5bbXNwGLGI2Slk3LqXfZtvzFKr+ElT9I6j6Tw9evDtdtNb6QXOrqZ00u80j8TnJxpd+Wb+Aqv8CVf8cqj5bZdrP0hVce4QamqYdesG4J2iqIkd/orzNbaj6a1D1r6HqTRbmSXW+4bPxL1niJSF6dffN9o4l/ibce2Dqv6G53+rMqqD9Ntesn7d3Gf3t57Jx1bojJ6UvuQcAUHSuAg==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAASwAAAAAAAAA=eF4BQAC//4/xUUTB8IG9PbKy9g2NZ72cV8517ICIPW/aJvAxXXC9zxQBZpALXD24ggxSbQuAPX1IW0TUOG+9kFo2CqbpNr0q7B5uAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAASwAAAAAAAAA=eF4BQAC//4/xUUTB8IG9PbKy9g2NZ72cV8517ICIPW/aJvAxXXC9zxQBZpALXD24ggxSbQuAPX1IW0TUOG+9kFo2CqbpNr0q7B5uAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAIQAAAAAAAAA=eF778x8E3tszgMEHe0YoDeP/gsr/g9Iw+f9QPgDaISFWAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAiwEAAAAAAAA=eF4BgAF//i84pr3mc2C9brYpTuy24DzfN6Z1ap1QvavNZ0yFxrQ8TxTM/g7JvrxmNkAoVNHQvA2IlR/RwWa911nQnbkIZL2CmqRD+XFSvX2x3iR0u1A99bZmlxBe8byV3irpTwlSPR1tpt+xymM9rynjCqMXYz1C6dCNGt9iPd/VipAZmUy9z5B5zMNxQr1i8Idh/MncPOYNyCRfky49iuxLFb4qWb01v/g8VhxSva2qmkkRmCU9+jr8wK4mUD0TPwzHkMRAvU/+ZhtJZVQ9++HrVblrH71vh8cOLh/wPFAf6qBtmuU8NJrkemXzQD0Dc77MA60zPcJK1muN5WQ9EYHHH4n8ZD3TVTWsFKRbPd0n1AZljkO9o7AuTrhGOr2Vu4IwYgVXPStDiEvWvDq9vHkARifyPL3jLX+RQoNDPbIfIi+Jr0Y9sEuhB5LuDr0mIT/5/PRDvbe/tsk9oi29fEPTIGBIUj1X4A5Qv1A3vdagsLAsUBA9KjCZNnF8Qz0JOIfckHlSvUUPufA=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCpCKHNDJpACIdGJE=AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAkwAAAAAAAAA=eF7b+eZZWuG3CbucnE4u5bq3dtcS99m6R/Xm7doJFS+ccnLLqSuZVu3auo0/787bZZ26aK4cw3KrIKj6L5VOdV57llmZQ8VjoeoVZm19c5G5ddekCxW7li6duGtfrvTezmv9VhOg5s+Gin+fkP1w5bTpVk1Q83nfmuY3fmnedRGq/hvUfHGo+Ceo+iCo+QAsHGplAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBD/YAAicBMA==AQAAAAAAAAAAgAAAAAAAADAAAAAAAAAAOQAAAAAAAAA=eF7rZ8uc4HZSwXZR/s9p65U0bPdN7vpXm6Jga+Tg1zZ5Pr9tZPAMzmx/RVsG3csf1OOFbAHw5hLWAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77v/LjJd+kAHsAGZgEqg==AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAHAAAAAAAAAA=eF5jYMAHPtjjlcaQh/GJ1YdLPy51mDQAp2EONQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQLlGaF0uxQmg1KAwAC8AAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jAQAADQAN
+ </AppendedData>
+</VTKFile>