[CTest/TRM] Reused BodyForce as a test for the total initial stress

input

ProcessLib/ThermoRichardsMechanics/Tests.cmake

@@ -448,6 +448,22 @@ AddTest(
     expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_ts_10_t_864000.000000.vtu sigma sigma 5e-2 1e-8
 )
 
+AddTest(
+    NAME ThermoRichardsMechanics_total_initial_stress_dont_apply_body_force_for_deformation
+    PATH ThermoRichardsMechanics/BodyForce
+    RUNTIME 1
+    EXECUTABLE ogs
+    EXECUTABLE_ARGS square_total_stress_test.xml
+    WRAPPER time
+    TESTER vtkdiff
+    DIFF_DATA
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu displacement displacement 1e-10 1e-10
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu pressure pressure 5e-2 1e-8
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu temperature temperature 1e-10 1e-10
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu epsilon epsilon 1e-10 1e-10
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu sigma sigma 5e-2 1e-8
+)
+
 if(OGS_USE_MFRONT)
     OgsTest(PROJECTFILE ThermoRichardsMechanics/MultiMaterialEhlers/square_1e1_2_matIDs.prj RUNTIME 1)
     OgsTest(PROJECTFILE ThermoRichardsMechanics/MultiMaterialEhlers/square_1e1_2_matIDs_restart.prj RUNTIME 1)

Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj

@@ -344,6 +344,14 @@
             <eigen>
                 <solver_type>SparseLU</solver_type>
             </eigen>
+            <petsc>
+                <parameters>
+                    -ksp_type bcgs
+                    -pc_type jacobi
+                    -ksp_rtol 1.e-20 -ksp_atol 1.e-10
+                    -ksp_max_it 4000
+                </parameters>
+            </petsc>
         </linear_solver>
     </linear_solvers>
 </OpenGeoSysProject>

Tests/Data/ThermoRichardsMechanics/BodyForce/square_total_stress_test.xml

+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProjectDiff base_file="square.prj">
+    <remove sel="/*/processes/process[1]/initial_stress"/>
+    <add sel="/*/processes/process[1]">
+        <initial_stress type = "total">total_stress0</initial_stress>
+    </add>
+    <add sel="/*/parameters">
+        <parameter>
+            <name>total_stress0</name>
+            <type>Constant</type>
+            <values>-4.5e6 -2.5e6 -6.0e6 0.0</values>
+        </parameter>
+    </add>
+    <replace msel="/*/time_loop/output/prefix/text()">
+        dont_apply_body_force_for_deformation_total_stess0_test
+    </replace>
+</OpenGeoSysProjectDiff>