tfix error in max permeability value - sphere - GPU-based 3D discrete element m… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit 686048360a86c7bfc8b4774e76ec8e2f63b197a2 | |
parent 3bb51c0797c7460da3e886a6f3d8051cfb6989d0 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Thu, 6 Nov 2014 14:43:27 +0100 | |
fix error in max permeability value | |
Diffstat: | |
M src/darcy.cuh | 4 +++- | |
1 file changed, 3 insertions(+), 1 deletion(-) | |
--- | |
diff --git a/src/darcy.cuh b/src/darcy.cuh | |
t@@ -379,6 +379,7 @@ __global__ void findDarcyPorosities( | |
// Make sure that the porosity is in the interval [0.0;1.0] | |
phi = fmin(0.9, fmax(0.1, void_volume/cell_volume)); | |
+ //phi = fmin(0.99, fmax(0.01, void_volume/cell_volume)); | |
//phi = void_volume/cell_volume; | |
Float dphi = phi - phi_0; | |
t@@ -642,7 +643,8 @@ __global__ void findDarcyPermeabilities( | |
// limit permeability [m*m] | |
// K_gravel = 3.0e-2 m/s => k_gravel = 2.7e-9 m*m | |
- k = fmin(2.7e-8, k); | |
+ //k = fmin(2.7e-9, k); | |
+ k = fmin(2.7e-10, k); | |
__syncthreads(); | |
dev_darcy_k[cellidx] = k; |