tFinish visualization - cngf-pf - continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
HTML git clone git://src.adamsgaard.dk/cngf-pf
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit ffb275d4bbd029b7884129340ebfcbddc74b2852
DIR parent 6cb6fc35f0aa16cadba133b4bd63f7376c7cdc58
HTML Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 5 Apr 2019 20:38:56 +0200
Finish visualization
Diffstat:
A 1d_fd_simple_shear.gp | 28 ++++++++++++++++++++++++++++
M 1d_fd_simple_shear.png | 0
M 1d_fd_simple_shear_damsgaard2013.h | 4 ++--
M Makefile | 17 +++++++++++++++++
M arrays.c | 7 +++++++
M arrays.h | 1 +
R 1d_fd_simple_shear.jl -> julia/1d_… | 0
C 1d_fd_simple_shear.png -> julia/1d… | 0
M main.c | 2 +-
9 files changed, 56 insertions(+), 3 deletions(-)
---
DIR diff --git a/1d_fd_simple_shear.gp b/1d_fd_simple_shear.gp
t@@ -0,0 +1,28 @@
+#!/usr/bin/env gnuplot
+
+#set terminal epslatex color size 8.6 cm, 7.6 cm
+#set output "plots/sheardisp.tex"
+#set terminal pdfcairo color size 8.6 cm, 7.6 cm
+#set output "plots/sheardisp.pdf"
+#set terminal pngcairo color size 8.6 cm, 7.6 cm
+set terminal pngcairo color size 18.6 cm, 17.6 cm
+set output "1d_fd_simple_shear.png"
+set xlabel "Horizontal velocity, v_x [m]"
+set ylabel "Vertical position, x_3 [m]" offset 2
+
+#set grid
+
+set yrange [0.0:0.73]
+
+set key bottom right #samplen 0.9
+
+set style line 1 linetype 1 linewidth 3 pointtype 1 pointsize 1
+set style line 2 linetype 2 linewidth 3 pointtype 2 pointsize 1
+
+plot "1d_fd_simple_shear_P10kPa.txt" u 2:1 w lp lt 1 lw 2 title "P_{wall} = 10 kPa", \
+ "1d_fd_simple_shear_P20kPa.txt" u 2:1 w lp lt 2 lw 2 title "P_{wall} = 20 kPa", \
+ "1d_fd_simple_shear_P40kPa.txt" u 2:1 w lp lt 3 lw 2 title "P_{wall} = 40 kPa", \
+ "1d_fd_simple_shear_P80kPa.txt" u 2:1 w lp lt 4 lw 2 title "P_{wall} = 80 kPa", \
+ "1d_fd_simple_shear_P120kPa.txt" u 2:1 w lp lt 5 lw 2 title "P_{wall} = 120 kPa"
+
+set xtics norotate # Restore defaults
DIR diff --git a/1d_fd_simple_shear.png b/1d_fd_simple_shear.png
Binary files differ.
DIR diff --git a/1d_fd_simple_shear_damsgaard2013.h b/1d_fd_simple_shear_damsgaard2013.h
t@@ -16,11 +16,11 @@ struct simulation init_sim(void)
sim.G = 9.81;
- sim.P_wall = 10e3; /* larger normal stress deepens the shear depth */
+ sim.P_wall = 120e3; /* larger normal stress deepens the shear depth */
sim.mu_wall = 0.40;
sim.v_x_bot = 0.0;
- sim.nz = 10;
+ sim.nz = 100;
/* lower values of A mean that the velocity curve can have sharper curves,
* e.g. at the transition from μ ≈ μ_s */
DIR diff --git a/Makefile b/Makefile
t@@ -17,11 +17,28 @@ default: 1d_fd_simple_shear
1d_fd_simple_shear: $(OBJ) $(HDR)
$(CC) $(LDFLAGS) $(OBJ) -o $@
+1d_fd_simple_shear.png: 1d_fd_simple_shear 1d_fd_simple_shear.gp
+ sed -i 's/P_wall = .*;/P_wall = 10e3;/' 1d_fd_simple_shear_damsgaard2013.h
+ rm main.o && make $< && ./$< > $<_P10kPa.txt
+ sed -i 's/P_wall = .*;/P_wall = 20e3;/' 1d_fd_simple_shear_damsgaard2013.h
+ rm main.o && make $< && ./$< > $<_P20kPa.txt
+ sed -i 's/P_wall = .*;/P_wall = 40e3;/' 1d_fd_simple_shear_damsgaard2013.h
+ rm main.o && make $< && ./$< > $<_P40kPa.txt
+ sed -i 's/P_wall = .*;/P_wall = 60e3;/' 1d_fd_simple_shear_damsgaard2013.h
+ rm main.o && make $< && ./$< > $<_P60kPa.txt
+ sed -i 's/P_wall = .*;/P_wall = 80e3;/' 1d_fd_simple_shear_damsgaard2013.h
+ rm main.o && make $< && ./$< > $<_P80kPa.txt
+ sed -i 's/P_wall = .*;/P_wall = 120e3;/' 1d_fd_simple_shear_damsgaard2013.h
+ rm main.o && make $< && ./$< > $<_P120kPa.txt
+ gnuplot $<.gp > $@
+
.PHONY: watch
watch:
echo $(SRC) $(HDR) | tr ' ' '\n' | entr -s 'make && ./1d_fd_simple_shear'
clean:
$(RM) *.png
+ $(RM) *.txt
$(RM) *.o
$(RM) 1d_fd_simple_shear
+ $(RM) 1d_fd_simple_shear.png
DIR diff --git a/arrays.c b/arrays.c
t@@ -109,6 +109,13 @@ void print_arrays(const double* a, const double* b, const int n)
printf("%.17g\t%.17g\n", a[i], b[i]);
}
+void print_arrays_2nd_normalized(const double* a, const double* b, const int n)
+{
+ double max_b = max(b, n);
+ for (int i=0; i<n; ++i)
+ printf("%.17g\t%.17g\n", a[i], b[i]/max_b);
+}
+
void copy_values(const double* in, double* out, const int n)
{
for (int i=0; i<n; ++i)
DIR diff --git a/arrays.h b/arrays.h
t@@ -27,6 +27,7 @@ double min(const double* a, const int n);
void print_array(const double* a, const int n);
void print_arrays(const double* a, const double* b, const int n);
+void print_arrays_2nd_normalized(const double* a, const double* b, const int n);
void copy_values(const double* in, double* out, const int n);
DIR diff --git a/1d_fd_simple_shear.jl b/julia/1d_fd_simple_shear.jl
DIR diff --git a/1d_fd_simple_shear.png b/julia/1d_fd_simple_shear.png
Binary files differ.
DIR diff --git a/main.c b/main.c
t@@ -30,7 +30,7 @@ int main(int argc, char** argv) {
compute_shear_strain_rate_plastic(&sim);
compute_shear_velocity(&sim);
- print_arrays(sim.z, sim.v_x, sim.nz);
+ print_arrays_2nd_normalized(sim.z, sim.v_x, sim.nz);
free_arrays(&sim);
return 0;