From 807fed956b7cc8d6e69bace559cee9ca60f7f7c9 Mon Sep 17 00:00:00 2001
From: Jan Thiedau <jan.thiedau@bgr.de>
Date: Fri, 12 Oct 2018 13:10:05 +0200
Subject: [PATCH] BGRa Verification examples by Peter Vogel
---
ProcessLib/ThermoMechanics/Tests.cmake | 14 +
.../Verification/m2_1D1bt/m2_1D1bt.gml | 3 +
.../Verification/m2_1D1bt/m2_1D1bt.prj | 287 +++++++++++++++
.../Verification/m2_1D1bt/m2_1D1bt.vtu | 3 +
.../m2_1D1bt_pcs_0_ts_200_t_1.500000.vtu | 3 +
.../m2_1D1bt_pcs_0_ts_400_t_3.000000.vtu | 3 +
.../Verification/m2_1D2bt/m2_1D2bt.gml | 3 +
.../Verification/m2_1D2bt/m2_1D2bt.prj | 287 +++++++++++++++
.../Verification/m2_1D2bt/m2_1D2bt.vtu | 3 +
.../m2_1D2bt_pcs_0_ts_200_t_15.000000.vtu | 3 +
.../m2_1D2bt_pcs_0_ts_400_t_30.000000.vtu | 3 +
.../Verification/m2_1Dcreep/m2_1Dcreep.gml | 3 +
.../Verification/m2_1Dcreep/m2_1Dcreep.prj | 274 +++++++++++++++
.../Verification/m2_1Dcreep/m2_1Dcreep.vtu | 3 +
.../m2_1Dcreep_pcs_0_ts_200_t_10.000000.vtu | 3 +
.../m2_1Dcreep_pcs_0_ts_400_t_20.000000.vtu | 3 +
.../m2_1Dlozenge/m2_1Dlozenge.gml | 3 +
.../m2_1Dlozenge/m2_1Dlozenge.prj | 314 +++++++++++++++++
.../m2_1Dlozenge/m2_1Dlozenge.vtu | 3 +
.../m2_1Dlozenge_pcs_0_ts_101_t_10.000000.vtu | 3 +
.../m2_1Dlozenge_pcs_0_ts_1_t_0.000000.vtu | 3 +
.../m2_1Dlozenge_pcs_0_ts_51_t_5.000000.vtu | 3 +
.../m2_1Dlozengebt/m2_1Dlozengebt.gml | 3 +
.../m2_1Dlozengebt/m2_1Dlozengebt.prj | 314 +++++++++++++++++
.../m2_1Dlozengebt/m2_1Dlozengebt.vtu | 3 +
...m2_1Dlozengebt_pcs_0_ts_350_t_0.700000.vtu | 3 +
...m2_1Dlozengebt_pcs_0_ts_500_t_1.000000.vtu | 3 +
.../Verification/m2_1Drelax/m2_1Drelax.gml | 3 +
.../Verification/m2_1Drelax/m2_1Drelax.prj | 278 +++++++++++++++
.../Verification/m2_1Drelax/m2_1Drelax.vtu | 3 +
.../m2_1Drelax_pcs_0_ts_100_t_0.100000.vtu | 3 +
.../m2_1Drelax_pcs_0_ts_200_t_1.100000.vtu | 3 +
.../Verification/m2_2Dload/m2_2Dload.gml | 3 +
.../Verification/m2_2Dload/m2_2Dload.prj | 289 +++++++++++++++
.../Verification/m2_2Dload/m2_2Dload.vtu | 3 +
.../m2_2Dload_pcs_0_ts_11_t_1.000000.vtu | 3 +
.../m2_2Dload_pcs_0_ts_6_t_0.500000.vtu | 3 +
.../Verification/m2_2Dload/m2_2Dload_ym45.gml | 3 +
.../Verification/m2_2Dload/m2_2Dload_ym45.prj | 288 +++++++++++++++
.../Verification/m2_2Dload/m2_2Dload_ym45.vtu | 3 +
...m2_2Dload_ym45_pcs_0_ts_101_t_1.000000.vtu | 3 +
.../m2_2Dload_ym45_pcs_0_ts_1_t_0.000000.vtu | 3 +
.../m2_2Dload_ym45_pcs_0_ts_51_t_0.500000.vtu | 3 +
.../Verification/m2_2Dloadbt/m2_2Dloadbt.gml | 3 +
.../Verification/m2_2Dloadbt/m2_2Dloadbt.prj | 307 ++++++++++++++++
.../Verification/m2_2Dloadbt/m2_2Dloadbt.vtu | 3 +
.../m2_2Dloadbt_pcs_0_ts_300_t_0.600000.vtu | 3 +
.../m2_2Dloadbt_pcs_0_ts_500_t_1.000000.vtu | 3 +
.../Verification/m2_3Dload/m2_3Dload.gml | 3 +
.../Verification/m2_3Dload/m2_3Dload.prj | 312 +++++++++++++++++
.../Verification/m2_3Dload/m2_3Dload.vtu | 3 +
.../m2_3Dload_pcs_0_ts_101_t_1.000000.vtu | 3 +
.../m2_3Dload_pcs_0_ts_51_t_0.500000.vtu | 3 +
.../Verification/m2_3Dloadbt/m2_3Dloadbt.gml | 3 +
.../Verification/m2_3Dloadbt/m2_3Dloadbt.prj | 330 ++++++++++++++++++
.../Verification/m2_3Dloadbt/m2_3Dloadbt.vtu | 3 +
.../m2_3Dloadbt_pcs_0_ts_250_t_0.500000.vtu | 3 +
.../m2_3Dloadbt_pcs_0_ts_500_t_1.000000.vtu | 3 +
58 files changed, 3432 insertions(+)
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt_pcs_0_ts_200_t_1.500000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt_pcs_0_ts_400_t_3.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt_pcs_0_ts_200_t_15.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt_pcs_0_ts_400_t_30.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep_pcs_0_ts_200_t_10.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep_pcs_0_ts_400_t_20.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_101_t_10.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_1_t_0.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_51_t_5.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt_pcs_0_ts_350_t_0.700000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt_pcs_0_ts_500_t_1.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax_pcs_0_ts_100_t_0.100000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax_pcs_0_ts_200_t_1.100000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_pcs_0_ts_11_t_1.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_pcs_0_ts_6_t_0.500000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_101_t_1.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_1_t_0.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_51_t_0.500000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt_pcs_0_ts_300_t_0.600000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt_pcs_0_ts_500_t_1.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload_pcs_0_ts_101_t_1.000000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload_pcs_0_ts_51_t_0.500000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.gml
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.prj
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt_pcs_0_ts_250_t_0.500000.vtu
create mode 100644 Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt_pcs_0_ts_500_t_1.000000.vtu
diff --git a/ProcessLib/ThermoMechanics/Tests.cmake b/ProcessLib/ThermoMechanics/Tests.cmake
index 8a66d3f29c9..eb2f809a2ac 100644
--- a/ProcessLib/ThermoMechanics/Tests.cmake
+++ b/ProcessLib/ThermoMechanics/Tests.cmake
@@ -1,3 +1,17 @@
+if (NOT OGS_USE_MPI)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.prj)
+ OgsTest(PROJECTFILE ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.prj)
+endif()
+
AddTest(
NAME ThermoMechanics_3D_ThermoElastic_Stress_Analysis
PATH ThermoMechanics
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.gml
new file mode 100644
index 00000000000..d18b8bdd2df
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8278ccbb5454158303f25bd40f01a6d3b633e694e688d302d4c5f2b29a7cca1b
+size 1778
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.prj
new file mode 100644
index 00000000000..01934d23140
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.prj
@@ -0,0 +1,287 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_1D1bt.vtu</mesh>
+ <geometry>m2_1D1bt.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>3</t_end>
+ <timesteps>
+ <pair>
+ <repeat>405</repeat>
+ <delta_t>7.5e-3</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_1D1bt</prefix>
+ <timesteps>
+ <pair>
+ <repeat>2</repeat>
+ <each_steps>200</each_steps>
+ </pair>
+ </timesteps>
+ <fixed_output_times>1.5 3.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>2.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>2.5e7</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.27</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Dirichlet_Tension_Scaling</name>
+ <type>Constant</type>
+ <value>1.0E-4</value>
+ </parameter>
+ <parameter>
+ <name>Dirichlet_Tension</name>
+ <type>CurveScaled</type>
+ <curve>Dirichlet_top_temporal</curve>
+ <parameter>Dirichlet_Tension_Scaling</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>Dirichlet_top_temporal</name>
+ <coords>0.0 5.0</coords>
+ <values>0.0 5.0</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1D1bt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1D1bt</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1D1bt</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1D1bt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet_Tension</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1D1bt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1D1bt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-12 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_1D1bt_pcs_0_ts_200_t_1.500000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>3e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D1bt_pcs_0_ts_200_t_1.500000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>3.0e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D1bt_pcs_0_ts_200_t_1.500000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D1bt_pcs_0_ts_400_t_3.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D1bt_pcs_0_ts_400_t_3.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>7e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D1bt_pcs_0_ts_400_t_3.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.vtu
new file mode 100644
index 00000000000..86f6c6b4549
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:829bf97d41b20ae6c3b1f74c6595d44d8d2344d2aa11704943ca3850f02cfcad
+size 2738
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt_pcs_0_ts_200_t_1.500000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt_pcs_0_ts_200_t_1.500000.vtu
new file mode 100644
index 00000000000..cb629a00fb2
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt_pcs_0_ts_200_t_1.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:5458e3ce5f698bf01b2931aa484e37a8949bbf7b7d72c7a288fc405f1f25fbb1
+size 9951
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt_pcs_0_ts_400_t_3.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt_pcs_0_ts_400_t_3.000000.vtu
new file mode 100644
index 00000000000..071c679b893
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D1bt/m2_1D1bt_pcs_0_ts_400_t_3.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f1b2955d0d083a5e6860551bdfc2aa664c3a9be0182abeb49661e84d14169ede
+size 9983
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.gml
new file mode 100644
index 00000000000..506ebf1081b
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:486715f4dc5a627fa74ff5dce1f8acc40cd690e8df71ab170d8553521dffaad8
+size 1778
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.prj
new file mode 100644
index 00000000000..007681d7a39
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.prj
@@ -0,0 +1,287 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_1D2bt.vtu</mesh>
+ <geometry>m2_1D2bt.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>30</t_end>
+ <timesteps>
+ <pair>
+ <repeat>405</repeat>
+ <delta_t>0.075</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_1D2bt</prefix>
+ <timesteps>
+ <pair>
+ <repeat>2</repeat>
+ <each_steps>200</each_steps>
+ </pair>
+ </timesteps>
+ <fixed_output_times>15.0 30.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.27</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Top</name>
+ <type>CurveScaled</type>
+ <curve>Neumann_top_temporal</curve>
+ <parameter>Neumann_Top_Scalar</parameter>
+ </parameter>
+ <parameter>
+ <name>Neumann_Top_Scalar</name>
+ <type>Constant</type>
+ <value>-1.0</value>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>Neumann_top_temporal</name>
+ <coords>0.0 31.0</coords>
+ <values>0.0 31.0</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1D2bt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1D2bt</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1D2bt</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1D2bt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Neumann_Top</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1D2bt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1D2bt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_1D2bt_pcs_0_ts_200_t_15.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D2bt_pcs_0_ts_200_t_15.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D2bt_pcs_0_ts_200_t_15.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D2bt_pcs_0_ts_400_t_30.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D2bt_pcs_0_ts_400_t_30.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1D2bt_pcs_0_ts_400_t_30.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.vtu
new file mode 100644
index 00000000000..86f6c6b4549
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:829bf97d41b20ae6c3b1f74c6595d44d8d2344d2aa11704943ca3850f02cfcad
+size 2738
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt_pcs_0_ts_200_t_15.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt_pcs_0_ts_200_t_15.000000.vtu
new file mode 100644
index 00000000000..9bd1ffee703
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt_pcs_0_ts_200_t_15.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c1fd2341bdf1856c39fba11e51fb9cf8bb51760effdc01ad60f8dfe146c37cad
+size 9719
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt_pcs_0_ts_400_t_30.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt_pcs_0_ts_400_t_30.000000.vtu
new file mode 100644
index 00000000000..57f2d9946d9
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1D2bt/m2_1D2bt_pcs_0_ts_400_t_30.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3cfd002bcfba0620325e5c3186c0a73411fb5c5dd06257305dc4953fc4d7af37
+size 9423
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.gml
new file mode 100644
index 00000000000..57e5a731ac8
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:d925934348a2b9a36c9e3c7771ef8157e0f15b53cdb8e988ee84a0f34afdf97f
+size 1780
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.prj
new file mode 100644
index 00000000000..877184fb63a
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.prj
@@ -0,0 +1,274 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_1Dcreep.vtu</mesh>
+ <geometry>m2_1Dcreep.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>20</t_end>
+ <timesteps>
+ <pair>
+ <repeat>405</repeat>
+ <delta_t>0.05</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_1Dcreep</prefix>
+ <timesteps>
+ <pair>
+ <repeat>2</repeat>
+ <each_steps>200</each_steps>
+ </pair>
+ </timesteps>
+ <fixed_output_times>10.0 20.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.27</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Top</name>
+ <type>Constant</type>
+ <value>-20</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1Dcreep</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dcreep</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dcreep</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dcreep</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Neumann_Top</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1Dcreep</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dcreep</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_1Dcreep_pcs_0_ts_200_t_10.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dcreep_pcs_0_ts_200_t_10.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dcreep_pcs_0_ts_200_t_10.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dcreep_pcs_0_ts_400_t_20.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dcreep_pcs_0_ts_400_t_20.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dcreep_pcs_0_ts_400_t_20.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.vtu
new file mode 100644
index 00000000000..86f6c6b4549
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:829bf97d41b20ae6c3b1f74c6595d44d8d2344d2aa11704943ca3850f02cfcad
+size 2738
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep_pcs_0_ts_200_t_10.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep_pcs_0_ts_200_t_10.000000.vtu
new file mode 100644
index 00000000000..4749ee1b89c
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep_pcs_0_ts_200_t_10.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:d79dd328efc2cdf347919686742239097cdfa1c2c7fada3579011ba40c16d914
+size 9771
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep_pcs_0_ts_400_t_20.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep_pcs_0_ts_400_t_20.000000.vtu
new file mode 100644
index 00000000000..40743d8d8cb
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dcreep/m2_1Dcreep_pcs_0_ts_400_t_20.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:7ec7bcd4296faa8db906b1d455fbf1dfc692d7e8010893f89eedc2ebeaae947a
+size 9803
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.gml
new file mode 100644
index 00000000000..090c3e0cfef
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3ed5adf7251cb9f434f049faa3cf17c9626098566aa34486f8d0dd90f4e29811
+size 2134
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.prj
new file mode 100644
index 00000000000..d651d2f2705
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.prj
@@ -0,0 +1,314 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_1Dlozenge.vtu</mesh>
+ <geometry>m2_1Dlozenge.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>10</t_end>
+ <timesteps>
+ <pair><repeat>1</repeat><delta_t>1.0e-10</delta_t></pair>
+ <pair><repeat>105</repeat><delta_t>1.0e-1</delta_t></pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_1Dlozenge</prefix>
+ <timesteps>
+ <pair><repeat>1</repeat><each_steps>1</each_steps></pair>
+ <pair><repeat>2</repeat><each_steps>50</each_steps></pair>
+ </timesteps>
+ <fixed_output_times>0.5 1.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.25</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension_Neg</name>
+ <type>Constant</type>
+ <value>-20</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension_Pos</name>
+ <type>Constant</type>
+ <value>20</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>POINTSxzero</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>POINTSyzero</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension_Pos</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>SURFACE7</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension_Pos</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>SURFACE5</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension_Neg</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>SURFACE6</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension_Neg</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>SURFACE5</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozenge</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_1_t_0.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_1_t_0.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>5e-9</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_1_t_0.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_51_t_5.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_51_t_5.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_51_t_5.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_101_t_10.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_101_t_10.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozenge_pcs_0_ts_101_t_10.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.vtu
new file mode 100644
index 00000000000..4fadd9ff7cc
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8d7638b101ba8b55845b28497aef50e4420c96d63bca5007ffb7539d563468e1
+size 40304
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_101_t_10.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_101_t_10.000000.vtu
new file mode 100644
index 00000000000..8bbd67335e0
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_101_t_10.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:baeaa307830f82a0b2b17654b7d0c39355eb429399e4db9340a0a28e8e05cda3
+size 144741
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_1_t_0.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_1_t_0.000000.vtu
new file mode 100644
index 00000000000..e7a0d3e5322
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_1_t_0.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f184c9f9fd93c78b03f603c9e0f3ac199e8cc833df6b4789d06b7da77176691c
+size 154621
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_51_t_5.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_51_t_5.000000.vtu
new file mode 100644
index 00000000000..a56135730ad
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozenge/m2_1Dlozenge_pcs_0_ts_51_t_5.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:40bfcb0f2a935ed856ad32bab214b974d8e47e336cb0f81a6451be6847f6b8e5
+size 143441
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.gml
new file mode 100644
index 00000000000..9c6585a0c03
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0a56a80ae19e61e081d76b6b8c1379ebf98011550f6245651f0e5ed32d253161
+size 2135
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.prj
new file mode 100644
index 00000000000..fbe7c1b5112
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.prj
@@ -0,0 +1,314 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_1Dlozengebt.vtu</mesh>
+ <geometry>m2_1Dlozengebt.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair><repeat>505</repeat><delta_t>0.002</delta_t></pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_1Dlozengebt</prefix>
+ <timesteps>
+ <pair><repeat>1</repeat><each_steps>350</each_steps></pair>
+ <pair><repeat>1</repeat><each_steps>150</each_steps></pair>
+ </timesteps>
+ <fixed_output_times>0.7 1.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.25</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension_Neg_Scal</name>
+ <type>Constant</type>
+ <value>-80</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension_Pos_Scal</name>
+ <type>Constant</type>
+ <value>80</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension_Neg</name>
+ <type>CurveScaled</type>
+ <curve>Neumann_Curve</curve>
+ <parameter>Neumann_Tension_Neg_Scal</parameter>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension_Pos</name>
+ <type>CurveScaled</type>
+ <curve>Neumann_Curve</curve>
+ <parameter>Neumann_Tension_Pos_Scal</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>Neumann_Curve</name>
+ <coords>0.0 1 2</coords>
+ <values>0.0 1 2</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>POINTSxzero</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>POINTSyzero</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension_Pos</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>SURFACE7</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension_Pos</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>SURFACE5</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension_Neg</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>SURFACE6</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension_Neg</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>SURFACE5</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Dlozengebt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_1Dlozengebt_pcs_0_ts_350_t_0.700000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozengebt_pcs_0_ts_350_t_0.700000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>5.0e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozengebt_pcs_0_ts_350_t_0.700000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozengebt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>5e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozengebt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Dlozengebt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.vtu
new file mode 100644
index 00000000000..4fadd9ff7cc
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8d7638b101ba8b55845b28497aef50e4420c96d63bca5007ffb7539d563468e1
+size 40304
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt_pcs_0_ts_350_t_0.700000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt_pcs_0_ts_350_t_0.700000.vtu
new file mode 100644
index 00000000000..6aa6d704dba
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt_pcs_0_ts_350_t_0.700000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:2ca9ff43784544bfa17f046f7907df54208259f8e3a40f5861aa54e1c66f03ee
+size 151709
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt_pcs_0_ts_500_t_1.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt_pcs_0_ts_500_t_1.000000.vtu
new file mode 100644
index 00000000000..1fe09ea133b
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Dlozengebt/m2_1Dlozengebt_pcs_0_ts_500_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f7780d364415ae8905fef1cd2f3805acae30766ecbe4e62a4dae9e84235b1e47
+size 148157
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.gml
new file mode 100644
index 00000000000..acb2bd0b426
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0085d0991a57dec4b1e4ce8ca2759a457c9f49459804c5dbc50a433ffe9b7696
+size 1780
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.prj
new file mode 100644
index 00000000000..574809a1b0a
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.prj
@@ -0,0 +1,278 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_1Drelax.vtu</mesh>
+ <geometry>m2_1Drelax.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1.1</t_end>
+ <timesteps>
+ <pair>
+ <repeat>100</repeat>
+ <delta_t>1.0e-3</delta_t>
+ </pair>
+ <pair>
+ <repeat>105</repeat>
+ <delta_t>1.0e-2</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_1Drelax</prefix>
+ <timesteps>
+ <pair>
+ <repeat>2</repeat>
+ <each_steps>100</each_steps>
+ </pair>
+ </timesteps>
+ <fixed_output_times>0.1 1.1</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.27</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Dirichlet_Tension</name>
+ <type>Constant</type>
+ <value>1.2E-3</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1Drelax</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Drelax</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Drelax</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Drelax</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet_Tension</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_1Drelax</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_1Drelax</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_1Drelax_pcs_0_ts_100_t_0.100000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Drelax_pcs_0_ts_100_t_0.100000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Drelax_pcs_0_ts_100_t_0.100000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Drelax_pcs_0_ts_200_t_1.100000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Drelax_pcs_0_ts_200_t_1.100000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_1Drelax_pcs_0_ts_200_t_1.100000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.vtu
new file mode 100644
index 00000000000..86f6c6b4549
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:829bf97d41b20ae6c3b1f74c6595d44d8d2344d2aa11704943ca3850f02cfcad
+size 2738
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax_pcs_0_ts_100_t_0.100000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax_pcs_0_ts_100_t_0.100000.vtu
new file mode 100644
index 00000000000..02275841052
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax_pcs_0_ts_100_t_0.100000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:45e2ce8d50175a9e9b1520d6e0e907abee206b8f8afc1292288d8e903acac39c
+size 9631
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax_pcs_0_ts_200_t_1.100000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax_pcs_0_ts_200_t_1.100000.vtu
new file mode 100644
index 00000000000..aa6205574aa
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_1Drelax/m2_1Drelax_pcs_0_ts_200_t_1.100000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:e6bfdf68d17f9b5867d8c958b659e5df68dcfb2a749aa775f13d19594040b0b9
+size 9815
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.gml
new file mode 100644
index 00000000000..56f8f4109c1
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:af04549ddfb506a89a7a6f880c4c4d21110cc8318f2c3772af70ca36edbc6553
+size 973
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.prj
new file mode 100644
index 00000000000..c5d4726b77d
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.prj
@@ -0,0 +1,289 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_2Dload.vtu</mesh>
+ <geometry>m2_2Dload.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair><repeat>1</repeat><delta_t>1e-10</delta_t></pair>
+ <pair><repeat>11</repeat><delta_t>0.1</delta_t></pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_2Dload</prefix>
+ <timesteps>
+ <pair><repeat>1</repeat><each_steps>6</each_steps></pair>
+ <pair><repeat>1</repeat><each_steps>5</each_steps></pair>
+ </timesteps>
+ <fixed_output_times>0.5 1.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.25</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension</name>
+ <type>Constant</type>
+ <value>50</value>
+ </parameter>
+ <parameter>
+ <name>Minus_Neumann_Tension</name>
+ <type>Constant</type>
+ <value>-50</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_2Dload</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dload</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dload</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dload</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dload</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dload</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_2Dload</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dload</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_2Dload_pcs_0_ts_6_t_0.500000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_pcs_0_ts_6_t_0.500000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_pcs_0_ts_6_t_0.500000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_pcs_0_ts_11_t_1.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_pcs_0_ts_11_t_1.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_pcs_0_ts_11_t_1.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>2.0e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.vtu
new file mode 100644
index 00000000000..0be925e4015
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c88f661790c07a35d3a2c6addc9790501055326596bffa980b9764c49526d33d
+size 21976
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_pcs_0_ts_11_t_1.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_pcs_0_ts_11_t_1.000000.vtu
new file mode 100644
index 00000000000..3f1e549ab72
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_pcs_0_ts_11_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:fe4d6048836e641e13ab13a4b8669dfea0a834802527be70d3cff9a2de264802
+size 93125
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_pcs_0_ts_6_t_0.500000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_pcs_0_ts_6_t_0.500000.vtu
new file mode 100644
index 00000000000..90c406aba1a
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_pcs_0_ts_6_t_0.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:048902ffb41ea98f9ae89c1cb2d2a02c419d5abbb24f5b941675f836a544442f
+size 93037
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.gml
new file mode 100644
index 00000000000..2d3d3ba37ed
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:cae792fbd173713d8490e7ae657502672aac71f7ad492d6acc689e0879c07e58
+size 1450
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.prj
new file mode 100644
index 00000000000..a74382338aa
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.prj
@@ -0,0 +1,288 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_2Dload_ym45.vtu</mesh>
+ <geometry>m2_2Dload_ym45.gml</geometry>
+ <processes>
+ <process>
+ <name>SD</name>
+ <type>SMALL_DEFORMATION</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <solid_density>rho_sr</solid_density>
+ <reference_temperature>273.15</reference_temperature>
+ <process_variables>
+ <process_variable>displacement</process_variable>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="SD">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair><repeat>1</repeat><delta_t>1e-10</delta_t></pair>
+ <pair><repeat>105</repeat><delta_t>0.01</delta_t></pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_2Dload_ym45</prefix>
+ <timesteps>
+ <pair><repeat>1</repeat><each_steps>1</each_steps></pair>
+ <pair><repeat>2</repeat><each_steps>50</each_steps></pair>
+ </timesteps>
+ <variables>
+ <variable>displacement</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+
+ <parameters>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.25</value>
+ </parameter>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension</name>
+ <type>Constant</type>
+ <value>35.35533905932738</value>
+ </parameter>
+ <parameter>
+ <name>Minus_Neumann_Tension</name>
+ <type>Constant</type>
+ <value>-35.35533905932738</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <!-- komplette Fixierung -->
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE5</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>plate_ym45</geometrical_set>
+ <geometry>SURFACE5</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Neumann_Tension</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_1_t_0.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_1_t_0.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_1_t_0.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_51_t_0.500000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_51_t_0.500000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>2.0e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_51_t_0.500000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_101_t_1.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_101_t_1.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dload_ym45_pcs_0_ts_101_t_1.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>2.0e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.vtu
new file mode 100644
index 00000000000..ad6cd0de0f6
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:6a3f7a0e51b80144c37c39208979146b371aae471a9d37f04e363ae71d87bd8e
+size 40304
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_101_t_1.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_101_t_1.000000.vtu
new file mode 100644
index 00000000000..7740f443fd0
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_101_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:6e07523b28723d3af3a76c9c9dac72015b60c26e977d2a30291fab6f00b4ba32
+size 180058
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_1_t_0.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_1_t_0.000000.vtu
new file mode 100644
index 00000000000..a6abfcd2f76
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_1_t_0.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:b15553f2172b5a087ef66442eac60f3352ab90d57db7c71d56d947cf1f7622da
+size 173707
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_51_t_0.500000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_51_t_0.500000.vtu
new file mode 100644
index 00000000000..c29822af19d
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dload/m2_2Dload_ym45_pcs_0_ts_51_t_0.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:20e68cb898ab435a73aa600af78b0e26aa0297ca0cf507504a8206bfe247e842
+size 178736
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.gml
new file mode 100644
index 00000000000..2f6eb25c699
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:eed89ce16025ed882a3dab9665b9124f64947920ee707946fa53bbfa6b41a0bb
+size 975
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.prj
new file mode 100644
index 00000000000..97c45678331
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.prj
@@ -0,0 +1,307 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_2Dloadbt.vtu</mesh>
+ <geometry>m2_2Dloadbt.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair><repeat>505</repeat><delta_t>0.002</delta_t></pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_2Dloadbt</prefix>
+ <timesteps>
+ <pair><repeat>1</repeat><each_steps>300</each_steps></pair>
+ <pair><repeat>1</repeat><each_steps>200</each_steps></pair>
+ </timesteps>
+ <fixed_output_times>0.6 1.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.25</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension_Scaling</name>
+ <type>Constant</type>
+ <value>50</value>
+ </parameter>
+ <parameter>
+ <name>Minus_Neumann_Tension_Scaling</name>
+ <type>Constant</type>
+ <value>-50</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension</name>
+ <type>CurveScaled</type>
+ <curve>NeumannCurve</curve>
+ <parameter>Neumann_Tension_Scaling</parameter>
+ </parameter>
+ <parameter>
+ <name>Minus_Neumann_Tension</name>
+ <type>CurveScaled</type>
+ <curve>NeumannCurve</curve>
+ <parameter>Minus_Neumann_Tension_Scaling</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>NeumannCurve</name>
+ <coords>0.0 2.0</coords>
+ <values>0.0 2.0</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_2Dloadbt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dloadbt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dloadbt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dloadbt</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dloadbt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dloadbt</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_2Dloadbt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_2Dloadbt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_2Dloadbt_pcs_0_ts_300_t_0.600000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dloadbt_pcs_0_ts_300_t_0.600000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dloadbt_pcs_0_ts_300_t_0.600000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dloadbt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dloadbt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_2Dloadbt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.vtu
new file mode 100644
index 00000000000..0be925e4015
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c88f661790c07a35d3a2c6addc9790501055326596bffa980b9764c49526d33d
+size 21976
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt_pcs_0_ts_300_t_0.600000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt_pcs_0_ts_300_t_0.600000.vtu
new file mode 100644
index 00000000000..e5a6e2658d9
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt_pcs_0_ts_300_t_0.600000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:807103fc682921dca955653bd034d4b34ff19890a50ec88ff19079b3a809d132
+size 93605
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt_pcs_0_ts_500_t_1.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt_pcs_0_ts_500_t_1.000000.vtu
new file mode 100644
index 00000000000..2cef084597a
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_2Dloadbt/m2_2Dloadbt_pcs_0_ts_500_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:106dc15c5e573682c86e04066f40a9eaaed17595a94ec6d787338ab0ac8e54b2
+size 93249
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.gml
new file mode 100644
index 00000000000..cf4b6d9962e
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:ceaba747fee60794999df330f3ea671932eb3c5788822da3775a25b41f4bf645
+size 5458
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.prj
new file mode 100644
index 00000000000..e2020c988d7
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.prj
@@ -0,0 +1,312 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_3Dload.vtu</mesh>
+ <geometry>m2_3Dload.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair><repeat>1</repeat><delta_t>1e-10</delta_t></pair>
+ <pair><repeat>105</repeat><delta_t>0.01</delta_t></pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_3Dload</prefix>
+ <timesteps>
+ <pair><repeat>1</repeat><each_steps>51</each_steps></pair>
+ <pair><repeat>1</repeat><each_steps>50</each_steps></pair>
+ </timesteps>
+ <fixed_output_times>0.5 1.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.25</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension</name>
+ <type>Constant</type>
+ <value>40</value>
+ </parameter>
+ <parameter>
+ <name>Minus_Neumann_Tension</name>
+ <type>Constant</type>
+ <value>-40</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <!-- Fixation -->
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <!-- Loads -->
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE7</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE11</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE12</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE13</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE16</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dload</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_3Dload_pcs_0_ts_51_t_0.500000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dload_pcs_0_ts_51_t_0.500000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dload_pcs_0_ts_51_t_0.500000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>2.0e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dload_pcs_0_ts_101_t_1.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dload_pcs_0_ts_101_t_1.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>3.0e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dload_pcs_0_ts_101_t_1.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.vtu
new file mode 100644
index 00000000000..0be925e4015
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c88f661790c07a35d3a2c6addc9790501055326596bffa980b9764c49526d33d
+size 21976
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload_pcs_0_ts_101_t_1.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload_pcs_0_ts_101_t_1.000000.vtu
new file mode 100644
index 00000000000..1a12e6ee259
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload_pcs_0_ts_101_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:5b8cd40013587d698017440bb817bb7a2a2a013e3f0ec87e4a29a3c1d3930691
+size 80201
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload_pcs_0_ts_51_t_0.500000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload_pcs_0_ts_51_t_0.500000.vtu
new file mode 100644
index 00000000000..a3755d183b7
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dload/m2_3Dload_pcs_0_ts_51_t_0.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:ca69a893858e4f9bbffe6b46bd3c34640d44fd1673ab3f38bab2ffc394307bc0
+size 79541
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.gml b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.gml
new file mode 100644
index 00000000000..205c91ae25f
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3cebdff3e254db4b36709c4ee993c7c05641e792cff6898f5e216b037293d24e
+size 5460
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.prj b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.prj
new file mode 100644
index 00000000000..aa2a10837ee
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.prj
@@ -0,0 +1,330 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>m2_3Dloadbt.vtu</mesh>
+ <geometry>m2_3Dloadbt.gml</geometry>
+ <processes>
+ <process>
+ <name>ThermoMechanics</name>
+ <type>THERMO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <constitutive_relation>
+ <type>CreepBGRa</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ <a>A</a>
+ <n>n</n>
+ <sigma0>sigma_f</sigma0>
+ <q>Q</q>
+ <nonlinear_solver>
+ <maximum_iterations>1000</maximum_iterations>
+ <error_tolerance>1e-8</error_tolerance>
+ </nonlinear_solver>
+ </constitutive_relation>
+ <reference_solid_density>rho_sr</reference_solid_density>
+ <linear_thermal_expansion_coefficient>alpha</linear_thermal_expansion_coefficient>
+ <specific_heat_capacity>cs</specific_heat_capacity>
+ <thermal_conductivity>lambda</thermal_conductivity>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <temperature>temperature</temperature>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="ThermoMechanics">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1e-10</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair><repeat>505</repeat><delta_t>0.002</delta_t></pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>m2_3Dloadbt</prefix>
+ <timesteps>
+ <pair><repeat>1</repeat><each_steps>250</each_steps></pair>
+ <pair><repeat>1</repeat><each_steps>500</each_steps></pair>
+ </timesteps>
+ <fixed_output_times>0.5 1.0</fixed_output_times>
+ <variables>
+ <variable>displacement</variable>
+ <variable>temperature</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>A</name>
+ <type>Constant</type>
+ <value>0.18</value>
+ </parameter>
+ <parameter>
+ <name>n</name>
+ <type>Constant</type>
+ <value>5.0</value>
+ </parameter>
+ <parameter>
+ <name>sigma_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>Q</name>
+ <type>Constant</type>
+ <value>54000</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>25000</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.25</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>cs</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>lambda</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>temperature_bc</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>Dirichlet0</name>
+ <type>Constant</type>
+ <value>0.</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension_Scaling</name>
+ <type>Constant</type>
+ <value>40</value>
+ </parameter>
+ <parameter>
+ <name>Minus_Neumann_Tension_Scaling</name>
+ <type>Constant</type>
+ <value>-40</value>
+ </parameter>
+ <parameter>
+ <name>Neumann_Tension</name>
+ <type>CurveScaled</type>
+ <curve>NeumannCurve</curve>
+ <parameter>Neumann_Tension_Scaling</parameter>
+ </parameter>
+ <parameter>
+ <name>Minus_Neumann_Tension</name>
+ <type>CurveScaled</type>
+ <curve>NeumannCurve</curve>
+ <parameter>Minus_Neumann_Tension_Scaling</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>NeumannCurve</name>
+ <coords>0.0 2.0</coords>
+ <values>0.0 2.0</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <!-- Fixation -->
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE2</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE3</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>Dirichlet0</parameter>
+ </boundary_condition>
+ <!-- Loads -->
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE7</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE11</geometry>
+ <type>Neumann</type>
+ <component>1</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE12</geometry>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE13</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Minus_Neumann_Tension</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE16</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>Neumann_Tension</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE1</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>m2_3Dloadbt</geometrical_set>
+ <geometry>SURFACE4</geometry>
+ <type>Dirichlet</type>
+ <parameter>temperature_bc</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>
+ <lis>-i bicgstab -p jacobi -tol 1e-11 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-10</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <file>m2_3Dloadbt_pcs_0_ts_250_t_0.500000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dloadbt_pcs_0_ts_250_t_0.500000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dloadbt_pcs_0_ts_250_t_0.500000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dloadbt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>displacement</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dloadbt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-12</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>m2_3Dloadbt_pcs_0_ts_500_t_1.000000.vtu</file>
+ <field>sigma</field>
+ <absolute_tolerance>1.0e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.vtu
new file mode 100644
index 00000000000..0be925e4015
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c88f661790c07a35d3a2c6addc9790501055326596bffa980b9764c49526d33d
+size 21976
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt_pcs_0_ts_250_t_0.500000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt_pcs_0_ts_250_t_0.500000.vtu
new file mode 100644
index 00000000000..749a9e2f0c6
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt_pcs_0_ts_250_t_0.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0fdd909b17405126c25b44f4ad02d0d825e1611a4996ce5c7bc4c274b6acfe4b
+size 78857
diff --git a/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt_pcs_0_ts_500_t_1.000000.vtu b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt_pcs_0_ts_500_t_1.000000.vtu
new file mode 100644
index 00000000000..d93aa7ce5d3
--- /dev/null
+++ b/Tests/Data/ThermoMechanics/CreepBGRa/Verification/m2_3Dloadbt/m2_3Dloadbt_pcs_0_ts_500_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:19d9b1de2b1d33d1fb891e030eca86c08afef0f0e02df2c87e5ca53e2a78b0e1
+size 78889
--
GitLab