README.md (3529B)
1 # 1d_fd_simple_shear 2 [](https://gitlab.com/admesg/1d_fd_simple_shear/commits/master) 3 4 This project contains a 1d implementation of a continuum 5 model for granular flows and pore-pressure dynamics with 6 transient porosity and strength. It is an extension of 7 [1d_fd_simple_shear](https://src.adamsgaard.dk/1d_fd_simple_shear). 8 9 WARNING: The code is currently under development, and the results may 10 in no way be accurate or reliable. 11 12 ## How to run 13 The implementation in C requires a C99-compatible compiler (e.g. `gcc` or 14 `clang`). Run `make test` for integration tests. 15 16 ### Advanced usage 17 The majority of simulation parameters can be adjusted from the command line: 18 19 ``` 20 usage: 1d_fd_simple_shear_transient [OPTIONS] [NAME] 21 runs a simulation and outputs state in files prefixed with NAME. 22 If NAME is not specified, intermediate output are not written. 23 24 Optional arguments: 25 -v, --version show version information 26 -h, --help show this message 27 -N, --normalize normalize output velocity 28 -G, --gravity VAL gravity magnitude [m/s^2] (default 9.81) 29 -P, --normal-stress VAL normal stress on top [Pa] (default 120000) 30 -m, --stress-ratio VAL applied stress ratio [-] (default 0.45) 31 -s, --set-shear-velocity VAL set top shear velocity to this value [m/s] 32 (default nan), overrides --stress-ratio 33 -l, --limit-shear-velocity VAL limit top shear velocity to this value [m/s] 34 (default nan), overrides --stress-ratio and 35 --set-shear-velocity 36 -V, --velocity-bottom VAL base velocity at bottom [m/s] (default 0) 37 -A, --nonlocal-amplitude VAL amplitude of nonlocality [-] (default 0.4) 38 -b, --rate-dependence VAL rate dependence beyond yield [-] (default 0.9377) 39 -f, --friction-coefficient VAL grain friction coefficient [-] (default 0.404026) 40 -C, --cohesion VAL grain cohesion [Pa] (default 0) 41 -p, --porosity VAL porosity fraction [-] (default 0.25) 42 -d, --grain-size VAL representative grain size [m] (default 0.04) 43 -r, --density VAL grain material density [kg/m^3] (default 2600) 44 -n, --resolution VAL number of cells in domain [-] (default 100) 45 -o, --origo VAL coordinate system origo [m] (default 0) 46 -L, --length VAL domain length [m] (default 0.7) 47 48 Optional arguments only relevant with transient (fluid) simulation: 49 -F, --fluid enable pore fluid computations 50 -c, --fluid-compressibility VAL fluid compressibility [Pa^-1] (default 3.9e-10) 51 -i, --fluid-viscosity VAL fluid viscosity [Pa*s] (default 0.001787) 52 -R, --fluid-density VAL fluid density [kg/m^3] (default 1000) 53 -k, --fluid-permeability VAL fluid intrinsic permeability [m^2] (default 1.9e-15) 54 -O, --fluid-pressure-top VAL fluid pressure at +z edge [Pa] (default 0) 55 -a, --fluid-pressure-ampl VAL amplitude of pressure variations [Pa] (default 0) 56 -q, --fluid-pressure-freq VAL frequency of pressure variations [s^-1] (default 1) 57 -H, --fluid-pressure-phase VAL fluid pressure at +z edge [Pa] (default 0) 58 -t, --time VAL simulation start time [s] (default 0) 59 -T, --time-end VAL simulation end time [s] (default 1) 60 -I, --file-interval VAL interval between output files [s] (default 0.1) 61 62 ```