[Problem] New parameters

src/sand-wedge-data/mygeometry.hh

@@ -85,7 +85,7 @@ namespace {
 
 double const lengthScale = reference::s;
 
-double const depth = 0.10;
+double const depth = 0.60 * lengthScale;
 
 LocalVector const A = rotate(reference::A);
 LocalVector const B = rotate(reference::B);

src/sand-wedge-data/parset.cfg

@@ -9,32 +9,32 @@ writeVTK        = false
 finalTime       = 1800  # [s]
 
 [body]
-bulkModulus     = 1e5   # [Pa]
-poissonRatio    = 0.3   # [1]      0.2 - 0.3
+bulkModulus     = 0.5e5 # [Pa]
+poissonRatio    = 0.3   # [1]
 [body.elastic]
 density         = 900   # [kg/m^3]
-shearViscosity  = 1e3   # [Pas]    0
-bulkViscosity   = 1e3   # [Pas]    0
+shearViscosity  = 1e3   # [Pas]
+bulkViscosity   = 1e3   # [Pas]
 [body.viscoelastic]
 density         = 1000  # [kg/m^3]
 shearViscosity  = 1e4   # [Pas]
 bulkViscosity   = 1e4   # [Pas]
 
 [boundary.friction]
-smallestDiameter= 1e-3  # [m]
+smallestDiameter= 2e-3  # [m]
 C               = 10    # [Pa]
-mu0             = 0.7   # [1]
+mu0             = 0.7   # [ ]
 V0              = 5e-5  # [m/s]
-L               = 2e-5  # [m]      ?
-initialAlpha    = 0.916290731874155 # [ ]      ?
+L               = 2.5e-5# [m]
+initialAlpha    = 0     # [ ]
 stateModel      = AgeingLaw
 frictionModel   = Truncated
 [boundary.friction.weakening]
-a               = 0.015 # [1]      ?
-b               = 0.030 # [1]      ?
+a               = 0.002 # [ ]
+b               = 0.014 # [ ]
 [boundary.friction.strengthening]
-a               = 0.030 # [1]      ?
-b               = 0.015 # [1]      ?
+a               = 0.025 # [ ]
+b               = 0.005 # [ ]
 
 [timeSteps]
 refinementTolerance = 1e-5
@@ -42,22 +42,22 @@ number = 100000
 scheme = newmark
 
 [u0.solver]
-tolerance         = 1e-10
+tolerance         = 1e-8
 maximumIterations = 100000
 verbosity         = quiet
 
 [a0.solver]
-tolerance         = 1e-10
+tolerance         = 1e-8
 maximumIterations = 100000
 verbosity         = quiet
 
 [v.solver]
-tolerance         = 1e-10
+tolerance         = 1e-8
 maximumIterations = 100000
 verbosity         = quiet
 
 [v.fpi]
-tolerance         = 1e-10
+tolerance         = 1e-8
 maximumIterations = 10000
 lambda            = 0.5