tset coefficients - 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 b520cff0a79ba8565828b298c2d1daa3e35fe21b
 (DIR) parent 567405123e189cf6c95a5d790b6b8e0bf2ff2205
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 21 Aug 2014 11:03:34 +0200
       
       set coefficients
       
       Diffstat:
         M python/diffusivity-test.py          |       2 ++
         M python/permeability-results-c=1.py  |      30 ++++++++++++++++--------------
       
       2 files changed, 18 insertions(+), 14 deletions(-)
       ---
 (DIR) diff --git a/python/diffusivity-test.py b/python/diffusivity-test.py
       t@@ -63,6 +63,8 @@ for sigma0 in sigma0_list:
            sim.setDEMstepsPerCFDstep(10)
            sim.setMaxIterations(2e5)
            sim.initTemporal(total = 5.0, file_dt = 0.01, epsilon=0.07)
       +    sim.c_grad_p[0] = c_grad_p
       +    sim.c_phi[0] = c_phi
        
            # Fix lowermost particles
            dz = sim.L[2]/sim.num[2]
 (DIR) diff --git a/python/permeability-results-c=1.py b/python/permeability-results-c=1.py
       t@@ -41,26 +41,28 @@ for sid in sids:
        
        fig = plt.figure()
        
       -plt.subplot(3,1,1)
       -plt.xlabel('Pressure gradient $\\Delta p/\\Delta z$ [Pa m$^{-1}$]')
       +#plt.subplot(3,1,1)
       +plt.xlabel('Pressure gradient $\\Delta p/\\Delta z$ [kPa m$^{-1}$]')
        plt.ylabel('Hydraulic conductivity $K$ [ms$^{-1}$]')
       -plt.plot(dpdz, K, '+')
       +plt.ticklabel_format(style='sci', axis='y', scilimits=(0,0))
       +dpdz /= 1000.0
       +plt.plot(dpdz, K, 'o-k')
        plt.grid()
        
       -plt.subplot(3,1,2)
       -plt.xlabel('Pressure gradient $\\Delta p/\\Delta z$ [Pa m$^{-1}$]')
       -plt.ylabel('Hydraulic flux $Q$ [m$^3$s$^{-1}$]')
       -plt.plot(dpdz, Q, '+')
       -plt.grid()
       +#plt.subplot(3,1,2)
       +#plt.xlabel('Pressure gradient $\\Delta p/\\Delta z$ [Pa m$^{-1}$]')
       +#plt.ylabel('Hydraulic flux $Q$ [m$^3$s$^{-1}$]')
       +#plt.plot(dpdz, Q, '+')
       +#plt.grid()
        
       -plt.subplot(3,1,3)
       -plt.xlabel('Pressure gradient $\\Delta p/\\Delta z$ [Pa m$^{-1}$]')
       -plt.ylabel('Mean porosity $\\bar{\\phi}$ [-]')
       -plt.plot(dpdz, phi_bar, '+')
       -plt.grid()
       +#plt.subplot(3,1,3)
       +#plt.xlabel('Pressure gradient $\\Delta p/\\Delta z$ [Pa m$^{-1}$]')
       +#plt.ylabel('Mean porosity $\\bar{\\phi}$ [-]')
       +#plt.plot(dpdz, phi_bar, '+')
       +#plt.grid()
        
        plt.tight_layout()
       -filename = 'permeability-dpdz-vs-K.png'
       +filename = 'permeability-dpdz-vs-K.pdf'
        plt.savefig(filename)
        print(os.getcwd() + '/' + filename)
        plt.savefig(filename)