tcomplete test to verify horizontal pressure gradient - 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 7c54905221253a4bee645c39829c788eb3121eac
(DIR) parent 7a1131484a0e2ac6b805cfdd87c07039e4b4bbfb
(HTM) Author: Anders Damsgaard Christensen <adc@geo.au.dk>
Date: Fri, 12 Aug 2016 13:59:05 -0700
complete test to verify horizontal pressure gradient
Diffstat:
M tests/cfd_tests_darcy.py | 15 ++++-----------
1 file changed, 4 insertions(+), 11 deletions(-)
---
(DIR) diff --git a/tests/cfd_tests_darcy.py b/tests/cfd_tests_darcy.py
t@@ -216,11 +216,11 @@ py = sphere.sim(sid = orig.sid, fluid = True)
print("# Pressure gradient along X")
orig.cleanup()
orig.p_f[:,:,:] = 0.0
-orig.p_f[0,:,-1] = 2.0
-orig.p_f[-1,:,-1] = 0.0
+orig.p_f[0,:,:] = 2.0
orig.setFluidXFixedPressure()
orig.setFluidYNoFlow()
-orig.setFluidTopFixedPressure()
+#orig.setFluidTopFixedPressure()
+orig.setFluidTopNoFlow()
orig.setFluidBottomNoFlow()
orig.initTemporal(total = 0.5, file_dt = 0.01, dt = 1.0e-6)
#orig.time_dt[0] *= 0.01
t@@ -230,14 +230,7 @@ orig.initTemporal(total = 0.5, file_dt = 0.01, dt = 1.0e-6)
orig.run(verbose=False)
py.readlast(verbose = False)
-# Fluid flow direction, opposite of gradient (i.e. towards -z)
-if ((py.v_f[:,:,:,2] < 0.0).all()):
- print("Flow field (Z):\t\t" + passed())
-else:
- print("Flow field (Z):\t\t" + failed())
- raise Exception("Failed")
-
-# Fluid flow direction, opposite of gradient (i.e. towards -z)
+# Fluid flow direction, opposite of gradient (i.e. towards +x)
if ((py.v_f[:,:,:,0] > 0.0).all()):
print("Flow field (X):\t\t" + passed())
else: