[T/THM] Added a test case for 0D volumetric source terms

ProcessLib/ThermoHydroMechanics/Tests.cmake

@@ -125,6 +125,23 @@ AddTest(
     expected_pointheatsource_linear-mesh_ts_10_t_50000.000000.vtu pointheatsource_linear-mesh_ts_10_t_50000.000000.vtu epsilon epsilon 1e-5 1e-5
     expected_pointheatsource_linear-mesh_ts_10_t_50000.000000.vtu pointheatsource_linear-mesh_ts_10_t_50000.000000.vtu sigma sigma 1e-5 1e-5
 )
+# ThermoHydroMechanics; Small deformation, linear poroelastic, point heat source consolidation, variation with a volumetric source term
+AddTest(
+    NAME ThermoHydroMechanics_point_heat_injection_vol_st
+    PATH ThermoHydroMechanics/Linear/Point_injection/with-volumetric-source-term
+    RUNTIME 45
+    EXECUTABLE ogs
+    EXECUTABLE_ARGS pointheatsource_quadratic-mesh.xml
+    WRAPPER time
+    TESTER vtkdiff
+    REQUIREMENTS NOT OGS_USE_MPI
+    DIFF_DATA
+    expected_pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu displacement displacement 1e-5 1e-5
+    expected_pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu pressure pressure 1e-5 1e-5
+    expected_pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu temperature temperature 1e-5 1e-5
+    expected_pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu epsilon epsilon 1e-5 1e-5
+    expected_pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu sigma sigma 1e-5 1e-5
+)
 # ThermoHydroMechanics; Small deformation, linear elastic, porosity=0, anisotropic thermal expansion
 AddTest(
     NAME ThermoHydroMechanics_cube_ortho-thermal-expansion-phi0

Tests/Data/ThermoHydroMechanics/Linear/Point_injection/with-volumetric-source-term/expected_pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu

+../expected_pointheatsource_quadratic-mesh_ts_10_t_50000.000000.vtu
\ No newline at end of file

Tests/Data/ThermoHydroMechanics/Linear/Point_injection/with-volumetric-source-term/pointheatsource_quadratic-mesh.xml

+<?xml version='1.0' encoding='ISO-8859-1'?>
+<OpenGeoSysProjectDiff base_file="../pointheatsource_quadratic-mesh.prj">
+    <replace sel="//source_term/type/text()">
+        Volumetric
+    </replace>
+
+    <!-- axial symmetry would multiply with an integral measure of 0 at the center point, thereby disabling the heat source -->
+    <remove sel="/*/meshes/mesh[text()='quarter_circle_geometry_center.vtu']/@axially_symmetric" />
+</OpenGeoSysProjectDiff>

Tests/Data/ThermoHydroMechanics/Linear/Point_injection/with-volumetric-source-term/quarter_002_2nd.vtu

+../quarter_002_2nd.vtu
\ No newline at end of file

Tests/Data/ThermoHydroMechanics/Linear/Point_injection/with-volumetric-source-term/quarter_circle_geometry_bottom.vtu

+../quarter_circle_geometry_bottom.vtu
\ No newline at end of file

Tests/Data/ThermoHydroMechanics/Linear/Point_injection/with-volumetric-source-term/quarter_circle_geometry_center.vtu

+../quarter_circle_geometry_center.vtu
\ No newline at end of file

Tests/Data/ThermoHydroMechanics/Linear/Point_injection/with-volumetric-source-term/quarter_circle_geometry_left.vtu

+../quarter_circle_geometry_left.vtu
\ No newline at end of file

Tests/Data/ThermoHydroMechanics/Linear/Point_injection/with-volumetric-source-term/quarter_circle_geometry_out.vtu

+../quarter_circle_geometry_out.vtu
\ No newline at end of file