diff --git a/ProcessLib/ThermoHydroMechanics/Tests.cmake b/ProcessLib/ThermoHydroMechanics/Tests.cmake
index a92655671d5a6a36c3cef9a10df4edda53cf9842..da8dedc413dca038c3d63a9fb6a6794c475007fd 100644
--- a/ProcessLib/ThermoHydroMechanics/Tests.cmake
+++ b/ProcessLib/ThermoHydroMechanics/Tests.cmake
@@ -3,6 +3,7 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE ThermoHydroMechanics/Linear/verification/thm2_1Dfixd/thm2_1Dfixd.prj RUNTIME 60)
OgsTest(PROJECTFILE ThermoHydroMechanics/A2/A2.prj RUNTIME 23)
OgsTest(PROJECTFILE ThermoHydroMechanics/A2/A2_heating.prj RUNTIME 23)
+ OgsTest(PROJECTFILE ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj RUNTIME 15)
endif()
AddTest(
diff --git a/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj
new file mode 100644
index 0000000000000000000000000000000000000000..dfeccc53982bba5a447c89a80335b3d0a705c75b
--- /dev/null
+++ b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj
@@ -0,0 +1,325 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<!--units: kg, m, s, K and derived units: N, Pa=N/m^2, Da=m^2 -->
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>mesh.vtu</mesh>
+ <mesh>mesh_left.vtu</mesh>
+ <mesh>mesh_right.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>HeatConduction</name>
+ <type>THERMO_HYDRO_MECHANICS</type>
+ <integration_order>2</integration_order>
+ <dimension>2</dimension>
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <temperature>temperature</temperature>
+ <pressure>pressure</pressure>
+ <displacement>displacement</displacement>
+ </process_variables>
+ <specific_body_force>0 0</specific_body_force>
+ <secondary_variables>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="velocity" output_name="velocity"/>
+ </secondary_variables>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <!-- Pore water phase -->
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>0.54</value> <!-- W/m/K-->
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>4186</value> <!-- J/kg/K-->
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1000</value> <!--kg/m^3-->
+ </property>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1.0e-3</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <!-- Pore ice phase -->
+ <type>FrozenLiquid</type>
+ <properties>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>2.37</value> <!-- W/m/K-->
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>2052</value> <!-- J/kg/K-->
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>900</value> <!--kg/m^3-->
+ </property>
+ <property>
+ <name>specific_latent_heat</name>
+ <type>Constant</type>
+ <value>334000</value> <!-- J/kg-->
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <!-- Porous solid phase -->
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>0</value> <!-- W/m/K-->
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>0</value> <!-- J/kg/K-->
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>0</value> <!--kg/m^3-->
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>1.0</value> <!--no porous solid present-->
+ </property>
+ <property>
+ <name>volume_fraction</name>
+ <type>TemperatureDependentFraction</type>
+ <steepness>5</steepness> <!-- 1/K -->
+ <characteristic_temperature>273.15</characteristic_temperature> <!-- K -->
+ </property>
+ <property>
+ <name>density</name>
+ <type>VolumeFractionAverage</type>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>VolumeFractionAverage</type>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1.e-11</value>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <parameters>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e10</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.3</value>
+ </parameter>
+ <parameter>
+ <name>T0</name>
+ <type>Constant</type>
+ <value>272.25</value>
+ </parameter>
+ <parameter>
+ <name>TB</name>
+ <type>Constant</type>
+ <value>373.15</value>
+ </parameter>
+ <parameter>
+ <name>U0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>zero</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>mesh_right</mesh>
+ <type>Dirichlet</type>
+ <parameter>T0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>mesh_left</mesh>
+ <type>Dirichlet</type>
+ <parameter>TB</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>zero</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>mesh</mesh>
+ <type>Dirichlet</type>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>1</order>
+ <initial_condition>U0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>mesh</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>mesh</mesh>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <time_loop>
+ <processes>
+ <process ref="HeatConduction">
+ <nonlinear_solver>NLS</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentResidual</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-10 1e-16 1e-16 1e-16</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>427000000</t_end> <!--427e6 s x 30780000 s/a ~ 13.5 a-->
+ <timesteps>
+ <pair>
+ <repeat>427</repeat>
+ <delta_t>1e6</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>Stefan_problem</prefix>
+ <timesteps>
+ <pair>
+ <repeat>10000</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable> temperature </variable>
+ <variable> heat_flux </variable>
+ </variables>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ </output>
+ </time_loop>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>NLS</name>
+ <type>Newton</type>
+ <!-- Many iterations are needed for the first time step. -->
+ <max_iter>500</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ <scaling>true</scaling>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <!--primary field-->
+ <vtkdiff>
+ <file>Stefan_problem_ts_46_t_46000000.000000.vtu</file>
+ <field>temperature</field>
+ <absolute_tolerance>2e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>Stefan_problem_ts_427_t_427000000.000000.vtu</file>
+ <field>temperature</field>
+ <absolute_tolerance>7e-11</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <!--secondary field-->
+ <vtkdiff>
+ <file>Stefan_problem_ts_46_t_46000000.000000.vtu</file>
+ <field>HeatFlowRate</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>Stefan_problem_ts_427_t_427000000.000000.vtu</file>
+ <field>HeatFlowRate</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem_ts_427_t_427000000.000000.vtu b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem_ts_427_t_427000000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..447b73addfab9928629991707045ed50ee01db12
--- /dev/null
+++ b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem_ts_427_t_427000000.000000.vtu
@@ -0,0 +1,39 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="166" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45" RangeMax="103" offset="176" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="400" format="appended" RangeMin="0" RangeMax="0" offset="256" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="400" format="appended" RangeMin="0" RangeMax="0" offset="348" />
+ </FieldData>
+ <Piece NumberOfPoints="202" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-2.4389855422" RangeMax="3.0246246598" offset="440" />
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-0" RangeMax="-0" offset="1272" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1352" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="201" offset="1440" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="1956" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="2040" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="2112" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="272.25" RangeMax="373.15" offset="2196" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="272.25" RangeMax="373.15" offset="3676" />
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="5156" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="5236" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="10.049875621" offset="5516" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="6400" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="7196" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="7504" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAKYAAAAAAAAAYwAAAAAAAAA=eF6FzDEKgDAMheG7ZO4kTr2KSIgaJWCTktZBxLvb1UXH9374LhCtvDlVMcVsbSG501kgDtcrmi/sELsASokhQpEtEUqG9hxpYkdbcbaUTVlrA/o7/BCci+wtfCLj/QAlSj0gAQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9Y1NbTUTU80Mbc0TzQ0AwApWwReAQAAAAAAAAAAgAAAAAAAAAAyAAAAAAAAIwAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAAAAAAAAAAAAAAB4NjIAAAE=AQAAAAAAAAAAgAAAAAAAAAAyAAAAAAAAIwAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAAAAAAAAAAAAAAB4NjIAAAE=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAATwIAAAAAAAA=eF5dlDtoFFEUhq8ovovxTQziRMV3SNy8ddfczWOzJkHGTQyIDwasbGQrbS9oYbkISrAasBa2SuXClILNlsFqEYuU6W2E+b8jrtP8/GfuPec/r/vq+88XL0sHvSu+4UoB8UqnwNG6eH5V3Hvx1rqwWZM9HRPmt4VhUP/9lHi2KB5Ocm5I3E0RryG7OyQeLYl3QbcHPMD/BfyOEacsvvMI3Q/wz/32vHiyKN6bRUcdP9P4mSHO9Uo/OnTeEWbL8Dn+rwl357k/if+b5E3d3Dg4rf/hBP/Piqf49Y/RVcWOPndB+G6N+7F4Th3dFfgq/BL8Hn6HDdFZ69fdm+Ge6aZfjnrFlscc57Fn1ify234mTDfQT3/a1pdfXwuImS9fJt4pzpH3t4ey96i7G8HvXXiJ+5PkeZ94nA9W78vope7udGFPK1sF+jNbsn95Lxz8QB5PK6+TgR+H3+zNf/+7Jym6//aZeNvrnT675eF2yroHt3kMA+gaFzaZI6trSr29+R0RWp3ckHCTOjerOt9aQB97Gc4JY5u7Y8JsRefa7EVM/wN722Ye3YTpgLM3Cfcjy4c5drn6m9gcHgHR6472x7U5ztDbo6/O8rwo3qR+UYO4vD8R71UKzw2X0enh1X6+C28l+GEPwnmhv6E4deJlq9xj/l1XfQ2z5M/7F9m7sA9dFfhx/D8RtumTZ79DJLT31u1HB3ULzJW7BWfvA+/CKO+iuybeK3HO6jlBHfgfLWHHT9f/p3ejQP9ZexIS7Ul4+0nnapvofd752NCe/AEPtuvjAQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAAGgAAAAAAAAA=eF5jYACDBoZRepQepUfpUXqUpoAGALHjZQE=AQAAAAAAAAAAgAAAAAAAAKAMAAAAAAAAIAAAAAAAAAA=eF7txbENAAAIwKCe7ukmPuECC3Um27Zt235+ARk3ygE=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAAYAEAAAAAAAA=eF4txddCCAAAAMAiIytlr1RSRkb2LE2UIquQnRFJ2TtEZtnKlhFZyd579Vc9uHu5gID/At3CLR3kVm7tNm7rYLdze3dwR3dyiDs71GHu4q7u5u7u4Z7u5d7u477u53D3d4QjHeUBjvZAxzjWgzzYQzzUcR7m4R7hkY73KI/2GI/1OI/3BE/0JE/2FE91ghM9zUlOdopTneZ0T/cMz3SGMz3LWc72bM9xjud6nud7gRc613le5MVe4nwv9TIv9wqv9CqvdoHXeK3Xeb0LvcEbXeRNLvZml7jUW7zV27zdO7zTu7zbe7zX+7zfB1zmgz7kwy73ER91hY/5uE/4pE/5tCtd5TM+63M+7wu+6Eu+7GrX+Iqv+pqv+4Zv+pZvu9Z3fNf3fN91fuCHrvcjP/YTP/UzN/i5G/3CL/3Kr/3Gb/3O7/3BH/3Jn/3FX/3N3/3DP/3Lv/3Hf/3PTW4GwEFPTg==AQAAAAAAAAAAgAAAAAAAAEAZAAAAAAAAHQAAAAAAAAA=eF7twTEBAAAAwqD1T+1rCKAAAAAAAAB4AxlAAAE=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAAFQAAAAAAAAA=eF5jYBgFo2AUjIJRMAooBwAGUAABAQAAAAAAAAAAgAAAAAAAAEAZAAAAAAAAHQAAAAAAAAA=eF7twTEBAAAAwqD1T+1rCKAAAAAAAAB4AxlAAAE=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAANQQAAAAAAAA=eF6F03s013ccx/FyGX42I6M0TnF+sxnHUeP8JpdeataOFcbpKCTs9Eus8vnekHB+btXZpG1axmmhJdH82DguczaFYSycuR1M5Laj2pzTdlxymY73X79/9vzn9cfr74dSuV5AHJxa25cSEAdVVetAlmMcgnb85V1kGQfF/cqsGlkcjl6ImP9tQcLdOYuOqWkJxqm1/mu9Er507Ora3iTBvf1a9ruVEp4lzPQfuSFhv8LmSWKWhFRvR7OiRAk7Gk5e6zglocUkp2YhSMKmir46u/clWGffGA9zluA4XMJybSSUFY2V9RtLWMuO6LFYE3G6/IRN5FMRoW57RtXDIiwy3tmv/auInoGCy6G1IlYSGpd/LBZhV/dDm/VVEX76vMPlNBH/Ot0P38yJsM2VX08MF1GV27N1xVfEy/+U6GZ6rG/Im3nmDiIOGofoqbeLUN4aTPIzEDHUq3ZYmBcQcC7c6860gOnSrpnIPgHBf7v7vNEswKVI8fnc9wLu8DkTzYUC4s0+VhZdEbBqPOl+KUVASXvyxYQzAuzLa6L4YwJ0822fSwcFNMWeDkl3E3BvzPFm/tsCTHnd2QYLAZfaznn/qS/g1utrbVYLPL7QHs84NsMj8Lz3p7f7eRxQ7Jx+3sLjuMKhMLiaR3zDZHPTtzz2qZKDFDk8Yv3GQ6vTeBQ/Lez34HmsRi83dUfyWOx4YH0mgEeZZc6s2T4esifHzdt28YhwHspLt17/q9xiPjThMdjoE2+5mcdI5pGK+TkO8gHPtZExDoOnwnw7uzmcL7ue3tLIoWVLVGZbBYcl7XKXvgIO3d5T7PEVDuWGsm2GKg62XxdoOTMO4ak1BsoIDrKsWa3CjzgEuCjvTXpxCBx6bLZrN4fVnP4HGTYc6ht21j7awiHMMjr/gPaLf3Vv1TMGB/OH0W9NMix5tI7e7GVwtXINtW1hMLmdUaiuZghPGojeW8zwjYsn1/cVw8XQ+rPsAoORq4+FaTzDhO+8TX0Ug1g66ak8ypCVEWm11YchpUJ+tnMPQ9UHOhOZ9gxxpe8ZfmfJ8LuxWuuQKcPoyV5VqhHDQ6+p8lYZw0vTh5uD9Bg2vUiHQUlOHMlJMjkJJCdO5CSAnBSQkxXVhpMkciInJ6PkxImcMHJiQk4ayMmiesPJNnIiJyeF5GSenJwgJ0Hk5DVy0kVOlsiJXMPJIjlxICe15MSInBiTE38NJ4PkxI+cPNJw4kZOKslJKjkxMNlwUkFOdpMTGTn5hZz8TE6MyImKnOSRk880nBwiJ1HkREVO/MlJKjmpIyevxmw4MejccFJPTszISQw5WSUnw+RkgpzYk5M/yEkKOWklJ8saTu7+j5PD5CSInOhc3XDyk4YTf3JiR04WyImCnLxCToLJSS45SSMneuRkhJx8Qk7SyIlITsrISSw56dZwMkZO9DWc/AeAJ+ZZAQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAANQQAAAAAAAA=eF6F03s013ccx/FyGX42I6M0TnF+sxnHUeP8JpdeataOFcbpKCTs9Eus8vnekHB+btXZpG1axmmhJdH82DguczaFYSycuR1M5Laj2pzTdlxymY73X79/9vzn9cfr74dSuV5AHJxa25cSEAdVVetAlmMcgnb85V1kGQfF/cqsGlkcjl6ImP9tQcLdOYuOqWkJxqm1/mu9Er507Ora3iTBvf1a9ruVEp4lzPQfuSFhv8LmSWKWhFRvR7OiRAk7Gk5e6zglocUkp2YhSMKmir46u/clWGffGA9zluA4XMJybSSUFY2V9RtLWMuO6LFYE3G6/IRN5FMRoW57RtXDIiwy3tmv/auInoGCy6G1IlYSGpd/LBZhV/dDm/VVEX76vMPlNBH/Ot0P38yJsM2VX08MF1GV27N1xVfEy/+U6GZ6rG/Im3nmDiIOGofoqbeLUN4aTPIzEDHUq3ZYmBcQcC7c6860gOnSrpnIPgHBf7v7vNEswKVI8fnc9wLu8DkTzYUC4s0+VhZdEbBqPOl+KUVASXvyxYQzAuzLa6L4YwJ0822fSwcFNMWeDkl3E3BvzPFm/tsCTHnd2QYLAZfaznn/qS/g1utrbVYLPL7QHs84NsMj8Lz3p7f7eRxQ7Jx+3sLjuMKhMLiaR3zDZHPTtzz2qZKDFDk8Yv3GQ6vTeBQ/Lez34HmsRi83dUfyWOx4YH0mgEeZZc6s2T4esifHzdt28YhwHspLt17/q9xiPjThMdjoE2+5mcdI5pGK+TkO8gHPtZExDoOnwnw7uzmcL7ue3tLIoWVLVGZbBYcl7XKXvgIO3d5T7PEVDuWGsm2GKg62XxdoOTMO4ak1BsoIDrKsWa3CjzgEuCjvTXpxCBx6bLZrN4fVnP4HGTYc6ht21j7awiHMMjr/gPaLf3Vv1TMGB/OH0W9NMix5tI7e7GVwtXINtW1hMLmdUaiuZghPGojeW8zwjYsn1/cVw8XQ+rPsAoORq4+FaTzDhO+8TX0Ug1g66ak8ypCVEWm11YchpUJ+tnMPQ9UHOhOZ9gxxpe8ZfmfJ8LuxWuuQKcPoyV5VqhHDQ6+p8lYZw0vTh5uD9Bg2vUiHQUlOHMlJMjkJJCdO5CSAnBSQkxXVhpMkciInJ6PkxImcMHJiQk4ayMmiesPJNnIiJyeF5GSenJwgJ0Hk5DVy0kVOlsiJXMPJIjlxICe15MSInBiTE38NJ4PkxI+cPNJw4kZOKslJKjkxMNlwUkFOdpMTGTn5hZz8TE6MyImKnOSRk880nBwiJ1HkREVO/MlJKjmpIyevxmw4MejccFJPTszISQw5WSUnw+RkgpzYk5M/yEkKOWklJ8saTu7+j5PD5CSInOhc3XDyk4YTf3JiR04WyImCnLxCToLJSS45SSMneuRkhJx8Qk7SyIlITsrISSw56dZwMkZO9DWc/AeAJ+ZZAQAAAAAAAAAAgAAAAAAAAKAMAAAAAAAAGgAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAAA+DAygAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAArwAAAAAAAAA=eF4txdciAgAAAMCQWSFbJCuzsjNDiuyMyo7y///gwd3LBQL/WtzqNgfd7g53usvd7nHIYUfc6z73O+oBD3rIwx7xqMc87pgnPOm4p5zwtGc86znPO+kFL3rJy17xqlNOO+M1r3vDm97ytnec9a73vO8DH/rIOR/7xKfO+8wFF33uC5d86Stf+8a3vvO9y37wo5/87IqrrvnFr37zuz/86S/X/e0fN9z0r/8AXBwTVw==AQAAAAAAAAAAgAAAAAAAAPASAAAAAAAAdQIAAAAAAAA=eF511stK1WEUhvGd2cnOWXYwraw7sLlejV5T3UmXIOiokTNBCAQHShKKIrF7WmzW+/it4Y/4nr+I8U4mN9/3b9P7sXGz78i/rk/v5+Df78u5g8E7h/Ltren9Grx/JN/bnd7xoHsyeOd08D2/B++fDb7zz6B7PuheDLqXg+7VoHs96E42E3jHzvfckvP+nJzvvC3n5uV07XTvyOneldO9J+fuy+na6T6Q012Q030o5x7J6drpPpbTfSKn+1TOPZPzjp3vsf97fv25nO+0cy/kdO107fxci3K6du6lnK6drp3uKzldO7ckp2una6f7Wk7Xzr2R07XTtdN9K6dr597J6drp2ukuy+naufdyuna6drorcrp2blVO107XTveDnK6d+yjn79HO99h538532rlPcrp2una6drp2bk1O107XTtdO1859lvP7stO107XTtXNfBl432yN9d6XvyPvuyn+/L+cOBu8cyvvuyveP5H13Zfdk8M6pvO+ufP9M3ndXds/lfXdl91Led1d2r+X/bzO97678ntpd+f6cvO+u7Na+yq69767s1u7Kbu2u7Na+yq69767s1u7Kbu2u7Na+yq69767s1u7Kbu2u7NaOynfsfXfNnP/3al/ld9q52lHZtffdlT9X7avs2rnaUdm1992V3dpX2bVztaOya++7K7u1r7Jr52pHZdfed1d2a19l187Vjsquve+u7Na+yq6dqx2VXXvfXdldkffdld3aUdm1992V3dpX2bVztYtmzt+jve+ufN/ed1d2axdl1953V3btfXdld01O1953V3btfXdlt3ZR/r7sfXdl1953V3Zrd838L+h5uNA=AQAAAAAAAAAAgAAAAAAAAIAMAAAAAAAAMgIAAAAAAAA=eF5dyuUXCAYARXHdHRPTzXQzXZvujqlpwzBs2nQM00w302Ziasb8Vz6478t9X37nvHMLFPi8gjgPv9dfCOery18YF6jLXwQXqstfFBepy18MF6vLXxyXqMtfApeqy18Sf1CXvxQuU5e/NC5Xl78MrlCXvyz+qC5/OVypLn95XKUufwVcrS5/RfxJXf5KuEZd/sq4Vl3+KrhOXf6q+LO6/F/gL+ryV8P16vJXxw3q8tfAjery18RN6vJ/iZvV5a+FW9Tlr41b1eWvg9vU5a+Lv6rLXw+3q8tfH3eoy98Ad6rL3xB3qcvfCHery98Y96jL3wT3qsvfFPepy98M96vL3xwPqMvfAn9Tl/8rPKguf0s8pC5/KzysLn9r/F1d/jZ4RF3+tnhUXf52eExd/vZ4XF3+DnhCXf6OeFJd/k54Sl3+znhaXf4ueEZd/q74h7r83fCsuvzd8Zy6/F/jeXX5e+AFdfl74kV1+XvhJXX5e+Nldfn74BV1+fviVXX5++E1dfn743V1+QfgDXX5B+JNdfkH4S11+b/B2+ryf4t/qss/GO+oyz8E76rLPxTvqcs/DO+ryz8cH6jLPwIfqss/Eh+pyz8KH6vLPxr/Upd/DD5Rl38s/q0u/zh8qi7/eHymLv8EfK4u/0R8oS7/JPxHXf7J+FJd/in4Sl3+qfhaXf5p+EZd/un4Vl3+Gfivuvzf4Tt1+Wfif+ryz8L36vLPxg/q8s/B/9Xln4sf1X0C7M6dCQ==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxgAAAAAAAAA=eF4txRF0QgEAAMB6LwiCIAiCIAiCIPgQBEEQBEEQBEEQBEEQBEEQBINBEARBEARBMBgEQRAMBoNBEARBEATdyUVCb1HHHHfCSaecdsZZ55x3wYGLLrnsiquuue6Gm2657Y677rnvgYceeeyJp575w5+ee+GlV15746133vvL3z746JPP/vGv//zvi6+++e6Hnw6F30UcdcxxJ5x0ymlnnHXOeRccuOiSy6646prrbrjpltvuuOue+x546JHHnnjqmV/7BioOAQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF7j5KQ9AACx7gOF
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem_ts_46_t_46000000.000000.vtu b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem_ts_46_t_46000000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..5756a72b1de5113b1099adbc09c1c7b8b8b697ce
--- /dev/null
+++ b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem_ts_46_t_46000000.000000.vtu
@@ -0,0 +1,39 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="166" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45" RangeMax="103" offset="176" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="400" format="appended" RangeMin="0" RangeMax="0" offset="256" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="400" format="appended" RangeMin="0" RangeMax="0" offset="348" />
+ </FieldData>
+ <Piece NumberOfPoints="202" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.010612610776" RangeMax="8.2790548572" offset="440" />
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-0" RangeMax="-0" offset="1468" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1548" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="201" offset="1636" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="2152" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="2236" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="2308" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="272.25" RangeMax="373.15" offset="2392" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="272.25" RangeMax="373.15" offset="3792" />
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="5192" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="5272" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="10.049875621" offset="5552" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="6436" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="7232" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="7540" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAKYAAAAAAAAAYwAAAAAAAAA=eF6FzDEKgDAMheG7ZO4kTr2KSIgaJWCTktZBxLvb1UXH9374LhCtvDlVMcVsbSG501kgDtcrmi/sELsASokhQpEtEUqG9hxpYkdbcbaUTVlrA/o7/BCci+wtfCLj/QAlSj0gAQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9Y1NbTUTU80Mbc0TzQ0AwApWwReAQAAAAAAAAAAgAAAAAAAAAAyAAAAAAAAIwAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAAAAAAAAAAAAAAB4NjIAAAE=AQAAAAAAAAAAgAAAAAAAAAAyAAAAAAAAIwAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAAAAAAAAAAAAAAB4NjIAAAE=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAA4gIAAAAAAAA=eF5VlF1Ik2EUx1/BIo1kWJCUH6sEzcQ0NQ1bezNTCz82W2WmttLUNG1U2pCSJzW8iVgJgWD6GtWNIF54EaXsJYT0ykEQXdVSCOlqlZCQZHD+/wXu5sf/fD3nOe9zFptYUBJ8ZNU1+R2cAffbBHoRaB4HtVIylXHx5A5QFYNGCeKspbSnQ6sUxuex7mEw6AItVeE4MHSW8THheOok9mGjTuY52ayn014OmuxH28c43s/nhNbtPJd1tLVp2HOZx35UgfhV/5jQXH8B+7du5v1BXtRl2Ic5t0AB/VPw/9Vhn6gBkxxCPaVRaCR0SrzndS/yFi6CW3B/86Yb8V1NyJ/HPANbPTb4WxFvz4KeOyPaWKD99gEwmvePOoF7xdVCR1qh0/OgvRFSJ/imGXo77uf7fFW0vgnz1oMu+A3UN8pqoJeaJV4VYc6eXtTxnfPgnHnMfW6gTfyWynroX9BaYgTqddSJ9g6hnnfyijB++JrQMdglcSEf3o01Bu+ryIp+154hrnqwSfyu1g70sdIorKjG/LxP64Wf5t22nSPj9y3TD/zODXuSiXtpsSTflbaK7/t/L06BKp8sxn20I9ABvnstDdrkfoX3Qu2FXQ+/82gweJp+vmMjfE4y7faNfm0X6x+lvYz2cN1U2gtZr5w6g3Hsx3Ke/WUzj3upssg97CMD7OaejL6EX/fC7mafd08yzkn/usxPpbVCj+M7BbsU9A98F/XBgbzcBnA35qYP3YJ/EXvlmOiD/wb6NiKwN+6EBtTbhnemuQqFunkB5/X3CN0/ryOvrh31y7Gnpv8S8iph991jv4t2sasp/O+pPpyX+Zx78xHvNdDDvXUdQv21FvT9NUfilgeRt5zbhj5W8H8xO9KJPlY6hKE4zMdSVS2cfFuLut/Rx2o8+stvR/1Q5B3hF7cSv3WJe/y+UujoqRA+8XOPcjCP5Xfo+9gA7LNO7N/m393Ch49bZl5Njsme/AOSMC9uAQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAAGgAAAAAAAAA=eF5jYACDBoZRepQepUfpUXqUpoAGALHjZQE=AQAAAAAAAAAAgAAAAAAAAKAMAAAAAAAAIAAAAAAAAAA=eF7txbENAAAIwKCe7ukmPuECC3Um27Zt235+ARk3ygE=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAAYAEAAAAAAAA=eF4txddCCAAAAMAiIytlr1RSRkb2LE2UIquQnRFJ2TtEZtnKlhFZyd579Vc9uHu5gID/At3CLR3kVm7tNm7rYLdze3dwR3dyiDs71GHu4q7u5u7u4Z7u5d7u477u53D3d4QjHeUBjvZAxzjWgzzYQzzUcR7m4R7hkY73KI/2GI/1OI/3BE/0JE/2FE91ghM9zUlOdopTneZ0T/cMz3SGMz3LWc72bM9xjud6nud7gRc613le5MVe4nwv9TIv9wqv9CqvdoHXeK3Xeb0LvcEbXeRNLvZml7jUW7zV27zdO7zTu7zbe7zX+7zfB1zmgz7kwy73ER91hY/5uE/4pE/5tCtd5TM+63M+7wu+6Eu+7GrX+Iqv+pqv+4Zv+pZvu9Z3fNf3fN91fuCHrvcjP/YTP/UzN/i5G/3CL/3Kr/3Gb/3O7/3BH/3Jn/3FX/3N3/3DP/3Lv/3Hf/3PTW4GwEFPTg==AQAAAAAAAAAAgAAAAAAAAEAZAAAAAAAAHQAAAAAAAAA=eF7twTEBAAAAwqD1T+1rCKAAAAAAAAB4AxlAAAE=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAAFQAAAAAAAAA=eF5jYBgFo2AUjIJRMAooBwAGUAABAQAAAAAAAAAAgAAAAAAAAEAZAAAAAAAAHQAAAAAAAAA=eF7twTEBAAAAwqD1T+1rCKAAAAAAAAB4AxlAAAE=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAA+AMAAAAAAAA=eF5101tQFHQUx3HWRG4RhQTB6oBgTKODgCIYTHkAFVA0zUwwZchyKeLSf2+gIBiyChEraDHg3QYDRUFGQ0niUiFxCbFMgRXQZUIp2GGiQQGRDtPvaWc6L9+Hc+a8fWQynreTaHv+nKiIBUl0M+NTw+YRNaW4rju/s05Nt9KfOWi1arItPFvWt0NNdMauIdpDTZ1pE185Tako4J7bKvtfVHR8+cD17SdUtK/nWu1YvIqiWiXdQ2+qyLtaujPkRRU1tx2Nf0GvpCUNtzvCrihpmCodJQeVlDJ4Iydgm5L8w5rnmSxRkotzuNgyS0lNdL8r8K6CNHrtxtYLCso6dU39OENBHiVRR1ojFWQ62ZMf7q2g8VZ3SbqFguzup86N0cspVXvrpHmNnGxD2ubHFMjJ+dCxSo2Q09n1FxfGrJfTTYkwWC+SU/1va1akmsnpy5ir/Zf/EBS/dDK/6idBfXnnLmiKBRXXWSe6agR90f/UkCsTVJj8j2tzqCDXTovy3xcL8rUctbpiJ+jUps8TImwEdWVLvFZaC8p5kp1psBL0sWlv47iloPLQT3z8uc/Mf00+YSHIc6g91pHrM5XQUWIu6MirL+Wt4l4yKewYMhOUssMx4yS3I0TpGMHd5tQoe4Vbrhvb0ztH0HfBYQXnuH2u78/bzU18OJATzq2rSE92405S9IEpU0H6rClZN3dLUNadau6aO5NFx7iHG4r2pnNXe1Rt3sXNnv/DaDi3aSjzdV/upv17Zi/gln6wwc+aq9XGX52YLWi49HzcIFd3d4S6uH/mznJp4R6XBj2+zm0Nyfi2nHsmt2ft19xd+sOlBdws6cP6HG6Pn+XBz7h9I9OGJG7v3JzhBK5+hTRZxq1Iq8mMmtlLLCy3cr/pkE6/xR1x17wXxl19Y6VD8Mx9cMmiN7g2HolaP27cwNGApVwXT/ViD+4yH+3W17gyt7JqN+4Gz5Z1ztzOJw7mUq7PO389sue2Rdc/suVWDi+bZcPtXzjtacUN2j0aZ8Z1uj12+TmuyczM/IOTSDhph5NUOGmDEys4WQ4nTXDiDicaOImFk7Vw4gInNXDiBicP4CQBTrzg5GU4+R5O9sLJPjhxh5OnE/85+RtOnocTNZzYGDkphpMWOKmFkzw4+QhOuuHkNJxkw0kBnLjAiQ+cnIYTHZzkwkksnFQYOfGGE184KYSTKjhJg5N2OImEk0twUgsnD4yc1MLJuJGTd+EkFE7y4SQQTg7Ayc9wstHIyaH/cTIIJwVw8iOcFMFJNJzsh5MuONHBiQ5O7sHJRTjpNXJigJNAOOmBE0s4kcGJFE484eRDIyddcOIFJ41wUgYnOjjxhxN7Iyf/AoPtp1Q=AQAAAAAAAAAAgAAAAAAAAFAGAAAAAAAA+AMAAAAAAAA=eF5101tQFHQUx3HWRG4RhQTB6oBgTKODgCIYTHkAFVA0zUwwZchyKeLSf2+gIBiyChEraDHg3QYDRUFGQ0niUiFxCbFMgRXQZUIp2GGiQQGRDtPvaWc6L9+Hc+a8fWQynreTaHv+nKiIBUl0M+NTw+YRNaW4rju/s05Nt9KfOWi1arItPFvWt0NNdMauIdpDTZ1pE185Tako4J7bKvtfVHR8+cD17SdUtK/nWu1YvIqiWiXdQ2+qyLtaujPkRRU1tx2Nf0GvpCUNtzvCrihpmCodJQeVlDJ4Iydgm5L8w5rnmSxRkotzuNgyS0lNdL8r8K6CNHrtxtYLCso6dU39OENBHiVRR1ojFWQ62ZMf7q2g8VZ3SbqFguzup86N0cspVXvrpHmNnGxD2ubHFMjJ+dCxSo2Q09n1FxfGrJfTTYkwWC+SU/1va1akmsnpy5ir/Zf/EBS/dDK/6idBfXnnLmiKBRXXWSe6agR90f/UkCsTVJj8j2tzqCDXTovy3xcL8rUctbpiJ+jUps8TImwEdWVLvFZaC8p5kp1psBL0sWlv47iloPLQT3z8uc/Mf00+YSHIc6g91pHrM5XQUWIu6MirL+Wt4l4yKewYMhOUssMx4yS3I0TpGMHd5tQoe4Vbrhvb0ztH0HfBYQXnuH2u78/bzU18OJATzq2rSE92405S9IEpU0H6rClZN3dLUNadau6aO5NFx7iHG4r2pnNXe1Rt3sXNnv/DaDi3aSjzdV/upv17Zi/gln6wwc+aq9XGX52YLWi49HzcIFd3d4S6uH/mznJp4R6XBj2+zm0Nyfi2nHsmt2ft19xd+sOlBdws6cP6HG6Pn+XBz7h9I9OGJG7v3JzhBK5+hTRZxq1Iq8mMmtlLLCy3cr/pkE6/xR1x17wXxl19Y6VD8Mx9cMmiN7g2HolaP27cwNGApVwXT/ViD+4yH+3W17gyt7JqN+4Gz5Z1ztzOJw7mUq7PO389sue2Rdc/suVWDi+bZcPtXzjtacUN2j0aZ8Z1uj12+TmuyczM/IOTSDhph5NUOGmDEys4WQ4nTXDiDicaOImFk7Vw4gInNXDiBicP4CQBTrzg5GU4+R5O9sLJPjhxh5OnE/85+RtOnocTNZzYGDkphpMWOKmFkzw4+QhOuuHkNJxkw0kBnLjAiQ+cnIYTHZzkwkksnFQYOfGGE184KYSTKjhJg5N2OImEk0twUgsnD4yc1MLJuJGTd+EkFE7y4SQQTg7Ayc9wstHIyaH/cTIIJwVw8iOcFMFJNJzsh5MuONHBiQ5O7sHJRTjpNXJigJNAOOmBE0s4kcGJFE484eRDIyddcOIFJ41wUgYnOjjxhxN7Iyf/AoPtp1Q=AQAAAAAAAAAAgAAAAAAAAKAMAAAAAAAAGgAAAAAAAAA=eF7twTEBAAAAwqD1T20ND6AAAAA+DAygAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAArwAAAAAAAAA=eF4txdciAgAAAMCQWSFbJCuzsjNDiuyMyo7y///gwd3LBQL/WtzqNgfd7g53usvd7nHIYUfc6z73O+oBD3rIwx7xqMc87pgnPOm4p5zwtGc86znPO+kFL3rJy17xqlNOO+M1r3vDm97ytnec9a73vO8DH/rIOR/7xKfO+8wFF33uC5d86Stf+8a3vvO9y37wo5/87IqrrvnFr37zuz/86S/X/e0fN9z0r/8AXBwTVw==AQAAAAAAAAAAgAAAAAAAAPASAAAAAAAAdQIAAAAAAAA=eF511stK1WEUhvGd2cnOWXYwraw7sLlejV5T3UmXIOiokTNBCAQHShKKIrF7WmzW+/it4Y/4nr+I8U4mN9/3b9P7sXGz78i/rk/v5+Df78u5g8E7h/Ltren9Grx/JN/bnd7xoHsyeOd08D2/B++fDb7zz6B7PuheDLqXg+7VoHs96E42E3jHzvfckvP+nJzvvC3n5uV07XTvyOneldO9J+fuy+na6T6Q012Q030o5x7J6drpPpbTfSKn+1TOPZPzjp3vsf97fv25nO+0cy/kdO107fxci3K6du6lnK6drp3uKzldO7ckp2una6f7Wk7Xzr2R07XTtdN9K6dr597J6drp2ukuy+naufdyuna6drorcrp2blVO107XTveDnK6d+yjn79HO99h538532rlPcrp2una6drp2bk1O107XTtdO1859lvP7stO107XTtXNfBl432yN9d6XvyPvuyn+/L+cOBu8cyvvuyveP5H13Zfdk8M6pvO+ufP9M3ndXds/lfXdl91Led1d2r+X/bzO97678ntpd+f6cvO+u7Na+yq69767s1u7Kbu2u7Na+yq69767s1u7Kbu2u7Na+yq69767s1u7Kbu2u7NaOynfsfXfNnP/3al/ld9q52lHZtffdlT9X7avs2rnaUdm1992V3dpX2bVztaOya++7K7u1r7Jr52pHZdfed1d2a19l187Vjsquve+u7Na+yq6dqx2VXXvfXdldkffdld3aUdm1992V3dpX2bVztYtmzt+jve+ufN/ed1d2axdl1953V3btfXdld01O1953V3btfXdlt3ZR/r7sfXdl1953V3Zrd838L+h5uNA=AQAAAAAAAAAAgAAAAAAAAIAMAAAAAAAAMgIAAAAAAAA=eF5dyuUXCAYARXHdHRPTzXQzXZvujqlpwzBs2nQM00w302Ziasb8Vz6478t9X37nvHMLFPi8gjgPv9dfCOery18YF6jLXwQXqstfFBepy18MF6vLXxyXqMtfApeqy18Sf1CXvxQuU5e/NC5Xl78MrlCXvyz+qC5/OVypLn95XKUufwVcrS5/RfxJXf5KuEZd/sq4Vl3+KrhOXf6q+LO6/F/gL+ryV8P16vJXxw3q8tfAjery18RN6vJ/iZvV5a+FW9Tlr41b1eWvg9vU5a+Lv6rLXw+3q8tfH3eoy98Ad6rL3xB3qcvfCHery98Y96jL3wT3qsvfFPepy98M96vL3xwPqMvfAn9Tl/8rPKguf0s8pC5/KzysLn9r/F1d/jZ4RF3+tnhUXf52eExd/vZ4XF3+DnhCXf6OeFJd/k54Sl3+znhaXf4ueEZd/q74h7r83fCsuvzd8Zy6/F/jeXX5e+AFdfl74kV1+XvhJXX5e+Nldfn74BV1+fviVXX5++E1dfn743V1+QfgDXX5B+JNdfkH4S11+b/B2+ryf4t/qss/GO+oyz8E76rLPxTvqcs/DO+ryz8cH6jLPwIfqss/Eh+pyz8KH6vLPxr/Upd/DD5Rl38s/q0u/zh8qi7/eHymLv8EfK4u/0R8oS7/JPxHXf7J+FJd/in4Sl3+qfhaXf5p+EZd/un4Vl3+Gfivuvzf4Tt1+Wfif+ryz8L36vLPxg/q8s/B/9Xln4sf1X0C7M6dCQ==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxgAAAAAAAAA=eF4txRF0QgEAAMB6LwiCIAiCIAiCIPgQBEEQBEEQBEEQBEEQBEEQBINBEARBEARBMBgEQRAMBoNBEARBEATdyUVCb1HHHHfCSaecdsZZ55x3wYGLLrnsiquuue6Gm2657Y677rnvgYceeeyJp575w5+ee+GlV15746133vvL3z746JPP/vGv//zvi6+++e6Hnw6F30UcdcxxJ5x0ymlnnHXOeRccuOiSy6646prrbrjpltvuuOue+x546JHHnnjqmV/7BioOAQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF7j5KQ9AACx7gOF
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh.vtu b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..2781f6a0b77a5da3e56047df730d4cfb0ad72e9a
--- /dev/null
+++ b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="202" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="201" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="2168" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="10.049875621" offset="3248" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="9724" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="14004" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="15084" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _UAYAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAEwAAAAAAAAATQAAAAAAAABOAAAAAAAAAE8AAAAAAAAAUAAAAAAAAABRAAAAAAAAAFIAAAAAAAAAUwAAAAAAAABUAAAAAAAAAFUAAAAAAAAAVgAAAAAAAABXAAAAAAAAAFgAAAAAAAAAWQAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABdAAAAAAAAAF4AAAAAAAAAXwAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYgAAAAAAAABjAAAAAAAAAGQAAAAAAAAAZQAAAAAAAABmAAAAAAAAAGcAAAAAAAAAaAAAAAAAAABpAAAAAAAAAGoAAAAAAAAAawAAAAAAAABsAAAAAAAAAG0AAAAAAAAAbgAAAAAAAABvAAAAAAAAAHAAAAAAAAAAcQAAAAAAAAByAAAAAAAAAHMAAAAAAAAAdAAAAAAAAAB1AAAAAAAAAHYAAAAAAAAAdwAAAAAAAAB4AAAAAAAAAHkAAAAAAAAAegAAAAAAAAB7AAAAAAAAAHwAAAAAAAAAfQAAAAAAAAB+AAAAAAAAAH8AAAAAAAAAgAAAAAAAAACBAAAAAAAAAIIAAAAAAAAAgwAAAAAAAACEAAAAAAAAAIUAAAAAAAAAhgAAAAAAAACHAAAAAAAAAIgAAAAAAAAAiQAAAAAAAACKAAAAAAAAAIsAAAAAAAAAjAAAAAAAAACNAAAAAAAAAI4AAAAAAAAAjwAAAAAAAACQAAAAAAAAAJEAAAAAAAAAkgAAAAAAAACTAAAAAAAAAJQAAAAAAAAAlQAAAAAAAACWAAAAAAAAAJcAAAAAAAAAmAAAAAAAAACZAAAAAAAAAJoAAAAAAAAAmwAAAAAAAACcAAAAAAAAAJ0AAAAAAAAAngAAAAAAAACfAAAAAAAAAKAAAAAAAAAAoQAAAAAAAACiAAAAAAAAAKMAAAAAAAAApAAAAAAAAAClAAAAAAAAAKYAAAAAAAAApwAAAAAAAACoAAAAAAAAAKkAAAAAAAAAqgAAAAAAAACrAAAAAAAAAKwAAAAAAAAArQAAAAAAAACuAAAAAAAAAK8AAAAAAAAAsAAAAAAAAACxAAAAAAAAALIAAAAAAAAAswAAAAAAAAC0AAAAAAAAALUAAAAAAAAAtgAAAAAAAAC3AAAAAAAAALgAAAAAAAAAuQAAAAAAAAC6AAAAAAAAALsAAAAAAAAAvAAAAAAAAAC9AAAAAAAAAL4AAAAAAAAAvwAAAAAAAADAAAAAAAAAAMEAAAAAAAAAwgAAAAAAAADDAAAAAAAAAMQAAAAAAAAAxQAAAAAAAADGAAAAAAAAAMcAAAAAAAAAyAAAAAAAAADJAAAAAAAAAA==IAMAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAEwAAAAAAAAATQAAAAAAAABOAAAAAAAAAE8AAAAAAAAAUAAAAAAAAABRAAAAAAAAAFIAAAAAAAAAUwAAAAAAAABUAAAAAAAAAFUAAAAAAAAAVgAAAAAAAABXAAAAAAAAAFgAAAAAAAAAWQAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABdAAAAAAAAAF4AAAAAAAAAXwAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYgAAAAAAAABjAAAAAAAAAA==8BIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACamZmZmZm5PwAAAAAAAAAAAAAAAAAAAACamZmZmZnJPwAAAAAAAAAAAAAAAAAAAAA0MzMzMzPTPwAAAAAAAAAAAAAAAAAAAACamZmZmZnZPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADgPwAAAAAAAAAAAAAAAAAAAAA0MzMzMzPjPwAAAAAAAAAAAAAAAAAAAABnZmZmZmbmPwAAAAAAAAAAAAAAAAAAAACamZmZmZnpPwAAAAAAAAAAAAAAAAAAAADNzMzMzMzsPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAACamZmZmZnxPwAAAAAAAAAAAAAAAAAAAAA0MzMzMzPzPwAAAAAAAAAAAAAAAAAAAADNzMzMzMz0PwAAAAAAAAAAAAAAAAAAAABnZmZmZmb2PwAAAAAAAAAAAAAAAAAAAAAAAAAAAAD4PwAAAAAAAAAAAAAAAAAAAACamZmZmZn5PwAAAAAAAAAAAAAAAAAAAAA0MzMzMzP7PwAAAAAAAAAAAAAAAAAAAADNzMzMzMz8PwAAAAAAAAAAAAAAAAAAAABnZmZmZmb+PwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwAQAAAAAAAAAAAAAAAAAAAAACamZmZmZkBQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYCQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMDQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAEQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwEQAAAAAAAAAAAAAAAAAAAAACamZmZmZkFQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYGQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMHQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwIQAAAAAAAAAAAAAAAAAAAAACamZmZmZkJQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYKQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMLQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAMQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwMQAAAAAAAAAAAAAAAAAAAAACamZmZmZkNQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYOQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMPQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAQQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYQQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwQQAAAAAAAAAAAAAAAAAAAAAAzMzMzMzMRQAAAAAAAAAAAAAAAAAAAAACamZmZmZkRQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAASQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYSQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwSQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMTQAAAAAAAAAAAAAAAAAAAAACamZmZmZkTQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAUQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYUQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwUQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMVQAAAAAAAAAAAAAAAAAAAAACamZmZmZkVQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAWQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYWQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwWQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMXQAAAAAAAAAAAAAAAAAAAAACamZmZmZkXQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAYQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYYQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwYQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMZQAAAAAAAAAAAAAAAAAAAAACamZmZmZkZQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAaQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYaQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwaQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMbQAAAAAAAAAAAAAAAAAAAAACamZmZmZkbQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAcQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYcQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwcQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMdQAAAAAAAAAAAAAAAAAAAAACamZmZmZkdQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAeQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYeQAAAAAAAAAAAAAAAAAAAAADNzMzMzMweQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMfQAAAAAAAAAAAAAAAAAAAAACamZmZmZkfQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAgQAAAAAAAAAAAAAAAAAAAAAAzMzMzMzMgQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYgQAAAAAAAAAAAAAAAAAAAAACamZmZmZkgQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwgQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAhQAAAAAAAAAAAAAAAAAAAAAAzMzMzMzMhQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYhQAAAAAAAAAAAAAAAAAAAAACamZmZmZkhQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwhQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAiQAAAAAAAAAAAAAAAAAAAAAAzMzMzMzMiQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYiQAAAAAAAAAAAAAAAAAAAAACamZmZmZkiQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwiQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAjQAAAAAAAAAAAAAAAAAAAAAA0MzMzMzMjQAAAAAAAAAAAAAAAAAAAAABnZmZmZmYjQAAAAAAAAAAAAAAAAAAAAACamZmZmZkjQAAAAAAAAAAAAAAAAAAAAADNzMzMzMwjQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAkQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAACamZmZmZm5PwAAAAAAAPA/AAAAAAAAAACamZmZmZnJPwAAAAAAAPA/AAAAAAAAAAA0MzMzMzPTPwAAAAAAAPA/AAAAAAAAAACamZmZmZnZPwAAAAAAAPA/AAAAAAAAAAAAAAAAAADgPwAAAAAAAPA/AAAAAAAAAAA0MzMzMzPjPwAAAAAAAPA/AAAAAAAAAABnZmZmZmbmPwAAAAAAAPA/AAAAAAAAAACamZmZmZnpPwAAAAAAAPA/AAAAAAAAAADNzMzMzMzsPwAAAAAAAPA/AAAAAAAAAAAAAAAAAADwPwAAAAAAAPA/AAAAAAAAAACamZmZmZnxPwAAAAAAAPA/AAAAAAAAAAA0MzMzMzPzPwAAAAAAAPA/AAAAAAAAAADNzMzMzMz0PwAAAAAAAPA/AAAAAAAAAABnZmZmZmb2PwAAAAAAAPA/AAAAAAAAAAAAAAAAAAD4PwAAAAAAAPA/AAAAAAAAAACamZmZmZn5PwAAAAAAAPA/AAAAAAAAAAA0MzMzMzP7PwAAAAAAAPA/AAAAAAAAAADNzMzMzMz8PwAAAAAAAPA/AAAAAAAAAABnZmZmZmb+PwAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwAQAAAAAAAAPA/AAAAAAAAAACamZmZmZkBQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYCQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMDQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAEQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwEQAAAAAAAAPA/AAAAAAAAAACamZmZmZkFQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYGQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMHQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAIQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwIQAAAAAAAAPA/AAAAAAAAAACamZmZmZkJQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYKQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMLQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAMQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwMQAAAAAAAAPA/AAAAAAAAAACamZmZmZkNQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYOQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMPQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAQQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYQQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwQQAAAAAAAAPA/AAAAAAAAAAAzMzMzMzMRQAAAAAAAAPA/AAAAAAAAAACamZmZmZkRQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAASQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYSQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwSQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMTQAAAAAAAAPA/AAAAAAAAAACamZmZmZkTQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAUQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYUQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwUQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMVQAAAAAAAAPA/AAAAAAAAAACamZmZmZkVQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAWQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYWQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwWQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMXQAAAAAAAAPA/AAAAAAAAAACamZmZmZkXQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAYQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYYQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwYQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMZQAAAAAAAAPA/AAAAAAAAAACamZmZmZkZQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAaQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYaQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwaQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMbQAAAAAAAAPA/AAAAAAAAAACamZmZmZkbQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAcQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYcQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwcQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMdQAAAAAAAAPA/AAAAAAAAAACamZmZmZkdQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAeQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYeQAAAAAAAAPA/AAAAAAAAAADNzMzMzMweQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMfQAAAAAAAAPA/AAAAAAAAAACamZmZmZkfQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAgQAAAAAAAAPA/AAAAAAAAAAAzMzMzMzMgQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYgQAAAAAAAAPA/AAAAAAAAAACamZmZmZkgQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwgQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAhQAAAAAAAAPA/AAAAAAAAAAAzMzMzMzMhQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYhQAAAAAAAAPA/AAAAAAAAAACamZmZmZkhQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwhQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAiQAAAAAAAAPA/AAAAAAAAAAAzMzMzMzMiQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYiQAAAAAAAAPA/AAAAAAAAAACamZmZmZkiQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwiQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAjQAAAAAAAAPA/AAAAAAAAAAA0MzMzMzMjQAAAAAAAAPA/AAAAAAAAAABnZmZmZmYjQAAAAAAAAPA/AAAAAAAAAACamZmZmZkjQAAAAAAAAPA/AAAAAAAAAADNzMzMzMwjQAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAkQAAAAAAAAPA/AAAAAAAAAAA=gAwAAAAAAAAAAAAAAAAAAAEAAAAAAAAAZgAAAAAAAABlAAAAAAAAAAEAAAAAAAAAAgAAAAAAAABnAAAAAAAAAGYAAAAAAAAAAgAAAAAAAAADAAAAAAAAAGgAAAAAAAAAZwAAAAAAAAADAAAAAAAAAAQAAAAAAAAAaQAAAAAAAABoAAAAAAAAAAQAAAAAAAAABQAAAAAAAABqAAAAAAAAAGkAAAAAAAAABQAAAAAAAAAGAAAAAAAAAGsAAAAAAAAAagAAAAAAAAAGAAAAAAAAAAcAAAAAAAAAbAAAAAAAAABrAAAAAAAAAAcAAAAAAAAACAAAAAAAAABtAAAAAAAAAGwAAAAAAAAACAAAAAAAAAAJAAAAAAAAAG4AAAAAAAAAbQAAAAAAAAAJAAAAAAAAAAoAAAAAAAAAbwAAAAAAAABuAAAAAAAAAAoAAAAAAAAACwAAAAAAAABwAAAAAAAAAG8AAAAAAAAACwAAAAAAAAAMAAAAAAAAAHEAAAAAAAAAcAAAAAAAAAAMAAAAAAAAAA0AAAAAAAAAcgAAAAAAAABxAAAAAAAAAA0AAAAAAAAADgAAAAAAAABzAAAAAAAAAHIAAAAAAAAADgAAAAAAAAAPAAAAAAAAAHQAAAAAAAAAcwAAAAAAAAAPAAAAAAAAABAAAAAAAAAAdQAAAAAAAAB0AAAAAAAAABAAAAAAAAAAEQAAAAAAAAB2AAAAAAAAAHUAAAAAAAAAEQAAAAAAAAASAAAAAAAAAHcAAAAAAAAAdgAAAAAAAAASAAAAAAAAABMAAAAAAAAAeAAAAAAAAAB3AAAAAAAAABMAAAAAAAAAFAAAAAAAAAB5AAAAAAAAAHgAAAAAAAAAFAAAAAAAAAAVAAAAAAAAAHoAAAAAAAAAeQAAAAAAAAAVAAAAAAAAABYAAAAAAAAAewAAAAAAAAB6AAAAAAAAABYAAAAAAAAAFwAAAAAAAAB8AAAAAAAAAHsAAAAAAAAAFwAAAAAAAAAYAAAAAAAAAH0AAAAAAAAAfAAAAAAAAAAYAAAAAAAAABkAAAAAAAAAfgAAAAAAAAB9AAAAAAAAABkAAAAAAAAAGgAAAAAAAAB/AAAAAAAAAH4AAAAAAAAAGgAAAAAAAAAbAAAAAAAAAIAAAAAAAAAAfwAAAAAAAAAbAAAAAAAAABwAAAAAAAAAgQAAAAAAAACAAAAAAAAAABwAAAAAAAAAHQAAAAAAAACCAAAAAAAAAIEAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAIMAAAAAAAAAggAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAhAAAAAAAAACDAAAAAAAAAB8AAAAAAAAAIAAAAAAAAACFAAAAAAAAAIQAAAAAAAAAIAAAAAAAAAAhAAAAAAAAAIYAAAAAAAAAhQAAAAAAAAAhAAAAAAAAACIAAAAAAAAAhwAAAAAAAACGAAAAAAAAACIAAAAAAAAAIwAAAAAAAACIAAAAAAAAAIcAAAAAAAAAIwAAAAAAAAAkAAAAAAAAAIkAAAAAAAAAiAAAAAAAAAAkAAAAAAAAACUAAAAAAAAAigAAAAAAAACJAAAAAAAAACUAAAAAAAAAJgAAAAAAAACLAAAAAAAAAIoAAAAAAAAAJgAAAAAAAAAnAAAAAAAAAIwAAAAAAAAAiwAAAAAAAAAnAAAAAAAAACgAAAAAAAAAjQAAAAAAAACMAAAAAAAAACgAAAAAAAAAKQAAAAAAAACOAAAAAAAAAI0AAAAAAAAAKQAAAAAAAAAqAAAAAAAAAI8AAAAAAAAAjgAAAAAAAAAqAAAAAAAAACsAAAAAAAAAkAAAAAAAAACPAAAAAAAAACsAAAAAAAAALAAAAAAAAACRAAAAAAAAAJAAAAAAAAAALAAAAAAAAAAtAAAAAAAAAJIAAAAAAAAAkQAAAAAAAAAtAAAAAAAAAC4AAAAAAAAAkwAAAAAAAACSAAAAAAAAAC4AAAAAAAAALwAAAAAAAACUAAAAAAAAAJMAAAAAAAAALwAAAAAAAAAwAAAAAAAAAJUAAAAAAAAAlAAAAAAAAAAwAAAAAAAAADEAAAAAAAAAlgAAAAAAAACVAAAAAAAAADEAAAAAAAAAMgAAAAAAAACXAAAAAAAAAJYAAAAAAAAAMgAAAAAAAAAzAAAAAAAAAJgAAAAAAAAAlwAAAAAAAAAzAAAAAAAAADQAAAAAAAAAmQAAAAAAAACYAAAAAAAAADQAAAAAAAAANQAAAAAAAACaAAAAAAAAAJkAAAAAAAAANQAAAAAAAAA2AAAAAAAAAJsAAAAAAAAAmgAAAAAAAAA2AAAAAAAAADcAAAAAAAAAnAAAAAAAAACbAAAAAAAAADcAAAAAAAAAOAAAAAAAAACdAAAAAAAAAJwAAAAAAAAAOAAAAAAAAAA5AAAAAAAAAJ4AAAAAAAAAnQAAAAAAAAA5AAAAAAAAADoAAAAAAAAAnwAAAAAAAACeAAAAAAAAADoAAAAAAAAAOwAAAAAAAACgAAAAAAAAAJ8AAAAAAAAAOwAAAAAAAAA8AAAAAAAAAKEAAAAAAAAAoAAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAogAAAAAAAAChAAAAAAAAAD0AAAAAAAAAPgAAAAAAAACjAAAAAAAAAKIAAAAAAAAAPgAAAAAAAAA/AAAAAAAAAKQAAAAAAAAAowAAAAAAAAA/AAAAAAAAAEAAAAAAAAAApQAAAAAAAACkAAAAAAAAAEAAAAAAAAAAQQAAAAAAAACmAAAAAAAAAKUAAAAAAAAAQQAAAAAAAABCAAAAAAAAAKcAAAAAAAAApgAAAAAAAABCAAAAAAAAAEMAAAAAAAAAqAAAAAAAAACnAAAAAAAAAEMAAAAAAAAARAAAAAAAAACpAAAAAAAAAKgAAAAAAAAARAAAAAAAAABFAAAAAAAAAKoAAAAAAAAAqQAAAAAAAABFAAAAAAAAAEYAAAAAAAAAqwAAAAAAAACqAAAAAAAAAEYAAAAAAAAARwAAAAAAAACsAAAAAAAAAKsAAAAAAAAARwAAAAAAAABIAAAAAAAAAK0AAAAAAAAArAAAAAAAAABIAAAAAAAAAEkAAAAAAAAArgAAAAAAAACtAAAAAAAAAEkAAAAAAAAASgAAAAAAAACvAAAAAAAAAK4AAAAAAAAASgAAAAAAAABLAAAAAAAAALAAAAAAAAAArwAAAAAAAABLAAAAAAAAAEwAAAAAAAAAsQAAAAAAAACwAAAAAAAAAEwAAAAAAAAATQAAAAAAAACyAAAAAAAAALEAAAAAAAAATQAAAAAAAABOAAAAAAAAALMAAAAAAAAAsgAAAAAAAABOAAAAAAAAAE8AAAAAAAAAtAAAAAAAAACzAAAAAAAAAE8AAAAAAAAAUAAAAAAAAAC1AAAAAAAAALQAAAAAAAAAUAAAAAAAAABRAAAAAAAAALYAAAAAAAAAtQAAAAAAAABRAAAAAAAAAFIAAAAAAAAAtwAAAAAAAAC2AAAAAAAAAFIAAAAAAAAAUwAAAAAAAAC4AAAAAAAAALcAAAAAAAAAUwAAAAAAAABUAAAAAAAAALkAAAAAAAAAuAAAAAAAAABUAAAAAAAAAFUAAAAAAAAAugAAAAAAAAC5AAAAAAAAAFUAAAAAAAAAVgAAAAAAAAC7AAAAAAAAALoAAAAAAAAAVgAAAAAAAABXAAAAAAAAALwAAAAAAAAAuwAAAAAAAABXAAAAAAAAAFgAAAAAAAAAvQAAAAAAAAC8AAAAAAAAAFgAAAAAAAAAWQAAAAAAAAC+AAAAAAAAAL0AAAAAAAAAWQAAAAAAAABaAAAAAAAAAL8AAAAAAAAAvgAAAAAAAABaAAAAAAAAAFsAAAAAAAAAwAAAAAAAAAC/AAAAAAAAAFsAAAAAAAAAXAAAAAAAAADBAAAAAAAAAMAAAAAAAAAAXAAAAAAAAABdAAAAAAAAAMIAAAAAAAAAwQAAAAAAAABdAAAAAAAAAF4AAAAAAAAAwwAAAAAAAADCAAAAAAAAAF4AAAAAAAAAXwAAAAAAAADEAAAAAAAAAMMAAAAAAAAAXwAAAAAAAABgAAAAAAAAAMUAAAAAAAAAxAAAAAAAAABgAAAAAAAAAGEAAAAAAAAAxgAAAAAAAADFAAAAAAAAAGEAAAAAAAAAYgAAAAAAAADHAAAAAAAAAMYAAAAAAAAAYgAAAAAAAABjAAAAAAAAAMgAAAAAAAAAxwAAAAAAAABjAAAAAAAAAGQAAAAAAAAAyQAAAAAAAADIAAAAAAAAAA==IAMAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAAAAQAAAAAAAAQBAAAAAAAACAEAAAAAAAAMAQAAAAAAABABAAAAAAAAFAEAAAAAAAAYAQAAAAAAABwBAAAAAAAAIAEAAAAAAAAkAQAAAAAAACgBAAAAAAAALAEAAAAAAAAwAQAAAAAAADQBAAAAAAAAOAEAAAAAAAA8AQAAAAAAAEABAAAAAAAARAEAAAAAAABIAQAAAAAAAEwBAAAAAAAAUAEAAAAAAABUAQAAAAAAAFgBAAAAAAAAXAEAAAAAAABgAQAAAAAAAGQBAAAAAAAAaAEAAAAAAABsAQAAAAAAAHABAAAAAAAAdAEAAAAAAAB4AQAAAAAAAHwBAAAAAAAAgAEAAAAAAACEAQAAAAAAAIgBAAAAAAAAjAEAAAAAAACQAQAAAAAAAA==ZAAAAAAAAAAJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJ
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh_left.vtu b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh_left.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..5db5370868791affd4c3332930f1fd343c8af8a0
--- /dev/null
+++ b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh_left.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="2" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="101" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="32" />
+ <DataArray type="UInt64" Name="bulk_face_ids" format="appended" RangeMin="3" RangeMax="3" offset="56" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="80" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="156" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="188" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="212" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _EAAAAAAAAAAAAAAAAAAAAGUAAAAAAAAACAAAAAAAAAAAAAAAAAAAAA==CAAAAAAAAAADAAAAAAAAAA==MAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAA=EAAAAAAAAAABAAAAAAAAAAAAAAAAAAAACAAAAAAAAAACAAAAAAAAAA==AQAAAAAAAAAD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh_right.vtu b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh_right.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..915396a9e71d2fefeb6cd22d22469164ddfd895c
--- /dev/null
+++ b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/mesh_right.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="2" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="100" RangeMax="201" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="99" RangeMax="99" offset="32" />
+ <DataArray type="UInt64" Name="bulk_face_ids" format="appended" RangeMin="1" RangeMax="1" offset="56" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="10" RangeMax="10.049875621" offset="80" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="156" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="188" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="212" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _EAAAAAAAAABkAAAAAAAAAMkAAAAAAAAACAAAAAAAAABjAAAAAAAAAA==CAAAAAAAAAABAAAAAAAAAA==MAAAAAAAAAAAAAAAAAAkQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAkQAAAAAAAAPA/AAAAAAAAAAA=EAAAAAAAAAAAAAAAAAAAAAEAAAAAAAAACAAAAAAAAAACAAAAAAAAAA==AQAAAAAAAAAD
+ </AppendedData>
+</VTKFile>