tdisplay values at a porosity of 0.3 - 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 48431f5c15a7f588b18bda59e3226a21515da310
   DIR parent e00575a0c51ed0d98599781875be8f9e83f66bc1
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 12 Nov 2014 14:16:08 +0100
       
       display values at a porosity of 0.3
       
       Diffstat:
         M python/sphere.py                    |       4 ++++
       
       1 file changed, 4 insertions(+), 0 deletions(-)
       ---
   DIR diff --git a/python/sphere.py b/python/sphere.py
       t@@ -3138,11 +3138,15 @@ class sim:
                    self.k_c[0] = k_c
                    phi = numpy.array([0.1, 0.9])
                    k = self.k_c * phi**3/(1.0 - phi**2)
       +            k_phi30 = self.k_c * 0.3**3/(1.0 - 0.3**2)
                    K = k * self.rho*numpy.abs(self.g[2])/self.mu
       +            K_phi30 = k_phi30 = self.rho*numpy.abs(self.g[2])/self.mu
                    print('Hydraulic permeability limits for porosity phi = [0.1, 0.9]:')
                    print('\tk = ' + str(k) + ' m*m')
       +            print('\tk(phi = 0.3) = ' + str(k_phi30) + ' m*m')
                    print('Hydraulic conductivity limits for porosity phi = [0.1, 0.9]:')
                    print('\tK = ' + str(K) + ' m/s')
       +            print('\tK(phi = 0.3) = ' + str(K_phi30) + ' m/s')
                else:
                    raise Exception('setPermeabilityPrefactor() only relevant for the '
                            + 'Darcy solver (cfd_solver = 1)')