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 03c9864eb8e1eb7f7fba6090654205d4d343271d
parent 23626aeafc710f8cb2e87cfa89d65ea0dc823f27
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Fri, 12 Apr 2019 09:26:46 +0200

Add comparison to Tulaczyk 2006 figure

Diffstat:
A1d_fd_simple_shear_rheology_tulaczyk.gp | 28++++++++++++++++++++++++++++
A1d_fd_simple_shear_rheology_tulaczyk.png | 0
MMakefile | 16+++++++++++++++-
MREADME.md | 10++++++++++
Atulaczyk2006-fig1.png | 0
5 files changed, 53 insertions(+), 1 deletion(-)

diff --git a/1d_fd_simple_shear_rheology_tulaczyk.gp b/1d_fd_simple_shear_rheology_tulaczyk.gp @@ -0,0 +1,28 @@ +#!/usr/bin/env gnuplot + +## resemble Tulaczyk 2006 plot: + +set terminal pngcairo color size 18.6 cm, 17.6 cm +set output "1d_fd_simple_shear_rheology_tulaczyk.png" +set xlabel "Normalized shear stress [-]" +set ylabel "Normalized shear strain rate [-]" offset 2 + +set xrange [0.1:10.0] +set yrange [1.0:1000.0] +set key top left #samplen 0.9 + +set logscale xy + +set style line 1 linetype 1 linewidth 3 pointtype 1 pointsize 1 + +# divide friction by threshold friction value (mu_s) +# divide shear velocity by 1 m/h +plot "1d_fd_simple_shear_rheology_b0.01_tulaczyk.txt" u ($1/0.502):($2/3.5e-7) w l lw 2 t "b = 0.01", \ + "1d_fd_simple_shear_rheology_b0.10_tulaczyk.txt" u ($1/0.502):($2/3.5e-8) w l lw 2 t "b = 0.1", \ + "1d_fd_simple_shear_rheology_b0.20_tulaczyk.txt" u ($1/0.502):($2/1.7e-8) w l lw 2 t "b = 0.2", \ + "1d_fd_simple_shear_rheology_b0.40_tulaczyk.txt" u ($1/0.502):($2/8.6e-9) w l lw 2 t "b = 0.4", \ + "1d_fd_simple_shear_rheology_b0.60_tulaczyk.txt" u ($1/0.502):($2/5.7e-9) w l lw 2 t "b = 0.6", \ + "1d_fd_simple_shear_rheology_b0.80_tulaczyk.txt" u ($1/0.502):($2/4.3e-9) w l lw 2 t "b = 0.8", \ + "1d_fd_simple_shear_rheology_b1.00_tulaczyk.txt" u ($1/0.502):($2/3.4e-9) w l lw 2 t "b = 1.0" + +set xtics norotate # Restore defaults diff --git a/1d_fd_simple_shear_rheology_tulaczyk.png b/1d_fd_simple_shear_rheology_tulaczyk.png Binary files differ. diff --git a/Makefile b/Makefile @@ -7,7 +7,8 @@ HDR=$(wildcard *.h) default: 1d_fd_simple_shear.png \ 1d_fd_simple_shear_rheology.png \ 1d_fd_simple_shear_rheology_kamb.png \ - 1d_fd_simple_shear_rheology_iverson.png + 1d_fd_simple_shear_rheology_iverson.png \ + 1d_fd_simple_shear_rheology_tulaczyk.png 1d_fd_simple_shear: $(OBJ) $(HDR) $(CC) $(LDFLAGS) $(OBJ) -o $@ @@ -60,6 +61,19 @@ default: 1d_fd_simple_shear.png \ done; done' gnuplot $<_rheology_iverson.gp > $@ +# shear velocity rate from 0.1 m/h to 100 m/h +1d_fd_simple_shear_rheology_tulaczyk.png: 1d_fd_simple_shear 1d_fd_simple_shear_rheology_tulaczyk.gp + /bin/bash -c '\ + for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 1.00); do \ + out="$<_rheology_b$${b}_tulaczyk.txt"; \ + rm -f "$$out"; \ + for t in $$(seq 0.1 0.002 0.9); do \ + printf "$$t\t" >> "$$out"; \ + ./$< -f 0.5 -P 10e3 --stress-ratio $$t -b $$b | \ + tail -n 1 | cut -f2 >> "$$out"; \ + done; done' + gnuplot $<_rheology_tulaczyk.gp > $@ + .PHONY: watch watch: echo $(SRC) $(HDR) | tr ' ' '\n' | entr -s 'make && ./1d_fd_simple_shear' diff --git a/README.md b/README.md @@ -56,3 +56,13 @@ Glass beads have *b* = 0.94. #### Continuum model ![1d_fd_simple_shear_rheology_iverson.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/1d_fd_simple_shear_rheology_iverson.png) + +--- + +#### Tulaczyk 2006, Whillans Ice Plain + +![tulaczyk2006-fig1.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/tulaczyk2006-fig1.png) + +#### Continuum model + +![1d_fd_simple_shear_rheology_tulaczyk.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/1d_fd_simple_shear_rheology_tulaczyk.png) diff --git a/tulaczyk2006-fig1.png b/tulaczyk2006-fig1.png Binary files differ.