twall thickness should depend on particle size not grid size - 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 a1f59ee7e5e8915d203ddc5d7203700eb8884dd2
 (DIR) parent 4208da212cc52b671a382207ec692867a4ba660a
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon, 31 Mar 2014 11:00:28 +0200
       
       wall thickness should depend on particle size not grid size
       
       Diffstat:
         M python/sphere.py                    |       3 ++-
       
       1 file changed, 2 insertions(+), 1 deletion(-)
       ---
 (DIR) diff --git a/python/sphere.py b/python/sphere.py
       t@@ -2051,6 +2051,7 @@ class sim:
                xmax = numpy.max(self.x[:,dim] + self.radius)
        
                cellsize = self.L[0] / self.num[0]
       +        d_max = numpy.max(self.radius)*2.0
        
                self.num[dim] = numpy.ceil(((xmax-xmin)*adjust + xmin)/cellsize)
                self.L[dim] = (xmax-xmin)*adjust + xmin
       t@@ -2062,7 +2063,7 @@ class sim:
                    self.w_x[idx] = numpy.array([xmin])
                #self.w_m[idx] = numpy.array([self.rho[0]*self.np*math.pi \
                #        *(cellsize/2.0)**3])
       -        self.w_m[idx] = numpy.array([self.rho*self.L[0]*self.L[1]*cellsize])
       +        self.w_m[idx] = numpy.array([self.rho*self.L[0]*self.L[1]*d_max])
        
            def consolidate(self, normal_stress = 10e3):
                '''