cngf-pf-exp1

commit 7640b6b7f669a0d19a397a1e7ffdb38c7b42822b
parent 933f8ef343f0090357c78e6c22e57ddf5fb6c8c5
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 12 Dec 2019 14:27:02 +0100

Calculate arrays of max. deformation depths, update 1d_fd_simple_shear

Diffstat:
M
fig-skin_depth/Makefile
 | 
25
++++++++++++++++++++++++-
M
fig-skin_depth/fig.gp
 | 
1
-

2 files changed, 24 insertions(+), 2 deletions(-)
diff --git a/fig-skin_depth/Makefile b/fig-skin_depth/Makefile
@@ -1,11 +1,34 @@
 FIG != basename $(PWD)
+BIN = ../1d_fd_simple_shear/max_depth_simple_shear
 
 default: ../$(FIG).pdf
 
-../$(FIG).pdf: fig.gp
+../$(FIG).pdf: fig.gp max_depth_ampl10e3.txt
 	gnuplot fig.gp > $@
 
+max_depth_ampl10e3.txt: $(BIN)
+	for A_f in 10e3 100e3; do\
+		out=max_depth_ampl$${A_f}.txt;\
+		rm -f $$out;\
+		for f in 1e-9 1e-8 1e-7 1e-6 1e-5 1e-4 1e-3; do\
+			for k in 1e-19 1e-18 1e-17 1e-16 1e-15 1e-14; do\
+				./$(BIN) \
+					--porosity 0.25 \
+					--fluid-compressibility 3.9e-10 \
+					--fluid-viscosity 1.787e-3 \
+					--fluid-permeability $$k \
+					--normal-stress 200e3 \
+					--fluid-pressure-top 100e3 \
+					--fluid-pressure-ampl $$A_f \
+					--fluid-pressure-freq $$f | \
+					awk '{print $$1}' | tr '\n' '\t' >> $$out;\
+			done;\
+			echo >> $$out;\
+		done;\
+	done
+
 clean:
 	rm -f ../$(FIG).pdf
+	rm -f max_depth_ampl*.txt
 
 .PHONY: default clean
diff --git a/fig-skin_depth/fig.gp b/fig-skin_depth/fig.gp
@@ -25,7 +25,6 @@ f_max = 2.0/3600.0
 
 # f: frequency, k: permeability
 skindepth(f,k) = (k/(phi*mu_f*beta_f*3.141592654*f))**(0.5)  # f: 1/s
-#skindepth(f,k) = (k/(phi*mu_f*beta_f*3.141592654*(f/3600)))**(0.5)  # f: 1/h
 
 set xlabel "Forcing frequency [1/s]"
 set ylabel "Permeability [m^2]"