tfix lowermost particles so everything doesnt slide horizontally - 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 8443f37adc22d395277cf77eae9fc4ff30424426
   DIR parent b74bec6ccdc448acc13f019a5fdec30d8db66aea
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 13 Aug 2014 14:02:52 +0200
       
       fix lowermost particles so everything doesnt slide horizontally
       
       Diffstat:
         M python/diffusivity-starter.py       |       7 +++++++
         M python/permeability-starter.py      |       4 ++++
       
       2 files changed, 11 insertions(+), 0 deletions(-)
       ---
   DIR diff --git a/python/diffusivity-starter.py b/python/diffusivity-starter.py
       t@@ -39,13 +39,20 @@ for sigma0_str in sys.argv[4:]:
            sim.cleanup()
            sim.adjustUpperWall()
            sim.zeroKinematics()
       +
            sim.consolidate(normal_stress = 10.0e3)
       +
            sim.initFluid(mu = 17.87e-4, p = 1.0e5, hydrostatic = True)
            sim.setFluidBottomNoFlow()
            sim.setFluidTopFixedPressure()
            sim.setDEMstepsPerCFDstep(10)
            sim.setMaxIterations(2e5)
            sim.initTemporal(total = 5.0, file_dt = 0.01, epsilon=0.07)
       +
       +    # Fix lowermost particles
       +    I = numpy.nonzero(sim.x[:,2] < 1.5*dz)
       +    sim.fixvel[I] = 1
       +    
            sim.run(dry=True)
            sim.run(device=0)
            #sim.writeVTKall()
   DIR diff --git a/python/permeability-starter.py b/python/permeability-starter.py
       t@@ -37,6 +37,10 @@ for dp_str in sys.argv[4:]:
                z = dz*iz
                sim.p_f[:,:,iz] = p_bottom + dp/sim.L[2] * z
        
       +    # Fix lowermost particles
       +    I = numpy.nonzero(sim.x[:,2] < 1.5*dz)
       +    sim.fixvel[I] = 1
       +    
            sim.setFluidTopFixedPressure()
            sim.setFluidBottomFixedPressure()
            sim.p_f[:,:,-1] = p_top