tUpdate readme with usage and diurnal animation. Reverse subplot order. - cngf-… | |
git clone git://src.adamsgaard.dk/cngf-pf | |
Log | |
Files | |
Refs | |
README | |
LICENSE | |
--- | |
commit 4dc0432607d8119b7a3b27c9911c65b771116458 | |
parent d87c0c129c509e0d5fa76b99bd1a0a31a7e66475 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Mon, 29 Apr 2019 16:58:52 +0200 | |
Update readme with usage and diurnal animation. Reverse subplot order. | |
Diffstat: | |
M 1d_fd_simple_shear_fluid.gp | 19 ++++++++++--------- | |
M Makefile | 3 ++- | |
M README.md | 46 +++++++++++++++++++++++++++++… | |
A diurnal.gif | 0 | |
4 files changed, 57 insertions(+), 11 deletions(-) | |
--- | |
diff --git a/1d_fd_simple_shear_fluid.gp b/1d_fd_simple_shear_fluid.gp | |
t@@ -7,28 +7,29 @@ set terminal pngcairo color size 60 cm, 17.6 cm | |
set multiplot layout 1,3 | |
# set multiplot layout 1,4 | |
-set ylabel "Vertical position, z [m]" offset 2 | |
set yrange [0.0:2.0] | |
set key bottom right #samplen 0.9 | |
-set xlabel "Horizontal velocity, v_x [m/s]" | |
-set xrange [0.0:0.025] | |
-plot filename u 2:1 w l lw 2 lc "red" t "" | |
+set xlabel "Water pressure, p_f [Pa]" | |
+set ylabel "Vertical position, z [m]" offset 2 | |
+set xrange [p_min:p_max] | |
+plot filename u 4:1 w l lw 2 lc "blue" t "" | |
set xlabel "Effective normal stress [Pa]" | |
set ylabel "" | |
set xrange [0:200e3] | |
plot filename u 3:1 w l lw 2 lc "black" t "" | |
-set xlabel "Water pressure, p_f [Pa]" | |
-set ylabel "" | |
-set xrange [p_min:p_max] | |
-plot filename u 4:1 w l lw 2 lc "blue" t "" | |
- | |
# set xlabel "Friction, mu [-]" | |
# set xrange [-0.05:1.5] | |
# plot filename u 5:1 w l lw 2 lc "gray" t "" | |
+set xlabel "Horizontal velocity, v_x [m/s]" | |
+set ylabel "" | |
+set xrange [0.0:0.025] | |
+plot filename u 2:1 w l lw 2 lc "red" t "" | |
+ | |
+ | |
unset multiplot | |
diff --git a/Makefile b/Makefile | |
t@@ -24,7 +24,8 @@ diurnal.mp4: 1d_fd_simple_shear 1d_fd_simple_shear_fluid.gp | |
ffmpeg -i diurnal.output%05d.txt.png -y diurnal.mp4 | |
diurnal.gif: diurnal.mp4 | |
- gm convert -delay 5 -loop 0 diurnal.output*.txt.png diurnal.gif | |
+ convert diurnal.output*.txt.png \ | |
+ -delay 1 -loop 0 -fuzz 10% -layers Optimize $@ | |
1d_fd_simple_shear: $(OBJ) $(HDR) | |
$(CC) $(LDFLAGS) $(OBJ) -o $@ | |
diff --git a/README.md b/README.md | |
t@@ -12,6 +12,44 @@ directory and an output PNG figure is generated. | |
Alternatively, an implementation in [Julia](https://julialang.org) resides in | |
**julia/**. | |
+### Advanced usage | |
+The majority of simulation parameters can be adjusted from the command line: | |
+ | |
+``` | |
+usage: 1d_fd_simple_shear [OPTIONS] [NAME] | |
+runs a simulation and outputs state in files prefixed with NAME. | |
+optional arguments: | |
+ -N, --normalize normalize output velocity | |
+ -G, --gravity VAL gravity magnitude [m/s^2] | |
+ -P, --normal-stress VAL normal stress on top [Pa] | |
+ -m, --stress-ratio VAL applied stress ratio [-] | |
+ -V, --velocity-bottom VAL base velocity at bottom [m/s] | |
+ -A, --nonlocal-amplitude VAL amplitude of nonlocality [-] | |
+ -b, --rate-dependence VAL rate dependence beyond yield [-] | |
+ -f, --friction-coefficient VAL grain friction coefficient [-] | |
+ -p, --porosity VAL porosity fraction [-] | |
+ -d, --grain-size VAL representative grain size [m] | |
+ -r, --density VAL grain material density [kg/m^3] | |
+ -n, --resolution VAL number of cells in domain [-] | |
+ -o, --origo VAL coordinate system origo [m] | |
+ -L, --length VAL domain length [m] | |
+ -F, --fluid enable pore fluid computations | |
+ -c, --fluid-compressibility VAL fluid compressibility [Pa^-1] | |
+ -i, --fluid-viscosity VAL fluid viscosity [Pa*s] | |
+ -R, --fluid-density VAL fluid density [kg/m^3] | |
+ -k, --fluid-permeability VAL fluid intrinsic permeability [m^2] | |
+ -O, --fluid-pressure-top VAL fluid pressure at +z edge [Pa] | |
+ -a, --fluid-pressure-ampl VAL amplitude of pressure variations [Pa] | |
+ -q, --fluid-pressure-freq VAL frequency of pressure variations [s^-1] | |
+ -H, --fluid-pressure-phase VAL fluid pressure at +z edge [Pa] | |
+ -t, --time VAL simulation start time [s] | |
+ -T, --time-end VAL simulation end time [s] | |
+ -D, --time-step VAL computational time step length [s] | |
+ -I, --file-interval VAL interval between output files [s] | |
+ -v, --version show version information | |
+ -h, --help show this message | |
+``` | |
+ | |
## Results | |
### Strain distribution | |
t@@ -26,7 +64,6 @@ probably be improved further. | |
| ![damsgaard2013-fig8.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/m… | |
### Stress and strain rate | |
- | |
The rheology is of Bingham type, where no deformation occurs beneath the | |
Mohr-Coulomb yield limit. Above it, deformation is highly non-linearly viscous… | |
The model has a parameter *b* for rate dependence beyond yield. Glass beads | |
t@@ -41,3 +78,10 @@ have *b* = 0.94 ([Forterre and Pouliquen | |
| ![iverson2010-fig2a.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/ma… | |
| Whillans Ice Plain ([Tulaczyk 2006](https://doi.org/10.3189/1727565067818286… | |
| ![tulaczyk2006-fig1.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/ma… | |
+ | |
+### Variable water pressure | |
+The model is expanded from the Henann and Kamrin 2013 model by including a | |
+diffusive porewater pressure parameterization. Below is an example of diurnal | |
+water-pressure variations that gradually propagate into the bed. | |
+ | |
+![diurnal.gif](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/diurnal… | |
diff --git a/diurnal.gif b/diurnal.gif | |
Binary files differ. |