manus_continuum_granular1_exp

Experiments for first paper with continuum granular model
git clone git://src.adamsgaard.dk/.manus_continuum_granular1_exp
Log | Files | Refs | Submodules | README | LICENSE

commit 5b1ea93b8d8fc4dec57c856931ee0c2eca03cc2d
parent 814e735794e22caa86e886221534455fa5fee840
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 10 Oct 2019 18:53:24 +0200

Perturb from a normal stress of 200 kPa

Diffstat:
Mfig-pulse_square/Makefile | 4++--
Mfig-pulse_square/fig.gp | 8++++----
Mfig-pulse_triangle/Makefile | 4++--
Mfig-pulse_triangle/fig.gp | 8++++----
4 files changed, 12 insertions(+), 12 deletions(-)

diff --git a/fig-pulse_square/Makefile b/fig-pulse_square/Makefile @@ -10,12 +10,12 @@ PF_SEQ=1 2 3 4 5 6 7 8 9 sim-4.output00000.txt: $(BIN) Makefile /bin/sh -c 'for i in $(PF_SEQ); do \ ./$< --length 8.0 \ - --normal-stress 100e3 \ + --normal-stress 200e3 \ --set-shear-velocity $$( echo "10000/(365*24*60*60)" | bc -l ) \ --fluid \ --fluid-permeability 2e-17 \ --fluid-pressure-top 0.0 \ - --fluid-pressure-ampl $$( printf "%s0e3" "$$i" ) \ + --fluid-pressure-ampl "$$( echo "2*$$i" | bc -l )0e3" \ --fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \ --fluid-pressure-pulse-time $$( echo "$$i*3600*24" | bc -l ) \ --fluid-pressure-pulse-shape square \ diff --git a/fig-pulse_square/fig.gp b/fig-pulse_square/fig.gp @@ -10,7 +10,7 @@ set multiplot layout 3,1 \ set key top right samplen 1.0 unset xlabel unset xtics -set yrange [0:320] +set yrange [0:380] set ylabel "[kPa]" set ytics 100 plot \ @@ -26,8 +26,8 @@ unset ytics #set y2tics 0.2 #set y2label "Effective friction [-]" #plot "" u ($0/1008*7):($4) w l lw 2 lt 7 t "" -set yrange [0:80] -set y2range [0:80] +set yrange [0:100] +set y2range [0:100] set y2tics 20 set y2label "Shear stress [kPa]" plot "" u ($0/1008*7):($4*$2/1000) w l lw 2 lt 1 t "" @@ -36,7 +36,7 @@ set xtics unset y2label unset y2tics set ytics 0.5 -set yrange [2.5:0.0] +set yrange [3.0:0.0] set xlabel "Time [d]" set ylabel "Depth of max. strain rate [m]" plot "max_strainrate_depth.txt" u ($0/1008*7):(8.0 - $1) w l lw 2 lt 4 t "" diff --git a/fig-pulse_triangle/Makefile b/fig-pulse_triangle/Makefile @@ -10,12 +10,12 @@ PF_SEQ=1 2 3 4 5 6 7 8 9 sim-4.output00000.txt: $(BIN) Makefile /bin/sh -c 'for i in $(PF_SEQ); do \ ./$< --length 8.0 \ - --normal-stress 100e3 \ + --normal-stress 200e3 \ --set-shear-velocity $$( echo "10000/(365*24*60*60)" | bc -l ) \ --fluid \ --fluid-permeability 2e-17 \ --fluid-pressure-top 0.0 \ - --fluid-pressure-ampl $$( printf "%s0e3" "$$i" ) \ + --fluid-pressure-ampl "$$( echo "2*$$i" | bc -l )0e3" \ --fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \ --fluid-pressure-pulse-time $$( echo "$$i*3600*24" | bc -l ) \ --fluid-pressure-pulse-shape triangle \ diff --git a/fig-pulse_triangle/fig.gp b/fig-pulse_triangle/fig.gp @@ -10,7 +10,7 @@ set multiplot layout 3,1 \ set key top right samplen 1.0 unset xlabel unset xtics -set yrange [0:320] +set yrange [0:380] set ylabel "[kPa]" set ytics 100 plot \ @@ -26,8 +26,8 @@ unset ytics #set y2tics 0.2 #set y2label "Effective friction [-]" #plot "" u ($0/1008*7):($4) w l lw 2 lt 7 t "" -set yrange [0:80] -set y2range [0:80] +set yrange [0:100] +set y2range [0:100] set y2tics 20 set y2label "Shear stress [kPa]" plot "" u ($0/1008*7):($4*$2/1000) w l lw 2 lt 1 t "" @@ -36,7 +36,7 @@ set xtics unset y2label unset y2tics set ytics 0.5 -set yrange [2.5:0.0] +set yrange [3.0:0.0] set xlabel "Time [d]" set ylabel "Depth of max. strain rate [m]" plot "max_strainrate_depth.txt" u ($0/1008*7):(8.0 - $1) w l lw 2 lt 4 t ""