tdisable vertical pressure forces at the top wall - 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 2fcfb8e66a446d872d8774f9d149a72649ac1e0f
 (DIR) parent 62db27fa140c1c6318255997a5d8cde1cda89f92
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon, 17 Nov 2014 10:28:23 +0100
       
       disable vertical pressure forces at the top wall
       
       Diffstat:
         M src/darcy.cuh                       |      10 +++++++---
       
       1 file changed, 7 insertions(+), 3 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -541,9 +541,9 @@ __global__ void findDarcyPressureForce(
                Float p_zp = dev_darcy_p[d_idx(i_x,i_y,i_z+1)];
        
                // Add Neumann BC at top wall
       -        if (i_z >= wall0_iz - 1)
       +        /*if (i_z >= wall0_iz - 1)
                    p_zp = p;
       -            //p_zp = p_zn;
       +            //p_zp = p_zn;*/
        
                // find particle volume (radius in x.w)
                const Float V = 4.0/3.0*M_PI*x.w*x.w*x.w;
       t@@ -557,12 +557,16 @@ __global__ void findDarcyPressureForce(
                // find pressure gradient force plus buoyancy force.
                // buoyancy force = weight of displaced fluid
                // f_b = -rho_f*V*g
       -        const Float3 f_p = -1.0*grad_p*V/(1.0 - phi);
       +        Float3 f_p = -1.0*grad_p*V/(1.0 - phi);
                    //- devC_params.rho_f*V*MAKE_FLOAT3(
                            //devC_params.g[0],
                            //devC_params.g[1],
                            //devC_params.g[2]);
        
       +        // Add Neumann BC at top wall
       +        if (i_z >= wall0_iz - 1)
       +            f_p.z = 0.0;
       +
                /*printf("%d,%d,%d findPF:\n"
                        "\tphi    = %f\n"
                        "\tp      = %f\n"