From be48dae6e6f4e461864ef38d31271294ac428245 Mon Sep 17 00:00:00 2001
From: Rafael Ravedutti <rafaelravedutti@gmail.com>
Date: Mon, 6 Nov 2023 13:24:26 +0100
Subject: [PATCH] Adjust gravity calculation
Signed-off-by: Rafael Ravedutti <rafaelravedutti@gmail.com>
---
examples/dem.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/examples/dem.py b/examples/dem.py
index 0d86063..be9afe4 100644
--- a/examples/dem.py
+++ b/examples/dem.py
@@ -89,7 +89,7 @@ def euler(i):
def gravity(i):
volume = (4.0 / 3.0) * pi * radius[i] * radius[i] * radius[i]
- force[i][2] += -(densityParticle_SI - densityFluid_SI) * volume * gravity_SI
+ force[i][2] = force[i][2] - (densityParticle_SI - densityFluid_SI) * volume * gravity_SI
cmd = sys.argv[0]
@@ -112,7 +112,7 @@ restitutionCoefficient = 0.1
collisionTime_SI = 5e-4
poissonsRatio = 0.22
timeSteps = 10000
-visSpacing = 1
+visSpacing = 100
denseBottomLayer = False
bottomLayerOffsetFactor = 1.0
kappa = 2.0 * (1.0 - poissonsRatio) / (2.0 - poissonsRatio) # from Thornton et al
--
GitLab