tsplit permeability experiments into two-per-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 0908b7b6fe945c097dfa569c738cf3d0b7fff89b
 (DIR) parent 7f0969d3d4369bfee18b6723b4f1440cf0f42536
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Fri,  4 Jul 2014 15:02:40 +0200
       
       split permeability experiments into two-per-test
       
       Diffstat:
         A python/permeability1.py             |      37 +++++++++++++++++++++++++++++++
         A python/permeability2.py             |      37 +++++++++++++++++++++++++++++++
         A python/permeability3.py             |      37 +++++++++++++++++++++++++++++++
       
       3 files changed, 111 insertions(+), 0 deletions(-)
       ---
 (DIR) diff --git a/python/permeability1.py b/python/permeability1.py
       t@@ -0,0 +1,37 @@
       +#!/usr/bin/env python
       +import sphere
       +import numpy
       +
       +
       +for dp in [1.0e3, 2.0e3]:
       +    # Read initial configuration
       +    sim = sphere.sim('diffusivity-relax')
       +    sim.readlast()
       +
       +    sim.sid = 'permeability-dp=' + str(dp)
       +    sim.cleanup()
       +
       +    sim.g[2] = 0.0
       +    sim.nw[0] = 0
       +    sim.initGrid()
       +    sim.zeroKinematics()
       +    sim.initFluid(mu = 17.87e-4, p = 1.0e5, hydrostatic=True)
       +
       +    # Initialize to linear hydraulic gradient
       +    p_bottom = 10.0
       +    p_top = p_bottom + dp
       +    dz = sim.L[2]/sim.num[2]
       +    for iz in range(sim.num[2]-1):
       +        #z = dz*iz + 0.5*dz # cell-center z-coordinate
       +        z = dz*iz
       +        sim.p_f[:,:,iz] = p_bottom + dp/sim.L[2] * z
       +
       +    sim.setFluidTopFixedPressure()
       +    sim.setFluidBottomFixedPressure()
       +    sim.p_f[:,:,-1] = p_top
       +    sim.setDEMstepsPerCFDstep(10)
       +    sim.initTemporal(total = 2.0, file_dt = 0.01, epsilon=0.07)
       +
       +    sim.run(dry=True)
       +    sim.run()
       +    sim.writeVTKall()
 (DIR) diff --git a/python/permeability2.py b/python/permeability2.py
       t@@ -0,0 +1,37 @@
       +#!/usr/bin/env python
       +import sphere
       +import numpy
       +
       +
       +for dp in [4.0e3, 10.0e3]:
       +    # Read initial configuration
       +    sim = sphere.sim('diffusivity-relax')
       +    sim.readlast()
       +
       +    sim.sid = 'permeability-dp=' + str(dp)
       +    sim.cleanup()
       +
       +    sim.g[2] = 0.0
       +    sim.nw[0] = 0
       +    sim.initGrid()
       +    sim.zeroKinematics()
       +    sim.initFluid(mu = 17.87e-4, p = 1.0e5, hydrostatic=True)
       +
       +    # Initialize to linear hydraulic gradient
       +    p_bottom = 10.0
       +    p_top = p_bottom + dp
       +    dz = sim.L[2]/sim.num[2]
       +    for iz in range(sim.num[2]-1):
       +        #z = dz*iz + 0.5*dz # cell-center z-coordinate
       +        z = dz*iz
       +        sim.p_f[:,:,iz] = p_bottom + dp/sim.L[2] * z
       +
       +    sim.setFluidTopFixedPressure()
       +    sim.setFluidBottomFixedPressure()
       +    sim.p_f[:,:,-1] = p_top
       +    sim.setDEMstepsPerCFDstep(10)
       +    sim.initTemporal(total = 2.0, file_dt = 0.01, epsilon=0.07)
       +
       +    sim.run(dry=True)
       +    sim.run()
       +    sim.writeVTKall()
 (DIR) diff --git a/python/permeability3.py b/python/permeability3.py
       t@@ -0,0 +1,37 @@
       +#!/usr/bin/env python
       +import sphere
       +import numpy
       +
       +
       +for dp in [20.0e3, 40.0e3]:
       +    # Read initial configuration
       +    sim = sphere.sim('diffusivity-relax')
       +    sim.readlast()
       +
       +    sim.sid = 'permeability-dp=' + str(dp)
       +    sim.cleanup()
       +
       +    sim.g[2] = 0.0
       +    sim.nw[0] = 0
       +    sim.initGrid()
       +    sim.zeroKinematics()
       +    sim.initFluid(mu = 17.87e-4, p = 1.0e5, hydrostatic=True)
       +
       +    # Initialize to linear hydraulic gradient
       +    p_bottom = 10.0
       +    p_top = p_bottom + dp
       +    dz = sim.L[2]/sim.num[2]
       +    for iz in range(sim.num[2]-1):
       +        #z = dz*iz + 0.5*dz # cell-center z-coordinate
       +        z = dz*iz
       +        sim.p_f[:,:,iz] = p_bottom + dp/sim.L[2] * z
       +
       +    sim.setFluidTopFixedPressure()
       +    sim.setFluidBottomFixedPressure()
       +    sim.p_f[:,:,-1] = p_top
       +    sim.setDEMstepsPerCFDstep(10)
       +    sim.initTemporal(total = 2.0, file_dt = 0.01, epsilon=0.07)
       +
       +    sim.run(dry=True)
       +    sim.run()
       +    sim.writeVTKall()