diff --git a/data/tools/config.ini b/data/tools/config.ini index 40b64b13c922d79bb6319d5c1907f75384e49378..8f7921b398f15dfd65b3199f64f5ceebfd445658 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 4502ed0014d43706de34bb82430699e0ccdf02ef..0cf531266f58f52893ee935590f9e11dfad7419c 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 fb03fce1796be8915c109756f0c236bcbc49c45c..9f1fc8e05f81c9b35305bb3886816ac83efb5241 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 # [ ]