1d_fd_simple_shear_transient

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

commit 8a887a0c13f67c89ab1f8ed8e4214d722b464107
parent fac32c7a7febb4dcc3aea5968ee0255046f60ccc
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Fri,  8 Nov 2019 11:17:28 +0100

Update README for fork

Diffstat:
MREADME.md | 48++++++------------------------------------------
1 file changed, 6 insertions(+), 42 deletions(-)

diff --git a/README.md b/README.md @@ -1,20 +1,20 @@ # 1d_fd_simple_shear [![pipeline status](https://gitlab.com/admesg/1d_fd_simple_shear/badges/master/pipeline.svg)](https://gitlab.com/admesg/1d_fd_simple_shear/commits/master) -This project contains a 1d implementation of the [Henann and Kamrin -2013](https://doi.org/10.1073/pnas.1219153110) model under simple shear, with -various extensions such as diffusive fluid coupling and cohesion. +This project contains a 1d implementation of a continuum +model for granular flows and pore-pressure dynamics with +transient porosity and strength. It is an extension of +[1d_fd_simple_shear](https://src.adamsgaard.dk/1d_fd_simple_shear). ## How to run The implementation in C requires a C99-compatible compiler (e.g. `gcc` or -`clang`). Visualization requires `gnuplot`. To run, simply run `make` in this -directory and an output PNG figure is generated. +`clang`). Run `make test` for integration tests. ### Advanced usage The majority of simulation parameters can be adjusted from the command line: ``` -usage: 1d_fd_simple_shear [OPTIONS] [NAME] +usage: 1d_fd_simple_shear_transient [OPTIONS] [NAME] runs a simulation and outputs state in files prefixed with NAME. If NAME is not specified, intermediate output are not written. @@ -57,39 +57,3 @@ Optional arguments only relevant with transient (fluid) simulation: -I, --file-interval VAL interval between output files [s] (default 0.1) ``` - -## Results - -### Strain distribution -See below for a comparison with the discrete-element method results from -[Damsgaard et al. 2013](https://doi.org/10.1002/2013JF002830). The continuum -model is not tuned rigorously against the DEM, and the correspondence could -probably be improved further. - -| Discrete-element model | Continuum Model | -| ----------------------- | --------------- | -| Damsgaard et al. 2013 | | -| ![damsgaard2013-fig8.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/doc/damsgaard2013-fig8.png) | ![1d_fd_simple_shear.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/examples/1d_fd_simple_shear.png) | - -### 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 -have *b* = 0.94 ([Forterre and Pouliquen -2003](https://doi.org/10.1017/S0022112003004555)). - -| Real material (laboratory or field study) | Continuum Model | -| -------------------------------------------- | --------------- | -| Upstream-B ([Kamb 1991](https://doi.org/10.1029/91jb00946)): | | -| ![kamb1991-fig1.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/doc/kamb1991-fig1.png) | ![1d_fd_simple_shear_rheology_kamb.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/examples/1d_fd_simple_shear_rheology_kamb.png) | -| Various subglacial tills ([Iverson 2010](https://doi.org/10.3189/002214311796406220)): | | -| ![iverson2010-fig2a.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/doc/iverson2010-fig2a.png) | ![1d_fd_simple_shear_rheology_iverson.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/examples/1d_fd_simple_shear_rheology_iverson.png) | -| Whillans Ice Plain ([Tulaczyk 2006](https://doi.org/10.3189/172756506781828601)): | | -| ![tulaczyk2006-fig1.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/doc/tulaczyk2006-fig1.png) | ![1d_fd_simple_shear_rheology_tulaczyk.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/examples/1d_fd_simple_shear_rheology_tulaczyk.png) | - -### 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/examples/diurnal.gif)