Bond atoms .. missing on proc

Submitted by boumediene on Thu, 07/18/2013 - 11:26

Hi!
I have a little trouble with the command "bond_type harmonic"
if I use this command (bond_type harmonic) under LIGGGHTS works without any problem,
when cfdem coupling, I got this message:
"
Courant Number mean: 0.00147829 max: 7.37258
Starting up LIGGGHTS
Executing command: 'run 10 pre no post no'
run 10 pre no post no
Step Atoms KinEng 1 Volume
91101 61 0.18215281 0.00059623725 0.00011248163
ERROR on proc 7: Bond atoms 53 54 missing on proc 7 at step 91104 (neigh_bond.cpp:60)
"

My file data is resumed here:
-------------------------------------------------
region reg block 0 0.066 0 0.01 0.02 0.18 units box
create_box 1 reg bonds 1 8
neighbor 0.002 bin
neigh_modify delay 0

bond_style harmonic
bond_coeff 1 0.3 0.001
region bc block 0.003 0.063 0.005 0.0055 0.155 0.1555 units box
lattice sc 0.001
create_atoms 1 region bc
group bonded region bc
set group all density 2000 diameter 0.001

special_bonds lj/coul 0 1 1 extra 8
mass 1 1.0
...
pair_style gran/hertz/history
pair_coeff * *
fix bdr all bond/create 1 1 1 0.001 1
....
-------------------------------------------------------------------------
I used before the command "bond_type gran"
I had the same message despite SEVERAL tests, so I decided to use the "type harmonic" to have an end of my problem simulation.

another point, when i using a pressure gradient , it will give me a message from the firsts iterations.

Thanks in advance