TST: Improve convergence with higher friction coefficient.

pmgbergen · Jan 22, 2025 · be4d693 · be4d693
1 parent 99e8d35
commit be4d693
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 15 deletions.
diff --git a/src/porepy/applications/test_utils/reference_dense_arrays.py b/src/porepy/applications/test_utils/reference_dense_arrays.py
@@ -7961,31 +7961,31 @@
     "test_anisotropic_damage_values": {
         "friction_damage": np.array(
             [
-                [1.08964648],
-                [1.08990014],
-                [1.08042241],
-                [1.08088030],
-                [1.07850799],
+                [1.09138658],
+                [1.09153955],
+                [1.08421782],
+                [1.08452635],
+                [1.08158935],
             ]
         ),
         # Dilation damage is the same as friction damage due to identical
         # parameters.
         "dilation_damage": np.array(
             [
-                [1.08964648],
-                [1.08990014],
-                [1.08042241],
-                [1.0808803],
-                [1.07850799],
+                [1.09138658],
+                [1.09153955],
+                [1.08421782],
+                [1.08452635],
+                [1.08158935],
             ]
         ),
         "damage_history": np.array(
             [
-                [0.21859257],
-                [0.21294128],
-                [0.43575474],
-                [0.42439987],
-                [0.48393954],
+                [0.18014312],
+                [0.17679804],
+                [0.34352722],
+                [0.33621384],
+                [0.40694277],
             ]
         ),
     }

diff --git a/tests/models/test_fracture_damage.py b/tests/models/test_fracture_damage.py
@@ -106,6 +106,8 @@ class AnisotropicDamage3D(
 def test_anisotropic_damage_values():
     params_local = copy.deepcopy(damage.model_params)
     solid_params_local = {**damage.solid_params, **additional_solid_params["both"]}
+    # Set higher friction coefficient to improve convergence.
+    solid_params_local["friction_coefficient"] = 0.5
     params_local["material_constants"] = {
         "solid": FractureDamageSolidConstants(**solid_params_local)
     }