tAdded fluid shear test - 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 980a99ca05ed04ff9cb51ed7fecef8bd6fa6f554
 (DIR) parent 874be72bf7a0fa969d777e4ad5e24061e6c7047a
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Tue,  1 Apr 2014 09:21:12 +0200
       
       Added fluid shear test
       
       Diffstat:
         A python/fluidshear.py                |      42 +++++++++++++++++++++++++++++++
       
       1 file changed, 42 insertions(+), 0 deletions(-)
       ---
 (DIR) diff --git a/python/fluidshear.py b/python/fluidshear.py
       t@@ -0,0 +1,42 @@
       +#!/usr/bin/env python
       +
       +import sphere
       +
       +sid = 'fluidshear'
       +
       +## Initialization from loose packing to a gravitationally collapsed state
       +## without fluids
       +sim = sphere.sim(sid + '-init', np = 2400, fluid = False)
       +#sim.cleanup()
       +sim.radius[:] = 0.05
       +sim.initRandomGridPos(gridnum = [12, 12, 9000])
       +sim.initTemporal(total = 5.0, file_dt = 0.05)
       +sim.g[2] = -9.81
       +sim.run(dry=True)
       +#sim.run()
       +#sim.writeVTKall()
       +
       +## Consolidation from a top wall without fluids
       +sim.readlast()
       +sim.sid = sid + '-cons'
       +sim.adjustUpperWall()
       +sim.consolidate(normal_stress = 10.0e3)
       +#sim.initFluid(mu = 17.87e-4, p = 1.0e5, hydrostatic = True)  # mu = water at 0 deg C
       +sim.initTemporal(total = 1.0, file_dt = 0.01)
       +sim.run(dry=True)
       +#sim.run()
       +#sim.writeVTKall()
       +sim.visualize('walls')
       +
       +## Shear with fluids
       +sim.readlast()
       +sim.sid = sid + '-shear'
       +sim.shear()
       +sim.fluid = True
       +sim.initFluid(mu = 17.87e-4, p = 1.0e5, hydrostatic = True)
       +sim.bc_bot[0] = 1  # Neumann BC
       +sim.initTemporal(total = 1.0, file_dt = 0.01)
       +sim.run(dry=True)
       +sim.run()
       +sim.writeVTKall()
       +sim.visualize('shear')