Contact with AllowBiggerThanPeriod
Hello,
I've noticed a weird behaviour that I cannot explain. Below there is a simple script reproducing the problem.
I have two spheres in a periodic cell. One sphere (blue in my script) is bigger than the period, so I am using the condition AllowBiggerThan
The problem comes when the width of the cell is exactly the same as the diameter of the big blue sphere, or slightly smaller. In these cases I have no interaction between the two sphere: one spheres goes through the other, and then the big sphere reaches the end of the cell, is remapped on the other side, and only at the second time finally hits the other ball.
And this only happens when the width of the cell is between 0.6 and 1 time the diameter of the sphere, and when I generated the big blue sphere first and the small red sphere after.
I am aware that it is a really specific case, but still I do not understand why it happens.. am I missing something? And why does the order in which I'm generating the particles affect the test?
Thank you very much
Giulia
#!/usr/bin/python
# -*- coding: utf-8 -*-
DbigSphere = 1.
radius_bigSphere= DbigSphere/2
radius_sph = DbigSphere/10
# Periodic cell
#cell0 = DbigSphere*
#cell0 = DbigSphere*0.605 # ok - Anything smaller than that works fine
cell0 = DbigSphere*1 # does NOT work
cell1 = DbigSphere*3.
cell2 = DbigSphere*3
O.periodic = True
O.cell.
# Big blue sphere
center_bigSphere = Vector3(
bigSphere = utils.sphere(
# Red small sphere
center_sph = Vector3(
sphere = utils.sphere(
#O.bodies.
O.bodies.
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop ([Ig2_Sphere_
ForceEngine (force=
NewtonIntegrator (damping=0.5),
PyRunner (command=
]
flagRFI = 0
def recordFirstInte
global flagRFI,iterFirst
try:
i=O.interacti
print'i.isReal ?',i.isReal,
if i.isReal:
if flagRFI==0:
iterFirst=
flagRFI=1
O.pause()
except StandardError:
i=0
O.dt=0.
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Giulia Macaro
- Solved:
- Last query:
- Last reply: