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 9b77bca8b4d73c3a5c3a4ce6997229ead220a699
parent e844d141252ab91ed58152d7174be6997aaa2f84
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 11 Apr 2019 14:05:07 +0200

Change b values for rheology plot

Diffstat:
M1d_fd_simple_shear_rheology.gp | 8+++-----
M1d_fd_simple_shear_rheology.png | 0
MMakefile | 2+-
MREADME.md | 2++
4 files changed, 6 insertions(+), 6 deletions(-)

diff --git a/1d_fd_simple_shear_rheology.gp b/1d_fd_simple_shear_rheology.gp @@ -11,14 +11,12 @@ 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.20.txt" u 1:2 w lp lw 2 t "b = 0.2", \ - "1d_fd_simple_shear_rheology_b0.30.txt" u 1:2 w lp lw 2 t "b = 0.3", \ +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_b0.10.txt" u 1:2 w lp lw 2 t "b = 0.1", \ + "1d_fd_simple_shear_rheology_b0.20.txt" u 1:2 w lp lw 2 t "b = 0.2", \ "1d_fd_simple_shear_rheology_b0.40.txt" u 1:2 w lp lw 2 t "b = 0.4", \ - "1d_fd_simple_shear_rheology_b0.50.txt" u 1:2 w lp lw 2 t "b = 0.5", \ "1d_fd_simple_shear_rheology_b0.60.txt" u 1:2 w lp lw 2 t "b = 0.6", \ - "1d_fd_simple_shear_rheology_b0.70.txt" u 1:2 w lp lw 2 t "b = 0.7", \ "1d_fd_simple_shear_rheology_b0.80.txt" u 1:2 w lp lw 2 t "b = 0.8", \ - "1d_fd_simple_shear_rheology_b0.90.txt" u 1:2 w lp lw 2 t "b = 0.9", \ "1d_fd_simple_shear_rheology_b1.00.txt" u 1:2 w lp lw 2 t "b = 1.0" set xtics norotate # Restore defaults diff --git a/1d_fd_simple_shear_rheology.png b/1d_fd_simple_shear_rheology.png Binary files differ. diff --git a/Makefile b/Makefile @@ -21,7 +21,7 @@ default: 1d_fd_simple_shear.png \ 1d_fd_simple_shear_rheology.png: 1d_fd_simple_shear 1d_fd_simple_shear_rheology.gp /bin/bash -c '\ - for b in $$(seq 0.20 0.10 1.00); do \ + for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 1.00); do \ out="$<_rheology_b$$b.txt"; \ rm "$$out"; \ for t in $$(seq 0.01 0.01 0.8); do \ diff --git a/README.md b/README.md @@ -32,4 +32,6 @@ Alternatively, an implementation in [Julia](https://julialang.org) resides in **Continuum model: Stress vs. strain rate** +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)