tadd explicit call for periodic boundaries - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
HTML git clone git://src.adamsgaard.dk/sphere
DIR Log
DIR Files
DIR Refs
DIR LICENSE
---
DIR commit 58668c68b1fce324d35d1c86991240dabd48591e
DIR parent ad617bec05a87782ca2fc5a55b50c8765bb216ff
HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Fri, 20 Jun 2014 21:40:33 +0200
add explicit call for periodic boundaries
Diffstat:
M python/cube-init.py | 2 ++
M python/fluidshear.py | 1 +
M python/segregation.py | 2 ++
M python/shear-test.py | 9 +++++++++
4 files changed, 14 insertions(+), 0 deletions(-)
---
DIR diff --git a/python/cube-init.py b/python/cube-init.py
t@@ -5,6 +5,8 @@ init = sphere.sim('cube-init', np=1e2)
init.generateRadii(psd='uni', radius_mean=0.01, radius_variance=0.002)
+init.periodicBoundariesXY()
+
# Initialize positions in random grid (also sets world size)
init.initRandomGridPos(gridnum=(6, 6, 1e12))
DIR diff --git a/python/fluidshear.py b/python/fluidshear.py
t@@ -9,6 +9,7 @@ sid = 'fluidshear'
sim = sphere.sim(sid + '-init', np = 24000, fluid = False)
#sim.cleanup()
sim.radius[:] = 0.05
+sim.periodicBoundariesXY()
sim.initRandomGridPos(gridnum = [20, 20, 9000])
sim.initTemporal(total = 10.0, file_dt = 0.05)
sim.g[2] = -9.81
DIR diff --git a/python/segregation.py b/python/segregation.py
t@@ -32,6 +32,8 @@ init.generateRadii(radius_mean = 0.08)
# Use default params
init.defaultParams(gamma_n = 100.0, mu_s = 0.4, mu_d = 0.4)
+init.periodicBoundariesXY()
+
# Initialize positions in random grid (also sets world size)
hcells = np**(1.0/3.0)*0.6
init.initRandomGridPos(gridnum = numpy.array([hcells, hcells, 1e9]))
DIR diff --git a/python/shear-test.py b/python/shear-test.py
t@@ -34,6 +34,9 @@ init.defaultParams(gamma_n = 100.0, mu_s = 0.6, mu_d = 0.6)
# Add gravity
init.g[2] = -9.81
+# Periodic x and y boundaries
+init.periodicBoundariesXY()
+
# Initialize positions in random grid (also sets world size)
hcells = np**(1.0/3.0)
init.initRandomGridPos(gridnum = numpy.array([hcells, hcells, 1e9]), periodic = 1, contactmodel = 2)
t@@ -71,6 +74,9 @@ for devs in devslist:
lastf = status(sim_id + "-init")
cons.readbin("../output/" + sim_id + "-init.output{:0=5}.bin".format(lastf), verbose=False)
+ # Periodic x and y boundaries
+ cons.periodicBoundariesXY()
+
# Setup consolidation experiment
cons.consolidate(deviatoric_stress = devs, periodic = init.periodic)
t@@ -117,6 +123,9 @@ for devs in devslist:
lastf = status(sim_id + "-cons-devs{}".format(devs))
shear.readbin("../output/" + sim_id + "-cons-devs{}.output{:0=5}.bin".format(devs, lastf), verbose = False)
+ # Periodic x and y boundaries
+ shear.periodicBoundariesXY()
+
# Setup shear experiment
shear.shear(shear_strain_rate = 0.05, periodic = init.periodic)