Particle Interactions

Asked by Jeff Steadman

I'm trying to calibrate rotational particle interaction values. In setting up the problem I noticed with the absence of applied forces, gravity and, damping that my particles were moving. Upon investigation I found that the particle interaction values are causing the particles to move. As I adjust the interaction values the particles move in different manners. How is this possible? Since no forces are being applied the particles should never move. My code is shown below.

I realize the interaction values in the code are not reasonable but, the values should not cause the particles to move in the absence of forces.

Thanks.

------------------------------------------------------------------------------------------------------------------------------------------------------------------
from esys.lsm import *
from esys.lsm.util import *
from esys.lsm.geometry import *
from POVsnaps import POVsnaps

sim = LsmMpi(numWorkerProcesses=8, mpiDimList=[4,2,1])
sim.initNeighbourSearch(
   particleType="RotSphere",
   gridSpacing=.625,
   verletDist=0.05
)

sim.setNumTimeSteps(100)
sim.setTimeStepSize(1)

domain = BoundingBox(Vec3(-20,-20,-20), Vec3(20,20,20))
sim.setSpatialDomain(domain)

r = .25 #radius

cube = CubicBoxPacker (
   radius = r,
   bBox = BoundingBox(
      Vec3(-0.50,1.75,-0.00),
      Vec3(4.50,2.25,0.00)
   ),
   circDimList = [False, False, False],
   tolerance = 1.0000e-05
)
cube.generate()
cube_particles = cube.getSimpleSphereCollection()

sim.createParticles(cube_particles)

sim.createConnections(
   ConnectionFinder(
      maxDist = 0.5,
      bondTag = 0,
      pList = cube_particles
   )
)

pp_bonds = sim.createInteractionGroup (
   RotBondPrms(
      name="pp_bonds",
      normalK=0.0005,
      shearK=0.0005,
      torsionK=0.0005,
      bendingK=0.0005,
      normalBrkForce=1000,
      shearBrkForce=1000,
      torsionBrkForce=1000,
      bendingBrkForce=1000,
      tag=0,
      scaling = True
   )
)

povcami = POVsnaps(sim=sim, interval=1)
povcami.configure(lookAt=Vec3(2,2,0), camPosn=(2,2,5), zoomFactor=0.2)
sim.addPostTimeStepRunnable(povcami)

sim.run()

Question information

Language:
English Edit question
Status:
Answered
For:
ESyS-Particle Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Jesús Fernando Vergara Núñez (jesufer1) said :
#1

Well, I notice that you have created connections with a maximum distance of 0.5 m between the particles. Any two particles within a distance between -1x10^5 and 0.5, and no 0, also creates a connection: there is a distance different from zero which counts in the computation of the bonding force.

I suggest you to choose a lower maxDist, because a value of 0.5 is greater than the distance between two diagonal particles, because they are suppose to have no connections, and because the distance of separation is 0.5*1.4142 - 0.5 = 0.2071 > 0.5.

And try to run the simulation with your new data, and reach 100000 timesteps or more to see what happens.

Revision history for this message
Jesús Fernando Vergara Núñez (jesufer1) said :
#2

Note: 0.2071 is the maximum distance in a cubic lattice arrangement.

Can you help with this problem?

Provide an answer of your own, or ask Jeff Steadman for more information if necessary.

To post a message you must log in.