1d_fd_simple_shear

Continuum model for granular flows with pore-pressure dynamics
git clone git://src.adamsgaard.dk/1d_fd_simple_shear
Log | Files | Refs | README | LICENSE

commit a408ece29c21ca2b524f4f1335b95354038effbc
parent 9b77bca8b4d73c3a5c3a4ce6997229ead220a699
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 11 Apr 2019 16:01:15 +0200

Add Kamb1991 figure and comparison

Diffstat:
M1d_fd_simple_shear_rheology.gp | 23+++++++++++++++++++++++
MREADME.md | 9+++++++++
Akamb1991-fig1.png | 0
3 files changed, 32 insertions(+), 0 deletions(-)

diff --git a/1d_fd_simple_shear_rheology.gp b/1d_fd_simple_shear_rheology.gp @@ -20,3 +20,26 @@ plot "1d_fd_simple_shear_rheology_b0.01.txt" u 1:2 w lp lw 2 t "b = 0.01", \ "1d_fd_simple_shear_rheology_b1.00.txt" u 1:2 w lp lw 2 t "b = 1.0" set xtics norotate # Restore defaults + +## + +set terminal pngcairo color size 18.6 cm, 17.6 cm +set output "1d_fd_simple_shear_rheology_md.png" +set xlabel "Applied friction, mu = tau/N [-]" +set ylabel "Shear velocity, v [m/d]" offset 2 + +set xrange [0.0:0.7] +set yrange [0.0:2.5] +set key top left #samplen 0.9 + +set style line 1 linetype 1 linewidth 3 pointtype 1 pointsize 1 + +plot "1d_fd_simple_shear_rheology_b0.01.txt" u 1:($2*3600*24) w l lw 2 t "b = 0.01", \ + "1d_fd_simple_shear_rheology_b0.10.txt" u 1:($2*3600*24) w l lw 2 t "b = 0.1", \ + "1d_fd_simple_shear_rheology_b0.20.txt" u 1:($2*3600*24) w l lw 2 t "b = 0.2", \ + "1d_fd_simple_shear_rheology_b0.40.txt" u 1:($2*3600*24) w l lw 2 t "b = 0.4", \ + "1d_fd_simple_shear_rheology_b0.60.txt" u 1:($2*3600*24) w l lw 2 t "b = 0.6", \ + "1d_fd_simple_shear_rheology_b0.80.txt" u 1:($2*3600*24) w l lw 2 t "b = 0.8", \ + "1d_fd_simple_shear_rheology_b1.00.txt" u 1:($2*3600*24) w l lw 2 t "b = 1.0" + +set xtics norotate # Restore defaults diff --git a/README.md b/README.md @@ -35,3 +35,12 @@ Alternatively, an implementation in [Julia](https://julialang.org) resides in Glass beads have *b* = 0.94. ![1d_fd_simple_shear_rheology.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/1d_fd_simple_shear_rheology.png) + +Vertical axis rescaled to Kamb 1991: + +![1d_fd_simple_shear_rheology_md.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/1d_fd_simple_shear_rheology_md.png) + +**Kamb 1991, subglacial team from Upstream-B site** + +Vertical scale as in Kamb 1991: +![kamb1991-fig1.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/kamb1991-fig1.png) diff --git a/kamb1991-fig1.png b/kamb1991-fig1.png Binary files differ.