hansen-zoet-exp

comparison of lab and model results of sediment advection
git clone git://src.adamsgaard.dk/hansen-zoet-exp # fast
git clone https://src.adamsgaard.dk/hansen-zoet-exp.git # slow
Log | Files | Refs | Submodules Back to index

commit 3c1823c1dfda25deb485bfd1290a023790c29ade
parent 23cb53072e37107bb08f145f3075fc6ed4c6864a
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon, 13 Jun 2022 13:17:11 +0200

update acme.dump

Diffstat:
Macme.dump | 267++++++++++---------------------------------------------------------------------
1 file changed, 31 insertions(+), 236 deletions(-)

diff --git a/acme.dump b/acme.dump @@ -1,239 +1,34 @@ -/home/ad/doc/todo +/home/ad/code/hansen-zoet-exp /lib/font/bit/lucsans/euro.8.font /lib/font/bit/lucm/unicode.9.font - 0.0000000 58.9843750 -w Newcol Kill Putall Dump Exit Dump /home/ad/code/hansen-zoet-exp/acme.dump -c 0 New Cut Paste Snarf Sort Zerox Delcol + 0.0000000 58.9062500 +w m Newcol Kill Putall Dump Exit Dump /home/ad/code/hansen-zoet-exp/acme.dump +c 0 New Cut Paste Snarf Sort Zerox Delcol cngf-pf(1) shear_flux(1) max_de c 1 New Cut Paste Snarf Sort Zerox Delcol -f 0 1 339 339 1.5804598 - 1 374 3853 0 0 /home/ad/code/hansen-zoet-exp/Makefile Del Snarf Undo | Look m mc˙git a acme.dump && git c -m 'update acme.dump'˙git a $%˙git c -m 'Makefile: fix pressure-amplitude values (Pa to kPa)'˙git push˙git pull˙git status˙git diff˙m strain_distribution_ampl1Pa_freq1.157e-5s.txt˙rm strain_distribution_ampl1Pa_freq1.157e-5s.txt˙max_depth_simple_shear -D 1.5e-5 -q 1.157e-5 -a 100e3 -f 1 2 0 0 1.5804598 - 2 52 5643 1 0 /home/ad/code/hansen-zoet-exp/ Del Snarf Get | Look -F 1 1 0 0 3.1609195 8513 - 12 33 8513 0 0 /man/cngf-pf(1) Del Snarf | Look -CNGF-PF(1) BSD General Commands Manual CNGF-PF(1) - -NAME - cngf-pf — simulate granular flows with pore-pressure dynamics - -SYNOPSIS - cngf-pf [-A grain-nonlocal-ampl] [-a fluid-pressure-ampl] - [-b grain-rate-dependence] [-C fluid-compressibility] - [-c grain-cohesion] [-D fluid-diffusivity] [-d grain-size] - [-e end-time] [-F] [-f applied-shear-friction] [-g gravity-accel] - [-H fluid-pressure-phase] [-h] [-I file-interval] - [-i fluid-viscosity] [-j time-step] [-K dilatancy-constant] - [-k fluid-permeability] [-L length] [-l applied-shear-vel-limit] - [-m grain-friction] [-N] [-n normal-stress] - [-O fluid-pressure-top] [-o origo] [-P grain-compressibility] - [-p grain-porosity] [-q fluid-pressure-freq] [-R fluid-density] - [-r grain-density] [-S fluid-pressure-pulse-shape] - [-s applied-shear-vel] [-T] [-t curr-time] [-U resolution] - [-u fluid-pulse-time] [-V bot-vel] [-v] [-Y max-porosity] - [-y min-porosity] [name] - -DESCRIPTION - The cngf-pf utility simulates the coupled continuum dynamics of granular - deformation and pore pressure. The simulation space is one-dimensional - and undergoes confined, simple shear at a pretedermined velocity or fric‐ - tion ratio. The effective stress is modulated by pore-pressure dynamics - when fluid dynamics are enabled. - - The arguments are as follows: - - -A grain-nonlocal-ampl - Amplitude of granular non-locality [-] (default 0.4). - - -a fluid-pressure-ampl - Amplitude of fluid-pressure perturbations [Pa] (default 0). Only - relevant with fluid dynamics enabled (-F). - - -b grain-rate-dependence - Granular rate dependence beyond yield [-], (default 0.9377). - - -C fluid-compressibility - Fluid adiabatic compressibility [Pa^-1] (default 3.9e-10). Only - relevant with fluid dynamics enabled (-F). - - -c grain-cohesion - Granular material cohesion [Pa] (default 0). - - -D fluid-diffusivity - Fluid diffusion coefficient [m^2/s] (default -1). Overrides - fluid permeability (-k), grain compressibility (-P), fluid com‐ - pressibility (-C), and fluid viscosity (-i). Do not use for - transient simulations (-T). Disabled when set to a negative - value. - - -d grain-size - Granular material representative grain size [m] (default 0.04). - - -e end-time - Simulation end time [s] (default 1). - - -F Enable pore-space fluid dynamics. - - -f applied-shear-friction - Applied ratio between shear stress and normal stress at the top - boundary [-] (default 0.45), overridden by -s and -l if these are - set. - - -g gravity-accel - Gravity magnitude [m/s^2] (default 9.81). - - -H fluid-pressure-phase - Phase of fluid-pressure perturbations [s] (default 0). Only rel‐ - evant with fluid dynamics enabled (-F). - - -h Show usage information. - - -I file-interval - Simulation time interval between writing output to disk [s] (de‐ - fault 1.0). - - -i fluid-viscosity - Fluid dynamic viscosity [Pa*s] (1.787e-3). Only relevant with - fluid dynamics enabled (-F). - - -j time-step - Override numerical time step length [s]. Only relevant with - fluid dynamics enabled (-F). Outside of stability testing, is - not recommended to specify the time step length manually. - - -K dilatancy-constant - Factor relating dilatancy angle to the volume fraction [-] (de‐ - fault 4.09). Only relevant with transient granular dynamics en‐ - abled (-T). - - -k fluid-permeability - Darcian intrinsic permeability of granular material [m^2] (de‐ - fault 1.9e-15). Only relevant with fluid dynamics enabled (-F), - and without transient dynamics. In transient simulations (-T) - the permeability value is overridden by the Kozeny-Carman rela‐ - tionship (k = d^2/180 * phi^3 * (1 - phi)). - - -L length - Simulation domain length [m] (default 1). - - -l applied-shear-vel-limit - Upper limit to modeled shear velocity [m/s] (default nan), over‐ - rides -f and -s. - - -m grain-friction - Grain friction coefficient [-] (default 0.404026). - - -N Normalize output velocities to the range [0;1]. - - -n normal-stress - Normal stress applied at the top boundary [Pa] (default 120e3). - - -O fluid-pressure-top - Fluid pressure applied at the top boundary [Pa] (default 0). - Only relevant with fluid dynamics enabled (-F). - - -o origo - Simulation domain origo [m] (default 0). - - -P grain-compressibility - Granular material compressibility [Pa^-1] (default 1e-8). - - -p grain-porosity - Granular material porosity [-] (default 0.25). - - -q fluid-pressure-freq - Frequency of fluid-pressure perturbations [s^-1] (default 1). - Only relevant with fluid dynamics enabled (-F). - - -R fluid-density - Fluid density [kg/m^3] (default 1000). Only relevant with fluid - dynamics enabled (-F). - - -r grain-density - Granular material density [kg/m^3] (default 2600). - - -S fluid-pressure-pulse-shape - Shape of pulse perturbations in fluid pressure at the top bound‐ - ary. Valid options are triangular (default) or square. Only - relevant with fluid dynamics enabled (-F). - - -s applied-shear-vel - Shear the material under constant velocity [m/s] instead of a - friction value. Overrides -f and is overridden by -l. - - -T Enable transient granular properties. - - -t curr-time - Simulation start time [s] (default 0). - - -U resolution - Numerical resolution of simulation domain [-] (default nan). If - not set, the cell size equals the grain size. - - -u fluid-pulse-time - Time of fluid-pressure pulse peak [s] (default nan). Only rele‐ - vant with fluid dynamics enabled (-F). - - -V bot-vel - shear velocity at base [m/s] (default 0). - - -v Show version information. - - -Y max-porosity - Granular material porosity [-] at inertia number I=1 in transient - simulations (-T) (default 0.55). - - -y min-porosity - Minimum granular material porosity [-] in transient simulations - (-T) (default 0.20). - - The final simulation state is written to stdout, see OUTPUT FORMAT for a - description of the data format. - -OUTPUT FORMAT - The output consists of the following tab-delimited fields, with one row - per cell in the simulation domain: - - 1. position [m] - 2. shear velocity [m/s] - 3. effective normal stress [Pa] - 4. fluid pressure [Pa] - 5. friction [-] - 6. shear strain rate [-] - 7. porosity [-] - 8. inertia number [-] - 9. shear stress [Pa] - 10. cumulative shear displacement [m] - -FILES - If name is set, the simulation state is written to files in the current - directory named in the format <name>.outputDDDDD.txt. Here, DDDDD is a - integer counter starting from 00000, and incrementing by 1 for each new - file. Each file contains the state at a single point in simulation time. - The output format is specified in the section OUTPUT FORMAT, and output - files are generated with the interval specified with -I file-interval. - -EXIT STATUS - cngf-pf exits 0 on success, and >0 if a runtime error occurs: - - 0 successful exit - 1 unspecified error - 2 pledge(2) error - 10 transient solver error - 11 fluid solver error - 12 granular solver error - 20 time step error - -EXAMPLES - Plot critical-state shear velocity as a function of depth: - - $ cngf-pf | gnuplot -e 'set term pdf; p "-" u 2:1 w lp' > out.pdf - -SEE ALSO - gnuplot(1), max_depth_simple_shear(1), shear_flux(1) - -AUTHORS - Anders Damsgaard <anders@adamsgaard.dk> - -BSD June 10, 2022 BSD -F 1 2 0 0 23.8026820 0 - 4 70 0 0 1 /home/ad/code/hansen-zoet-exp/+Errors Del Snarf Undo Put | Look Clear +f 0 1 339 339 1.2108262 + 1 369 3761 0 0 /home/ad/code/hansen-zoet-exp/Makefile Del Snarf | Look m mc˙git a acme.dump && git c -m 'update acme.dump'˙git a $%˙git c -m 'Makefile: fix pressure-amplitude values (Pa to kPa)'˙git push˙git pull˙git status˙git diff˙m strain_distribution_ampl1Pa_freq1.157e-5s.txt˙rm strain_distribution_ampl1Pa_freq1.157e-5s.txt˙max_depth_simple_shear -D 1.5e-5 -q 1.157e-5 -a 100e3 +f 1 2 0 0 1.2108262 + 2 52 6626 1 0 /home/ad/code/hansen-zoet-exp/ Del Snarf Get | Look +F 1 1 391 391 16.3817664 391 + 4 70 391 0 1 /home/ad/code/hansen-zoet-exp/+Errors Del Snarf Undo Put | Look Clear +sh -c '\ +for freq in 1.157e-5 ; do \ + for ampl in 1 34 67 100; do \ + ./cngf-pf/cngf-pf \ + -d 530.7e-6 \ + -D 1.5e-5 \ + -m 0.6 \ + -c 0 \ + -s 3.169e-6 \ + -L 0.5 \ + -n 20e3 \ + -A 0.33 \ + -b 0.002 \ + -F \ + -e 345600 \ + -I 1800 \ + -q ${freq} \ + -a ${ampl}e3 \ + strain_distribution_ampl${ampl}kPa_freq${freq}s \ + >strain_distribution_ampl${ampl}kPa_freq${freq}s.txt; \ +done; done'