Hello,
I want to fill up a cylinder with particles. I try the following code, but get an error "Expected integer parameter in input script or data file". Since LIGGGHTS does not tell which line has the error I tried changing every float value to integer to see what's wrong, but I could not find the problem.
I would be grateful if someone helps me out! The code's attached herewith.
Attachment | Size |
---|---|
cylinder.txt | 1.5 KB |
mschramm | Thu, 11/28/2019 - 20:52
Echo both
Does it happen right away?
Using the -echo both usually helps narrow down where the error is.
Nothing in your input script jumps out at me.
Saurabh Das | Fri, 11/29/2019 - 00:03
Happens after 'neigh_modify delay 0.0001'
Thanks for the tip. The error comes up after the line 'neigh_modify delay 0.0001' statement, it seems. How do I rectify it?
---------------------------log file----------------------------------------
LIGGGHTS (Version LIGGGHTS-PUBLIC 3.7.0, compiled 2017-07-18-23:47:28 by root, git commit unknown)
#Stirred Mill
atom_style granular
boundary f f f
newton off
communicate single vel yes
units si
#region
region outer_cylinder cylinder y 0 0 0.18 0 0.148 units box
create_box 1 outer_cylinder
Created orthogonal box = (-0.18 0 -0.18) to (0.18 0.148 0.18)
1 by 1 by 1 MPI processor grid
neighbor 0.05 bin
neigh_modify delay 0.0001
ERROR: Expected integer parameter in input script or data file (/build/liggghts-LzNUbp/liggghts-3.7.0+repack1/src/neighbor.cpp:1914)
---------------------------------log file-------------------------------------------------------
mschramm | Fri, 11/29/2019 - 05:07
Delay needs to be an int
neigh_modify delay 0.0001 is incorrect.
Delay tells liggghts how many steps to wait until the neighbor list is updated.
neigh_modify delay 0 should get the code working.
Saurabh Das | Fri, 11/29/2019 - 07:49
"Numeric index is out of bounds" error after pa
The previous error got sorted out, but now this error appears, after pair style description.
-----------------------------------------------------------------------
#Stirred Mill
atom_style granular
boundary f f f
newton off
communicate single vel yes
units si
#region
region outer_cylinder cylinder y 0 0 0.18 0 0.148 units box
create_box 1 outer_cylinder
Created orthogonal box = (-0.18 0 -0.18) to (0.18 0.148 0.18)
1 by 1 by 1 MPI processor grid
neighbor 0.05 bin
neigh_modify delay 0
#material properties
fix m1 all property/global youngsModulus peratomtype 5e6
fix m1 all property/global poissonsRatio peratomtype 0.45
fix m1 all property/global coefficientRestitution peratomtypepair 1 0.9
fix m1 all property/global coefficientFriction peratomtypepair 1 0.05
fix m1 all property/global characteristicVelocity scalar 2e-2
#pair style
pair_style gran model hooke tangential history
pair_coeff **
ERROR: Numeric index is out of bounds (/build/liggghts-LzNUbp/liggghts-3.7.0+repack1/src/force.cpp:785)
--------------------------------------------------------------------
mschramm | Fri, 11/29/2019 - 14:51
Fix names
You cannot use the name m1 for all of your fix names. Try,
fix m1 all property/global youngsModulus peratomtype 5e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.9
fix m4 all property/global coefficientFriction peratomtypepair 1 0.05
fix m5 all property/global characteristicVelocity scalar 2e-2
Saurabh Das | Mon, 12/02/2019 - 08:58
Still getting the same errors
Still getting the same errors, I am afraid.
mschramm | Mon, 12/02/2019 - 16:57
pair_coeff is incorrect
pair_coeff ** should be pair_coeff * *
there should be a space between the *.
This then leads to another error "Unable to calculate region volume."
You can get it to run by using "all_in off", but there is an additional problem that makes particles "explode".