timprove permeability estimation function - 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 2f4bac35a58c5e192e3aafb331545f19badbb523
 (DIR) parent ce939231ee80ea8b4ac5b23cf76f2c534d2a6dcd
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon,  2 Feb 2015 15:02:01 +0100
       
       improve permeability estimation function
       
       Diffstat:
         M python/sphere.py                    |      15 ++++++++++-----
       
       1 file changed, 10 insertions(+), 5 deletions(-)
       ---
 (DIR) diff --git a/python/sphere.py b/python/sphere.py
       t@@ -2872,15 +2872,20 @@ class sim:
                '''
                self.dx = self.L/self.num
        
       -    def hydraulicConductivity(self):
       +    def hydraulicConductivity(self, phi=0.35):
                '''
                Determine the hydraulic conductivity (K) [m/s] from the permeability
       -        prefactor. This value is stored in `self.K_c`. This function only works
       -        for the Darcy solver (`self.cfd_solver == 1`)
       +        prefactor and a chosen porosity.  This value is stored in `self.K_c`.
       +        This function only works for the Darcy solver (`self.cfd_solver == 1`)
       +
       +        :param phi: The porosity to use in the Kozeny-Carman relationship
       +        :type phi: float
       +        :returns: The hydraulic conductivity [m/s]
       +        :return type: float
                '''
                if self.cfd_solver[0] == 1:
       -            self.K_c = self.k_c*self.rho_f*numpy.abs(self.g[2])/self.mu
       -            return self.K_c
       +            k = self.k_c * phi**3/(1.0 - phi**2)
       +            self.K_c = k*self.rho_f*numpy.abs(self.g[2])/self.mu
                else:
                    raise Exception('This function only works for the Darcy solver')