1d_fd_simple_shear_transient

commit c1aa4dce15eba37cc67913aa232337e9205feb70
parent 7a535a7e52e53c1d5967573ab4b97a062198c120
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Fri,  5 Apr 2019 08:55:34 +0200

Improve comments

Diffstat:
M
1d_fd_simple_shear.jl
 | 
42
++++++++++++++++++++++++++----------------

1 file changed, 26 insertions(+), 16 deletions(-)
diff --git a/1d_fd_simple_shear.jl b/1d_fd_simple_shear.jl
@@ -12,43 +12,53 @@ G = 9.81
 
 ## Wall parameters
 
+### Effective normal stress exerted by top wall
 # A larger normal stress deepens the deformational depth
 P_wall_ = [10, 20, 40, 60, 80, 120] .* 1e3  # normal stress [Pa]
 
-v_x_bot = 0.0   # bottom velocity along x [m/s]
-μ_wall = 0.40  # stress ratio at top wall
+### bottom velocity along x [m/s]
+v_x_bot = 0.0
 
-## Nodes along z
+# stress ratio at top wall
+μ_wall = 0.40
+
+### Nodes along z
 nz = 100
 
 ## Material properties
 
+### nonlocal amplitude [-]
 # lower values of A mean that the velocity curve can have sharper curves, e.g.
 # at the transition from μ ≈ μ_s
-A = 0.48        # nonlocal amplitude [-]
-
+A = 0.48        # Henann and Kamrin 2016
 
+### rate dependence beyond yield [-]
 # lower values of b mean larger shear velocity for a given stress ratio above
 # yield
-b = 0.9377      # rate dependence beyond yield [-]
+b = 0.9377      # Henann and Kamrin 2016
 
-μ_s = 0.3819    # static yield friction coefficient [-]
+### bulk and critical state static yield friction coefficient [-]
+#μ_s = 0.3819               # Henann and Kamrin 2016
 μ_s = atan(deg2rad(22.0))  # Damsgaard et al 2013
-ϕ = 0.38        # porosity [-]
 
+### porosity [-]
+#ϕ = 0.38        # Henann and Kamrin 2016
+ϕ = 0.28        # Damsgaard et al 2013
+
+# representative grain size [m]
 # lower values of d mean that the shear velocity curve can have sharper curves,
 # e.g. at the transition from μ ≈ μ_s
-#const d = 1e-3        # representative grain size [m]
-d = 0.04        # representative grain size [m] - Damsgaard et al 2013
+#d = 1e-3        # Henann and Kamrin 2016
+d = 0.04        # Damsgaard et al 2013
 
-ρ_s = 2.485e3   # grain material density [kg/m^3]
+### grain material density [kg/m^3]
+#ρ_s = 2.485e3   # Henann and Kamrin 2016
+ρ_s = 2.6e3     # Damsgaard et al 2013
 
-## Domain size
+## Spatial settings
 origo_z = 0.0
-#const L_z = 20.0*d
-L_z = 0.7
-
-## Spatial coordiantes
+#L_z = 20.0*d    # Henann and Kamrin 2016
+L_z = 0.7       # Damsgaard et al 2013
 z = collect(range(origo_z, L_z, length=nz))
 Δz = z[2] - z[1]