From cdeb7bf9a511e1751528fb5b514b8dd55bcdf59f Mon Sep 17 00:00:00 2001
From: podlesny <podlesny@zedat.fu-berlin.de>
Date: Sat, 6 Mar 2021 16:55:59 +0100
Subject: [PATCH] .
---
data/tools/config.ini | 2 +-
data/tools/main.py | 6 +++---
src/foam/foam.cfg | 2 +-
3 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/data/tools/config.ini b/data/tools/config.ini
index 40b64b13..8f7921b3 100644
--- a/data/tools/config.ini
+++ b/data/tools/config.ini
@@ -3,6 +3,6 @@
#/home/joscha/Desktop/strikeslip/viscosity-1e3/
[directories]
-simulation = /home/joscha/software/dune/build-release/dune-tectonic/src/foam/output/test/
+simulation = /home/joscha/software/dune/build-release/dune-tectonic/src/foam/output/tresca/
experiment = ~/group/publications/2016-RosenauCorbiDominguezRudolfRitterPipping
output = generated
diff --git a/data/tools/main.py b/data/tools/main.py
index 4502ed00..0cf53126 100644
--- a/data/tools/main.py
+++ b/data/tools/main.py
@@ -34,7 +34,7 @@ def build_patch(coords, percentage):
NBODIES = 2
FINAL_TIME = 15 # s
-FINAL_VELOCITY = 1e-5 # m/s
+FINAL_VELOCITY = 2e-4 # m/s
THRESHOLD_VELOCITY = 0.5*FINAL_VELOCITY # 1000e-6 + FINAL_VELOCITY
TANGENTIAL_COORDS = 1
@@ -72,9 +72,9 @@ for body_ID in range(NBODIES):
continue
coords = np.array(h5file[body + '/coordinates'])
- patch = build_patch(coords, 0.1)
+ patch = build_patch(coords, 0.05)
- friction_stats(h5file, body_ID, FINAL_TIME, patch, [])
+ friction_stats(h5file, body_ID, FINAL_TIME, [], [0, 3])
plt.show()
diff --git a/src/foam/foam.cfg b/src/foam/foam.cfg
index fb03fce1..9f1fc8e0 100644
--- a/src/foam/foam.cfg
+++ b/src/foam/foam.cfg
@@ -34,7 +34,7 @@ bulkViscosity = 0.0 # [Pas]
[boundary.friction]
C = 0.0 # [Pa]
-mu0 = 0.6 # [ ]
+mu0 = 0.3 # [ ]
V0 = 1e-6 # [m/s]
L = 1e-5 # [m]
initialAlpha = 0.0 #-10.0 # [ ]
--
GitLab