tMade time step more conservative, validated new integration scheme - 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 e99f67a88ac006ee17fb3945b30a5349cb3e7fc5
 (DIR) parent a8f9bb4ed81e686acb8e1a51b4ba60fedcf03dae
 (HTM) Author: Anders Damsgaard <adc@geo.au.dk>
       Date:   Wed, 14 Nov 2012 10:21:42 +0100
       
       Made time step more conservative, validated new integration scheme
       
       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@@ -743,7 +743,8 @@ class Spherebin:
        
            # Computational time step (O'Sullivan et al, 2003)
            #self.time_dt[0] = 0.17 * math.sqrt((4.0/3.0 * math.pi * r_min**3 * self.rho[0]) / numpy.amax([self.k_n[:], self.k_t[:]]) )
       -    self.time_dt[0] = 0.12 * math.sqrt((4.0/3.0 * math.pi * r_min**3 * self.rho[0]) / numpy.amax([self.k_n[:], self.k_t[:]]) )
       +    # Computational time step (Zhang and Campbell, 1992)
       +    self.time_dt[0] = 0.075 * math.sqrt((4.0/3.0 * math.pi * r_min**3 * self.rho[0]) / numpy.amax([self.k_n[:], self.k_t[:]]) )
            
            # Time at start
            self.time_current[0] = current