tbegin implementing fluid oscillation experiment - hansen-zoet-exp - comparison… | |
git clone git://src.adamsgaard.dk/hansen-zoet-exp | |
Log | |
Files | |
Refs | |
Submodules | |
--- | |
commit 34849286277708d5f92ba66797431c8c1444e27e | |
parent 29fdaa065b7f780d0dce0d988065804be087579e | |
Author: Anders Damsgaard <[email protected]> | |
Date: Tue, 7 Jun 2022 15:39:24 +0200 | |
begin implementing fluid oscillation experiment | |
Diffstat: | |
M Makefile | 51 +++++++++++++++++++++++------… | |
1 file changed, 39 insertions(+), 12 deletions(-) | |
--- | |
diff --git a/Makefile b/Makefile | |
t@@ -12,10 +12,19 @@ stressvals_ != range -l -N 40 0.01 100 | |
stressvals_ += 8 22 50 | |
stressvals != echo ${stressvals_} | sort -n | |
velvals = 3.169e-6 9.506e-6 | |
-ampvals = range -N 4 0.1 1.0 | |
-freq = 1.157e-05 # daily | |
+fluid_ampl_vals != range -N 4 1 10 | |
+fluid_freq_vals = 1.157e-5 # daily | |
+grain_size = 530.7e-6 | |
+grain_friction = 0.6 | |
+grain_cohesion = 0 | |
+grain_nonlocal_ampl = 0.33 | |
+grain_rate_dependence = 0.002 | |
+# grain_nonlocal_ampl doubles -> scale factor doubles, stress exponent is cons… | |
+# sediment flux is proportional to cuberoot(N) | |
+# grain_nonlocal_ampl = 0.66 -> a = 2.62 | |
+# grain_nonlocal_ampl = 0.33 -> a = 1.79 | |
-default: ${PDFS} | |
+default: ${PDFS} strain_distribution_ampl1Pa_freq1.157e-5s.txt | |
${BIN}: ${REPO}/cngf-pf.c | |
CFLAGS="-march=native -O3 -pipe" make -C ${REPO} | |
t@@ -28,22 +37,40 @@ strain_distribution_N50kPa_v3.169e-6ms.txt: ${BIN} Makefile | |
for v in ${velvals}; do \ | |
for P in ${stressvals}; do \ | |
./${BIN} \ | |
- -d 530.7e-6 \ | |
- -m 0.6 \ | |
- -c 0e3 \ | |
+ -d ${grain_size} \ | |
+ -m ${grain_friction} \ | |
+ -c ${grain_cohesion} \ | |
-s $${v} \ | |
-L 0.11 \ | |
-n $${P}e3 \ | |
- -A 0.33 \ | |
- -b 0.002 \ | |
+ -A ${grain_nonlocal_ampl} \ | |
+ -b ${grain_rate_dependence} \ | |
> strain_distribution_N$${P}kPa_v$${v}ms.txt; \ | |
done; done' # also close fit with lower A and higher b | |
- # A doubles -> scale factor doubles, stress exponent is constant at 1/3 | |
- # q \propto cuberoot(N) | |
- # A = 0.66 -> a = 2.62 | |
- # A = 0.33 -> a = 1.79 | |
# higher values of L do not make a difference in flux | |
+strain_distribution_ampl1Pa_freq1.157e-5s.txt: ${BIN} | |
+ sh -c '\ | |
+ for freq in ${fluid_freq_vals}; do \ | |
+ for ampl in ${fluid_ampl_vals}; do \ | |
+ ./${BIN} \ | |
+ -d ${grain_size} \ | |
+ -D 1.5e-5 \ | |
+ -m ${grain_friction} \ | |
+ -c ${grain_cohesion} \ | |
+ -s 3.169e-6 \ | |
+ -L 0.5 \ | |
+ -n 20e3 \ | |
+ -A ${grain_nonlocal_ampl} \ | |
+ -b ${grain_rate_dependence} \ | |
+ -F \ | |
+ -e 345600 \ | |
+ -I 1800 \ | |
+ -q $${freq} \ | |
+ -a $${ampl} \ | |
+ > strain_distribution_ampl$${ampl}Pa_freq$${freq}s.txt… | |
+ done; done' | |
+ | |
mohr_coulomb_v3.169e-6ms.txt: strain_distribution_N50kPa_v3.169e-6ms.txt | |
sh -c '\ | |
rm -f mohr_coulomb_v*ms.txt; \ |