tpermeability BC suitable for shear stress experiments - 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 ad51c8045076b389f7ee142b3b4bcdb3c27da4a5
 (DIR) parent d3bb81a40174dfb91af6e86033020917a20d00f0
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon,  9 Feb 2015 15:29:23 +0100
       
       permeability BC suitable for shear stress experiments
       
       Diffstat:
         M src/darcy.cuh                       |       8 ++++----
       
       1 file changed, 4 insertions(+), 4 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -976,8 +976,8 @@ __global__ void firstDarcySolution(
                // Dirichlet BC at dynamic top wall. wall0_iz will be larger than the
                // grid if the wall isn't dynamic
                if ((bc_bot == 0 && z == 0) || (bc_top == 0 && z == nz-1)
       -                || (z >= wall0_iz-1 && bc_top == 0))
       -                //|| (z >= wall0_iz && bc_top == 0))
       +                //|| (z >= wall0_iz-1 && bc_top == 0))
       +                || (z >= wall0_iz && bc_top == 0))
                    dp_expl = 0.0;
        
        #ifdef REPORT_FORCING_TERMS
       t@@ -1129,8 +1129,8 @@ __global__ void updateDarcySolution(
                // Dirichlet BC at dynamic top wall. wall0_iz will be larger than the
                // grid if the wall isn't dynamic
                if ((bc_bot == 0 && z == 0) || (bc_top == 0 && z == nz-1)
       -                || (z >= wall0_iz-1 && bc_top == 0))
       -                //|| (z >= wall0_iz && bc_top == 0))
       +                //|| (z >= wall0_iz-1 && bc_top == 0))
       +                || (z >= wall0_iz && bc_top == 0))
                    dp_impl = 0.0;
                    //p_new = p;