tchange world size, particle number, add wall viscous damping - sphere - GPU-ba… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit dd60c1b8ba232d30b61afe79541a10399bc5f257 | |
parent cb145bb36bca6fad2bdcdbe633d5c135eb3bf66a | |
Author: Anders Damsgaard Christensen <[email protected]> | |
Date: Wed, 15 Jun 2016 18:10:29 -0700 | |
change world size, particle number, add wall viscous damping | |
Diffstat: | |
M python/alejandro-init-cons.py | 8 ++++---- | |
1 file changed, 4 insertions(+), 4 deletions(-) | |
--- | |
diff --git a/python/alejandro-init-cons.py b/python/alejandro-init-cons.py | |
t@@ -12,7 +12,7 @@ plots = True | |
device = 0 | |
# Number of particles | |
-np = 2e3 | |
+np = 1e4 | |
# Common simulation id | |
sim_id = "alejandro" | |
t@@ -28,8 +28,9 @@ init = sphere.sim(np=np, nd=3, nw=0, sid=sim_id + "-init") | |
# Save radii | |
init.generateRadii(mean=0.01) | |
-# Use default params | |
+# Add viscous damping to quickly dissipate kinetic energy | |
init.defaultParams(k_n=1.16e7, k_t=1.16e7, gamma_n=100.0, mu_s=0.5, mu_d=0.5) | |
+init.gamma_wn[0] = 10000.0 | |
# Add gravity | |
init.g[2] = -9.81 | |
t@@ -38,8 +39,7 @@ init.g[2] = -9.81 | |
init.periodicBoundariesX() | |
# Initialize positions in random grid (also sets world size) | |
-hcells = np**(1.0/3.0) | |
-init.initRandomGridPos(gridnum=[hcells, hcells, 1e9]) | |
+init.initRandomGridPos(gridnum=[12, 12, 1e9]) | |
# Set duration of simulation | |
init.initTemporal(total=10.0, epsilon=0.07) |