Tutorial problems

Submitted by mytxii on Thu, 04/16/2015 - 14:45

Hey guys,

I try to do the hopper tutorial. Therefore I take following code :

# Import mesh from cad:
fix cad1 all mesh/surface file meshes/hopper.stl type 2 scale 0.001
# Use the imported mesh as granular wall
fix geometry all wall/gran/hertz/history mesh n_meshes 1 meshes cad1
# Create stopper for funnel
fix stopper all wall/gran/hertz/history primitive type 2 zplane0.0
# Define the physics
pair_style gran/hertz/history
pair_coeff * *

If I start the runscript it shows the following message:
LIGGGHTS (Version LIGGGHTS-PUBLIC 3.1.0, compiled 2015-04-16-11:55:31 by nils, git commit b3030a8fd1f722c77f69150c7b2d24b386eab848 based on LAMMPS 23 Nov 2013)
Created orthogonal box = (-0.138 -0.138 -0.0045) to (0.138 0.138 0.43)
1 by 1 by 1 MPI processor grid
Fix particledistribution/discrete (id pdd): distribution based on mass%:
pts: d=3.000000e-03 (max. bounding sphere) mass%=100.000000%
Fix particledistribution/discrete (id pdd): distribution based on number%:
pts: d=3.000000e-03 (max. bounding sphere) number%=100.000000%
Reading STL file 'meshes/hopper.stl'
ERROR: Invalid fix style: "wall/gran/hertz/history" (../modify.cpp:794)

I can't found anything at google or here.
What's wrong ?? May somebody can help :)

Thank you, guys.

ckloss's picture

ckloss | Fri, 04/17/2015 - 09:31

Hi,

as it is mentioned on the download page

http://www.cfdem.com/tutorials

you'll need to convert the systax from LIGGGHTS 2.X to 3.X - just some minor changes

Christoph

mytxii | Fri, 04/17/2015 - 11:26

Thank you, it helps a lot.

But I got another general question :

If there is an error message (e.g.: ERROR: Could not find fix group ID (../modify.cpp:730)), how can I localize the error ?
(Other errors contain some keywords like "fix m1 ...", .. but that is very universal I guess.)

Nils

cstoltz | Tue, 04/21/2015 - 15:33

As a crude method of debugging something like this, just include a line that contains the word 'stop' somewhere in the input deck. LIGGGHTS won't recognize it, so it will crash at that point and tell me it doesn't recognize 'stop'. Just keep moving it down the deck until you get a crash not related to the 'stop' line, and then you know exactly which line has failed and can then go figure out why. Of course this is really painful if it takes a long time for the crash condition to occur each time.

Like I said - crude, but works.

Regards,
Chris

richti83's picture

richti83 | Wed, 04/22/2015 - 10:16

I think it's more convenient to use echo both at the beginning of the script and look at console output at which line the parser stops.
The error you mentioned (Could not find fix group ID) says that the group ID does not exist, but you did not provide the script, so I guess you changed fix m1 all property/global to fix m1 foo property/global without defining group foo before.
Most of the error messages are self-explanatory when you know the syntax of the command (eg. fix ID GROUP type $EVERY $PARAMETER)
cheers,
Christian

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

anhvt2 | Mon, 08/10/2015 - 22:26

I also encounter the same problem; just like ckloss has said

change

fix geometry all wall/gran/hertz/history mesh n_meshes 1 meshes cad1

to

fix geometry all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad1

hth,
at

anhvt2 | Mon, 08/10/2015 - 22:34

i update a few changes from the original, basically just the syntax.

**********************************************************************************************************************
# bin-flow simulation
# drag from online

### init
# preliminaries

units si
atom_style sphere
atom_modify map array
boundary f f f
newton off
communicate single vel yes
# processor 2 2 3

# declare domain
region domain block -0.138 0.138 -0.138 0.138 -0.0045 0.43 units box
create_box 2 domain

# neighbor listing
neighbor 0.003 bin
neigh_modify every 1 delay 0 check no # add delay 0

### setup

# material and interaction properties req'd
fix m1 all property/global youngsModulus peratomtype 2.5e7 2.5e7
fix m2 all property/global poissonsRatio peratomtype 0.25 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.5 0.5 0.5 0.5
fix m4 all property/global coefficientFriction peratomtypepair 2 0.2 0.175 0.175 0.5

# particle insetion
fix pts all particletemplate/sphere 1 atom_type 1 density constant 1000 radius constant 0.0015
fix pdd all particledistribution/discrete 63243 1 pts 1.0
region factory cylinder z 0 0 0.10 0.4 0.41 units box
fix ins all insert/rate/region seed 123481 distributiontemplate pdd nparticles 150000 particlerate 150000 insert_every 1000 overlapcheck yes vel constant 0. 0. -1.0 region factory ntry_mc 10000

# import mesh from cad
fix cad1 all mesh/surface file hopper.stl type 2 scale 0.001

# use imported mesh as granular wall
fix geometry all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad1

# create stopper for funnel
fix stopper all wall/gran model hertz tangential history primitive type 2 zplane 0.0

# define the physics
# pair_style gran/hertz/history
pair_style gran model hertz tangential history

pair_coeff * *

### detailed settings
# integrator
fix integrate all nve/sphere

# gravity
fix grav all gravity 9.81 vector 0.0 0.0 -1.0

# time-step
timestep 0.00001

# thermodynamic output settings
thermo_style custom step atoms ke cpu
thermo 10000
thermo_modify norm no lost ignore

# check time-step and init dump file
fix ctg all check/timestep/gran 1 0.01 0.01
run 1
unfix ctg

# create imaging info
dump dumpstl all stl 10000 dump*.stl
dump dmp all custom 10000 dump.1 id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius

### execution and further settings

# fill the hopper
run 150000 upto

# pull the hopper
unfix stopper
run 300000

**********************************************************************************************************************

VOny | Sat, 03/25/2017 - 00:32

It turns out that several typo/ none spacing error as well as syntax change from 2.x.x. to 3.x.x make it unable to run, the updated script is really helpful and also put the file in the folder of input file also helps