problem using the hopper example with smaller particles
Hi,
Using the clump-hopper-test located at trunk/examples/
Is it possible to run this gravity simulation with small (50e-09) particles? Reducing the timestep even further doesnt seem to be the way to go here.
Your help always appreciated
Jesse
The script I have pasted below.
from yade import pack,export,qt,geom
import itertools
from numpy import *
scalefac = 1e-6 #be careful of timestep
kinEnergyMax = 100000
# Parameters
tc=0.001# collision time
en=.3 # normal restitution coefficient
es=.3 # tangential restitution coefficient
frictionAngle=
density=2700
# facets material
facetMat=
# default spheres material
dfltSpheresMat=
#O.dt=.05*tc # time step
O.dt=.05*tc * 1e-3# time step
#O.dt= 1e-09
Rs=0.05 # particle radius
# Create geometry
x0=0.; y0=0.; z0=0.; ab=.7; at=2.; h=1.; hl=h; al=at*3
Rs *= scalefac
ab *= scalefac
at *= scalefac
h *= scalefac
zb=z0; x0b=x0-ab/2.; y0b=y0-ab/2.; x1b=x0+ab/2.; y1b=y0+ab/2.
zt=z0+h; x0t=x0-at/2.; y0t=y0-at/2.; x1t=x0+at/2.; y1t=y0+at/2.
zl=z0-hl;
vibrationPlate = O.bodies.
dBunker=
wallMask=
# Create engines
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
DomainLimiter(
NewtonIntegrat
PyRunner(
PyRunner(
]
numSphereGen = 0
def addBodies():
global numSphereGen
# Create clumps...
clumpColor=(0.0, 0.5, 0.5)
for k,l in itertools.
#clpId,
clpId,
numSphereGen += len(sphId)
# ... and spheres
spheresColor=(0.4, 0.4, 0.4)
for k,l in itertools.
#sphAloneId=
sphAloneId=
numSphereGen += len(sphAloneId)
#O.bodies.
def state():
global numSphereGen
print "Iter %d: Total number of generated spheres %d, removed particles %d, current particles %d, kinEnergy %g"%(O.iter, numSphereGen, domLim.nDeleted, numSphereGen-
if (utils.
print "Kinetic energy is over a threshold value! Error!"
addBodies()
from yade import qt
qt.View()
O.saveTmp()
#O.run(50001)
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask rhaven for more information if necessary.