From 5c0f35f2d6a39ed1b57209c0b1251bb093150bec Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sat, 19 Jul 2014 00:19:31 +0200
Subject: [PATCH] [Problem] New parameters

---
 src/sand-wedge-data/mygeometry.hh |  2 +-
 src/sand-wedge-data/parset.cfg    | 32 +++++++++++++++----------------
 2 files changed, 17 insertions(+), 17 deletions(-)

diff --git a/src/sand-wedge-data/mygeometry.hh b/src/sand-wedge-data/mygeometry.hh
index 19e05533..fdfcdea4 100644
--- a/src/sand-wedge-data/mygeometry.hh
+++ b/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);
diff --git a/src/sand-wedge-data/parset.cfg b/src/sand-wedge-data/parset.cfg
index ab6670d0..45a42435 100644
--- a/src/sand-wedge-data/parset.cfg
+++ 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
 
-- 
GitLab