tneumann test passes with set2 - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
(HTM) git clone git://src.adamsgaard.dk/sphere
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) LICENSE
---
(DIR) commit 86fdfad131a636d5ce80131892b66d157250f0c5
(DIR) parent 2b64c71fd50ab41816d321bf665d7ba0b2e0475e
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Tue, 17 Jun 2014 13:46:51 +0200
neumann test passes with set2
Diffstat:
M tests/cfd_tests_neumann.py | 8 +++-----
1 file changed, 3 insertions(+), 5 deletions(-)
---
(DIR) diff --git a/tests/cfd_tests_neumann.py b/tests/cfd_tests_neumann.py
t@@ -15,7 +15,7 @@ cleanup(orig)
orig.defaultParams(mu_s = 0.4, mu_d = 0.4)
orig.defineWorldBoundaries([0.4, 0.4, 1], dx = 0.1)
orig.initFluid(mu = 8.9e-4)
-orig.initTemporal(total = 1.0, file_dt = 0.05, dt = 1.0e-4)
+orig.initTemporal(total = 0.5, file_dt = 0.05, dt = 1.0e-4)
py = sphere.sim(sid = orig.sid, fluid = True)
orig.bc_bot[0] = 1 # No-flow BC at bottom (Neumann)
#orig.run(dry=True)
t@@ -34,7 +34,6 @@ else:
print("Flow field:\t\t" + failed())
raise Exception("Failed")
-"""
print('''# Neumann bottom, Dirichlet top BC.
# Gravity, pressure gradients => transient flow''')
orig = sphere.sim("neumann", fluid = True)
t@@ -47,7 +46,7 @@ orig.g[2] = -10.0
orig.bc_bot[0] = 1 # No-flow BC at bottom (Neumann)
#orig.run(dry=True)
orig.run(verbose=False)
-#orig.writeVTKall()
+orig.writeVTKall()
py.readlast(verbose = False)
ideal_grad_p_z = numpy.linspace(
orig.p_f[0,0,0] + orig.L[2]*orig.rho_f*numpy.abs(orig.g[2]),
t@@ -62,5 +61,4 @@ else:
print("Flow field:\t\t" + failed())
raise Exception("Failed")
-#orig.cleanup()
-"""
+orig.cleanup()