commit 9c642d72b8e7cbc2355cc6ae865500327be0d5f1
parent f054a123409ad313c173719cad47e6e5a60b2796
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Tue, 2 Jul 2019 10:41:32 +0200
Add skindepth to stick-slip-depth figures
Diffstat:
6 files changed, 33 insertions(+), 10 deletions(-)
diff --git a/fig2/fig.gp b/fig2/fig.gp
@@ -7,6 +7,7 @@ set terminal pdfcairo color size 7.5 cm, 8.0 cm
# set yrange [0.0:0.01]
# set key top left #samplen 0.9
+
set multiplot layout 2,1 \
margins 0.19,0.95,0.23,0.95 \
spacing 0.03,0.03
diff --git a/fig3/Makefile b/fig3/Makefile
@@ -3,7 +3,7 @@ FIG = fig3
default: ../$(FIG).pdf
-unnamed.output00000.txt: $(BIN) Makefile
+unnamed.output00000.txt: $(BIN) Makefile ../skindepth.gp
/bin/sh -c '\
./$< --resolution 200 --length 8.0 --normal-stress 150e3 \
--fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 \
diff --git a/fig3/fig.gp b/fig3/fig.gp
@@ -20,24 +20,32 @@ filename(t) = sprintf('unnamed.output%05d.txt', t)
set pm3d map
set palette rgb 10,6,6
+load '../skindepth.gp'
+k = 2e-17
+f = 1.0/(3600.0*24.0)
+f(x) = 8.0 - skindepth(f,k)
+
unset colorbox
set xlabel "Water pressure [kPa]"
set ylabel "Vertical position, z [m]"
set xrange [0:100]
-plot for [t = t0:tend] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t "", \
+ f(x) title "" w l
set xlabel "Eff. normal stress [kPa]"
set ylabel ""
set xrange [0:200]
unset ytics
-plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t "", \
+ f(x) title "" w l
set xlabel "Shear velocity [m/s]"
set ylabel ""
set xrange [0.0:0.06]
set xtics (0,0.02,0.04,0.06)
unset ytics
-plot for [t = t0:tend] filename(t) u 2:1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend] filename(t) u 2:1:(t-t0) w l lc palette lw 1 t "", \
+ f(x) title "" w l
set xlabel "Shear strain rate [1/s]"
set ylabel ""
@@ -47,7 +55,8 @@ set xtics (0,0.001,0.002,0.003)
unset ytics
prev1 = prev2 = 0
shift(x) = (prev2 = prev1, prev1 = x)
-plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1:(t-t0) w l lc palette lw 1 t "", \
+ f(x) title "" w l
unset format x
unset multiplot
diff --git a/fig4/Makefile b/fig4/Makefile
@@ -3,7 +3,7 @@ FIG = fig4
default: ../$(FIG).pdf
-unnamed.output00000.txt: $(BIN) Makefile
+unnamed.output00000.txt: $(BIN) Makefile ../skindepth.gp
/bin/sh -c '\
./$< --resolution 200 --length 8.0 --normal-stress 150e3 \
--fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 \
diff --git a/fig4/fig.gp b/fig4/fig.gp
@@ -20,24 +20,32 @@ filename(t) = sprintf('unnamed.output%05d.txt', t)
set pm3d map
set palette rgb 10,6,6
+load '../skindepth.gp'
+k = 2e-17
+f = 1.0/(3600.0*24.0)
+f(x) = 8.0 - skindepth(f,k)
+
unset colorbox
set xlabel "Water pressure [kPa]"
set ylabel "Vertical position, z [m]"
set xrange [0:100]
-plot for [t = t0:tend] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t "", \
+ f(x) title "" w l
set xlabel "Eff. normal stress [kPa]"
set ylabel ""
set xrange [0:200]
unset ytics
-plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t "", \
+ f(x) title "" w l
set xlabel "Shear velocity [-]"
set ylabel ""
set xrange [0.0:1.0]
#set xtics (0,0.01,0.02)
unset ytics
-plot for [t = t0:tend] filename(t) u 2:1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend] filename(t) u 2:1:(t-t0) w l lc palette lw 1 t "", \
+ f(x) title "" w l
set xlabel "Shear strain rate [1/s]"
set ylabel ""
@@ -48,7 +56,8 @@ set xtics (0,0.05,0.10,0.15)
unset ytics
prev1 = prev2 = 0
shift(x) = (prev2 = prev1, prev1 = x)
-plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1:(t-t0) w l lc palette lw 1 t "", \
+ f(x) title "" w l
unset format x
unset multiplot
diff --git a/skindepth.gp b/skindepth.gp
@@ -0,0 +1,4 @@
+phi = 0.25
+beta_f = 4.5e-10
+mu_f = 0.001
+skindepth(f,k) = sqrt(k/(phi*mu_f*beta_f*3.141592654*f)) # f: 1/s
