Enforcing velocity of particle in a periodic cell
Dear sir/ madam,
As part of a larger setup, I have an issue in enforcing motion of particles in a periodic cell.
When running the simulation without the periodic cell, the particle moves at the specified velocity for infinite iterations.
However, when the periodic boundary conditions are activated, the simulation runs for some iterations and can then no longer proceed (without pausing or error prompt)
Is there some way to work around this issue or an explanation why this is happening?
In order to summarize the problem, I have attached a shorter script, I kindly ask you to comment and uncomment the periodic boundary section to see the effect.
Looking forward to your reply...
#######
from yade import qt, ymport, pack, plot, qt, wrapper
import numpy as np
Mat1 = O.materials.
rad = 0.001
O.bodies.
O.bodies.
vel = 0.01
L = 20*rad
w = 10*rad
h = 5*rad
# setup the periodic boundary (comment and uncomment)
O.periodic = True
O.cell.hSize = Matrix3(L, 0, 0, 0, h, 0, 0, 0, w)
O.engines = [
ForceResett
InsertionSo
Interaction
),
NewtonInteg
PyRunner(
]
O.dt = 0.8*utils.
qt.View()
O.pause()
def init(vel):
#######
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: