diff --git a/Documentation/ProjectFile/prj/processes/process/RICHARDS_MECHANICS/t_mass_lumping.md b/Documentation/ProjectFile/prj/processes/process/RICHARDS_MECHANICS/t_mass_lumping.md
new file mode 100644
index 0000000000000000000000000000000000000000..01fff2a4ff88f5086ca1f730f1d3dcacc33ffe12
--- /dev/null
+++ b/Documentation/ProjectFile/prj/processes/process/RICHARDS_MECHANICS/t_mass_lumping.md
@@ -0,0 +1 @@
+A tag for enabling diagonal lumping of the mass matrix.
diff --git a/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp b/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
index 72edbee59de593f14c3f991777bcc09eed099c71..f22c601ea16313d8bfea91ded6b448da5ba27156 100644
--- a/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
+++ b/ProcessLib/RichardsMechanics/CreateRichardsMechanicsProcess.cpp
@@ -188,6 +188,10 @@ std::unique_ptr<Process> createRichardsMechanicsProcess(
ProcessLib::RichardsFlow::createRichardsFlowMaterialProperties(
flow_material_config, material_ids, parameters);
+ auto const mass_lumping =
+ //! \ogs_file_param{prj__processes__process__RICHARDS_MECHANICS__mass_lumping}
+ config.getConfigParameter<bool>("mass_lumping", false);
+
RichardsMechanicsProcessData<DisplacementDim> process_data{
material_ids,
std::move(flow_material),
@@ -197,7 +201,8 @@ std::unique_ptr<Process> createRichardsMechanicsProcess(
solid_density,
solid_bulk_modulus,
temperature,
- specific_body_force};
+ specific_body_force,
+ mass_lumping};
SecondaryVariableCollection secondary_variables;
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index 6cfa2b4c248bffcb34084ef6fda7fa1e730ec158..5656d8d15382df0804820291219174c7c7524f28 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -281,6 +281,13 @@ void RichardsMechanicsLocalAssembler<
.noalias() += N_p.transpose() * S_L * rho_LR * alpha *
identity2.transpose() * B * w;
}
+
+ if (_process_data.apply_mass_lumping)
+ {
+ auto Mpp = M.template block<pressure_size, pressure_size>(
+ pressure_index, pressure_index);
+ Mpp = Mpp.colwise().sum().eval().asDiagonal();
+ }
}
template <typename ShapeFunctionDisplacement, typename ShapeFunctionPressure,
@@ -553,6 +560,11 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
dNdx_p.transpose() * rho_LR * k_rel * rho_Ki_over_mu * b * w;
}
+ if (_process_data.apply_mass_lumping)
+ {
+ storage_p = storage_p.colwise().sum().eval().asDiagonal();
+ }
+
// pressure equation, pressure part.
local_Jac
.template block<pressure_size, pressure_size>(pressure_index,
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsProcessData.h b/ProcessLib/RichardsMechanics/RichardsMechanicsProcessData.h
index 0ea1097941182905676b0d20e47af91675379b46..b12ca3cfaf43cb1dbcb6b097e694b7d67fa70036 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsProcessData.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsProcessData.h
@@ -48,7 +48,8 @@ struct RichardsMechanicsProcessData
ParameterLib::Parameter<double> const& solid_bulk_modulus_,
ParameterLib::Parameter<double> const& temperature_,
Eigen::Matrix<double, DisplacementDim, 1>
- specific_body_force_)
+ specific_body_force_,
+ bool const apply_mass_lumping_)
: material_ids(material_ids_),
flow_material{std::move(flow_material_)},
solid_materials{std::move(solid_materials_)},
@@ -57,7 +58,8 @@ struct RichardsMechanicsProcessData
solid_density(solid_density_),
solid_bulk_modulus(solid_bulk_modulus_),
temperature(temperature_),
- specific_body_force(std::move(specific_body_force_))
+ specific_body_force(std::move(specific_body_force_)),
+ apply_mass_lumping(apply_mass_lumping_)
{
}
@@ -105,6 +107,8 @@ struct RichardsMechanicsProcessData
MeshLib::PropertyVector<double>* element_saturation = nullptr;
MeshLib::PropertyVector<double>* pressure_interpolated = nullptr;
+ bool const apply_mass_lumping;
+
EIGEN_MAKE_ALIGNED_OPERATOR_NEW;
};
diff --git a/ProcessLib/RichardsMechanics/Tests.cmake b/ProcessLib/RichardsMechanics/Tests.cmake
index d318ed709b28111518ed306895169847b3da4cd1..39ec1743fd21c5479833ebe5b7d89ff64e590e28 100644
--- a/ProcessLib/RichardsMechanics/Tests.cmake
+++ b/ProcessLib/RichardsMechanics/Tests.cmake
@@ -130,6 +130,26 @@ AddTest(
GLOB Richards_2d_geometry_OBSERVATION_POINT_pcs_0_ts_*.vtu displacement displacement 2e-14 0
GLOB Richards_2d_geometry_OBSERVATION_POINT_pcs_0_ts_*.vtu pressure pressure 1e-7 1e-15
)
+AddTest(
+ NAME RichardsMechanics_RichardsFlow_2d_small_masslumping
+ PATH RichardsMechanics
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS RichardsFlow_2d_small_masslumping.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu displacement displacement 2e-14 0
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu sigma sigma 1e-8 0
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu epsilon epsilon 1e-15 0
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu pressure pressure 1e-7 1e-15
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu pressure_interpolated pressure_interpolated 1e-7 1e-15
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu saturation saturation 1e-11 1e-15
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu saturation_avg saturation_avg 1e-11 1e-15
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu velocity velocity 1e-15 1e-15
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu HydraulicFlow HydraulicFlow 1e-13 0
+ GLOB RichardsFlow_2d_small_masslumping_pcs_0_ts_*.vtu NodalForces NodalForces 1e-9 0
+)
AddTest(
NAME RichardsMechanics_RichardsFlow_2d_quasinewton
PATH RichardsMechanics
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj
new file mode 100644
index 0000000000000000000000000000000000000000..0422929797780966ad4e6c92c7410cbff8e83270
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj
@@ -0,0 +1,291 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>Richards_2d.vtu</mesh>
+ <geometry>Richards_2d.gml</geometry>
+ <processes>
+ <process>
+ <name>RM</name>
+ <type>RICHARDS_MECHANICS</type>
+ <integration_order>3</integration_order>
+ <dimension>2</dimension>
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <solid_bulk_modulus>K_SR</solid_bulk_modulus>
+ <fluid_bulk_modulus>K_LR</fluid_bulk_modulus>
+ <biot_coefficient>alpha</biot_coefficient>
+ <solid_density>rho_sr</solid_density>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <pressure>pressure</pressure>
+ </process_variables>
+ <material_property>
+ <fluid>
+ <density>
+ <type>Constant</type>
+ <value> 1 </value>
+ </density>
+ <viscosity>
+ <type>Constant</type>
+ <value> 1.e-3 </value>
+ </viscosity>
+ </fluid>
+ <porous_medium>
+ <porous_medium id="0">
+ <permeability>
+ <permeability_tensor_entries>k</permeability_tensor_entries>
+ <type>Constant</type>
+ </permeability>
+ <porosity>
+ <type>Constant</type>
+ <porosity_parameter>phi</porosity_parameter>
+ </porosity>
+ <storage>
+ <type>Constant</type>
+ <value> 0.0 </value>
+ </storage>
+ <capillary_pressure>
+ <type>vanGenuchten</type>
+ <m>0.789029535864979</m>
+ <pc_max>36333.30</pc_max>
+ <pd>3633.33</pd>
+ <smax>0.95</smax>
+ <sr>0.1689</sr>
+ </capillary_pressure>
+ <relative_permeability>
+ <type>WettingPhaseVanGenuchten</type>
+ <krel_min>1e-12</krel_min>
+ <m>0.789029535864979</m>
+ <smax>0.95</smax>
+ <sr>0.1689</sr>
+ </relative_permeability>
+ </porous_medium>
+ </porous_medium>
+ </material_property>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="sigma" output_name="sigma"/>
+ <secondary_variable type="static" internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable type="static" internal_name="velocity" output_name="velocity"/>
+ <secondary_variable type="static" internal_name="saturation" output_name="saturation"/>
+ </secondary_variables>
+ <specific_body_force>0 -9.81</specific_body_force>
+ <mass_lumping>true</mass_lumping>
+ <temperature> temp </temperature>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="RM">
+ <nonlinear_solver>nonlinear_solver</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-7 1e-11 1e-10</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>2000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>2</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>5</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>10</delta_t>
+ </pair>
+ <pair>
+ <repeat>50</repeat>
+ <delta_t>20</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>RichardsFlow_2d_small_masslumping</prefix>
+ <timesteps>
+ <pair>
+ <repeat>9</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ <pair>
+ <repeat>100</repeat>
+ <each_steps>10</each_steps>
+ </pair>
+ </timesteps>
+ <output_iteration_results>false</output_iteration_results>
+ <variables>
+ <variable>displacement</variable>
+ <variable>pressure</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity</variable>
+ <variable>saturation</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <!-- Mechanics -->
+ <parameter>
+ <name>K_SR</name>
+ <type>Constant</type>
+ <value>160e300</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>5e9</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.1</value>
+ </parameter>
+ <!-- Model parameters -->
+ <parameter>
+ <name>k</name>
+ <type>Constant</type>
+ <value>4.46e-13</value>
+ </parameter>
+ <parameter>
+ <name>alpha</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>phi</name>
+ <type>Constant</type>
+ <value>0.38</value>
+ </parameter>
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>K_LR</name>
+ <type>Constant</type>
+ <value>2.2e300</value>
+ </parameter>
+ <parameter>
+ <name>temp</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>dirichlet0</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ <parameter>
+ <name>p0</name>
+ <type>Constant</type>
+ <value>-5000</value>
+ </parameter>
+ <parameter>
+ <name>p_Dirichlet_bottom</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>2</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>Richards_2d_geometry</geometrical_set>
+ <geometry>LEFT</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>Richards_2d_geometry</geometrical_set>
+ <geometry>RIGHT</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>Richards_2d_geometry</geometrical_set>
+ <geometry>BOTTOM</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>Richards_2d_geometry</geometrical_set>
+ <geometry>TOP</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>p0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>Richards_2d_geometry</geometrical_set>
+ <geometry>BOTTOM</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_Dirichlet_bottom</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>Richards_2d_geometry</geometrical_set>
+ <geometry>TOP</geometry>
+ <type>Dirichlet</type>
+ <parameter>p0</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>nonlinear_solver</name>
+ <type>Newton</type>
+ <max_iter>100</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ <!--damping>0.9</damping-->
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ <scaling>true</scaling>
+ <!--
+ <solver_type>PardisoLU</solver_type>
+ -->
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_0_t_0.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..48c18f3419c8c95e9e96ff7ea7b76b324dcc3afd
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_0_t_0.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4e141d24f636f15fda9732dc31c9cd1aa17c6e1c895c000dd271241c62dcb725
+size 8735
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_104_t_2000.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_104_t_2000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..eb65da22434a58c00e6e4a11c865d592d89f8b56
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_104_t_2000.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:98f8f70779afa8f0c959622fe2be4d06d4c43e225037020d581aea22bf40c096
+size 41155
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_19_t_318.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_19_t_318.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..e372a4d8e48807eac3f330d07365b77fcc4f46aa
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_19_t_318.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3103596ab508ffdc503875d9205f15ed518fe943602ab94c569fa33d10e68589
+size 40255
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_1_t_1.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_1_t_1.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..2569e9392d837184d74da16df40daac6ec0f6ac6
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_1_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:5eb924b177d7865986dc8203559a034684b251481fc1aa54fa191a77bb1c9a29
+size 39503
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_29_t_518.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_29_t_518.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..73d86903ea5172d5fef6f891cb671c8c0a208fee
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_29_t_518.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4f041203838263b72e672519f357cfa59a962817576df9bb637a9ad35f674829
+size 40539
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_2_t_3.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_2_t_3.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..3a03f3a9cad161187282a8a9984dcf3686de3e54
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_2_t_3.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:cfe3c9e5844c035ed21aae97412563ae24b26bf28679e516f244248feec1cbc8
+size 39535
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_39_t_718.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_39_t_718.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..eded4fc2ba8ae775e3f0a2176a40ed9c2de88af0
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_39_t_718.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:6450e7e3b8b21638f0c841d43d344263a2a4cb0a9131aa9d8adcb1b2e71b53b4
+size 40687
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_3_t_8.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_3_t_8.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d589598045362b005aad2cb0f3ea423f6d354fee
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_3_t_8.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:ecea9f5bc566a259c99d86c2078afe7d67d5b6915479717d1325bd022cbdf93e
+size 39635
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_49_t_918.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_49_t_918.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..6c483c3c7f877b394d663a48bf6b160544d9df11
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_49_t_918.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:9f9e6e14db286908c5645ffd4c5ffe834a8e9e84a96dc6a19b2aad1ca8bb5a3d
+size 40703
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_4_t_18.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_4_t_18.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d9c948e44385813e2b32d3af29686a70ece043aa
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_4_t_18.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:7284e3f3f2b3314f3379844316306ee7ba1ef2ce8de5832023654af2628f6f2e
+size 39743
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_59_t_1118.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_59_t_1118.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..a94bd7d1118092580b4d266d0cf9d4107b81e546
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_59_t_1118.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:99c5b702a22e2006eb01a6d87ed2a161f6af5dff96adaa9a79544d838612d64b
+size 41007
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_5_t_38.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_5_t_38.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..eb5e2cc7595c731b03200662a4980bb9475ce627
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_5_t_38.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f2446df14b368485382f1cf28b0d9077318350b9c6151ead7b37b9bc54a7c0fd
+size 40007
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_69_t_1318.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_69_t_1318.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..df503293b833cc7c6b26a996f29369494e867fa9
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_69_t_1318.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:44835c4f4b4bc0de12edb89ed77f9c4a4ff582eb0e6686153e517f154f0d94bd
+size 40739
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_6_t_58.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_6_t_58.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9dbf29f9b28d909a374ab403ca43839d9df733a7
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_6_t_58.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:fb5ec11fbe3d41eed9ce5dd200336cc126e4adb271839d6337a2774df213cd5b
+size 39855
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_79_t_1518.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_79_t_1518.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..e0eef6e10423498f388a0fe428d8f2ab7bd1c7b2
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_79_t_1518.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:6aa5bd7973b7fa43d98c371cbbb727a35ecbf7c59a244f55269524063614fe16
+size 41099
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_7_t_78.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_7_t_78.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..6cf3ea7778a096eef547b5a559c7cf3e7e6ae30b
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_7_t_78.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:33a4b36e4bc7f660e3b38d393b06b3d89a48ad5137d5033209e3130c99c4eb1a
+size 39935
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_89_t_1718.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_89_t_1718.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..cdb77a343ac9ef1766a90e8260308325469a9e34
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_89_t_1718.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:99b1448c7fd315a546cef496557e9a3a9c30593c46581f8aa1d6a98e2bafb9fa
+size 41087
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_8_t_98.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_8_t_98.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..26102387039959151f09dfc808fba612fa9d67e1
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_8_t_98.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:39baad583ea33da47d23709c1e084b82384639eddda8f4e8f9442ce9c31d5bd9
+size 40043
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_99_t_1918.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_99_t_1918.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..4b1f6eed0c88a0be1e0b6599db4eb0d6253a640e
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_99_t_1918.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:715743b11ec79fbec93460658167dc9c959998601b680353bd0827145aec74d1
+size 41091
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_9_t_118.000000.vtu b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_9_t_118.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..73f68703a910f93b51d4b79e83d55100fd2ab79f
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping_pcs_0_ts_9_t_118.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8387630a2e20c59e64401aaba0f08eb63ba78e3708c63bbfae145932036b362b
+size 40151