tCorrected interaction force in velocity prediction - 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 2d996a75c458c8904aded32488f73e5f37dc9135
(DIR) parent eadf3fe7a64f77d94b1a24a69011372c629dff04
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Thu, 3 Apr 2014 10:47:38 +0200
Corrected interaction force in velocity prediction
Diffstat:
M src/navierstokes.cuh | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
---
(DIR) diff --git a/src/navierstokes.cuh b/src/navierstokes.cuh
t@@ -1691,9 +1691,9 @@ __global__ void findPredNSvelocities(
Float3 v_p = v
+ pressure_term
+ 1.0/devC_params.rho_f*div_phi_tau*devC_dt/phi
- + devC_dt
- *MAKE_FLOAT3(devC_params.g[0], devC_params.g[1], devC_params.g[2])
- - devC_dt*(f_i)
+ + MAKE_FLOAT3(devC_params.g[0], devC_params.g[1], devC_params.g[2])
+ *devC_dt
+ - devC_dt/(devC_params.rho_f*phi)*f_i
- v*dphi/phi
- div_phi_vi_v*devC_dt/phi;