granular-basin

commit 5e50d09f3c24e4b80cc61b4f64c53774901d1106
parent 845a5bc453c83e6d59a61d86b9511575b001085d
Author: esbenpalmstrom <esbenpalmstroem@gmail.com>
Date:   Wed,  1 Dec 2021 09:29:38 +0100

minor changes

Diffstat:
M
deform_basin.jl
 | 
13
++++++++-----
M
init_basin.jl
 | 
2
+-

2 files changed, 9 insertions(+), 6 deletions(-)
diff --git a/deform_basin.jl b/deform_basin.jl
@@ -10,7 +10,7 @@ id = "simulation500"   # folder name of simulation
 
 hw_ratio = 0.2          # height/width ratio of indenter
 grain_radius = 0.05     # grain radius of grains in indenter
-def_time = 2.0          # time spent deforming
+def_time = 2.0         # time spent deforming
 
 deformation_type = "shortening" # "diapir" or "shortening"
                                 # diapir will only introduce an indenter while
@@ -137,16 +137,19 @@ Granular.addWallLinearFrictionless!(sim, [0.,1.],
 
 global checked_done = false
 
+#sim.walls[1].vel = -boomerang_vel
+#sim.walls[2].vel = boomerang_vel
+
 while sim.time < sim.time_total
 
-    if sim.grains[left_edge_index].lin_vel[1] > boomerang_vel/2 && checked_done == false
-        sim.walls[1].vel = boomerang_vel
-        sim.walls[2].vel = -boomerang_vel
+    if sim.grains[right_edge_index].lin_vel[1] > boomerang_vel/3 && checked_done == false
+        sim.walls[1].vel = -boomerang_vel
+        sim.walls[2].vel = boomerang_vel
         global checked_done = true
     end
 
     #sim.walls[1].vel = sim.grains[left_edge_index].lin_vel[1]
-    #sim.walls[2].vel = -sim.grains[right_edge_index].lin_vel[1]
+    #sim.walls[2].vel = sim.grains[right_edge_index].lin_vel[1]
 
     Granular.run!(sim,single_step = true)
 
diff --git a/init_basin.jl b/init_basin.jl
@@ -12,7 +12,7 @@ t_stack = 0.5                   # duration for each stack to settle [s]
 
 g = [0.,-9.8]                   # vector for direction and magnitude of gravitational acceleration of grains
 
-ngrains = 500                  # total number of grains
+ngrains = 500                   # total number of grains
 aspect_ratio = 4                # should be x times as wide as it is tall
 
 mkpath("simulation$(ngrains)")