cngf-pf-exp1

commit 0f6b925eed1255799f8724643f772bd85f77a05e
parent 63709c7bdabc7953803e85e4b1ac902388f619de
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon,  1 Jul 2019 07:43:11 +0200

Update 1d_fd_simple_shear, increase domain size in stick-slip experiments

Diffstat:
M
fig2/Makefile
 | 
1
-
M
fig3/Makefile
 | 
3
+--
M
fig3/fig.gp
 | 
2
+-
M
fig4/Makefile
 | 
3
+--
M
fig4/fig.gp
 | 
2
+-
M
p_f_analytical.jl
 | 
4
++--

6 files changed, 6 insertions(+), 9 deletions(-)
diff --git a/fig2/Makefile b/fig2/Makefile
@@ -9,7 +9,6 @@ unnamed.output00000.txt: $(BIN)
 	--fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 \
 	--fluid-pressure-ampl 50e3 \
 	--fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \
-	--time-step 1e-1 \
 	--file-interval $$( echo "60*10" | bc -l ) \
 	--time-end $$( echo "3600*24*7" | bc -l )'
 
diff --git a/fig3/Makefile b/fig3/Makefile
@@ -5,11 +5,10 @@ default: ../$(FIG).pdf
 
 unnamed.output00000.txt: $(BIN) Makefile
 	/bin/sh -c '\
-	./$< --resolution 50 --length 2.0 --normal-stress 150e3 \
+	./$< --resolution 50 --length 8.0 --normal-stress 150e3 \
 	--fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 \
 	--fluid-pressure-ampl 50e3 \
 	--fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \
-	--time-step 1e-1 \
 	--file-interval $$( echo "3600" | bc -l ) \
 	--time-end $$( echo "3600*24*3 + 3600" | bc -l )'
 
diff --git a/fig3/fig.gp b/fig3/fig.gp
@@ -8,7 +8,7 @@ set multiplot layout 1,4 \
 	margins 0.09,0.97,0.23,0.95 \
 	spacing 0.03,0.03
 
-set yrange [0.0:2.0]
+set yrange [4.0:8.0]
 
 set key bottom right #samplen 0.9
 
diff --git a/fig4/Makefile b/fig4/Makefile
@@ -5,11 +5,10 @@ default: ../$(FIG).pdf
 
 unnamed.output00000.txt: $(BIN) Makefile
 	/bin/sh -c '\
-	./$< --resolution 50 --length 2.0 --normal-stress 150e3 \
+	./$< --resolution 50 --length 8.0 --normal-stress 150e3 \
 	--fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 \
 	--fluid-pressure-ampl 50e3 \
 	--fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \
-	--time-step 1e-1 \
 	--file-interval $$( echo "3600" | bc -l ) \
 	--time-end $$( echo "3600*24*3 + 3600" | bc -l ) \
 	--normalize'
diff --git a/fig4/fig.gp b/fig4/fig.gp
@@ -8,7 +8,7 @@ set multiplot layout 1,4 \
 	margins 0.09,0.97,0.23,0.95 \
 	spacing 0.03,0.03
 
-set yrange [0.0:2.0]
+set yrange [4.0:8.0]
 
 set key bottom right #samplen 0.9
 
diff --git a/p_f_analytical.jl b/p_f_analytical.jl
@@ -1,8 +1,8 @@
 #!/usr/bin/env julia
 import PyPlot
 
-n = 50                          # resolution
-z = range(0, stop=2, length=n)  # spatial grid
+n = 200                         # resolution
+z = range(0, stop=8, length=n)  # spatial grid
 k = 2e-17                       # permeability
 ϕ = 0.25                        # porosity
 μ_f = 1e-3                      # water viscosity